高斯投影坐标正反算编程报告

精品文档 1欢迎下载 高斯投影坐标正反算编程报告高斯投影坐标正反算编程报告 1 1 编程思想编程思想 进行高斯投影坐标正反算的编程需要牵涉到大量的公式 为了使程序条理更清楚 各块的数据复用性更强 这里采取了结构化的编程思想 程序由四大块组成 GeodesyHomework cpp 文件用于存放 main 函数 是整个程序的入口 通过结构 化的编程尽力使 main 函数变得简单 MyFunction h 和 MyFunction cpp 用于存放计算过程中进行角度弧度换算时所要用 到的一些自定的转换函数 Zhengsuan h 和 Zhengsuan cpp 用于存放 Zhengsuan 类 在 Zhengsuan 类中声明了 高斯投影坐标正算所要用到的所有变量 在类的构造函数中进行成员变量的初始化及 正算计算 通过 get 函数获得相应的正算结果 Fansuan h 和 Fansuan cpp 用于存放 Fansuan 类 类似于 Zhengsuan 类 Fansuan 类中声明了高斯投影坐标反算所要用到的所有变量 在类的构造函数中进行成员变量 的初始化及反算计算 通过 get 函数获得相应的反算结果 2 2 计算模型计算模型 高斯投影正算公式 6425 6 44223 4 2 2 5861 cossin 720 495 cos 24 cossin 2 lttBB N ltBsimB N lBB N Xx 5222425 5 3223 3 5814185 cos 120 1 cos 6 cos ltttB N ltB N lB N y 高斯投影反算公式 精品文档 2欢迎下载 22242 5 5 22 3 3 642 5 42222 3 2 8624285 cos120 21 cos6cos 459061 720 935 242 fffff ff ff ffff ff ff f ffff ff f ff f f ttt BN y t BN y BN y l ytty NM t ytt NM t y NM t BB 3 3 程序框图程序框图 4 4 计算结果计算结果 5 5 附录 程序代码附录 程序代码 主函数入口 GeodesyHomework cppGeodesyHomework cpp include MyFunction h include Zhengsuan h include Fansuan h 精品文档 3欢迎下载 include using namespace std void fansuan void zhengsuan void main zhengsuan fansuan printf n over void zhengsuan double myB myL cout 正算 endl cout 请输入大地纬度 B endl myB angleToDegree cout 请输入大地经度 L endl myL angleToDegree Zhengsuan myZhengsuan1 myB myL printf Radian B f L f n myZhengsuan1 getrB myZhengsuan1 getrL myZhengsuan1 printLocation 精品文档 4欢迎下载 void fansuan double myX myY cout 反算 endl cout 请输入国家统一坐标 X Y 例如 3378627 1819 20243953 4517 myX myY Fansuan myFansuan1 myX myY myFansuan1 printLocation 自定功能函数库 MyFunction hMyFunction h define PI 3 1415926 include using namespace std double angleToDegree int du int fen float miao double angleToDegree 将度分秒换算为度 double degreeToRadian double degree double degreeToRadian 将角度换算为弧度 MyFunction cppMyFunction cpp include MyFunction h double angleToDegree int du int fen float miao 精品文档 5欢迎下载 double result 0 result miao 3600 0 fen 60 0 du return result double angleToDegree int du fen float miao double result cout 请输入度 分 秒 例如 30 20 00 du fen miao result angleToDegree du fen miao return result double degreeToRadian double degree double result 0 result degree 57 295779513082321 return result double degreeToRadian double result degree degree angleToDegree result degreeToRadian degree 精品文档 6欢迎下载 return result 正算类 Zhengsuan hZhengsuan h Zhengsuan h interface for the Zhengsuan class if defined AFX ZHENGSUAN H 2655EA28 E810 44A3 8F14 56421A7B4466 INCLUDED define AFX ZHENGSUAN H 2655EA28 E810 44A3 8F14 56421A7B4466 INCLUDED if MSC VER 1000 pragma once endif MSC VER 1000 define rouSecond 206264 806247096355 include MyFunction h include include using namespace std class Zhengsuan public 精品文档 7欢迎下载 Zhengsuan Zhengsuan double fB double fL double getX double getY double getrB double getrL void printLocation virtual Zhengsuan private double x double y 大地坐标 double X double Y 国家统一坐标 double B double rB int Bsecond double L double rL 输入的大地纬度 B 大地经度 L rB rL 为对应弧度表示值 Bsecond 为换算成 秒数值 精品文档 8欢迎下载 int n 带号 n double L0 中央经线纬度 L0 double LDot 纬度差 L L0 int LDotSecond 换算成秒的纬度差 double l double N double a0 double a3 double a4 double a5 double a6 七个计算参数 endif defined AFX ZHENGSUAN H 2655EA28 E810 44A3 8F14 56421A7B4466 INCLUDED Zhengsuan cppZhengsuan cpp Zhengsuan cpp implementation of the Zhengsuan class 精品文档 9欢迎下载 include Zhengsuan h Construction Destruction Zhengsuan Zhengsuan Zhengsuan Zhengsuan double fB double fL B fB rB degreeToRadian fB L fL rL degreeToRadian fL Bsecond B 3600 初始化大地经度 L 大地纬度 B Bsecond 按弧度的大地纬度 rB n int L 6 1 初始化带号 n L0 6 n 3 中央经线经度 角度单位 精品文档 10欢迎下载 LDot L L0 经度差 LDotSecond LDot 3600 l LDot 3600 rouSecond 计算参数 l N 6399698 902 21562 267 108 973 0 612 cos rB cos rB cos rB cos rB cos rB cos rB 计算参数 N a0 32140 404 135 3302 0 7092 0 004 cos rB cos rB cos rB cos rB cos rB cos rB 计算参数 a0 a4 0 25 0 00252 cos rB cos rB cos rB cos rB 0 04166 计算参数 a4 a6 0 166 cos rB cos rB 0 084 cos rB cos rB 计算参数 a6 a3 0 3333333 0 001123 cos rB cos rB cos rB cos rB 0 1666667 计 算参数 a3 a5 0 0083 0 1667 0 1968 0 004 cos rB cos rB cos rB cos rB cos rB cos rB 计算参数 a5 x 6367558 4969 Bsecond rouSecond a0 0 5 a4 a6 l l l l l l N sin rB cos rB 正算 x y 1 a3 a5 l l l l l N cos rB 正算 y X x Y n 1000000 y 500000 国家统一坐标 精品文档 11欢迎下载 Zhengsuan Zhengsuan double Zhengsuan getX return X double Zhengsuan getY return Y void Zhengsuan printLocation printf 正算得国家统一坐标为 X 8 8f Y 8 8f n X Y double Zhengsuan getrB return rB double Zhengsuan getrL return rL 精品文档 12欢迎下载 反算类 Fansuan hFansuan h Fansuan h interface for the Fansuan class if defined AFX FANSUAN H 5B5E4668 EE81 463F 8D65 FFE2AEACE745 INCLUDED define AFX FANSUAN H 5B5E4668 EE81 463F 8D65 FFE2AEACE745 INCLUDED if MSC VER 1000 pragma once endif MSC VER 1000 define rouSecond 206264 806247096355 include include MyFunction h include using namespace std class Fansuan public Fansuan 精品文档 13欢迎下载 Fansuan double X double Y double getB double getL void printLocation virtual Fansuan private double x double y 高斯投影坐标 double X double Y int N 国家统一坐标 N 为带号 double B Bsecond double L 最后反算得到 B L double L0 中央经线经度 double l lsecond L L0 l L0 6 N 3 double Bf BfSecond BfDegree double beta betaSecond betaDegree double Z double Nf double b2 double b3 精品文档 14欢迎下载 double b4 double b5 计算的 8 个参数 endif defined AFX FANSUAN H 5B5E4668 EE81 463F 8D65 FFE2AEACE745 INCLUDED Fansuan cppFansuan cpp Fansuan cpp implementation of the Fansuan class include Fansuan h Construction Destruction Fansuan Fansuan Fansuan Fansuan double X double Y 精品文档 15欢迎下载 this X X this Y Y 初始化 x y N int Y 1000000 取出带号 L0 6 N 3 初始化该带号的中央经线经度 x X y Y 1000000 N 500000 beta x 6367558 4969 初始化 beta 弧度单位 betaSecond beta rouSecond 初始化 beta 秒单位 betaDegree betaSecond 3600 初始化 beta 整度数单位 Bf beta 50221746 293622 2350 22 cos beta cos beta cos beta cos b eta cos beta cos beta 1e 10 sin beta cos beta 初始化 Bf 弧度单位 BfSecond Bf rouSecond 初始化 Bf 秒单位 BfDegree BfSecond 3600 初始化 Bf 整度数单位 Nf 6399698 902 21562 267 108 973 0 612 cos Bf cos Bf cos Bf cos Bf