第一篇：攝影測量實習報告-單片空間后方交會

攝影測量實習報告

實習內容:單片空間后方交會編程

實習者:李友兵

學號:0810050121

指導老師:張金平老師

實習時間:2011.05.30——2011.06.03

一、實習目的與任務

此次攝影測量實習主要是要自主編程實現單像空間后方交會，通過已知的內方位元素和控制點像點坐標和地面坐標求解六個外方位元素，在此過程中深入理解單像空間后方交會的原理和對編程的熟悉和理解（我用的是C語言編程），和對時間的合理運用，對知識的綜合運用，培養理論的實際運用能力，任務是在一個星期內自主完成。

二、單片空間后方交會理論基礎

單像空間后方交會：是通過以像點平面坐標為觀測值，以控制點坐標為已知值，利用共線條件方程和最小二乘原理，運用間接平差方法，通過迭代求解6個外方位元素。程序設計的思路與流程

三、程序設計的思路與流程（編程框架和步驟）

1、根據內方位元素和已知的數據將控制點的框標坐標轉換為像平面坐標系坐標：x=x′-x0，y=y′-y0

2、確定未知數的初始值：角元素初始值κ0=ω0=Φ0=0；線元素初始值Zs0=H=mf，Xs0=(X1+X2+X3+X4)/4,Ys0=（Y1+Y2+Y3+Y4）/4

3、利用角元素初始值計算方向余弦值組成旋轉矩陣R a1=cosΦ*cosκ-sinΦ*sinω*sinκ,a2=-cosΦ*sinκ-sinΦ*sinω*cosκ，a3=-sinΦ*cosω

b1=cosω*sinκ,b2=cosω*cosκ,b3=-sinω

c1=sinΦ*cosκ+cosΦ*sinω*sinκ,c2=-sinΦ*sinκ+cosΦ*sinω*cosκ,c3=cosΦ*cosω

4、逐點計算控制點的像點坐標的近似值，共線方程為: x=(-f*(a1*(X-Xs)+b1*(Y-Ys)+c1*(Z-Zs)))/(a3*(X-Xs)+b3*(Y-Ys)+c3*(Z-Zs))

y=(-f*(a2*(X-Xs)+b2*(Y-Ys)+c2*(Z-Zs)))/(a3*(X-Xs)+b3*(Y-Ys)+c3*(Z-Zs))

(x)=(-f*(a1*(X-Xs)+b1*(Y-Ys)+c1*(Z-Zs)))/(a3*(X-Xs)+b3*(Y-Ys)+c3*(Z-Zs))(y)=(-f*(a2*(X-Xs)+b2*(Y-Ys)+c2*(Z-Zs)))/(a3*(X-Xs)+b3*(Y-Ys)+c3*(Z-Zs))

5、組成誤差方程：Vx=a11*dXs+a12*dYs+a13*dZs+a14*dΦ+a15*dω+a16*dκ+(x)-x Vy=a21*dXs+a22*dYs+a23*dZs+a24*dΦ+a25*dω+a26*dκ+(y)-y 系數求解(是共線方程分別外方位元素求導，是共線方程線性化的系數)：變量代換

A= a1*(X-Xs)+b1*(Y-Ys)+c1*(Z-Zs)B= a2*(X-Xs)+b2*(Y-Ys)+c2*(Z-Zs)C= a3*(X-Xs)+b3*(Y-Ys)+c3*(Z-Zs)a11=(a1*f+a3*x)/C,a12=(b1*f+b3*x)/C,a13=(c1*f+c3*f)/C,a14=y*sinω-((x*(x*cosκ-y*sinκ))/f+f*cosκ)cosω,a15=-f*sinκ-(x*(x*sinκ+y*cosκ)/f),a16=y a21=(a2*f+a3*x)/C,a22=(b2*f+b3*x)/C,a23=（c2*f+c3*x）/C,a24=-xsinω-((x*(x*cosκ-y*sinκ))/f-f*sinκ)cosω,a25=-f*cosκ-(y*(x*sinκ+y*cosκ)/f),a26 =-x 誤差方程的常系數（是像點坐標觀測值與計算的近似值的差值）： lx=x-(x),ly=y-(y)

6、組成法方程，解求外方位元素改正數X=(ATA)-1ATL(A為誤差方程的系數矩陣，L為誤差方程的常系數矩陣，通過步驟5求得，此處先求ATA再求矩陣的逆矩陣，解得的改正數加上相應的近似值得到外方位元素新的近似值)

7、檢查計算是否收斂：將求得的外方位元素改正數與規定的限差比較，大于限差繼續迭代，小于限差則終止。

四、各子函數詳細設計的關鍵技術參數

子函數（輸入函數、Input，矩陣求積Matrixmultiply，計算函數Resection，矩陣轉置Matrixtranspose，矩陣求逆Matrixinverse，輸出函數Output）主要用了Matrixtranspose，矩陣的行變列，列變行,參數為要轉置的矩陣，轉置后的矩陣，要轉置矩陣的行列數，Matrixinverse,求矩陣的代數余子式，參數有要求逆的矩陣和，逆矩陣的行數，Matrixmultiply，一矩陣的行乘以二矩陣的列，參數為一矩陣，二矩陣，所求矩陣，一的行，一的列，二的列

五、像片外方位元素解算結

六、實習體會

從周一開始進行整個實習框架進行構建，編寫了實習的思路既步驟；然后是梳理框架，要用到的子函數，子函數參數，怎么編寫，有不懂的東西網上參考資料；再是自己動手編程，發現編程要求很細膩，出不得一點差錯，程序編寫出來后，有很多的錯誤，要求逐一進行改正，修改；最后才得以運行。通過本次實習，我深刻理解了單片空間后方交會原理，進一步熟悉理解了C語言編程，認識到了時間的搭配的重要性和參考資料的必要性，當遇到難題是要迎難而上，達到突破，當完成是能感到一絲絲欣慰和成就感。

附錄:源代碼

#include “stdio.h” #include “math.h” #include “Matrixmultiply.c” #include “Matrixtranspose.c” #include “Matrixinverse.c” void main(){ int i,j,k,f=0;double x0=0.00018, y0=0.00026,fk=0.15324;

//內方位元素

double m=40000;//估算比例尺

double R[3][3],XG[6][1],AT[6][8],ATA[6][6],ATL[6][1];double Xs=0.0, Ys=0.0, Zs=0.0,Q=0.0,W=0.0,K=0.0;

double X,Y,Z,L[8][1],A[8][6];

