Get geocode from ECEF
mainly based on reverse_geocode and reverse-geocoder.
using Pkg
Pkg.add("https://github.com/htyeim/HTrg.jl")
using Geodesy
using HTrg
p = LLA(22.2710, 113.5767)
e = ECEFfromLLA(wgs84)(p)
@show code, city, dis, ci = get_geocode(e);
(code, city, dis, ci) = get_geocode(e) = ("CN", "Zhuhai", 1130.3958719970858, ("CHN", 156, "China", "Beijing", 9.59696e6, 1.39273e9))