地图学实习讲解.docVIP

地图学实习报告 实习一 正轴等角割圆锥投影源程序的设计 圆锥投影，特别是正轴等角割圆锥投影，在我国以及一些中等纬度的国家应用得非常广泛。例如，编制中华人民共和国地图（南海诸岛作插图）以及各省（区）地图都是采用这种投影。新中国成立之前，我国地形图也是采用该投影作为数学基础的。可见，等角割圆锥投影可适用编制我国不同比例尺、不同类型的地图。 一、投影公式 Φ0:长度比最小的纬度; M0:比例尺分母。 二、制图区域 河北省 制图区域中央经线：λ0=116000′ 标准纬线：Ф1=37030′,Ф2=41000′ 主比例尺：1：500万 经纬网密度：△Ф=△λ=10 三、要求 1、编写正轴等角圆锥投影的源程序。在编制源程序时，采用通用程序为好，为此， 地球椭球的数据可采用IAG-1975椭球：a=6378140m，b=6356755m。 2、输出成果要求完整、整齐，有条件时可采用汉字化、表格化。输出成果包括： U1,U2; 投影常数：α 积分常数：K 长度比为最小的纬线纬度：φ0 纬度：φ 经度：λ 直角坐标：X，Y 3.绘制该地图的经纬线网格。 四．实习实现代码 #include stdio.h #include iostream #include cmath using namespace std; const double PI=3.1415926; const double a1=6378140; const double b1=6356755; double get_e() { return(sqrt((a1*a1-b1*b1)/(a1*a1))); } double get_N(double m) { return(a/sqrt(1-get_e()*get_e()*sin(m*PI/180)*sin(m*PI/180))); } double get_r(double m) { return(get_N(m)*cos(m*PI/180)); } double get_U(double m) //得到U值的函数 { return(tan(45*PI/180+m*PI/360)/pow((1-e*sin(m*PI/360)/(1+e*sin(m*PI/360))),e/2); } double get_a(double m1, double m2) //得到a的函数 { double r1,r2; double U1,U2; r1=log10(get_r(m1)); r2=log10(get_r(m2)); U1=log10(get_U(m1)); U2=log10(get_U(m2)); /*double m; m=log10(r1);*/ return((r1-r2)/(U2-U1)); } double get_C(double m1,double m2,) //得到C的函数 { double U=pow(get_U(m1),get_a(m1,m2)); double no=get_r(m1)*U/get_a(m1,m2); return(( no)); } double get_weidu(double m1, double m2) { double wd=asin(get_a( m1, m2)); return(wd*180/PI); } double get_jiajiao(double m1, double m2,double y) { return(get_a(m1,m2)*y); } double get_p(double m1,double m2,double x) { return(get_C(m1,m2)/pow(get_U(x),get_a(m1,m2))); } double get_x(double ps,double m1, double m2,double x,double y) { return(ps-get_p(m1,m2,x,M)*cos(get_jiajiao(m1,m2,y))); } double get_y(double m1, double m2,int M,double x,double y) { return(get_p(m1,m2,x)*sin(get_jiajiao(m1,m2,y))); } //******************************************************************// //