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)