您所在位置: 网站首页 / 文档列表 / 图形图像 / 文档详情
高斯投影正算与反算的理论方法与实.pdf 立即下载
上传人:文库****品店 上传时间:2024-09-08 格式:PDF 页数:4 大小:107KB 金币:10 举报 版权申诉
预览加载中,请您耐心等待几秒...

高斯投影正算与反算的理论方法与实.pdf

高斯投影正算与反算的理论方法与实.pdf

预览

在线预览结束,喜欢就下载吧,查找使用更方便

10 金币

下载文档

如果您无法下载资料,请参考说明:

1、部分资料下载需要金币,请确保您的账户上有足够的金币

2、已购买过的文档,再次下载不重复扣费

3、资料包下载后请先用软件解压,在使用对应软件打开

高斯投影正算与反算的理论方法与实现代码高斯投影是正形投影的一种,同一坐标系中的高斯投影换带计算公式是根据正形投影原理推导出的两个高斯坐标系间的显函数式。在同一大地坐标系中(例如1954北京坐标系或1980西安坐标系),如果两个高斯坐标系只是主子午线的经度不同,那么显函数式前的系数可以根据坐标系使用的椭球元素和主子午线经度唯一确定。但如果两个高斯坐标系除了主子午线的经度不同以外,还存在其他线性系,则将线性变换公式代入换带计算的显函数式中,仍然可以得到严密的坐标变换公式。此时显函数式前的系数等价于使用两个坐标系主子午线的经度和线性变换参数联合求解得到的,可以唯一确定。//6度带宽54北京坐标系//高斯投影由大地坐标(Unit:Metres)反算经纬度(Unit:DD)voidGaussProjInvCal(doubleX,doubleY,double*longitude,double*latitude){intProjNo;intZoneWide;////带宽doublelongitude1,latitude1,longitude0,latitude0,X0,Y0,xval,yval;doublee1,e2,f,a,ee,NN,T,C,M,D,R,u,fai,iPI;iPI=0.0174532925199433;////3.1415926535898/180.0;a=6378245.0;f=1.0/298.3;//54年北京坐标系参数////a=6378140.0;f=1/298.257;//80年西安坐标系参数ZoneWide=6;////6度带宽ProjNo=(int)(X/1000000L);//查找带号longitude0=(ProjNo-1)*ZoneWide+ZoneWide/2;longitude0=longitude0*iPI;//中央经线X0=ProjNo*1000000L+500000L;Y0=0;xval=X-X0;yval=Y-Y0;//带内大地坐标e2=2*f-f*f;e1=(1.0-sqrt(1-e2))/(1.0+sqrt(1-e2));ee=e2/(1-e2);M=yval;u=M/(a*(1-e2/4-3*e2*e2/64-5*e2*e2*e2/256));fai=u+(3*e1/2-27*e1*e1*e1/32)*sin(2*u)+(21*e1*e1/16-55*e1*e1*e1*e1/32)*sin(4*u)+(151*e1*e1*e1/96)*sin(6*u)+(1097*e1*e1*e1*e1/512)*sin(8*u);C=ee*cos(fai)*cos(fai);T=tan(fai)*tan(fai);NN=a/sqrt(1.0-e2*sin(fai)*sin(fai));R=a*(1-e2)/sqrt((1-e2*sin(fai)*sin(fai))*(1-e2*sin(fai)*sin(fai))*(1-e2*sin(fai)*sin(fai)));D=xval/NN;//计算经度(Longitude)纬度(Latitude)longitude1=longitude0+(D-(1+2*T+C)*D*D*D/6+(5-2*C+28*T-3*C*C+8*ee+24*T*T)*D*D*D*D*D/120)/cos(fai);latitude1=fai-(NN*tan(fai)/R)*(D*D/2-(5+3*T+10*C-4*C*C-9*ee)*D*D*D*D/24+(61+90*T+298*C+45*T*T-256*ee-3*C*C)*D*D*D*D*D*D/720);//转换为度DD*longitude=longitude1/iPI;*latitude=latitude1/iPI;}//高斯投影由经纬度(Unit:DD)正算平面坐标(含带号,Unit:Metres)voidGaussProjCal(doublelongitude,doublelatitude,double*X,double*Y){intProjNo=0;intZoneWide;////带宽doublelongitude1,latitude1,longitude0,latitude0,X0,Y0,xval,yval;doublea,f,e2,ee,NN,T,C,A,M,iPI;iPI=0.0174532925199433;////3.1415926535898/180.0;ZoneWide=6;//
单篇购买
VIP会员(1亿+VIP文档免费下)

扫码即表示接受《下载须知》

高斯投影正算与反算的理论方法与实

文档大小:107KB

限时特价:扫码查看

• 请登录后再进行扫码购买
• 使用微信/支付宝扫码注册及付费下载,详阅 用户协议 隐私政策
• 如已在其他页面进行付款,请刷新当前页面重试
• 付费购买成功后,此文档可永久免费下载
年会员
99.0
¥199.0

6亿VIP文档任选,共次下载特权。

已优惠

微信/支付宝扫码完成支付,可开具发票

VIP尽享专属权益

VIP文档免费下载

赠送VIP文档免费下载次数

阅读免打扰

去除文档详情页间广告

专属身份标识

尊贵的VIP专属身份标识

高级客服

一对一高级客服服务

多端互通

电脑端/手机端权益通用

手机号注册 用户名注册
我已阅读并接受《用户协议》《隐私政策》
已有账号?立即登录
我已阅读并接受《用户协议》《隐私政策》
已有账号?立即登录
登录
手机号登录 微信扫码登录
微信扫一扫登录 账号密码登录

首次登录需关注“豆柴文库”公众号

新用户注册
VIP会员(1亿+VIP文档免费下)
年会员
99.0
¥199.0

6亿VIP文档任选,共次下载特权。

已优惠

微信/支付宝扫码完成支付,可开具发票

VIP尽享专属权益

VIP文档免费下载

赠送VIP文档免费下载次数

阅读免打扰

去除文档详情页间广告

专属身份标识

尊贵的VIP专属身份标识

高级客服

一对一高级客服服务

多端互通

电脑端/手机端权益通用