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

我不 我不知道主題的投影。如果WT文件轉(zhuǎn)換成SHP文件，首先要注意mapgis中的單位和arcgis中的單位是不一樣的。前者是mm，后者是m，需要投影變換。另外，arcgis中的經(jīng)緯度是用度表示的，mapgis多為度、分、秒。您可以直接在arcgis中添加xy數(shù)據(jù)并將其導(dǎo)出為shp文件。

下面以54坐標(biāo)系為例，介紹具體的計(jì)算方法。

完成經(jīng)緯度BL到平面直角坐標(biāo)XY的轉(zhuǎn)換。在EXCEL中，選擇輸入公式的起始單元格，例如第2行第1列(單元格A2)為起始單元格，每個(gè)單元格的格式如下:

細(xì)胞；單元格內(nèi)容；解釋

a2；輸入中央經(jīng)線的度、分、秒；如果是115度30分，輸入115.30；初始數(shù)據(jù)L0

B2；INT(A2)(INT(A2 * 100)-INT(A2)* 100)/60(A2 * 10000-INT(A2 * 100)* 100)/3600；將L0轉(zhuǎn)換成度數(shù)

C2；以十進(jìn)制形式輸入緯度值，如38° 14′20 ″,然后輸入38.1420；初始數(shù)據(jù)b

D2；以十進(jìn)制度數(shù)輸入經(jīng)度值；初始數(shù)據(jù)l

E2；INT(C2)(INT(C2 * 100)-INT(C2)* 100)/60(C2 * 10000-INT(C2 * 100)* 100)/3600；將b轉(zhuǎn)換成度數(shù)

F2；INT(D2)(INT(D2 * 100)-INT(D2)* 100)/60(D2 * 10000-INT(D2 * 100)* 100)/3600；將l轉(zhuǎn)換成度數(shù)

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個(gè)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

計(jì)算結(jié)果x

T2；((((L2-18)*L2-(58*L2-14)*K2 5)* O2/20 M2-L2)* O2/6 1)* N2 *(H2 * J2)

計(jì)算結(jié)果y

根據(jù)上表中的公式，將其輸入相應(yīng)的單元格后，就可以很容易地由經(jīng)緯度得到平面直角坐標(biāo)。當(dāng)所有的經(jīng)緯度都輸入后，用鼠標(biāo)下拉就可以得到所有的計(jì)算結(jié)果。

原則:

設(shè)地球上一點(diǎn)的經(jīng)度為λ，緯度為φ，

那么這個(gè)點(diǎn)的空間坐標(biāo)是

xcos(φ)* cos(λ)

ycos(phi)*sin(lambda)

zsin(φ)

設(shè)地球上兩點(diǎn)的空間坐標(biāo)分別為(x1，y1，Z1)和(x2，y2，z2)。

那么它們的夾角就是

Aacos(x1*x2 y1*y2 z1*z2)，其中a為角度。

那么兩地之間的距離是

A/180*pi*R，其中R是地球的平均半徑6371。

誤差不超過1%