坐标转换源代码.docVIP

#include iostream.h #include math.h #include stdio.h #define pi 3.1415926535897932 //圆周率 void xyz_xyz(double xyz[],double xyz1[],double canshu[]); void gstyz(double blh[],double xy[],double para2[]); void blh_xyz(double blh[],double xyz[],double para1[]); double ziwuhu(double B,double a,double e2); double xyz_blh(double blh[],double xyz[],double para2[]); double huahu(double b); main() { double blh[3],xyz[3],xyz1[3],xy[2],para1[2],para2[2],canshu[7]; int i; cout请输入大地坐标B,L,H，角度输入方式如下：endl; cout若输入的角度为30度30分30秒，对应的代码为：30.3030endl; cinblh[0]blh[1]blh[2]; coutendl; //将输入的角度形式转化为弧度 blh[0]=huahu(blh[0]); blh[1]=huahu(blh[1]); cout请输入WGS84椭球的椭球参数，长半轴以及扁率的倒数：endl; cinpara1[0]para1[1]; coutendl; cout请输入北京54椭球的椭球参数，长半轴以及扁率的倒数：endl; cinpara2[0]para2[1]; coutendl; //调用函数，实现大地坐标向空间直角坐标系中的转换 blh_xyz(blh,xyz,para1); cout对应的空间直角坐标系中的坐标为：endl; printf(x=%10.6f,y=%10.6f,z=%10.6f,xyz[0],xyz[1],xyz[2]); coutendl; cout请分别输入两坐标系之间的缩放参数，旋转参数和位移参数：endl; for(i=0;i7;i++) { cincanshu[i]; } coutendl; //将输入的角度转换为弧度的形式 canshu[1]=huahu(canshu[1]); canshu[2]=huahu(canshu[2]); canshu[3]=huahu(canshu[3]); //调用函数，实现两个空间直角坐标系之间的转换 xyz_xyz(xyz,xyz1,canshu); //调用函数，将大地坐标转换为直角坐标 xyz_blh(blh,xyz,para2); //调用函数，计算对应的高斯投影面上的坐标 gstyz(blh,xy,para2); //输出结果 printf(对应的高斯投影坐标系中的坐标为：\n); printf(x=%10.6f,y=%10.6f,xy[0],xy[1]); } //大地坐标转换为直角坐标 void blh_xyz(double blh[],double xyz[],double para1[]) { double b,e2,N; b=para1[0]*(1-1/para1[1]); e2=1-b*b/(para1[0]*para1[0]); N=para1[0]/sqrt(1-e2*sin(blh[0])*sin(blh[0])); xyz[0]=(N+blh[2])*cos(blh[0])*cos(blh[1]); xyz[1]=(N+blh[2])*cos(blh[0])*sin(blh[1]); xyz[2]=(N*(1-e2)+blh[2])*sin(blh[0]); } //计算子午弧长