| 140 | } |
| 141 | |
| 142 | Imath::V2d StandardRadialLensModel::undistort( Imath::V2d p ) |
| 143 | { |
| 144 | Imath::V2d dn( UVtoDN( p ) ); |
| 145 | |
| 146 | const double dnx2( dn.x*dn.x ); |
| 147 | const double dny2( dn.y*dn.y ); |
| 148 | const double dnx4( dnx2*dnx2 ); |
| 149 | const double dny4( dny2*dny2 ); |
| 150 | |
| 151 | dn.x = dn.x * (1. + m_cxx*dnx2 + m_cxy*dny2 + m_cxxx*dnx4 + m_cxxy*dnx2*dny2 + m_cxyy*dny4); |
| 152 | dn.y = dn.y * (1. + m_cyx*dnx2 + m_cyy*dny2 + m_cyxx*dnx4 + m_cyyx*dnx2*dny2 + m_cyyy*dny4); |
| 153 | |
| 154 | return DNtoUV( dn ); |
| 155 | } |
| 156 | |
| 157 | Imath::V2d StandardRadialLensModel::distort( Imath::V2d p ) |
| 158 | { |