攝影測量學作業(yè)一(編輯修改稿)
2024-11-16 22:27 本頁面
【文章內(nèi)容簡介】 sing namespace std。const int n=6。void inverse(double c[n][n])。templatevoid transpose(T1*mat1,T2*mat2,int a,int b)。templatevoid multi(T1*mat1,T2 * mat2,T2 * result,int a,int b,int c)。templatevoid input(T*mat,int a,int b)。templatevoid output(T*mat,char*s,int a,int b)。int main(){ ofstream outFile。(5)。double x0=, y0=。double fk=。//內(nèi)方位元素 double m=39689。//估算比例尺double B[4][5]={},R[3][3],XG[6][1],AT[6][8],ATA[6][6],ATL[6][1]。input(B,4,5)。//從文件中讀取控制點的影像坐標和地面坐標，存入數(shù)組B double Xs=, Ys=, Zs=,Q=,W=,K=。double X,Y,Z,L[8][1],A[8][6]。//確定未知數(shù)的出始值for(int i=0。iYs=Ys+B[i][3]。Zs=Zs+B[i][4]。} Xs=Xs/4。Ys=Ys/4。Zs=Zs/4+m*fk。int f=0。do//迭代計算{f++。//組成旋轉(zhuǎn)矩陣R[0][0]=cos(Q)*cos(K)sin(Q)*sin(W)*sin(K)。R[0][1]=cos(Q)*sin(K)sin(Q)*sin(W)*cos(K)。R[0][2]=sin(Q)*cos(W)。R[1][0]=cos(W)*sin(K)。R[1][1]=cos(W)*cos(K)。R[1][2]=sin(W)。R[2][0]=sin(Q)*cos(K)+cos(Q)*sin(W)*sin(K)。R[2][1]=sin(Q)*sin(K)+cos(Q)*sin(W)*cos(K)。R[2][2]=cos(Q)*cos(W)。//計算系數(shù)陣和常數(shù)項for(int i=0,k=0,j=0。i{X=R[0][0]*(B[i][2]Xs)+R[1][0]*(B[i][3]Ys)+R[2][0]*(B[i][4]Zs)。Y=R[0][1]*(B[i][2]Xs)+R[1][1]*(B[i][3]Ys)+R[2][1]*(B[i][4]Zs)。Z=R[0][2]*(B[i][2]Xs)+R[1][2]*(B[i][3]Ys)+R[2][2]*(B[i][4]Zs)。L[j][0]=B[i][0](x0fk*X/Z)。L[j+1][0]=B[i][1](y0fk*Y/Z)。j++。A[k][0]=(R[0][0]*fk+R[0][2]*(B[i][0]x0))/Z。A[k][1]=(R[1][0]*fk+R[1][2]*(B[i][0]x0))/Z。A[k][2]=(R[2][0]*fk+R[2][2]*(B[i][0]x0))/Z。A[k][3]=(B[i][1]y0)*sin(W)((B[i][0]x0)*((B[i][0]x0)*cos(K)(B[i][1]y0)*sin(K))/fk+fk*cos(K))*cos(W)。A[k][4]=fk*sin(K)(B[i][0]x0)*((B[i][0]x0)*sin(K)+(B[i][1]y0)*cos(K))/fk。A[k][5]=B[i][1]y0。A[k+1][0]=(R[0][1]*fk+R[0][2]*(B[i][1]y0))/Z。A[k+1][1]=(R[1][1]*fk+R[1][2]*(B[i][1]y0))/Z。A[k+1][2]=(R[2][1]
公司管理相關推薦
鄂ICP備17016276號-1