PHP 언어로하세요 :
// 고도를 0으로 가정
$ earthR = 6371; // km 단위 (= 3959 마일)
$ LatA = 37.418436;
$ LonA = -121.963477;
$ DistA = 0.265710701754;
$ LatB = 37.417243;
$ LonB = -121.961889;
$ DistB = 0.234592423446;
$ LatC = 37.418692;
$ LonC = -121.960194;
$ DistC = 0.0548954278262;
/ *
#Authoric Sphere 사용
# 타원체를 사용하는 경우이 단계는 약간 다릅니다
# 측지 위도 / 경도를 ECEF xyz로 변환
# 1. 위도 / 경도를 라디안으로 변환
# 2. Lat / Long (radians)을 ECEF로 변환
* /
$ xA = $ earthR * (cos (deg2rad ($ LatA)) * cos (deg2rad ($ LonA)));
$ yA = $ earthR * (cos (deg2rad ($ LatA)) * sin (deg2rad ($ LonA)));
$ zA = $ earthR * (sin (deg2rad ($ LatA)));
$ xB = $ earthR * (cos (deg2rad ($ LatB)) * cos (deg2rad ($ LonB)));
$ yB = $ earthR * (cos (deg2rad ($ LatB)) * sin (deg2rad ($ LonB)));
$ zB = $ earthR * (sin (deg2rad ($ LatB)));
$ xC = $ earthR * (cos (deg2rad ($ LatC)) * cos (deg2rad ($ LonC)));
$ yC = $ earthR * (cos (deg2rad ($ LatC)) * sin (deg2rad ($ LonC)));
$ zC = $ earthR * (sin (deg2rad ($ LatC)));
/ *
설치 :
sudo 배 설치 Math_Vector-0.7.0
sudo 배 설치 Math_Matrix-0.8.7
* /
// PEAR :: Math_Matrix 포함
// /usr/share/php/Math/Matrix.php
// include_path = ". : / usr / local / php / pear /"
require_once 'Math / Matrix.php';
require_once 'Math / Vector.php';
require_once 'Math / Vector3.php';
$ P1vector = 새로운 Math_Vector3 (배열 ($ xA, $ yA, $ zA));
$ P2vector = 새로운 Math_Vector3 (배열 ($ xB, $ yB, $ zB));
$ P3vector = 새로운 Math_Vector3 (배열 ($ xC, $ yC, $ zC));
#wikipedia에서 : http://en.wikipedia.org/wiki/Trilateration
원점에서 원 1을 얻는 #transform
x 축에서 원 2를 얻기위한 #transform
// CALC EX
$ P2minusP1 = Math_VectorOp :: substract ($ P2vector, $ P1vector);
$ 1 = 새로운 Math_Vector ($ P2minusP1);
$ P2minusP1_length = $ l-> 길이 ();
$ norm = 새로운 Math_Vector3 (배열 ($ P2minusP1_length, $ P2minusP1_length, $ P2minusP1_length));
$ d = $ norm; // calc D 저장
$ ex = Math_VectorOp :: divide ($ P2minusP1, $ norm);
// echo "ex :". $ ex-> toString (). "\ n";
$ ex_x = floatval ($ ex-> _ tuple-> getData () [0]);
$ ex_y = floatval ($ ex-> _ tuple-> getData () [1]);
$ ex_z = floatval ($ ex-> _ tuple-> getData () [2]);
$ ex = 새로운 Math_Vector3 (배열 ($ ex_x, $ ex_y, $ ex_z));
// CALC i
$ P3minusP1 = Math_VectorOp :: substract ($ P3vector, $ P1vector);
$ P3minusP1_x = floatval ($ P3minusP1-> _ tuple-> getData () [0]);
$ P3minusP1_y = floatval ($ P3minusP1-> _ tuple-> getData () [1]);
$ P3minusP1_z = floatval ($ P3minusP1-> _ tuple-> getData () [2]);
$ P3minusP1 = 새로운 Math_Vector3 (배열 ($ P3minusP1_x, $ P3minusP1_y, $ P3minusP1_z));
$ i = Math_VectorOp :: dotProduct ($ ex, $ P3minusP1);
// echo "i = $ i \ n";
// CALC 아이
$ iex = Math_VectorOp :: scale ($ i, $ ex);
// echo "iex =". $ iex-> toString (). "\ n";
$ P3P1iex = Math_VectorOp :: substract ($ P3minusP1, $ iex);
// echo "P3P1iex =". $ P3P1iex-> toString (). "\ n";
$ 1 = 새로운 Math_Vector ($ P3P1iex);
$ P3P1iex_length = $ l-> 길이 ();
$ norm = 새로운 Math_Vector3 (배열 ($ P3P1iex_length, $ P3P1iex_length, $ P3P1iex_length));
// echo "norm :". $ norm-> toString (). "\ n";
$ ey = Math_VectorOp :: divide ($ P3P1iex, $ norm);
// echo "ey =". $ ey-> toString (). "\ n";
$ ey_x = floatval ($ ey-> _ tuple-> getData () [0]);
$ ey_y = floatval ($ ey-> _ tuple-> getData () [1]);
$ ey_z = floatval ($ ey-> _ tuple-> getData () [2]);
$ ey = 새로운 Math_Vector3 (배열 ($ ey_x, $ ey_y, $ ey_z));
// CALC EZ
$ ez = Math_VectorOp :: crossProduct ($ ex, $ ey);
// echo "ez =". $ ez-> toString (). "\ n";
// CALC D
// 전에해라
$ d = floatval ($ d-> _ tuple-> getData () [0]);
// echo "d = $ d \ n";
// CALC J
$ j = Math_VectorOp :: dotProduct ($ ey, $ P3minusP1);
// echo "j = $ j \ n";
# 위키 백과에서
위의 값을 사용하는 #plug and chug
$ x = (pow ($ DistA, 2)-pow ($ DistB, 2) + pow ($ d, 2)) / (2 * $ d);
$ y = ((pow ($ DistA, 2)-pow ($ DistC, 2) + pow ($ i, 2) + pow ($ j, 2)) / (2 * $ j))-(($ i / $ j) * $ x);
# 여기에 표시된 하나의 사례 만
$ z = sqrt (pow ($ DistA, 2)-pow ($ x, 2)-pow ($ y, 2));
// echo "x = $ x-y = $ y-z = $ z \ n";
#triPt는 ECEF x, y, z의 평가 포인트가있는 배열입니다.
$ xex = Math_VectorOp :: scale ($ x, $ ex);
$ yey = Math_VectorOp :: scale ($ y, $ ey);
$ zez = Math_VectorOp :: scale ($ z, $ ez);
// CALC $ triPt = $ P1vector + $ xex + $ yey + $ zez;
$ triPt = Math_VectorOp :: add ($ P1vector, $ xex);
$ triPt = Math_VectorOp :: add ($ triPt, $ yey);
$ triPt = Math_VectorOp :: add ($ triPt, $ zez);
// echo "triPt =". $ triPt-> toString (). "\ n";
$ triPt_x = floatval ($ triPt-> _ tuple-> getData () [0]);
$ triPt_y = floatval ($ triPt-> _ tuple-> getData () [1]);
$ triPt_z = floatval ($ triPt-> _ tuple-> getData () [2]);
ECEF에서 위도 / 경도로 다시 변환
#도 단위로 변환
$ lat = rad2deg (asin ($ triPt_z / $ earthR));
$ lon = rad2deg (atan2 ($ triPt_y, $ triPt_x));
echo $ lat. ','. $ lon;