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)
である.