我不我不知道主题的投影。如果wt文件转换成shp文件,首先要注意mapgis中的单位和arcgis中的单位是不一样的。前者是mm,后者是m,需要投影变换。另外,arcgis中的经纬度是用度表示的,mapgis多为度、分、秒。您可以直接在arcgis中添加xy数据并将其导出为shp文件。
下面以54坐标系为例,介绍具体的计算方法。
完成经纬度bl到平面直角坐标xy的转换。在excel中,选择输入公式的起始单元格,例如第2行第1列(单元格a2)为起始单元格,每个单元格的格式如下:
细胞;单元格内容;解释
a2;输入中央经线的度、分、秒;如果是115度30分,输入115.30;初始数据l0
b2;int(a2)(int(a2*100)-int(a2)*100)/60(a2*10000-int(a2*100)*100)/3600;将l0转换成度数
c2;以十进制形式输入纬度值,如38°14′20″,然后输入38.1420;初始数据b
d2;以十进制度数输入经度值;初始数据l
e2;int(c2)(int(c2*100)-int(c2)*100)/60(c2*10000-int(c2*100)*100)/3600;将b转换成度数
f2;int(d2)(int(d2*100)-int(d2)*100)/60(d2*10000-int(d2*100)*100)/3600;将l转换成度数
g2;f2-b2;左-l0
h2;g2/57.2955130823;变成弧度
i2;tan(弧度(e2));谭(乙)
j2;cos(弧度(e2));cos(b)
k2;0.006738525415j2*j2
l2;i2*i2
m2;1个k2
n2;6399698.9018/m2
o2;h2*h2*j2*j2
p2;i2*j2
q2;p2*p2
r2;(32005.78006q2*(133.92133q2*0.7031))
s2;6367558.49686*e2/57.29577951308-p2*j2*r2((((l2-58)*l261)*
o2/30(4*k25)*m2-l2)*o2/121)*n2*i2*o2/2
计算结果x
t2;((((l2-18)*l2-(58*l2-14)*k25)*o2/20m2-l2)*o2/61)*n2*(h2*j2)
计算结果y
根据上表中的公式,将其输入相应的单元格后,就可以很容易地由经纬度得到平面直角坐标。当所有的经纬度都输入后,用鼠标下拉就可以得到所有的计算结果。
原则:
设地球上一点的经度为λ,纬度为φ,
那么这个点的空间坐标是
xcos(φ)*cos(λ)
ycos(phi)*sin(lambda)
zsin(φ)
设地球上两点的空间坐标分别为(x1,y1,z1)和(x2,y2,z2)。
那么它们的夹角就是
aacos(x1*x2y1*y2z1*z2),其中a为角度。
那么两地之间的距离是
a/180*pi*r,其中r是地球的平均半径6371。
误差不超过1%
原文标题:坐标xy与度分秒之间的关系 用mapgis转了shp文件以后导入mapgis无法编辑为什么?,如若转载,请注明出处:https://www.bjhtrmyl.com/tag/20893.html
免责声明:此资讯系转载自合作媒体或互联网其它网站,「亨通号」登载此文出于传递更多信息之目的,并不意味着赞同其观点或证实其描述,文章内容仅供参考。