我不 我不知道主題的投影。如果WT文件轉換成SHP文件，首先要注意mapgis中的單位和arcgis中的單位是不一樣的。前者是mm，后者是m，需要投影變換。另外，arcgis中的經緯度是用度表示的，m

我不 我不知道主題的投影。如果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%