Now, I am trying to implement the inverse projection according to the description provided on page 14 of the above paper.
However, the problem is that the value for the
spherical distance, z, is not getting properly calculated (it is about 0.5-2 units “off” from the value that results after using the forward ISEA formula for the same lat,long point). The
azimuth Az, however, is getting properly calculated, and the only problem it seems is with the value of z. Therefore, the latitude and longitude of a point gets incorrectly calculated from the inverse, and is off by 50-100 kms from the true location.
I am attaching the main part of the Python code to this email, and it uses the same constants defined in the forward projection code.
If anyone can provide some guidance as to where I am going wrong, or what the problem could be, I would greatly appreciate it. If any additional information is
needed, please let me know.