double B[4][5]={-0.08615,-0.06899,36589.41,25273.32,2195.17,-0.05340,0.08221,37631.08,31324.51,728.69,-0.01478,-0.07663,39100.97,24934.98,2386.80,0.01046,0.06443,40426.54,30319.81,757.31};

for(i=0;i<4;i++){

Xs=Xs+B[i][2];

Ys=Ys+B[i][3];

Zs=Zs+B[i][4];} Xs=Xs/4;Ys=Ys/4;Zs=m*fk;//求得外方位線元素的初始值

do//迭代計算

{

f++;//迭代次數

//組成旋轉矩陣

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);

//計算系數陣和常數項

for(i=0, k=0,j=0;i<=3;i++,k++,j++)

{

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]-(x0-fk*X/Z);

L[j+1][0]=B[i][1]-(y0-fk*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]*fk+R[2][2]*(B[i][1]-y0))/Z;

A[k+1][3]=-(B[i][0]-x0)*sin(W)-((B[i][1]-y0)*((B[i][0]-x0)*cos(K)-(B[i][1]-y0)*sin(K))/fk-fk*sin(K))*cos(W);

A[k+1][4]=-fk*cos(K)-(B[i][1]-y0)*((B[i][0]-x0)*sin(K)+(B[i][1]-y0)*cos(K))/fk;

A[k+1][5]=-(B[i][0]-x0);

k++;

}

Matrixtranspose(A,AT,8,6);//此為轉置函數的調用,求AT

Matrixmultiply(AT,A,ATA,6,8,6);//此為矩陣相乘函數的調用，求ATA

Matrixinverse(ATA,6);//此為求矩陣逆函數的調用，求ATA的逆

Matrixmultiply(AT,L,ATL,6,8,1);//此為矩陣相乘函數的調用，求ATL

Matrixmultiply(ATA,ATL,XG,6,6,1);//此為矩陣相乘函數的調用，求外方位元素改正數

Xs=Xs+XG[0][0];Ys=Ys+XG[1][0];Zs=Zs+XG[2][0];

Q=Q+XG[3][0];W=W+XG[4][0];K=K+XG[5][0];//初始值加外方位元素改正數進行迭代

}while(XG[3][0]>=0.00000001||XG[4][0]>=0.00000001||XG[5][0]>=0.00000001);//當限差滿足要求時要再一次進行旋轉矩陣的求解

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);printf(“迭代次數:%d”,f);

//屏幕輸出誤差方程系數陣、常數項、改正數

printf(“nn誤差方程系數矩陣A為:nn”);for(i=0;i<6;i++){

for(j=0;j<6;j++)

printf(“%13.5e ”,A[i][j]);

printf(“n”);}

printf(“n常數項L為:nn”);

for(i=0;i<8;i++){

for(j=0;j<1;j++)

printf(“%13.5e ”,L[i][j]);

printf(“n”);} printf(“n改正數XG為:nn”);

for(i=0;i<6;i++){

for(j=0;j<1;j++)

printf(“%13.5e ”,XG[i][j]);

printf(“n”);}

printf(“n相片的外方位元素為:nn”);

printf(“ Xs=%13.7e, Ys=%13.7e, Zs=%13.7e nn”,Xs,Ys,Zs);

printf(“ Q=%13.7e, W=%13.7e, K=%13.7e n”,Q,W,K);printf(“n旋轉矩陣R為:nn”);

for(i=0;i<3;i++){

for(j=0;j<3;j++)

printf(“%13.5e ”,R[i][j]);

printf(“n”);}

} //子函數

#include

#include

#include

int Matrixinverse(a,n)

int n;

double a[];

{ int *is,*js,i,j,k,l,u,v;

double d,p;

is=malloc(n*sizeof(int));

js=malloc(n*sizeof(int));

for(k=0;k<=n-1;k++)

{ d=0.0;

for(i=k;i<=n-1;i++)

for(j=k;j<=n-1;j++)

{ l=i*n+j;p=fabs(a[l]);

if(p>d){ d=p;is[k]=i;js[k]=j;}

}

if(d+1.0==1.0)

{ free(is);free(js);printf(“err**not invn”);

return(0);

}

if(is[k]!=k)

for(j=0;j<=n-1;j++)

{ u=k*n+j;v=is[k]*n+j;

p=a[u];a[u]=a[v];a[v]=p;

}

if(js[k]!=k)

for(i=0;i<=n-1;i++)

{ u=i*n+k;v=i*n+js[k];

p=a[u];a[u]=a[v];a[v]=p;

}

l=k*n+k;

a[l]=1.0/a[l];

for(j=0;j<=n-1;j++)

if(j!=k)

{ u=k*n+j;a[u]=a[u]*a[l];}

for(i=0;i<=n-1;i++)

if(i!=k)

for(j=0;j<=n-1;j++)

if(j!=k)

{ u=i*n+j;

a[u]=a[u]-a[i*n+k]*a[k*n+j];

}

for(i=0;i<=n-1;i++)

if(i!=k)

{ u=i*n+k;a[u]=-a[u]*a[l];}

}

for(k=n-1;k>=0;k--)

{ if(js[k]!=k)

for(j=0;j<=n-1;j++)

{ u=k*n+j;v=js[k]*n+j;

p=a[u];a[u]=a[v];a[v]=p;

}

if(is[k]!=k)

for(i=0;i<=n-1;i++)

{ u=i*n+k;v=i*n+is[k];

p=a[u];a[u]=a[v];a[v]=p;

}

}

free(is);free(js);

return(1);

} //子函數

