Guess the projection from a point with unknown CRS - akin to https://stevage.github.io/WhatTheProj
Install via ]add https://github.com/JuliaGeo/WhatTheProj.jl
.
julia> using WhatTheProj;
julia> cache = fill!(EPSGGuesserCache(#= crs =#));
julia> berlin_latlong = (13.4050, 52.5200)
(13.405, 52.52)
julia> berlin_webmerc = (1.492238e6, 6.8947e6)
(1.492238e6, 6.8947e6)
julia> epsg_codes, distances = WhatTheProj.guess(cache, berlin_webmerc #= unknown =#, berlin_latlong #= known =#)
([3785, 3857, 7362, 6820, 6822, 8316, 8123, 7572, 3205, 3256], [10.694386392423878, 10.694386392423878, 257337.7611506837, 426271.3327776441, 426271.3327776441, 464075.3776994989, 536543.033609665, 536701.7777924604, 678768.517736908, 708965.1169710964])