高斯-克吕格投影分带后是怎么拼接在同的？

匿名用户 6小时前 109 e1 高斯 高斯投影 高斯-克里格投影有6度分带和3度分带，分带投影后是像花瓣一样的一些列投影带，最后是怎么编程我们常见的完整的地图的呢？ 0 赞 0 踩 其他回答 两部分
第一部分：用圆柱投影是怎么画世界地图的
第二部分：用高斯-克里格投影坐标系下的坐标是怎么画地图的

第一部分：
哈哈哈哈哈哈，高斯-克吕格投影，是一种横轴等角圆柱投影，这所以加上6度带或者3度带的限制，是为了在工程应用等相关应用中，控制长度的变形

严格来说，任何一种投影方式都能画世界地图，或者任何区域的地图，只是变形程度的大小不同罢了

下图是高斯-克吕格投影

是的，没错，常用的世界地图中，有一种投影方式，用的就是切圆柱投影，
但是它不叫做高斯-克吕格投影，它叫做墨卡托投影
区别在于
1：高斯-克吕格投影中，那个圆柱是横着的，而在墨卡托投影中，那个圆柱是竖着的，也就是，所谓的“中央子午线”变成了赤道
2：墨卡托投影，不分带，它是直接把整个椭球投影到圆柱壳上，再把壳展开，如下的样子

也就是说，若要用圆柱投影画世界地图，我们根本就不分带！就不分花瓣，直接一次成型！原因在于：世界地图上，因为比例尺极小，所以精度什么的，不那么重要的，看的就是个形，你看那个墨卡托投影的世界地图，上方与下方的纬线网都变形成什么样子了！但是你在乎嘛？

第二部分：
假若你手里有四五个高斯-克里格投影的花瓣，怎么把这四五个花瓣拼成一个完整的地图呢？
步骤有以下几个
1：首先，把花瓣都重新贴在椭球上，也就是说，把每一个花瓣上的所有点，通过高斯反算，算出对应的经纬度（大地坐标）
2：选择一种你喜欢的投影方式，再重新投影不就完了？

这并不是我瞎BB的想当然的给你说，实质上，高斯-克里格投影坐标系下的地图，要缩编成比例尺更小、高斯-克里格投影无法表示的地图的时候，简单说来，就是以上两步走的，也只能是以上两步走的 热心网民 6小时前 0条评论 0 赞 0 踩 有一个变形公式----高斯正算公式，知道该点的经纬度以及分带的经度，代入公式计算就将其投影到了平面上，得到X，Y值~~~
下面是一段高斯正算 matlab算法：

function [x,y]=gauss_zhengsuan(B,L,L0)

l=L-L0;
l=deg2rad(l);
B=deg2rad(B);
a=6378137;
f=1/298.257223563;
b=a-a*f;
c=a^2/b;
e=sqrt(a^2-b^2)/a;
e1=sqrt(a^2-b^2)/b;
p=3600*180/pi;
Beta0=1-(3/4)*e1^2+(45/64)*e1^4-(175/256)*e1^6+(11025/16384)*e1^8;
Beta2=Beta0-1;
Beta4=(15/32)*e1^4-(175/384)*e1^6+(3675/8192)*e1^8;
Beta6=-(35/96)*e1^6+(735/2048)*e1^8;
Beta8=(315/1024)*e1^8;
Xb0=c*(Beta0*(B)+sin(B)*(Beta2*cos(B)+Beta4*cos(B)^3+Beta6*cos(B)^5+Beta8*cos(B)^7));
Np=a/sqrt(1-(e*sin(B))^2);
m0=l*cos(B);
t=tan(B);
in2=(e1)^2*cos(B)^2;
x=Xb0+(1/2)*Np*t*m0^2+(1/24)*(5-t^2+9*in2+4*in2^2)*Np*t*m0^4+(1/720)*(61-58*t^2+t^4+270*in2-330*in2*t^2)*Np*t*m0^6;
y=Np*m0+(1/6)*(1-t^2+in2)*Np*m0^3+(1/120)*(5-18*t^2+t^4+14*in2-58*t^2*in2)*Np*m0^5;

Ding 6小时前 0条评论 0 赞 0 踩 登录后可回答 提交回答