void Matrixmultiply(a,b,c,m,n,k)int m,n,k;double a[],b[],c[];{ int i,j,l,u;for(i=0;i

for(j=0;j

{

u=i*k+j;c[u]=0.0;

for(l=0;l

c[u]+=a[i*n+l]*b[l*k+j];

}

return;} //子函數

void Matrixtranspose(a,b,m,n)int m,n;double a[],b[];{ int i,j,u;for(i=0;i

for(j=0;j

{

u=j*m+i;b[u]=0.0;b[j*m+i]=a[i*n+j];

} } return;}

第二篇：攝影測量空間后前交會VC++

#include #include using namespace std;const int N=4;const int n=6;/*---------------矩陣相乘---------------------*/

void mult(double *m1,double *m2,double *result,int i_1,int j_12,int j_2){

int i,j,k;for(i=0;i

for(j=0;j

result[i*j_2+j]=0.0;

for(k=0;k

result[i*j_2+j]+=m1[i*j_12+k]*m2[j+k*j_2];

} return;} /*---------------矩陣求逆---------------------*/ void inverse(double c[n][n]){

int i,j,h,k;double p;double q[n][12];for(i=0;i

for(j=0;j

q[i][j]=c[i][j];

for(i=0;i

for(j=n;j<12;j++){

if(i+6==j)

q[i][j]=1;

else q[i][j]=0;

}

for(h=k=0;k

for(i=k+1;i

if(q[i][h]==0)

continue;

p=q[k][h]/q[i][h];

for(j=0;j<12;j++){

q[i][j]*=p;

q[i][j]-=q[k][j];

} } for(h=k=n-1;k>0;k--,h--)// 消去對角線以上的數據

for(i=k-1;i>=0;i--){

if(q[i][h]==0)

continue;

p=q[k][h]/q[i][h];

for(j=0;j<12;j++){

q[i][j]*=p;

q[i][j]-=q[k][j];

} } for(i=0;i

p=1.0/q[i][i];

for(j=0;j<12;j++)

q[i][j]*=p;} for(i=0;i

for(j=0;j

c[i][j]=q[i][j+6];} /*---------------矩陣轉置---------------------*/

void transpose(double *m1,double *m2,int m,int n){ //矩陣轉置

int i,j;

for(i=0;i

for(j=0;j

m2[j*m+i]=m1[i*n+j];

return;

} void main(){ double Xs,Ys,Zs,q,w,k;double Xsr,Ysr,Zsr,qr,wr,kr;// double a[3],b[3],c[3];double ar[3],br[3],cr[3];// double x0,y0,f;double x[N],y[N];double X[N],Y[N],Z[N];double x1[N],y1[N];double xr[N],yr[N];// double x1r[N],y1r[N];// double m;double L[2*N];double Lr[2*N];// double XX[6];double XXr[6];// double A[2*N][6];double B[2*N][6];// double X0[N],Y0[N],Z0[N],At[6][2*N],result1[6][6],result2[6][1];double X0r[N],Y0r[N],Z0r[N],Bt[6][2*N],result1r[6][6],result2r[6][1];// double R1[3][3],R2[3][3],XI1[3],XII1[3],XI2[3],XII2[3],XI3[3],XII3[3],XI4[3],XII4[3],XI5[3],XII5[3];double N1,N2,N3,N4,N5,BX,BY,BZ,G[5][3];double RS1[3][1],RSR1[3][1],RS2[3][1],RSR2[3][1],RS3[3][1],RSR3[3][1],RS4[3][1],RSR4[3][1],RS5[3][1],RSR5[3][1];int i,j,n=0,nr=0;double sum=0;double m1,n1,m2,n2,m3,n3,m4,n4,m5,n5,p1,q1,p2,q2,p3,q3,p4,q4,p5,q5;

/*---------------輸入點地面坐標---------------------*/ X[0]=5083.205;X[1]=5780.02;X[2]=5210.879;X[3]=5909.264;Y[0]=5852.099;Y[1]=5906.365;Y[2]=4258.446;Y[3]=4314.283;Z[0]=527.925;Z[1]=571.549;Z[2]=461.81;Z[3]=455.48;m1=0.051758;n1=0.081555;p1=-0.039953;q1=0.078463;m2=0.014618;n2=-0.000231;p2=-0.076016;q2=0.000036;m3=0.04988;n3=-0.000792;p3=-0.042201;q3=-0.001022;m4=0.086243;n4=-0.001346;p4=-0.007706;q4=-0.002112;m5=0.048135;n5=-0.079962;p5=-0.044438;q5=-0.079736;/*---------------輸入點像片坐標---------------------*/ x[0]=16.012;x[1]=88.56;x[2]=13.362;x[3]=82.24;y[0]=79.963;y[1]=81.134;y[2]=-79.37;y[3]=-80.027;

xr[0]=-73.93;xr[1]=-5.252;xr[2]=-79.122;xr[3]=-9.887;yr[0]=78.706;yr[1]=78.184;yr[2]=-78.879;yr[3]=-80.089;

/*-----------------設定外方位元素初始值--------------*/ x0=0;y0=0;f=152.00;m=10000;Xs=0;Ys=0;Zs=f*m/1000;Xsr=0;Ysr=0;Zsr=f*m/1000;// q=0;w=0;k=0;qr=0;wr=0;kr=0;// XX[3]=1;XXr[3]=1;// /*------------------迭代計算左片Xs,Ys,Zs,q,w,k--------------------------*/ while((XX[3]>6/206265 || XX[4]>6/206265 || XX[5]>6/206265)&&n<100){// /*----------------旋轉矩陣R-----------------------*/ a[0]=cos(q)*cos(k)-sin(q)*sin(w)*sin(k);a[1]=-cos(q)*sin(k)-sin(q)*sin(w)*cos(k);a[2]=-sin(q)*cos(w);

b[0]=cos(w)*sin(k);b[1]=cos(w)*cos(k);b[2]=-sin(w);c[0]=sin(q)*cos(k)+cos(q)*sin(w)*sin(k);c[1]=-sin(q)*sin(k)+cos(q)*sin(w)*cos(k);c[2]=cos(q)*cos(w);/*-----------------像點坐標計算值------------------*/ for(i=0;i

A[2*i][0]=((a[0]*f+a[2]*(x[i]-x0)))/Z0[i];

A[2*i][1]=((b[0]*f+b[2]*(x[i]-x0)))/Z0[i];

A[2*i][2]=((c[0]*f+c[2]*(x[i]-x0)))/Z0[i];

A[2*i][3]=(y[i]-y0)*sin(w)-((x[i]-x0)*((x[i]-x0)*cos(k)-y[i]*sin(k))/f+f*cos(k))*cos(w);

A[2*i][4]=-f*sin(k)-(x[i]-x0)*((x[i]-x0)*sin(k)+(y[i]-y0)*cos(k))/f;

A[2*i][5]=y[i]-y0;

L[2*i]=x[i]-x1[i];

A[1+2*i][0]=((a[1]*f+a[2]*(y[i]-y0)))/Z0[i];

A[1+2*i][1]=((b[1]*f+b[2]*(y[i]-y0)))/Z0[i];

A[1+2*i][2]=((c[1]*f+c[2]*(y[i]-y0)))/Z0[i];

A[1+2*i][3]=-(x[i]-x0)*sin(w)-((y[i]-y0)*((x[i]-x0)*cos(k)-(y[i]-y0)*sin(k))/f-f*sin(k))*cos(w);

A[1+2*i][4]=-f*cos(k)-(y[i]-y0)*((x[i]-x0)*sin(k)+(y[i]-y0)*cos(k))/f;

A[1+2*i][5]=-x[i]+x0;

L[1+2*i]=y[i]-y1[i];}

/*-------------------解法方程--------------------*/ transpose(&A[0][0],&At[0][0],2*N,6);//求A轉置矩陣

mult(&At[0][0],&A[0][0],&result1[0][0],6,2*N,6);//A矩陣 X A轉置矩陣-> result1 inverse(result1);// result1 求逆

mult(&At[0][0],L,&result2[0][0],6,2*N,1);// A轉置矩陣 X L-> result2 mult(&result1[0][0],&result2[0][0],&XX[0],6,6,1);// result1 X result2-> XX Xs+=XX[0];

} Ys+=XX[1];Zs+=XX[2];q+=XX[3];w+=XX[4];k+=XX[5];n++;/*------------------迭代計算有片Xsr,Ysr,Zsr,qr,wr,kr--------------------------*/

while((XXr[3]>6/206265 || XXr[4]>6/206265 || XXr[5]>6/206265)&&nr<100){//右點算Xs,Ys,Zs,q,w,k /*----------------旋轉矩陣R-----------------------*/

ar[0]=cos(qr)*cos(kr)-sin(qr)*sin(wr)*sin(kr);

ar[1]=-cos(qr)*sin(kr)-sin(qr)*sin(wr)*cos(kr);

ar[2]=-sin(qr)*cos(wr);

br[0]=cos(wr)*sin(kr);

br[1]=cos(wr)*cos(kr);

br[2]=-sin(wr);

cr[0]=sin(qr)*cos(kr)+cos(qr)*sin(wr)*sin(kr);

cr[1]=-sin(qr)*sin(kr)+cos(qr)*sin(wr)*cos(kr);

cr[2]=cos(qr)*cos(wr);

/*-----------------像點坐標計算值------------------*/

for(i=0;i

X0r[i]=ar[0]*(X[i]-Xsr)+br[0]*(Y[i]-Ysr)+cr[0]*(Z[i]-Zsr);

Y0r[i]=ar[1]*(X[i]-Xsr)+br[1]*(Y[i]-Ysr)+cr[1]*(Z[i]-Zsr);

Z0r[i]=ar[2]*(X[i]-Xsr)+br[2]*(Y[i]-Ysr)+cr[2]*(Z[i]-Zsr);

x1r[i]=x0-f*X0r[i]/Z0r[i];

y1r[i]=y0-f*Y0r[i]/Z0r[i];

}

/*-------------誤差方程中各偏導數的值--------------*/

for(i=0;i

B[2*i][0]=((ar[0]*f+ar[2]*(xr[i]-x0)))/Z0r[i];

B[2*i][1]=((br[0]*f+br[2]*(xr[i]-x0)))/Z0r[i];

B[2*i][2]=((cr[0]*f+cr[2]*(xr[i]-x0)))/Z0r[i];

B[2*i][3]=(yr[i]-y0)*sin(wr)-((xr[i]-x0)*((xr[i]-x0)*cos(kr)-yr[i]*sin(kr))/f+f*cos(kr))*cos(wr);

B[2*i][4]=-f*sin(kr)-(xr[i]-x0)*((xr[i]-x0)*sin(kr)+(yr[i]-y0)*cos(kr))/f;

B[2*i][5]=yr[i]-y0;

Lr[2*i]=xr[i]-x1r[i];

B[1+2*i][0]=((ar[1]*f+ar[2]*(yr[i]-y0)))/Z0r[i];

B[1+2*i][1]=((br[1]*f+br[2]*(yr[i]-y0)))/Z0r[i];

B[1+2*i][2]=((cr[1]*f+cr[2]*(yr[i]-y0)))/Z0r[i];

B[1+2*i][3]=-(xr[i]-x0)*sin(wr)-((yr[i]-y0)*((xr[i]-x0)*cos(kr)-(yr[i]-y0)*sin(kr))/f-f*sin(kr))*cos(wr);

}

} B[1+2*i][4]=-f*cos(kr)-(yr[i]-y0)*((xr[i]-x0)*sin(kr)+(yr[i]-y0)*cos(kr))/f;B[1+2*i][5]=-xr[i]+x0;Lr[1+2*i]=yr[i]-y1r[i];/*-------------------解法方程--------------------*/ transpose(&B[0][0],&Bt[0][0],2*N,6);//求A轉置矩陣

mult(&Bt[0][0],&B[0][0],&result1r[0][0],6,2*N,6);//A矩陣 X A轉置矩陣-> result1 inverse(result1r);// result1 求逆

mult(&Bt[0][0],Lr,&result2r[0][0],6,2*N,1);// A轉置矩陣 X L-> result2 mult(&result1r[0][0],&result2r[0][0],&XXr[0],6,6,1);// result1 X result2-> XX Xsr+=XXr[0];Ysr+=XXr[1];Zsr+=XXr[2];qr+=XXr[3];wr+=XXr[4];kr+=XXr[5];nr++;BX=Xsr-Xs;//BX,BY,BZ BY=Ysr-Ys;BZ=Zsr-Zs;for(i=0;i<3;i++){//R1,R2賦值

R1[0][i]=a[i];

R1[1][i]=b[i];

R1[2][i]=c[i];

R2[0][i]=ar[i];

R2[1][i]=br[i];

R2[2][i]=cr[i];} RS1[0][0]=m1;RS1[1][0]=n1;RS1[2][0]=-f/1000;// 給 將要與R1，R2相乘的 RS矩陣和RSR矩陣 賦值

RS2[0][0]=m2;RS2[1][0]=n2;RS2[2][0]=-f/1000;RS3[0][0]=m3;RS3[1][0]=n3;RS3[2][0]=-f/1000;RS4[0][0]=m4;RS4[1][0]=n4;RS4[2][0]=-f/1000;RS5[0][0]=m5;RS5[1][0]=n5;RS5[2][0]=-f/1000;

RSR1[0][0]=p1;RSR1[1][0]=q1;RSR1[2][0]=-f/1000;RSR2[0][0]=p2;RSR2[1][0]=q2;RSR2[2][0]=-f/1000;RSR3[0][0]=p3;RSR3[1][0]=q3;RSR3[2][0]=-f/1000;RSR4[0][0]=p4;RSR4[1][0]=q4;RSR4[2][0]=-f/1000;RSR5[0][0]=p5;RSR5[1][0]=q5;RSR5[2][0]=-f/1000;

mult(&R1[0][0],&RS1[0][0],&XI1[0],3,3,1);//R1 X RS mult(&R1[0][0],&RS2[0][0],&XI2[0],3,3,1);mult(&R1[0][0],&RS3[0][0],&XI3[0],3,3,1);mult(&R1[0][0],&RS4[0][0],&XI4[0],3,3,1);mult(&R1[0][0],&RS5[0][0],&XI5[0],3,3,1);mult(&R2[0][0],&RSR1[0][0],&XII1[0],3,3,1);//R2 X RSR mult(&R2[0][0],&RSR2[0][0],&XII2[0],3,3,1);mult(&R2[0][0],&RSR3[0][0],&XII3[0],3,3,1);mult(&R2[0][0],&RSR4[0][0],&XII4[0],3,3,1);mult(&R2[0][0],&RSR5[0][0],&XII5[0],3,3,1);N1=(BX*XII1[2]-BZ*XII1[0])/(XI1[0]*XII1[2]-XI1[2]*XII1[0]);//給N1 賦值 有五個點 N2=(BX*XII2[2]-BZ*XII2[0])/(XI2[0]*XII2[2]-XI2[2]*XII2[0]);N3=(BX*XII3[2]-BZ*XII3[0])/(XI3[0]*XII3[2]-XI3[2]*XII3[0]);N4=(BX*XII4[2]-BZ*XII4[0])/(XI4[0]*XII4[2]-XI4[2]*XII4[0]);N5=(BX*XII5[2]-BZ*XII5[0])/(XI5[0]*XII5[2]-XI5[2]*XII5[0]);G[0][0]=Xs+N1*XI1[0];//五個點XA,YA,ZA計算G[i][j] i:第幾個點，j：是x，y還是z G[0][1]=Ys+N1*XI1[1];G[0][2]=Zs+N1*XI1[2];G[1][0]=Xs+N2*XI1[0];G[1][1]=Ys+N2*XI1[1];G[1][2]=Zs+N2*XI1[2];G[2][0]=Xs+N3*XI1[0];G[2][1]=Ys+N3*XI1[1];G[2][2]=Zs+N3*XI1[2];G[3][0]=Xs+N4*XI1[0];G[3][1]=Ys+N4*XI1[1];G[3][2]=Zs+N4*XI1[2];G[4][0]=Xs+N5*XI1[0];G[4][1]=Ys+N5*XI1[1];G[4][2]=Zs+N5*XI1[2];/*----------------旋轉矩陣R-----------------------*/ cout<<“--------左--------”<

”<

cout<<“Ys

”<

cout<<“Zs

”<

cout<<“迭代次數為:”<

”<

cout<<“Ysr

”<

cout<<“Zsr

”<

5”<

%-.5f

%.5f

%-.5f

%.5f

%-.5f n”,XI1[0],XI2[0],XI3[0],XI4[0],XI5[0]);printf(“ Y1

%-.5f

%.5f

%-.5f

%.5f

%-.5f n”,XI1[1],XI2[1],XI3[1],XI4[1],XI5[1]);printf(“ Z1

%-.5f

%.5f

%-.5f

%.5f

%-.5f n”,XI1[2],XI2[2],XI3[2],XI4[2],XI5[2]);printf(“ X2

%-.5f

%.5f

%-.5f

%.5f

%-.5f n”,XII1[0],XII2[0],XII3[0],XII4[0],XII5[0]);printf(“ Y2

%-.5f

%.5f

%-.5f

%.5f

%-.5f n”,XII1[1],XII2[1],XII3[1],XII4[1],XII5[1]);printf(“ Z2

%-.5f

%.5f

%-.5f

%.5f

%-.5f n”,XII1[2],XII2[2],XII3[2],XII4[2],XII5[2]);

printf(“ N1

%-.5f

%.5f

%-.5f

%.5f

%-.5f nn”,N1,N2,N3,N4,N5);printf(“ XA

%-.5f

%.5f

%-.5f

%.5f

%-.5f n”,G[0][0],G[1][0],G[2][0],G[3][0],G[4][0]);printf(“ YA

%-.5f

%.5f

%-.5f

%.5f

%-.5f n”,G[0][1],G[1][1],G[2][1],G[3][1],G[4][1]);printf(“ ZA

%-.5f

%.5f

%-.5f

%.5f

%-.5f n”,G[0][2],G[1][2],G[2][2],G[3][2],G[4][2]);

}

第三篇：攝影測量實習報告

攝影測量實習報告

一、引言：

數字攝影測量是基于測量的基本原理，通過對所獲取的數字/數字化影像進行處理自動（半自動）提取被攝對象用數字方式表達的幾何與物理信息，從而獲得各種形式的數字產品和目視化產品，攝影測量實習報告。數字攝影測量實習實在學完《數字攝影測量學》課程之后，進行數字攝影測量操作基本技能強化的一個重要實踐環節。

二、實驗目的和要求：

1、了解數字攝影測量生產流程

2、掌握立體像對定向建模型過程

3、掌握數字攝影測量測圖方法

三、實驗內容

使用JX4G全數字攝影測量軟件，按照相應的規范和規程，進行地理信息數據采集，完成全數字測圖實習。利用JX4G數字化成圖軟件測繪地形圖，具體內容包括：像對內定向、像對相對定向、像點坐標測量、匹配生成核線、數字地面模型（DEM）、編輯修改等高線、地形圖測圖。

四、實驗步驟

1、建模1—1 新建工程

啟動Geoway軟件，點擊“文件”→“新建工程”，在彈出的對話框中新建名為2196的工程，并指定新建工程的存儲路徑在D盤的名為0933的文件夾中。

a)點擊“工程管理”菜單，彈出其下拉菜單；

b)選擇“創建工程目錄”項，出現輸入對話框如圖22：

圖22 航片創建工程

c)在輸入對話框中，輸入新建的工程目錄名稱2196；然后點擊“瀏覽”，選擇所建立的目錄的保存路徑；

d)點擊“確定”，完成建立工程目錄，則在所建立的工程目錄下生成21*.ini文件——該文件記錄有關工程目錄的配置信息，實習報告《攝影測量實習報告》。

a)單擊“工程管理”菜單，彈出其下拉菜單；

b)選擇“選擇工程目錄”項，出現選擇路徑窗口如圖23：

圖23 航片選擇工程

c)在選擇路徑窗口中，選擇所需的目錄文件夾；

d)點擊“確定”，完成工程目錄選擇。

輸入相機信息時，點擊“工程管理”→“輸入文件”→“輸入相機文件”，彈出相機信息輸入窗口如圖24：

圖24 航片相機信息輸入

相機信息要依據相機自身的有關報告輸入。

在“焦距”欄輸入相機焦距，單位為毫米。

鼠標左鍵雙擊列表框內部，在激活的文本條內輸入框標坐標x和y值。

根據校正記錄讀取的變形值，設置透鏡變形參數，用與框標相同的方法輸入。

相機列表欄中，點擊“新建”按鈕創建新的相機文件；

點擊“確認”按鈕則在該相機文件中保存了當前對話框中的全部信息。

在“工程管理”菜單的“輸入文件”子菜單項中

點擊選擇“輸入文件”→“輸入控制點文件”，彈出控制點輸入窗口如圖26：

圖26 航片控制點輸入

在該窗口中輸入控制點坐標：點號XYZ。

設置航空影像立體模型建立所需的影像信息，點擊“輸入文件”→“設置影像路徑”，設置航空影像放置目錄。操作如下所述。

對于各單模型建立作業時，在創建像對之前需要設置影像放置目錄。在工程目錄下，要求用來建立像對的影像文件必須放在同一個目錄下；然后，通過如下步驟進行影像目錄設置：

a)選擇“輸入文件”→“設置影像路徑”，彈出設置路徑窗口，如圖28：

圖28 航片設置影像路徑

b)在該窗口中，只能通過單擊“瀏覽”，在彈出的選擇路徑窗口選擇影像放置路徑；

c)點擊“確定”，完成影像目錄設置，系統會記錄該信息到當前操作的工程目錄的ini文件中。

