说明

之前一个帖子给出了正反算的C语言代码
链接如下：
link

我把它移植到了matlab中
支持向量输入

高斯投影正算

高斯投影正算matlab代码

function [X, Y] = GaussProjCal(Lon, Lat)
%%  高斯投影正算： 经纬度 算 直角坐标
% Longitude  Latitude 均为角度制
Lon = deg2rad(Lon);
Lat = deg2rad(Lat);ZoneWide = deg2rad(6); %带宽为6度
a=6378140.0; f=1/298.257;  %80年西安坐标系参数ProjNo = floor(Lon./ZoneWide);lon0 = ProjNo.*ZoneWide+ZoneWide./2;e2 = 2*f-f*f;
ee = e2*(1.0-e2);
NN = a./sqrt(1.0-e2.*sin(Lat).*sin(Lat));
T = tan(Lat).*tan(Lat);
C = ee.*cos(Lat).*cos(Lat);
A = (Lon-lon0).*cos(Lat);
M = a.*((1-e2/4-3*e2*e2/64-5*e2*e2*e2/256).*Lat-(3*e2/8+3*e2*e2/32+45*e2*e2*e2/1024).*sin(2.*Lat)...+(15*e2*e2/256+45*e2*e2*e2/1024).*sin(4.*Lat)-(35*e2*e2*e2/3072).*sin(6.*Lat));xval = NN.*(A+(1-T+C).*A.*A.*A./6+(5-18.*T+T.*T+72.*C-58.*ee).*A.*A.*A.*A.*A/120);
yval = M+NN.*tan(Lat).*(A.*A./2+(5-T+9.*C+4.*C.*C).*A.*A.*A.*A./24 ...+(61-58.*T+T.*T+600.*C-330.*ee).*A.*A.*A.*A.*A.*A./720);X0 = 1000000.*(ProjNo+1)+500000;
Y0 = 0; xval = xval+X0; yval = yval+Y0;
X = xval;
Y = yval;

高斯投影反算

高斯投影反算matlab代码。

function [Lon, Lat] = GaussProjInvCal(X,Y)
%%   高斯投影反算：  直角坐标 反算 经纬度
a=6378140.0; f=1/298.257; %80年西安坐标系参数
ZoneWide = deg2rad(6); %6度带宽
ProjNo = floor(X./1000000);longitude0 = (ProjNo-1).*ZoneWide+ZoneWide./2;X0 = ProjNo.*1000000+500000;
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）
Lon = 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);
Lon = rad2deg(Lon);
Lat = 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);
Lat = rad2deg(Lat);

高斯投影法正反算代码MATLAB版本相关推荐

高斯坐标正反算c语言,高斯投影坐标正反算编程报告
<高斯投影坐标正反算编程报告>由会员分享,可在线阅读,更多相关<高斯投影坐标正反算编程报告(11页珍藏版)>请在人人文库网上搜索. 1.高斯投影坐标正反算编程报告1. 编程思想 ...
高斯投影坐标正反算公式
§8.3高斯投影坐标正反算公式任何一种投影①坐标对应关系是最主要的:②如果是正形投影,除了满足正形投影的条件外(C-R偏微分方程),还有它本身的特殊条件. 8.3.1高斯投影坐标正算公式: 高斯投 ...
matlab写函数进行坐标正反算,坐标正反算的MATLAB GUI设计与应用
第 30 卷第 4 期 2018 年 10 月黄河水利职业技术学院学报 JournalofYellowRiverConservancyTechnicalInstitute Vol．30 No．4 ...
高斯投影坐标正反算算法（python）
# 高斯坐标正反算算法 # 在正算的时候,y的第一个数字表示是用几度带换算的,第二个数字和第三个数字表示带号 # 该程序是根据武汉大学出版社出版的第三版的<大地测量学基础>编写 # 用于课 ...
matlab写函数进行坐标正反算,基于matlab的坐标正反算
<基于matlab的坐标正反算>由会员分享,可在线阅读,更多相关<基于matlab的坐标正反算(6页珍藏版)>请在人人文库网上搜索. 1.基于 matlab 的坐标正反算测量程 ...
OTB官方评估代码matlab版本
OTB最初就是MATLAB版本评估,后来又推出了python版本,不过精度貌似没有matlab高 python版本配置可以参考我另一篇博客:OTB官方评估代码python版本 tracker_benc ...
c语言坐标正反算代码,C语言代码逆向的实战十篇（原创）
); } printf("\n"); system("pause"); } (之前从VS2008复制过来的时候代码有高亮的,现在换了VC6就木有了..) 进入正 ...
高斯投影正反算C语言程序代码,高斯投影正反算c代码
<高斯投影正反算c代码>由会员分享,可在线阅读,更多相关<高斯投影正反算c代码(11页珍藏版)>请在人人文库网上搜索. 1.高斯投影正反算程序设计一程序设计流程本程序的设计思路 ...
matlab高斯正反算程序6,基于matlab的高斯投影正反算与相邻带坐标换算程序设计...
第卷第期在月中国水运基于的高斯投影正反算与相邻带坐标换算程序设计徐翰 ,周强波 (核工业二三研究所 ,湖南长沙 ) 摘要 :地图投影方法众多 , ...

高斯投影法正反算代码MATLAB版本

高斯投影法正反算代码MATLAB版本

说明

高斯投影正算

高斯投影反算

高斯投影法正反算代码MATLAB版本相关推荐

最新文章

热门文章