DCL:MATH1:MAPLIB : 地図投影変換:メソッド(元サブルーチン)の説明
ランベルト正積方位図法の投影, 逆投影を行なう.
CALL MPFAZA(XLON,YLAT,X,Y)
CALL MPIAZA(X,Y,XLON,YLAT)
xlon,ylat (R) 経度, 緯度. x,y (R) x, y座標.
x = r sin(λ)
y = -r cos(λ)
ここで, r は
r = 2sin(θ/2)
である.