从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表中的差异。