Option Explicit
Public Const Pi = 3.14159265 'PI
Public Const R = 6.371229 * 1000000# '地球半徑

Function Distance(Lon1 As Double, Lat1 As Double, Lon2 As Double, lat2 As Double) As Double
'該函數(shù)用于已知兩個(gè)點(diǎn)的經(jīng)緯度信息時(shí)的距離計(jì)算
'輸入的經(jīng)緯度信息，單位為度
'輸出為距離，單位為m

Dim x As Double, y As Double

'PI * R * Cos(((Lat1 + Lat2) / 2) * PI / 180) -- 在兩個(gè)緯度他們中間緯度的經(jīng)度間距，即每度經(jīng)度多少米
x = (Lon2 - Lon1) * Pi * R * Cos(((Lat1 + lat2) / 2) * Pi / 180) / 180 '兩個(gè)經(jīng)度差換算的經(jīng)度距離

'( PI * R )/180 -- 每度緯度多少米
y = (lat2 - Lat1) * Pi * R / 180 '兩個(gè)緯度差換算得到的緯度距離

Distance = Sqr(x * x + y * y) '兩點(diǎn)間的距離，平方和后開根號

End Function