在“工程管理”菜單子項中，創建新像對。單擊選擇菜單“像對”→“建立新像對”，彈出輸入像對信息對話框如圖29：

圖29 輸入像對對話框

在像對信息對話框中，選擇像片、建立像對，并進行像對有關的設置。

第四篇：攝影測量實習報告

攝影測量實習報告

一、實習目的與要求

本次實習是在攝影測量的教學基礎上，理論實際相聯系的動手操作實習，是我們在學習測量專業的一個重要的實習環節。一方面是培養我們的實踐操作能力和運用軟件解算數據的能力，另一方面培養我們在今后遇到問題應該如何去解決的能力，通過實習發現自己在實踐動手方面的不足并想辦法解決，為以后的工作實踐打下扎實的基礎。使我們熟練地掌握攝影測量及遙感的原理，信息獲取的途徑，數字處理系統和應用處理方法。并進一步鞏固和深化理論知識，使理論與實踐相結合。切實加強我們大家的實踐動手能力，提高大家對這門新技術的認識和把握，全面培養我們的應用能了、創新能力和探索精神。

二、實習地點

桂林市雁山區大埠鄉

桂林理工大學博文管理學院機房

三、實習用具

小比例尺航片兩張、畫圖板一個、透明紙兩張、鉛筆、橡皮；電子計算機、ENVI遙感圖像處理系統、編程軟件（MATLAB、Visual Basic）

