sharpetronics/Map-Projections
1
0
Fork 0
mirror of https://github.com/csharpee/Map-Projections.git synced 2025-12-15 00:00:14 -05:00
Code Issues Projects Releases Wiki Activity

Goodbye, cool curls

Browse Source 
I just can't justify using this much computing power on an initial guess
for a bisection scheme. For Newton-Raphon it might have made sense. I'm
really not sure why that didn't converge.
This commit is contained in:
Justin Kunimune 2019-10-16 08:41:51 -04:00
parent 5cf6bcd6e2
commit 5240f3863d
1 changed files with 1 additions and 3 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

4
src/maps/Lenticular.java
View File

@ -303,11 +303,9 @@ public class Lenticular {
else if (y == 0)
return new double[] {0, x};
double logitLatGuess = Math.log(Math.hypot(x, y+Math.PI/2)/Math.hypot(x, y-Math.PI/2));
double lat = NumericalAnalysis.bisectionFind(
(ph)->(Math.pow(x, 2) + Math.pow(y - ph, 2) - 2*(y - ph)/Math.tan(ph)),
0, Math.PI/2*(Math.exp(logitLatGuess) - 1)/(Math.exp(logitLatGuess) + 1), // this guess is pretty complicated, but I'm proud of it
1e-4);
0, Math.PI/2, 1e-4);
if (Double.isNaN(lat))
return null;
return new double[] { lat, Math.atan2(x, -(y - lat - 1/Math.tan(lat)))/Math.sin(lat) };