def project(x, y, fromproj=KnownCoordinateSystems.geographic.world.WGS1984, toproj=KnownCoordinateSystems.geographic.world.WGS1984):
"""
Project geographic coordinates from one projection to another.
:param x: (*array_like*) X coordinate values for projection.
:param y: (*array_like*) Y coordinate values for projection.
:param fromproj: (*ProjectionInfo*) From projection. Default is longlat projection.
:param toproj: (*ProjectionInfo*) To projection. Default is longlat projection.
:returns: (*array_like*, *array_like*) Projected geographic coordinates.
"""
if isinstance(fromproj, str):
fromproj = ProjectionInfo.factory(fromproj)
if isinstance(toproj, str):
toproj = ProjectionInfo.factory(toproj)
if isinstance(x, (tuple, list)):
x = np.array(x)
if isinstance(y, (tuple, list)):
y = np.array(y)
if isinstance(x, NDArray):
outxy = Reproject.reproject(x.asarray(), y.asarray(), fromproj, toproj)
return NDArray(outxy[0]), NDArray(outxy[1])
else:
inpt = PointD(x, y)
outpt = Reproject.reprojectPoint(inpt, fromproj, toproj)
return outpt.X, outpt.Y
def reproject(a, x=None, y=None, fromproj=None, xp=None, yp=None, toproj=None, method='bilinear'):
"""
Project array
:param a: (*array_like*) Input array.
:param x: (*array_like*) Input x coordinates.
:param y: (*array_like*) Input y coordinates.
:param fromproj: (*ProjectionInfo*) Input projection.
:param xp: (*array_like*) Projected x coordinates.
:param yp: (*array_like*) Projected y coordinates.
:param toproj: (*ProjectionInfo*) Output projection.
:param method: Interpolation method: ``bilinear`` or ``neareast`` .
:returns: (*NDArray*) Projected array
"""
if x is None or y is None:
if isinstance(a, DimArray):
y = a.dimvalue(a.ndim - 2)
x = a.dimvalue(a.ndim - 1)
else:
raise ValueError('Input x/y coordinates are None')
if fromproj is None:
if isinstance(a, DimArray):
fromproj = a.proj
else:
fromproj = KnownCoordinateSystems.geographic.world.WGS1984
if toproj is None:
toproj = KnownCoordinateSystems.geographic.world.WGS1984
if method == 'bilinear':
method = ResampleMethods.Bilinear
else:
method = ResampleMethods.NearestNeighbor
if xp is None or yp is None:
pr = Reproject.reproject(a.asarray(), x.aslist(), y.aslist(), fromproj, toproj, method)
return NDArray(pr[0]), NDArray(pr[1]), NDArray(pr[2])
if isinstance(xp, (list, tuple)):
xp = NDArray(xp)
if isinstance(yp, (list, tuple)):
yp = NDArray(yp)
xp, yp = ArrayUtil.meshgrid(xp.asarray(), yp.asarray())
r = Reproject.reproject(a.asarray(), x.aslist(), y.aslist(), xp, yp, fromproj, toproj, method)
return NDArray(r)