四、實習任務與要求

掌握航片調繪的方法步驟

掌握使用編程軟件設計解算移動曲面法數字高程模型內插子程序

掌握使用編程軟件設計解算空間后方交會

掌握使用ENVI遙感圖像處理系統處理遙感影像

五、實習步驟

航片調繪

本次實習的遙感圖像調繪主要判讀航片測區地物屬性，在透明紙上勾出邊界，必要時進行清繪。

在進行野外調繪之前，將調繪航片平放在畫圖板上，然后再將比調繪圖稍大一些的透明紙蓋于調繪航片上，用膠帶粘好，連同調繪航片用夾子固定于畫圖板。

第一天先將測區走過一遍，確定繪圖邊界，確定調繪路線，并對測區的情況有一個大致的了解。比較實際測區和航片的差別，并知道那些地方的地物是發生了變化的，以便于以后的判讀調繪工作進行。

接下來的兩天根據預定路線進行測區航片的地物判讀和屬性的標注。一般按照由遠及近、從總貌到碎部、邊走邊判，遠看近判的原則進行

調繪時的注意事項：

即地物地貌的調繪要連續進行，避免調繪不連貫和遺漏。

當地理名稱注記過密時，可適當取舍。

調繪工作應按照國家標準的地形圖圖式進行，說明性質的注記應采用“簡注表”，不得任意命名。

