Skip to content

Instantly share code, notes, and snippets.

@aont
Created April 18, 2013 06:58
Show Gist options
  • Select an option

  • Save aont/5410729 to your computer and use it in GitHub Desktop.

Select an option

Save aont/5410729 to your computer and use it in GitHub Desktop.
日本測地系に変換
a=6377.39715500
e=sqrt(0.006674372230614))
h=0
phi=34
lambda=130
N=a/sqrt(1-e*e*sin(phi)*sin(phi))
X=(N+h)*cos(phi)*cos(lambda)
Y=(N+h)*cos(phi)*sin(lambda)
Z=(N*(1-e*e)+h)*sin(phi)
参考
http://vldb.gsi.go.jp/sokuchi/surveycalc/trans_alg/trans_alg.html
http://homepage1.nifty.com/aida/jr1huo_calsat32/Calsat32Keiido.htm
Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment