diff options
Diffstat (limited to 'src/lsq_rectilinear.mac')
| -rw-r--r-- | src/lsq_rectilinear.mac | 60 |
1 files changed, 0 insertions, 60 deletions
diff --git a/src/lsq_rectilinear.mac b/src/lsq_rectilinear.mac deleted file mode 100644 index 95ed273..0000000 --- a/src/lsq_rectilinear.mac +++ /dev/null @@ -1,60 +0,0 @@ -/* - * Copyright 2006-2009 Johannes Hofmann <Johannes.Hofmann@gmx.de> - */ - -/* - * rectilinear (pinhole) projection model with distortion correction - */ - -/* switch to cartesian coordinates */ -c_x : cos(m_nick) * cos(m_view - c_view)$ -c_y : cos(m_nick) * sin(m_view - c_view)$ -c_z : sin(m_nick)$ - -/* nick rotation */ -c_x_rot : cos(c_nick) * c_x + sin(c_nick) * c_z$ -c_y_rot : c_y$ -c_z_rot : -sin(c_nick) * c_x + cos(c_nick) * c_z$ - -/* pinhole projection */ -x : c_y_rot / c_x_rot$ -y : -c_z_rot / c_x_rot$ - -/* rotation around view axis */ -x_rot : y * sin(c_tilt) + x * cos(c_tilt) + x0$ -y_rot : y * cos(c_tilt) - x * sin(c_tilt)$ - -/* distortion correction */ -d : sqrt(x_rot ^ 2 + y_rot ^ 2)$ -dist_fact : d ^ 2 * k1 + d * k0$ -x_dist : x_rot * (1 + dist_fact) * scale$ -y_dist : y_rot * (1 + dist_fact) * scale$ - -/* - * Some mangling for C code generation - */ - -load("expr2c.mac")$ - -x_expand : trigexpand(x_dist)$ -y_expand : trigexpand(y_dist)$ - -args: "double c_view, double c_nick, double c_tilt, double scale, double k0, double k1, double x0, double m_view, double m_nick"$ - -expr2c("ProjectionRectilinear::mac_x", args, x_expand)$ -expr2c("ProjectionRectilinear::mac_y", args, y_expand)$ - -expr2c("ProjectionRectilinear::mac_x_dc_view", args, diff(x_expand, c_view))$ -expr2c("ProjectionRectilinear::mac_x_dc_nick", args, diff(x_expand, c_nick))$ -expr2c("ProjectionRectilinear::mac_x_dc_tilt", args, diff(x_expand, c_tilt))$ -expr2c("ProjectionRectilinear::mac_x_dscale", args, diff(x_expand, scale))$ -expr2c("ProjectionRectilinear::mac_x_dk0", args, diff(x_expand, k0))$ -expr2c("ProjectionRectilinear::mac_x_dk1", args, diff(x_expand, k1))$ -expr2c("ProjectionRectilinear::mac_x_dx0", args, diff(x_expand, x0))$ -expr2c("ProjectionRectilinear::mac_y_dc_view", args, diff(y_expand, c_view))$ -expr2c("ProjectionRectilinear::mac_y_dc_nick", args, diff(y_expand, c_nick))$ -expr2c("ProjectionRectilinear::mac_y_dc_tilt", args, diff(y_expand, c_tilt))$ -expr2c("ProjectionRectilinear::mac_y_dscale", args, diff(y_expand, scale))$ -expr2c("ProjectionRectilinear::mac_y_dk0", args, diff(y_expand, k0))$ -expr2c("ProjectionRectilinear::mac_y_dk1", args, diff(y_expand, k1))$ -expr2c("ProjectionRectilinear::mac_y_dx0", args, diff(y_expand, x0))$ |