Di tempat kerja saya, saya telah melakukan beberapa pemeliharaan pada perangkat lunak internal yang mengklaim untuk mengimplementasikan proyeksi Stereografis yang rumus-rumusnya disalin dari sistem lama sejak lama.
Dari apa yang saya lihat, tampaknya menerapkan "Oblique and Equatorial Stereographic" dari EPSG (kode 9809). Tetapi ada satu langkah yang tidak bisa saya identifikasi. Nilai "R" yang digunakan dalam kode kami dihitung dengan parameter berikut:
(phi_n, lambda_n) - null distortion coordinates
(PHI_nc, LAMBDA_nc) - conformal null distortion coordinates
(PHI_o, LAMBDA_o) - conformal projection origin coordinates
e - eccenttricity
r_eq - earth radius at equator
Mengingat ini, rumus untuk mendapatkan R adalah:
a = 1 + cos(PHI_nc) * cos(PHI_o) * cos(LAMBDA_nc - LAMBA_o)
b = sin(PHI_nc) * sin(PHI_o)
c = 2 * cos(PHI_nc) * [1 - (sin(phi_n) * e/2) ^ 2] ^ 1/2
R = r_eq * cos(phi_n) * (a + b) / c
Adakah yang tahu apa yang sedang terjadi di sini? Saya telah mencari dan membaca tentang proyeksi stereografi, tetapi saya tidak dapat menemukan hal seperti itu. Bahkan, saya belum melihat proyeksi yang menggunakan "koordinat nol distorsi" ini.