Interpolate topography to a plane

[albert-de-montserrat]

InterpolateTopographyOnPlane interpolates the topopgraphy of a GeoData object to a plane defined by x::AbstractRange and y::AbstractRange. Example:

LonLat = load("ExportToba_2.jld2", "TobaTopography_LATLON")
proj   = ProjectionPoint(; Lat = 2.19, Lon = 98.91)

surf   = interpolate_topography2plane(LonLat, proj, x, y)
120×120 Matrix{Float64}:
   0.0      0.0      0.0    …  273.912  517.351   560.975
   0.0      0.0      0.0       461.027  540.289   772.906
   0.0      0.0      0.0       589.679  844.387  1026.88
   0.0      0.0      0.0       938.434  986.609  1065.11
   0.0      0.0      0.0       805.56   908.139   919.15
   0.0      0.0      0.0    …  790.199  793.64    761.984
   ⋮                        ⋱                    
 205.9    179.577  236.505       0.0      0.0       0.0
 188.478  163.321  239.682  …    0.0      0.0       0.0
 178.089  163.979  214.634       0.0      0.0       0.0
 159.713  171.233  184.283       0.0      0.0       0.0
 164.447  172.294  205.19        0.0      0.0       0.0
 153.36   187.248  219.04        0.0      0.0       0.0