調繪要按照實地情況嚴格進行，不得偽造、篡改。

在調繪好的透明紙上，圖名注于調繪片正上方，調繪者姓名及調繪日期等信息在調繪航片的右下角。

編輯移動曲面法數字高程模型內插子程序

要求利用二次曲面擬合法：Z=Ax2+Bxy+Cy2+Dx+Ey+F，根據已知坐標（Xn,Yn,Zn）和待求點的平面坐標（Xn,Yn）,求出待求點的高程P。

解算步驟：

讀入已知點的坐標，建立以待定點為原點的局部坐標系；

建立誤差方程式：v=Xi2A+XiYiB+Yi2C+XiD+YiE+F-Zi

組成法方程，解算六個系數：X=(MTPM)?1MTPZ,其中pi=1/di2。

編輯空間后方交會程序

要求一直攝影機主距f四個控制點的像點坐標與相應的地面坐標，利用共線

方程的線性化形式，計算近似垂直攝影情況下像片的外方位元素。

解算步驟：

獲取已知數據：m,x0,y0,f,Xtp,Ytp,Ztp；

量測控制點像點坐標：x，y；

確定未知數初值；

組成誤差方程式：若P=I，X=(ATA)?1ATL；

解求外方位元素改正數、外方位元素的近似值；

檢查迭代是否收斂，是否需要重復計算。

使用ENVI系統處理遙感影像

主要要求學會使用ENVI系統對遙感影像進行監督分類和非監督分類

監督分類

制作分類模版：打開一幅遙感影像，在影像窗口打開Overlay-Region of interest，在Zoom窗口依次繪制可識別地物類別的區域；

監督分類：Classification-Supervised-Minimum Distance；

監督分類后處理：Classification-post Classification——Clump classes。

非監督分類

打開一幅遙感影像，單擊主菜單Classification-Unsupervised-Isodata，得到非監督分類的結果；

點擊Classification-post Classification-Combine classes合并相同或相似類別進行監督分類后處理

