楕円体上の距離で若干無茶をしてみた

GIS 保留

id:yellow_73:20080504#p3 の続き。スターありがとうございます。
ITRF3次元でやってみる。
Pn(n=1,2,...)に対してP0を与えます。
ただし、max(ST_distance_spheroid(Pn-P0)^2)<700[m] となっていました。

En=(d2(Pn,P0)-d1(Pn,P0))/d1(Pn,P0)
ただし
  d1(Pn,P0)=ST_distance_spheroid(Pn,P0)
  d2(Pn,P0)=sqrt*1;
    e2=f*(2.0-f)
    a=6378137, 1/f=298.257222101 (GSR1980)

なお、ST_distance_spheroid()以外は、いきがかり上:-) PHPで計算しました。
で、Enの分散は 2.41e-11 でした。ま、問題無いかなと。

*1:Xn-X0)^2+(Yn-Y0)^2+(Zn-Z0)^2) Xn=N0*cos(φn)*cos(λn) Yn=N0*cos(φn)*sin(λn) Zn=N0*(1.0-e2)*sin(φn) X0=N0*cos(φ0)*cos(λ0) Y0=N0*cos(φ0)*sin(λ0) Z0=N0*(1.0-e2)*sin(φ0) N0=a/sqrt(1.0-e2*sin^2(φ0