攝影測量作業(yè)3(編輯修改稿)
2025-04-21 02:40 本頁面
【文章內(nèi)容簡介】 double YS = XX(0, 1)。 double ZS = XX(0, 2)。 double w = XX(0, 4)。 double k = XX(0, 5)。 double q = XX(0, 3)。 double a1 = cos(q)*cos(k) sin(q)*sin(w)*sin(k)。 double a2 = cos(q)*sin(k) sin(q)*sin(w)*cos(k)。 double a3 = sin(q)*cos(w)。 double b1 = cos(w)*sin(k)。 double b2 = cos(w)*cos(k)。 double b3 = sin(w)。 double c1 = sin(q)*cos(k) + cos(q)*sin(w)*sin(k)。 double c2 = sin(q)*sin(k) + cos(q)*sin(w)*cos(k)。 double c3 = cos(q)*cos(w)。 f = f / 。 for (int i = 0。 i iRow。 i++) { double ZA = xyXYZ(i, 4)。 double XA = xyXYZ(i, 2)。 double YA = xyXYZ(i, 3)。 double x = (xyXYZ(i, 0)_wtof(strx0))/。 double y = (xyXYZ(i, 1)_wtof(stry0))/。 double _Z = a3*(XA XS) + b3*(YA YS) + c3*(ZA ZS)。 A(2*i, 0) = (a1*f + a3*x) / _Z。 A(2*i, 1) = (b1*f + b3*x) / _Z。 A(2*i, 2) = (c1*f + c3*x) / _Z。 A(2*i, 3) = y*sin(w) (x*(x*cos(k) y*sin(k)) / f + f*cos(k))*cos(w)。 A(2*i, 4) = f*sin(k) x*(x*sin(k) + y*cos(k)) / f。 A(2*i, 5) = y。 A(2*i+1, 0) = (a2*f + a3*y) / _Z。 A(2 * i + 1, 1) = (b2*f + b3*y) / _Z。 A(2 * i + 1, 2) = (c2*f + c3*y) / _Z。 A(2 * i + 1, 3) = x*sin(w) (y*(x*cos(k) y*sin(k)) / f f*sin(k))*cos(w)。 A(2 * i + 1, 4) = f*cos(k) y/ f*(x*sin(k) + y*cos(k))。 A(2 * i + 1, 5) = x。} return A。}CMatrix CKongJianHouFangJiaoHuiDlg::GetL(CMatrix xyXYZ, double f, CMatrix XX)//計(jì)算L矩陣{ int iRow = ()。 CMatr
醫(yī)療健康相關(guān)推薦
鄂ICP備17016276號-1