Trong công việc của mình, tôi đã thực hiện một số bảo trì trên phần mềm nội bộ tuyên bố thực hiện phép chiếu Stereographic có công thức được sao chép từ một hệ thống cũ từ lâu.
Từ những gì tôi đã thấy, có vẻ như triển khai "Xiên và Xích đạo" từ EPSG (mã 9809). Nhưng có một bước mà tôi không thể xác định được. Giá trị "R" được sử dụng trong mã của chúng tôi được tính toán với các tham số sau:
(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
Với điều này, công thức để có R là:
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
Có ai biết những gì đang xảy ra ở đây? Tôi đã tìm kiếm và đọc về các phép chiếu lập thể, nhưng tôi không thể tìm thấy bất cứ thứ gì như thế. Trên thực tế, tôi chưa thấy bất kỳ dự đoán nào sử dụng "tọa độ biến dạng null" này.