goole earth 基站扇區(qū)立體圖

昨天才發(fā)現(xiàn) 還有 goole earth 這種東西，就隨便做了一個。

接個圖給大家看看，不足之處請大家海涵。

樓下有問到是怎樣實(shí)現(xiàn)的，現(xiàn)在我說一下。


由于我對MAPIinfo BASIC并不熟悉，我是用EXCEL和VBA 實(shí)現(xiàn)的，畢竟VBA 相對來說 功能強(qiáng)大和友善。

也就是通過EXCEL的宏實(shí)現(xiàn)的。 把基站的相關(guān)信息與經(jīng)緯度做成EXCEL，然后通過宏的方式生成相關(guān)KML。

下面是一部分進(jìn)行繪圖計(jì)算的源碼。

Const Pi = 3.1415926535898
Private a, b, C, alpha, e, e2, w, V As Double
'a長軸半徑
'b短軸
'c極曲率半徑
'alpha扁率
'e第一偏心率
'e2第二偏心率
'W第一基本緯度函數(shù)
'V第二基本緯度函數(shù)
Private B1, L1, B2, L2 As Double
'B1點(diǎn)1的緯度
'L1點(diǎn)1的經(jīng)度
'B2點(diǎn)1的緯度
'L2點(diǎn)2的經(jīng)度
Private S As Double'''''大地線長度
Private A1, A2 As Double
'A1點(diǎn)1到點(diǎn)2的方位角
'A2點(diǎn)2到點(diǎn)1的方位角
Private ellipsetype As Integer '''''橢圓類型，對應(yīng)于列數(shù)
'1克拉索夫斯基橢球體
'21975年國際橢球體
'3WGS-84橢球體
'4海福特橢球
'5Bessel橢球體
'6自定義橢球體1
'7自定義橢球體2

'''''''''''''''公共參數(shù)定義完成

Private Sub getellipseparameter()
a = 6378245
b = 6356752.3142
C = a ^ 2 / b
alpha = (a - b) / a
e = Sqr(a ^ 2 - b ^ 2) / a
e2 = Sqr(a ^ 2 - b ^ 2) / b
End Sub
Private Function computerw()
w = Sqr(1 - e ^ 2 * (Sin(B1) ^ 2))
V = w * (a / b)
End Function
Function Computation(STARTLAT, STARTLONG, ANGLE1, DISTANCE As Double) As String '''''正算
ellipsetype = ellipsetypes
B1 = STARTLAT
L1 = STARTLONG
A1 = ANGLE1
S = DISTANCE
Call getellipseparameter
If B1 = 0 Then
If A1 = 90 Then
A2 = 270
B2 = 0 '''緯度相同
L2 = L1 + S / a * 180 / Pi
End If
If A1 = 270 Then
A2 = 90
B2 = 0 '''緯度相同
L2 = L1 - S / a * 180 / Pi
End If
Exit Function
End If

B1 = rad(B1)
L1 = rad(L1)
A1 = rad(A1)
Call computerw

Dim W1 As Double
E1 = e ''''第一偏心率

'// 計(jì)算起點(diǎn)的歸化緯度
W1 = w ''Sqr(1 - e1 * e1 * Sin(B1 ) * Sin(B1 ))
sinu1 = Sin(B1) * Sqr(1 - E1 * E1) / W1
cosu1 = Cos(B1) / W1
'// 計(jì)算輔助函數(shù)值
sinA0 = cosu1 * Sin(A1)
cotq1 = cosu1 * Cos(A1)
sin2q1 = 2 * cotq1 / (cotq1 ^ 2 + 1)
cos2q1 = (cotq1 ^ 2 - 1) / (cotq1 ^ 2 + 1)

'// 計(jì)算系數(shù)AA,BB,CC及AAlpha, BBeta的值。
cos2A0 = 1 - sinA0 ^ 2
e2 = Sqr(a ^ 2 - b ^ 2) / b
k2 = e2 * e2 * cos2A0

Dim aa, BB, cc, EE22, AAlpha, BBeta As Double
aa = b * (1 + k2 / 4 - 3 * k2 * k2 / 64 + 5 * k2 * k2 * k2 / 256)
BB = b * (k2 / 8 - k2 * k2 / 32 + 15 * k2 * k2 * k2 / 1024)
cc = b * (k2 * k2 / 128 - 3 * k2 * k2 * k2 / 512)
e2 = E1 * E1
AAlpha = (e2 / 2 + e2 * e2 / 8 + e2 * e2 * e2 / 16) - (e2 * e2 / 16 + e2 * e2 * e2 / 16) * cos2A0 + (3 * e2 * e2 * e2 / 128) * cos2A0 * cos2A0
BBeta = (e2 * e2 / 32 + e2 * e2 * e2 / 32) * cos2A0 - (e2 * e2 * e2 / 64) * cos2A0 * cos2A0


'// 計(jì)算球面長度
q0 = (S - (BB + cc * cos2q1) * sin2q1) / aa
sin2q1q0 = sin2q1 * Cos(2 * q0) + cos2q1 * Sin(2 * q0)
cos2q1q0 = cos2q1 * Cos(2 * q0) - sin2q1 * Sin(2 * q0)
q = q0 + (BB + 5 * cc * cos2q1q0) * sin2q1q0 / aa
'// 計(jì)算經(jīng)度差改正數(shù)



theta = (AAlpha * q + beta * (sin2q1q0 - sin2q1)) * sinA0

'// 計(jì)算終點(diǎn)大地坐標(biāo)及大地方位角
sinu2 = sinu1 * Cos(q) + cosu1 * Cos(A1) * Sin(q)
B2 = Atn(sinu2 / (Sqr(1 - E1 * E1) * Sqr(1 - sinu2 * sinu2))) * 180 / Pi
lamuda = Atn(Sin(A1) * Sin(q) / (cosu1 * Cos(q) - sinu1 * Sin(q) * Cos(A1))) * 180 / Pi

If (Sin(A1) > 0) Then
If (Sin(A1) * Sin(q) / (cosu1 * Cos(q) - sinu1 * Sin(q) * Cos(A1)) > 0) Then
lamuda = Abs(lamuda)
Else
lamuda = 180 - Abs(lamuda)
End If
Else
If (Sin(A1) * Sin(q) / (cosu1 * Cos(q) - sinu1 * Sin(q) * Cos(A1)) > 0) Then
lamuda = Abs(lamuda) - 180
Else
lamuda = -Abs(lamuda)
End If
End If
L2 = L1 * 180 / Pi + lamuda - theta * 180 / Pi

A2 = Atn(cosu1 * Sin(A1) / (cosu1 * Cos(q) * Cos(A1) - sinu1 * Sin(q))) * 180 / Pi
If (Sin(A1) > 0) Then
If (cosu1 * Sin(A1) / (cosu1 * Cos(q) * Cos(A1) - sinu1 * Sin(q)) > 0) Then
A2 = 180 + Abs(A2)
Else
A2 = 360 - Abs(A2)
End If
Else
If (cosu1 * Sin(A1) / (cosu1 * Cos(q) * Cos(A1) - sinu1 * Sin(q)) > 0) Then
A2 = Abs(A2)
Else
A2 = 180 - Abs(A2)
End If
End If

Computation = L2 & "," & B2
End Function
''''''下面是基本函數(shù)
Private Function rad(ByVal angle_d As Double) As Double
rad = angle_d * Pi / 180
End Function

[ 本帖最后由 whoami926 于 2008-3-5 13:54 編輯 ]