PHP 大地坐标转空间直角坐标

公式就不写了，直上代码

<?php

  $a = 6378137;            //长半轴
  $f = 1/298.257223563;    //WGS84 扁率
  $b = $a * (1 - $f);
  $e = sqrt($a * $a - $b * $b) / $a;   //椭球第一偏心率
  
  $lat = 25;
  $lon = 117;
  $h = 100;
  
  $N = $a / sqrt(1 - $e * $e * sin($lat * M_PI / 180) * sin($lat * M_PI / 180));
  $WGS84_X = ($N + $h) * cos($lat * M_PI / 180) * cos($lon * M_PI / 180);
  $WGS84_Y = ($N + $h) * cos($lat * M_PI / 180) * sin($lon * M_PI / 180);
  $WGS84_Z = ($N * (1 - ($e * $e)) + $h) * sin($lat * M_PI / 180);
 
?>