坐标xy与度分秒之间的关系 用mapgis转了shp文件以后导入mapgis无法编辑为什么?
我不 我不知道主题的投影。如果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.006738525415 J2 * J2
L2;I2*I2
M2;1个K2
N2;6399698.9018/M2
O2;H2*H2*J2*J2
P2;I2*J2
Q2;P2*P2
R2;(32005.78006 Q2 *(133.92133 Q2 * 0.7031))
S2;6367558.49686 * E2/57.29577951308-P2 * J2 * R2((((L2-58)* L2 61)*
O2/30(4 * K2 5)* M2-L2)* O2/12 1)* N2 * I2 * O2/2
计算结果x
T2;((((L2-18)*L2-(58*L2-14)*K2 5)* O2/20 M2-L2)* O2/6 1)* N2 *(H2 * J2)
计算结果y
根据上表中的公式,将其输入相应的单元格后,就可以很容易地由经纬度得到平面直角坐标。当所有的经纬度都输入后,用鼠标下拉就可以得到所有的计算结果。
原则:
设地球上一点的经度为λ,纬度为φ,
那么这个点的空间坐标是
xcos(φ)* cos(λ)
ycos(phi)*sin(lambda)
zsin(φ)
设地球上两点的空间坐标分别为(x1,y1,Z1)和(x2,y2,z2)。
那么它们的夹角就是
Aacos(x1*x2 y1*y2 z1*z2),其中a为角度。
那么两地之间的距离是
A/180*pi*R,其中R是地球的平均半径6371。
误差不超过1%
版权声明:本文内容由互联网用户自发贡献,本站不承担相关法律责任.如有侵权/违法内容,本站将立刻删除。