qVP - `MetricPics3D`: Copyright 2009 Jeffrey M. Setterholm -0.231167656030E+03 0.231366783505E+03 -0.186457185269E-01 0.75806847368E+06 -0.414215353840E+04 0.414214095644E+04 -0.140203072523E+00 0.13593333398E+08 0.411151538960E+04 -0.411152917160E+04 -0.141405217852E+00 -0.13493482416E+08 0.917132401014E-03 -0.917132401014E-03 0.000000000000E+00 -0.30097028572E+01 EcefXyz = ^ hEcefQpix(4,4) * QvpPixels 0.495243592757E+01 0.446393933738E+01 -0.507901863801E+01 0.44179693549E+08 0.496717882878E+01 0.172816310036E+01 -0.236844438905E+01 0.19674873031E+08 -0.453975226394E+00 -0.576299095244E+01 -0.129801622632E+01 -0.20323482240E+08 -0.449253399075E-05 0.833660048155E-03 -0.825980366500E-03 0.74669048553E+04 QvpPixels= ^inv(hEcefQpix(4,4)) * EcefXyz -0.382811111982E-02 0.710365977862E+00 -0.703822081965E+00 0.63625877009E+07 0.995637371276E+00 -0.629097677323E-01 -0.689099850761E-01 0.27586170134E+06 -0.932285926343E-01 -0.701015362614E+00 -0.707026089261E+00 -0.19971002107E+05 0.000000000000E+00 0.000000000000E+00 0.000000000000E+00 0.10000000000E+01 SocUvw = ^ hSocEcef(4,4) * EcefXyz Internal qVP geometry (+X:in,+Y:right,+Z:down): 8744 4492 -pixels in image 0.874400000000E+03 0.898400000000E+03 0.200000000000E+00 -extent/PixSize -0.852105098998E+03 -0.138612287979E+03 -0.138228702083E+04 -left eye loc. -0.852105098998E+03 0.794587340224E+02 -0.138228702083E+04 -right eye loc. 0.000000000000E+00 0.000000000000E+00 -Shift L-eye subimage (pixels) 0.000000000000E+00 0.000000000000E+00 -Shift R-eye subimage (pixels) (`Shift`:provides minor adjustment by editing these values in the .jpg file.) 2010.05.11.1819.09 -qVP creation date & time P0102011 4 24076 2006.04.23.1904:06.73 -left eye: Run,Camera,Frame,Date P0102011 4 24081 2006.04.23.1904:09.73 -right eye: Run,Camera,Frame,Date Source images-with-geodesy: courtesy of GeoSpan Corp. All distances are in meters. Pixel:(0.,0.) is the image upper-left corner ECEF: native-GPS WGS-84 Earth-Centered,Earth-Fixed coords. Numerical Example (crossed-eye): Select the midpoint of the stereo pixel field: QvpPixels:(iLeft,iRight,jDown,1)= 6558.00 2186.00 2246.00 1.00 EcefXyz = hEcefQpix(4,4) * QvpPixels ...remember to divide homog. outputs by the 4th component. EcefXyz:(X,Y,Z,1)= -252203.104 -4516504.272 4481715.144 1.000 meters WGS-84 LLH:(Lat,Lon,Ht,1.)= 44.92630381557 -93.19609622338 228.855 1.000 Elevation - with respect to the Ellipsoid, not the Geoid. & Reversing... (it works!...) EcefXyz:(X,Y,Z,1)= -252203.104 -4516504.272 4481715.144 1.000 meters QvpPixels= inv(hEcefQpix(4,4)) * EcefXyz QvpPixels:(iLeft,iRight,jDown,1)= 6558.00 2186.00 2246.00 1.00 Note: QvpPixels=(0.,0.,0.,1.) is the image upper-left-corner Second example: Move the left-eye pixel 9.99 right QvpPixels:(iLeft,iRight,jDown,1)= 6567.99 2186.00 2246.00 1.00 EcefXyz = hEcefQpix(4,4) * QvpPixels EcefXyz:(X,Y,Z,1)= -252201.758 -4516503.492 4481726.916 1.000 meters WGS-84 LLH:(Lat,Lon,Ht,1.)= 44.92638424583 -93.19607975097 236.564 1.000 Final example: If the center of the user-defined Gameing Area is: WGS-84 LLH:(Lat,Lon,Ht,1.)= 44.92280000000 -93.19670000000 230.000 1.000 EcefXyz:(X,Y,Z,1)= -252266.078 -4516776.976 4481440.247 1.000 meters QvpPixels:(iLeft,iRight,jDown,1)= 6425.00 2053.00 4203.00 1.00 Corresponding points are meant to be in the same pixel row. Expressing the same point w.r.t. the Stereo-optic-center (Soc) SocUvw = hSocEcef(4,4) * EcefXyz SocUvw:(U,V,W,1)= 852.105 29.577 1382.287 1.000 meters !--- Example completed: @ 2010.05.11.1819.47 L