测绘应用程序设计.doc

华北科技学院课程设计说明书课程名称： 测绘程序课程设计 专业班级： 测绘B08-3班 学 号： 200805064333 姓 名： 周 涛 指导教师： 赵 亚 红 坐标转换程序设计 一、设计目的和任务计算机程序设计综合课程设计是测绘工程专业一个综合性的实践教学环节，其目的在于促进学生复习和巩固计算机软件设计知识，加深对软件设计方法、 软件设计技术和设计思想的理解，并能运用所学软件设计知识和C语言技术进行综合软件设计，通过本课程设计能够进行简单软件系统的开发，掌握软件设计的方法和C语言程序设计的基本技术，提高学生的综合应用能力。本次课程设计共两周，在这两周内，要求同学们主要完成一项软件设计与调试任务，通过课程设计，每人提交运行通过的相应软件1套，程序代码1套，设计说明书1本。二、软件设计的一般步骤与方法通常计算机软件的设计以工程的管理方式、方法进行，包括需求分析、系统设计、编制程序，调试与测试。1. 需求分析 准确把握所要开发系统的功能需求。进行业务流的分析，根据系统所要完成的功能和所要处理的数据，画出反映系统逻辑模型的数据流图，对于重要的加工和数据项在数据字典中进行进一步的补充说明。2、总体设计 根据系统所完成的功能，准确划分各功能模块，明确各模块的功能职责、所要完成的数据处理与其他模块的数据传递关系，进行数据表的设计和对象类的设计，明确对象类间的关系。对各个模块进行初步的算法设计。3、详细设计 对各功能模块，准确划出其流程图，详细地描述算法。对于各对象类，详细地描述其内部功能实现。在具体的DBMS，如SQL SERVER中，完成数据的物理实施。4、编程调试 选用自己熟悉的编程工具，如C语言、VB或C#完成相应模块的设计。5、单元测试、组装测试 完成各模块的测试，成功后，由小组负责人将同一小组分别完成的各模块进行组装测试，成功后，提交软件成果。三.实验原理与过程1.问题需求分析 本次设计软件为高斯投影坐标转换的程序设计，包括坐标正反算，分为6度带和3度带，即大地坐标与测量坐标间的转换，同时还包括点所在不同带的带号，比以前所做的高斯投影较为复杂。此程序可以实现以下功能：1.1由大地坐标转换为6度带高斯平面坐标，其中包含点所在高斯坐标系的带号；1.2由高斯平面坐标转换为大地坐标；1.3统一转换为点所在国家统一坐标系中的坐标。2.实验要求2.1实验尽量要求可以应用实际，并且符合实际工作需要；2.2程序应该有一定容错能力，并且最好可以检验数据的小量错误。2.3数据流程图数据输出开始数据输入数据处理数据显示说明结束3.程序流程图 开始大地坐标转换为高斯平面坐标高斯平面坐标转换为大地坐标输入点大地经度和纬度选择要转换类型计算点所在6度带带号和高斯平面坐标，以及国家统一坐标系坐标输入高斯平面坐标及所在带号计算点大地经度和纬度，以及转化为国家统一坐标系坐标输出结果结束四、算例成果图五课程设计心得本次实验为测绘程序设计课程综合性设计，主要目的是巩固和检验学习的内容，能掌握测绘程序的编写及应用，掌握更多编程方面的知识。此次实验采用C+软件进行程序的编写，经过编写及实验，最终得出结果，并将结果保存附于实验报告中。这次课程设计在以前实验做过，但又不完全一样，内容更多了，所得出的结果更丰富，更符合实际的需要。通过这次综合性课程设计，我学到了很多，也对以前所学的知识能更好的掌握了，基本上达到了实验的目的。六程序代码及算例验证#include iostream.h#include iomanip.h#include math.h#define p 206264.806247096355double x,y,B,L,Bm,Lm;int Bd,Bf,Ld,Lf,n;void ZBZS(int Bd,int Bf,double Bm,int Ld,int Lf,double Lm);void ZBFS(double x,double y,int n);void main()int k;for(int i=0;i3;i+)coutPress 1:ZBZS; 2:ZBFSk; cout setiosflags(ios:fixed)setprecision(11); /坐标正算if(k=1) coutPlease input Bd,Bf,Bm,Ld,Lf,Lm:BdBfBmLdLfLm; /定6度带带号（逢余进一） if(Ld%6=0&(Lf*1.0/6+Lm*1.0/3600)=0) n=(Ld+Lf*1.0/6+Lm*1.0/3600)/6; else n=1+(Ld+Lf*1.0/60+Lm*1.0/3600)/6;/判断点是否在重叠区内 if(Ld*60+Lf+Lm/60-6*(n-1)*60)-15)|(Ld*60+Lf+Lm/60-6*n*60)-15) cout该点位于6度相邻两带重叠区，有两套坐标endl; n=1+(Ld+Lf*1.0/60+Lm*1.0/3600)/6; ZBZS(Bd,Bf,Bm,Ld,Lf,Lm); n=n-1; ZBZS(Bd,Bf,Bm,Ld,Lf,Lm); /在重叠区内 else /不在重叠区内 ZBZS(Bd,Bf,Bm,Ld,Lf,Lm); /用坐标反算检查结果 ZBFS(x,y,n); break;/坐标反算 else if(k=2) coutPlease input x,y,n:xyn; ZBFS(x,y,n); ZBZS(Bd,Bf,Bm,Ld,Lf,Lm); break; elsecoutPlease press 1 or 2!endl; continue;void ZBZS(int Bd,int Bf,double Bm,int Ld,int Lf,double Lm)double l,N,a0,a4,a6,a3,a5,t1;l=N=a0=a4=a6=a3=a5=t1=x=y=0.0;int L0=0;L0=6*n-3;l=(Ld*3600+Lf*60+Lm-L0*3600)/p;t1=cos(Bd*3600+Bf*60+Bm)/p)*cos(Bd*3600+Bf*60+Bm)/p);N=6399698.902-(21562.267-(108.973-0.612*t1)*t1)*t1;a0=32140.404-(135.3302-(0.7092-0.004*t1)*t1)*t1;a4=(0.25+0.00252*t1)*t1-0.04166;a6=(0.166*t1-0.084)*t1;a3=(0.3333333+0.001123*t1)*t1-0.1666667;a5=0.0083-(0.1667-(0.1968+0.004*t1)*t1)*t1;x=6367558.4969*(Bd*3600+Bf*60+Bm)/p-(a0-(0.5+(a4+a6*l*l)*l*l)*l*l*N)*sin(Bd*3600+Bf*60+Bm)/p)*cos(Bd*3600+Bf*60+Bm)/p);y=(1+(a3+a5*l*l)*l*l)*l*N*cos(Bd*3600+Bf*60+Bm)/p);coutZBZS:endl;coutl=lendl; coutN=Nendl; couta0=a0endl;couta4=a4endl;couta6=a6endl;couta3=a3endl;couta5=a5endl;cout该点在6度第n带的高斯坐标为:endl;coutx=xendl;couty=yendl;cout该点的国家统一坐标为：x,ny+500000endl;void ZBFS(double x,double y,int n)double beta,Bf1,l,Z,Nf,b2,b3,b4,b5,t2,t3;beta=Bf1=Z=Nf=b2=b3=b4=b5=t2=t3=l=0.0;int L0=0; beta=(x/6367558.4969)*p;t2=cos(beta/p)*cos(beta/p);Bf1=beta+(50221746+(293622+(2350+22*t2)*t2)*t2)*sin(beta/p)*cos(beta/p)*pow(10,-10)*p;t3=cos(Bf1/p)*cos(Bf1/p);Nf=6399698.902-(21562.267-(108.973-0.612*t3)*t3)*t3;Z=y/(Nf*cos(Bf1/p);b2=(0.5+0.003369*t3)*sin(Bf1/p)*cos(Bf1/p);b3=0.333333-(0.166667-0.001123*t3)*t3;b4=0.25+(0.16161+0.00562*t3)*t3;b5=0.2-(0.1667-0.0088*t3)*t3;B=Bf1-(1-(b4-0.12*Z*Z)*Z*Z)*Z*Z*b2*p;l=(1-(b3-b5*Z*Z)*Z*Z)*Z*p;L0=6*n-3;Bd=(int)(B/3600);Bf=(B/3600-Bd)*60;Bm=B-Bd*3600-Bf*60;Ld=L0+(int)(l/3600);Lf=(l/3600-(int)(l/3600)*60;Lm=l-(int)(l/3600)*3