DCL:MATH1:MAPLIB : 地図投影変換:サブルーチンの説明
正距円錐図法の投影, 逆投影を行なう.
CALL MPFCON(XLON,YLAT,X,Y)
CALL MPICON(X,Y,XLON,YLAT)
CALL MPSCON(YLAT1)
XLON,YLAT (R) 経度, 緯度. X,Y (R) x, y座標. YLAT1 (R) 標準緯線の緯度
x = r sin(kλ)
y = -r cos(kλ)
ここで, r, k は
k = cos(θ1)
r = θ-θ1 + tan(θ1)
である (θは余緯度, θ1は標準緯線の余緯度).