從kml文件計算地面覆蓋角落的經度/ lng

Question

我需要在php或javascript中找到在kml文件中給出的地面覆蓋圖的緯度/經度。

I.e.對於一個具體的例子，我需要得到：

<LatLonBox> 
    <north>60.406505416667</north> 
    <south>60.400570555556</south> 
    <east>5.3351572222222</east> 
    <west>5.3190577777778</west> 
    <rotation>3.7088732260919</rotation> 
    </LatLonBox> 

到頂點座標

SW: 60.400316388889;5.3194425 
SE: 60.400824722222;5.3355405555556 
NE: 60.406759444444;5.3347738888889 
NW: 60.406251388889;5.3186730555556 

我可以

$w=($nw_lng+$sw_lng)/2; 
$e=($ne_lng+$se_lng)/2; 
$n=($ne_lat+$nw_lat)/2; 
$s=($se_lat+$sw_lat)/2; 
$rot= rad2deg (atan (($nw_lng - $sw_lng)/($sw_lat - $nw_lat)/2 )); 

得到其他方式（大約至少，PHP代碼中給出）應該很容易回來，但我已經用了幾個小時沒有到達那裏。有小費嗎？

Answer 1

您需要使用spherical trigonometry，作爲spherical geometry的一部分才能獲得完整的準確性。然而，由於你只處理球體的一小部分，如果你記住一件事情，歐幾里德幾何圖形就會做。

隨着緯度的增加，經度線越來越接近。例如，在北極附近，緯線幾乎接觸。因此，要處理你的緯度差異，通過多因素cos（緯度）來減少它們。這會爲您的應用帶來足夠的準確度。

$n = 60.406505416667; 
$s = 60.400570555556; 
$e = 5.3351572222222; 
$w = 5.3190577777778; 
$rotn = 3.7088732260919; 

$a = ($e + $w)/2.0; 
$b = ($n + $s)/2.0; 
$squish = cos(deg2rad($b)); 
$x = $squish * ($e - $w)/2.0; 
$y = ($n - $s)/2.0; 

$ne = array(
    $a + ($x * cos(deg2rad($rotn)) - $y * sin(deg2rad($rotn))) /$squish, 
    $b + $x * sin(deg2rad($rotn)) + $y *cos(deg2rad($rotn)) 
    ); 
$nw = array(
    $a - ($x * cos(deg2rad($rotn)) + $y * sin(deg2rad($rotn))) /$squish, 
    $b - $x * sin(deg2rad($rotn)) + $y *cos(deg2rad($rotn)) 
    ); 
$sw = array(
    $a - ($x * cos(deg2rad($rotn)) - $y * sin(deg2rad($rotn))) /$squish, 
    $b - $x * sin(deg2rad($rotn)) - $y *cos(deg2rad($rotn)) 
    ); 
$se = array(
    $a + ($x * cos(deg2rad($rotn)) + $y * sin(deg2rad($rotn))) /$squish, 
    $b + $x * sin(deg2rad($rotn)) - $y *cos(deg2rad($rotn)) 
    ); 
print_r(array(
'sw'=>$sw, 
'se'=>$se, 
'ne'=>$ne, 
'nw'=>$nw, 
)); 

我的$squish變量是我提到的cos（lat）。有水平長度的相對部分去除。正弦表看起來像這樣：

NE: (a + x cos A - y sin A, b + x sin A + y cos A) 
NW: (a - x cos A - y sin A, b - x sin A + y cos A) 
SW: (a - x cos A + y sin A, b - x sin A - y cos A) 
SE: (a + x cos A + y sin A, b + x sin A - y cos A) 

也許tttppp可能會解釋tttppp表中的差異。