PHP рассчитать широту / долготу квадрата вокруг заданной точки на поверхности? - PullRequest
Новая
       35

PHP рассчитать широту / долготу квадрата вокруг заданной точки на поверхности?

1 голос
/ 10 октября 2011

Учитывая широту / долготу точки земной поверхности (и расстояния, скажем, в километрах), мне нужно вычислить координаты сверху-слева, сверху-справа, снизу-справа и слева-вниз. Независимо от ориентации квадрата.

Я нашел хороший класс PHP для вычисления расстояния между двумя точками на поверхности, но, к сожалению, нет функции относительно квадрата (или даже круга) вокруг точки.

1 Ответ

4 голосов
/ 10 октября 2011

Это может сработать для вас.Подход здесь заключается в том, чтобы найти точки назначения из заданной точки с севера, с востока на запад, с запада и с юга.Функция getSquareAroundPoint возвращает массив, определяющий эти четыре точки, а функция getMinMaxCoords возвращает минимальную и максимальную координаты этих точек, как того требует упомянутый вами REST API (он возвращается как массив массива, есливходная координата находится около 180-градусного меридиана.) Она находится в свободном доступе. 

// Distance is in km, alat and alon are in degrees
function getDestinationPoint($alat, $alon, $distance, $bearing){
 $pi=3.14159265358979;
 $alatRad=$alat*$pi/180;
 $alonRad=$alon*$pi/180;
 $bearing=$bearing*$pi/180;
 $alatRadSin=sin($alatRad);
 $alatRadCos=cos($alatRad);
 // Ratio of distance to earth's radius
 $angularDistance=$distance/6370.997;
 $angDistSin=sin($angularDistance);
 $angDistCos=cos($angularDistance);
 $xlatRad = asin( $alatRadSin*$angDistCos +
                                   $alatRadCos*$angDistSin*cos($bearing) );
 $xlonRad = $alonRad + atan2(
            sin($bearing)*$angDistSin*$alatRadCos,
            $angDistCos-$alatRadSin*sin($xlatRad));
 // Return latitude and longitude as two element array in degrees
 $xlat=$xlatRad*180/$pi;
 $xlon=$xlonRad*180/$pi;
 if($xlat>90)$xlat=90;
 if($xlat<-90)$xlat=-90;
 while($xlat>180)$xlat-=360;
 while($xlat<=-180)$xlat+=360;
 while($xlon>180)$xlon-=360;
 while($xlon<=-180)$xlon+=360;
 return array($xlat,$xlon);
}

// Distance is in km, lat and lon are in degrees
function getSquareAroundPoint($lat,$lon,$distance){
 return array(
  getDestinationPoint($lat,$lon,$distance,0), // Get north point
  getDestinationPoint($lat,$lon,$distance,90), // Get east point
  getDestinationPoint($lat,$lon,$distance,180), // Get south point
  getDestinationPoint($lat,$lon,$distance,270) // Get west point
 );
}

// Returns array containing an array with min lat, max lat, min lon, max lon
// If the square defining these points crosses the 180-degree meridian, two
// such arrays are returned.  Otherwise, one such array (within another array)
// is returned.
function getMinMaxCoords($lat,$lon,$distance){
 $s=getSquareAroundPoint($lat,$lon,$distance);
 if($s[3][1]>$s[1][1]){// if west longitude is greater than south longitude
  // Crossed the 180-degree meridian
  return array(
    array($s[2][0],$s[0][0],$s[3][1],180),
    array($s[2][0],$s[0][0],-180,$s[1][1])
  );
 } else {
  // Didn't cross the 180-degree meridian (usual case)
  return array(
    array($s[2][0],$s[0][0],$s[3][1],$s[1][1])
  );
 }
}

// Example: Gets extreme coordinates around point at (10.0,20.0)
print_r(getSquareAroundPoint(10.0,180,100));
print_r(getMinMaxCoords(10.0,180,100));
...