Jasper Yu says to YSITD function lat_to_point(lat){ sin = Math.sin(lat * 0.0174532925); d0 = Math.log((sin + 1) / (1 - sin)) * -0.5 d1 = 667544.214 d8 = 2097152 return d8 + d1 * d0 }