平差作业2.docx

课程设计报告张弛测绘一班课后作业2——航空摄影航带设计已知条件：需要测绘某地区的1:2000地形图，该地区的地面四角平面坐标分别为：左下角（0,0），右下角（0，7360），左上角（12880,0），右上角（12880，7360）。地面平均高程为750m。现采用1:8000比例尺的航空摄影，航摄相机主距f=150mm，像幅大小为23cm×23cm。 在上述条件下，飞机航摄飞行速度为200km/h，要求满足60%的影像航向重叠度和30%的影像旁向重叠度，第1航带中飞机飞行方向为自西向东，第1摄站摄影中心的平面坐标位置为（-368.000，184.000），第2航带中飞机飞行方向则为自东向西，以此类推。假设飞机为水平飞行、垂直摄影，并且不考虑地面起伏影响。（上述条件中若无特别注明，平面坐标的单位均为m。）请计算：1）飞机在每个航带航摄飞行过程中，相邻像片曝光拍摄的时间间隔应为多少？2）航带内相邻两像片航摄基线长度为多少？3）完成对上述整个地区进行航空摄影的航带设计，即：设计并计算对上述整个地区进行航摄时各个摄站摄影中心的三维坐标位置。要求：1）分析并描述摄影瞬间的解析几何关系；2）采用一种高级语言编程实现对上述所有问题的解算；3）编写并提交课程设计报告（2周时间内）。1. 设计任务2. 测区概况与摄影机参数3. 摄影解析几何关系描述4. 算法描述5. 程序流程6. 源程序7. 计算结果8. 结果分析9. 心得体会与建议解：设m为摄影比例尺，分别为像幅长度、宽度，p、q分别为横向旁向重叠度摄影基线=m*(1-p)=8000*23*0.01*(1-60%)=736米时间间隔T===13.25s航线间隔=m*(1-q)=8000*23*0.01*(1-30%)=1288米飞行高度H=f*m=0.15*8000=1200米C++程序编码：\*N代表摄区在（y方向）上的长度，M代表摄区在南北方向（x方向）上长度，A为第一摄站的x坐标，B为第一摄站的y坐标，S为摄影基线长，H为航线间隔。*\#includeiostreamusing namespace std;void main(){const int M=12880,N=7360,A=-368,B=184,H=1288,S=736; /*N代表摄区在（y方向）上的长度，M代表摄区在南北方向（x方向）上长度，A为第一摄站的x坐标，B为第一摄站的y坐标，S为摄影基线长，H为航线间隔。默认左下角坐标为（0,0）}*/int a,b,c,s,t,i,j,n;a=A;b=B;c=1;for(i=0;aM+H;i++){if(c==1) for(j=0;bN+S;j++){cout(a,b,1200);b+=736*c;}else for(j=0;b-S;j++){cout(a,b,1200);b+=736*c;}a+=1288;c=-c; coutendlendl;}}摄影瞬间解析几何关系：飞机上航摄像机与地面构成两个以飞行方向上某点为顶点的两个等腰三角形。航片上任一线段l与地面上相应线段L的比，以及摄影机主距f与摄影机物镜中心与摄区平均高程面之间的距离H之比均等于摄影比例尺1/m。同时摄影测量中还存在像片倾角、航线弯曲、像片旋偏等。 #include stdio.h #include malloc.h #include iostream#define MaxSize 50 typedef char ElemType;using namespace std;int main(){int i,x,y,X,Y,x1,y1,n;printf(请输入要航摄地区的地面右上角的坐标\n);printf(X=);scanf(%d,X);printf(Y=);scanf(%d,Y);printf(请输入第1摄站摄影中心的平面坐标\n);printf(x=);scanf(%d,x);printf(y=);scanf(%d,y);for(n=0;x+1288*n+920X;n++){if((n+1)%2==1){for(i=0;y+736*i+920Y;i++){x1=x+1288*n;y1=y+736*i;printf(第%d个航带第%d个摄站摄影中心的坐标为(%d,%d)\n,n+1,i+1,x1,y1);}x1=x+1288*n;y1=y+736*i;printf(第%d个航带第%d个摄站摄影中心的坐标为(%d,%d)\n,n+1,i+1,x1,y1);}else{for(i=9;i0;i--){x1=x+1288*n;y1=y+736*i;printf(第%d个航带第%d个摄站摄影中心的坐标为(%d,%d)\n,n+1,i+1,x1,y1);}x1=x+1288*n;y1=y+736*i;p