六、實習心得

攝影測量是一門專業的測繪學科，也是一門應用很廣的學科，隨著遙感技術的不斷發展，這門學科正從幾何學向信息科學發展。它的發展及運用對我們測繪來說是很有幫助的。而攝影測量實習則可以提高我們對攝影測量知識的理解，加強我們的實際運用能力。因此學校安排了三周的攝影測量實習，這對提升我們的攝影測量實際操作能力是很有幫助的。

此次實習分四個板塊，分別是全數字立體測圖；數字攝影測量的編程；遙感影像自動分類；像片的判讀與調繪等。全數字立體測圖是利用計算機代替解析測圖儀、用數字影像代替模擬像片、用數字光標代替光學光標，直接在計算機上進行數字化測圖的作業方法。這個實習要求我們學會使用ENVY軟件構建立體模型，制作測區的DEM、DOM和等高線圖，同時熟練使用交互式數字影像測圖系統在立體影像上量測不同類地物，并時行地物數據采集及編輯，生成數字測圖文件，按標準的制圖符號將之輸出為矢量地形圖。數字攝影測量的編程則要求我們學習使用Matlab進行攝影測量編程，掌握移動曲面法數字高程模型內插子程序的設計方法和空間后方交會程序的設計方法。遙感影像自動分類則是讓我們了解并掌握督與非督分類的過程和方法，并利用監督分類結果制作一幅影像地圖。像片的判讀與調繪則是讓我們利用學過的幾類常用遙感影像的判讀技術與方法，完成航空像片或彩紅外片的判讀和外業調繪工作，掌握全野外調繪的基本技能。

為使學生明確本次實習的總體任務及每一實習項目具體的作業程序、作業方法，指導教師在各項實習內容開展之前進行集中講解，做到任務明確、過程清晰；實習過程中，分組指導和定期集中討論相

結合，啟發學生解決作業中出現的實際問題。本次實習不僅使學生理論知識得到鞏固、操作能力得到加強，同時也使學生運用知識的能力得到提高。

在航片調繪實習過程中不免出現些錯誤和困難，但是我們都沒有因此而放棄。我個人覺得在實習過程中細心是非常必要的，例如在選擇界點時，一不小心就有可能將航片像點中的界點找錯。還有在航片調繪清繪時，如果我們不細心，在沒有記住航片中現在有所改變的地方，我們的成果就會因我們的粗心大意而失去意義。所以我認為養成一個細心嚴謹的態度是非常必要的，這將減少一些不必要的錯誤和損失。其實，我覺得本次實習沒有什么特別困難的地方，只要大家能夠做到認真細心，我們的實習就會很順利。

在全數字立體測圖實習過程中讓我深深體會到，理論指導實踐這一真理。在本次全數字立體測圖實習過程中，我發現我們要做的工作其實很簡單，只要點擊有關的按鈕，計算機就自動幫我們完成要做的工作了。但是，如果我們沒有扎實的基礎知識，就無法正常并順利地操作計算機去完成我們要的指令操作。當我們再次遇到類似的問題時就無法解決了。對于我們來說，如果只有理論知識，而實踐操作為零，那也只是紙上談兵；但是指揮操作，不懂理論知識，那也不能獨立完成工作。所以，我們要好好學習理論知識，這樣才能指導實踐，加強我們的動手能力。將來畢業了，才能是個合格的測繪工作人員。

攝影測量實習只有三周，但是在這三周中我們學到了很多東西。比如剛開始的像片的判讀與調繪，雖然我們剛開始時要天天去野外調

繪踏勘，可能有點辛苦，但也讓我們學會了在野外如何建立判讀標志，根據判讀任務擬訂判讀編輯指示或技術補充指示；這室內判讀收集調查判讀樣片資料，對那些難于識別和需要判讀而不常見的地物，拍攝地面立體照片，建立作業區的判讀標志庫，對復雜地物提出謄建議預測，根據預測做好人員分工和組織工作等。這些東西都是我們平時在課堂上所不能學到和理解到的。在地理要素的調繪中，讓我們學會了認真與仔細，因為調繪片上表示的各項地理要素一定要齊全，綜合取舍要合理；在圖面上的各種數字注記要齊全；位置要恰當。所有的這些都有助于我們養成做事認真負責的態度。

在這次實習中，我們學到了很多課堂上學不到的東西，也讓我們了解了攝影測量的應用及以后的發展狀況，讓我們知道了攝影測量的重要性及對我們將來工作的重要性。同時這次實習也讓我們對以后的測繪方式有了全新的了解與認識，增長了我們的知識及見解，也讓我們對測繪工作有了更深刻的理解。

第五篇：攝影測量實習報告

學 院 班

級 學生姓名

攝影測量實習報告

測繪與城市空間信息學院

專 業

學 號

指導教師

地理信息科學

一、實習目的.........................................................................................4

二、實習要求.........................................................................................4

三、實習過程.........................................................................................5

四、心得體會.....................................................................................12

一、實習目的

1.與實際工作結合，體驗實際工作。

在學長指導下進行外業數據采集，在實習過程中了解實際工作的要求，與已經參與工作的學長交流從而習得實際工作經驗，在實際的工作中結合課本理論，將理論與實踐相結合，深入理解攝影測量在實際工作中怎樣應用，將課本知識融會貫通。

2.對4D產品概念的認識與了解

4D產品的簡述及其應用

數字高程模型（Digital Elevation Model，縮寫DEM）是在某一投影平面（如高斯投影平面）上規則格網點的平面坐標（X，Y）及高程（Z）的數據集。DEM的格網間隔應與其高程精度相適配，并形成有規則的格網系列。根據不同的高程精度，可分為不同類型。為完整反映地表形態，還可增加離散高程點數據。

數字正射影像圖（Digital Orthophoto Map，縮寫DOM）是利用數字高程模型（DEM）對經掃描處理的數字化航空像片，經逐像元進行投影差改正、鑲嵌，按國家基本比例尺地形圖圖幅范圍剪裁生成的數字正射影像數據集。它是同時具有地圖幾何精度和影像特征的圖像，具有精度高、信息豐富、直觀真實等優點。

數字線劃地圖（Digital Elevation Model，縮寫DLG）是現有地形圖要素的矢量數據集，保存各要素間的空間關系和相關的屬性信息，全面地描述地表目標。

