摄影测量学单像空间后方交会程序设计作业.doc

usingSystem;usingSystem、Collections、Generic;usingSystem、Linq;usingSystem、Text;namespace单像空间后方交会{classProgram{staticvoidMain(string[]args){intx0,y0,i,j;doublef,m;Console、Write("请输入像片比例尺:");m=double、Parse(Console、ReadLine());Console、Write("请输入像片得内方位元素x0:");//均以毫米为单位x0=int、Parse(Console、ReadLine());Console、Write("请输入像片得内方位元素y0:");y0=int、Parse(Console、ReadLine());Console、Write("请输入摄影机主距f:");f=double、Parse(Console、ReadLine());Console、WriteLine();//输入坐标数据double[,]zuobiao=newdouble[4,5];for(i=0;i<4;i++){for(j=0;j<5;j++){if(j<3){Console、Write("请输入第{0}个点得第{1}个地面坐标:",i+1,j+1);zuobiao[i,j]=double、Parse(Console、ReadLine());}else{Console、Write("请输入第{0}个点得第{1}个像点坐标:",i+1,j-2);zuobiao[i,j]=double、Parse(Console、ReadLine());}}Console、WriteLine();}//归算像点坐标for(i=0;i<4;i++){for(j=3;j<5;j++){if(j==3)zuobiao[i,j]=zuobiao[i,j]-x0;elsezuobiao[i,j]=zuobiao[i,j]-y0;}}//计算与确定初值doublezs0=m*f,xs0=0,ys0=0;for(i=0;i<4;i++){xs0=xs0+zuobiao[i,0];ys0=ys0+zuobiao[i,1];}xs0=xs0/4;ys0=ys0/4;//逐点计算误差方程系数double[,]xishu=newdouble[8,6];