Stereographic projection

##
Stereographic projection

Up to the release jul23, the GILDAS implementation was:

|  |  | (3)
|  |  | (4)
|  |  | (5)
|  |  | (6)
|  |  | (7)

where are the celestial coordinates and the projection plane coordinates. is an intermediate quantity. is the distance from the reference point in the projection plane, and the inclination angle (not detailed here). After a few computations one can show that the angle is compatible with the AIPS definition (AIPS memo 46 Eq. 3):

| (8)

But, on the other hand, the coordinate modulus in the AIPS convention is (AIPS memo 46 section 2.1.2):

| (9)

which can be written as

| (10)

thanks to traditional trigonometric identities 6, i.e. twice the modulus used in GILDAS. In other words, the spatial scaling in the GILDAS and AIPS conventions differ by a factor 2. Of these two definitions, the GILDAS convention was the unsatisfying one because angles from the reference point are projected as twice less distance on the projection plane ( for small angles, i.e. for the GILDAS definition, while one would prefer for correct distance measurement in the projected plane). This has been considered as a bug in the implementation and this was fixed with the appropriate factor in the aug23 GILDAS release. Doing so, the stereographic projection implementation in GILDAS is now fully compatible with the AIPS one, and thus to the CG02 one.