數字柵格地圖（Digital Raster Graphic，縮寫DRG）是現有紙質地形圖經計算機處理后得到的柵格數據文件。每一幅地形圖在掃描數字化后，經幾何糾正，并進行內容更新和數據壓縮處理，彩色地形圖還應經色彩校正，使每幅圖像的色彩基本一致。數字柵格地圖在內容上、幾何精度和色彩上與國家基本比例尺地形圖保持一致。

3.學習攝影測量相關軟件（JX4，PHOTOMOD等）

熟悉4D產品的作業流程。結合平常上機練習積累經驗知識在JX-4G平臺下完成對城建校區及其周圍區域DLG的制作，在PHOTOMOD軟件中學習空三加密過程并完成DEM、DOM的制作。

二、實習要求

1.航帶設計要求

像素14500*10000 焦距50mm 像元尺寸6.5微米 高差達到像對行高1/6 如有漏曝光需要航帶補飛

2.外業要求：

外業利用GPS做像控，在學長指導下了解外業工作，實習過程中要求同學們注意安全，聽從學長安排，利用GPS接收機進行RTK動態測量采集像控點數據，之后將數據挑揀匯總，為內業處理做準備。要求動態測量坐標誤差不超過2cm 3.內業要求：

利用PHOTOMOD做空三加密，并制作DEM，制作DOM，挑選像對利用JX4制作DLG，主要掌握4D產品的制作過程，4.本次實習分為兩個攝影測量測區：

（1）平頂山新城區測區，（2）平頂山市區測區

1）測區與資料分析

測區分析：兩個測區內部有樓房、道路、花壇、林地、操場、線桿等地物，也都有地形起伏，河流、湖泊等地物地貌，具有較強的地形地貌代表性。

資料分析：攝影主距35.6830mm，掃描影像像素大小為0.05mm，攝影比例尺為1：1000，航攝面積較小，航片的清晰度比較好，適合實習操作。2）成果要求：

DEM透視圖圖片一幅：比例尺為1：1000，要求精度：等高線間隔 2m，DEM拼接中誤差差限為2，大于3倍中誤差的點所占的比率小于1％。DOM圖片一幅，比例尺為1：1000，檢查每個模型的接邊處，保證影象無或很少變形及扭曲等錯誤。DLG圖片，比例尺為1：1000。

上交成果過程截圖保存至實習報告，最終每人上交實習報告一份。

三、實習過程

1.外業像控點采集

A.早晨到達指定實習地點，按照兩人一個小組，分組領取各組儀器：

每組一臺GPS接收機，一個手持桿，一個手溥，4塊接收機電池，外業記錄紙，水筆 B.由組長架設基準站，設定參數： 波特率9600，頻道1，接口7 C.小組長負責連接手溥，并根據已知點檢驗精度 D.由學長帶隊進入測區開始測量 外業具體測量過程如下：

a)連接接收機和手溥，準備好紙，筆 b)打開接收機，設置移動站，打開手溥，按照機器號藍牙連接各自接收機 c)打開手溥軟件，新建工程，以日期命名，d)移動到需要觀測的地點，點擊點測量，保存點名（字母+數字格式），同時在外業觀測記錄紙上記錄點名。

E.當晚由組長匯總當天外業數據以備內業處理

2.JX4操作過程（PHOTOMOD軟件操作）

1.安裝、打開可以進行空三加密的軟件

2.點擊創建新工程，設置完成后，添加航帶，如圖所示

3.添加像片到列表，并完成格式轉變

4.設置相機參數，焦距：50.2 像元大小X=6.8 Y= 6.8 打開外方位元素列表

5.導入EO

6.通過外方位元素劃分航帶，顯示圖像

7.設置相機方位角，進行同步

8.UAV自動空

三、自動匹配連接點、測區平差、設置參數、進行計算，結果如下

9.自動計算點，配置格網屬性、提取點云

10.構建TIN、通過TIN創建DEM、保存濾波后的DEM

DLG圖

四、心得體會

本次攝影測量實習，不僅使學生理論知識得到鞏固、操作能力得到加強，同時也使學生運用知識的能力得到提高。

在實習過程中不免出現些錯誤和困難，但是我們都沒有因此而放棄。我個人覺得在實習過程中細心是非常必要的，例如在選擇同名點時，一不小心就有可能將同名像點找錯。還有在影像匹配后編輯時，如果我們不細心，在沒有保存我們成果的情況下就關閉了窗口，我們的成果就會因我們的粗心大意而失去。所以我認為養成一個細心嚴謹的態度是非常必要的，這將減少一些不必要的錯誤和損失。其實，我覺得本次實習沒有什么特別困難的地方，只要大家能夠做到認真細心，我們的實習就會很順利。

本次實習讓我深深體會到，理論指導實踐這一真理。在本次實習中，我發現我們要做的工作其實很簡單，只要點擊有關的按鈕，計算機就自動幫我們完成要做的工作了。但是，如果我們沒有扎實的基礎知識，就無法正常并順利地操作計算機去完成我們要的指令操作。當我們再次遇到類似的問題時就無法解決了。對于我們來說，如果只有理論知識，而實踐操作為零，那也只是紙上談兵；但是指揮操作，不懂理論知識，那也不能獨立完成工作。所以，我們要好好學習理論知識，這樣才能指導實踐，加強我們的動手能力。將來畢業了，才能是個合格的測繪工作人員。

時間飛逝，兩周的實習就這樣結束了。雖然時間不長，我從中學到了許多在課堂中無法學會的知識，理解了許多在課堂上無法理解的知識。很珍惜這些天的實習，讓我更加清晰的認識了攝影測量學這門學科，進一步了解了相關軟件的操作和使用，鍛煉了我的動手能力。最后，很感謝李軍杰老師給我們提供的實習機會，每天不辭辛苦的陪伴著我們，給我們技術上的指導，生活上的管理。雖然，本次實習在學校機房，但依然當不住寒冬的侵襲，老師這樣陪伴著我們，讓我很感動。我知道只有優秀的實習報告與成果才能回報老師的辛欣工作，與默默付出。我相信，我的實習總結會讓老師得到安慰，覺得一些付出都是值得的了。在今后的工作和學習中，這次實習會給我源源不絕的動力和力量，我相信我會更加自信的面對今后的生活和工作，更加努力的學習和工作。

通過這次實習我們的收獲很大，在很多情況下我們都得到了很多意外的收獲，獲益匪淺！不僅對書本上的理論知識有了通盤的理解，更重要的是從實踐中檢驗了它的真理，了解了它的適應范圍之廣和作用之大，為我們以后從事工作而需要它打下了堅實的理論基礎與實踐經驗。