geoutils.Raster.xy2ij

Raster.xy2ij(x, y, op=<class 'numpy.float32'>, precision=None, shift_area_or_point=None)

Get indexes (row,column) of coordinates (x,y).

By default, the indexes are shifted with the interpretation of pixel coordinates “AREA_OR_POINT” of the raster, to ensure that the indexes of points represent the right location. See parameter description of shift_area_or_point for more details.

This function is reversible with ij2xy for any pixel interpretation.

Parameters:

• x (Union[_SupportsArray[dtype[Any]], _NestedSequence[_SupportsArray[dtype[Any]]], bool, int, float, complex, str, bytes, _NestedSequence[Union[bool, int, float, complex, str, bytes]]]) – X coordinates.

• y (Union[_SupportsArray[dtype[Any]], _NestedSequence[_SupportsArray[dtype[Any]]], bool, int, float, complex, str, bytes, _NestedSequence[Union[bool, int, float, complex, str, bytes]]]) – Y coordinates.

• op (type) – Operator to compute index.

• precision (float | None) – Precision passed to rasterio.transform.rowcol().

• shift_area_or_point (bool | None) – Whether to shift with pixel interpretation, which shifts to center of pixel indexes if self.area_or_point is “Point” and maintains corner pixel indexes if it is “Area” or None. Defaults to True. Can be configured with the global setting geoutils.config[“shift_area_or_point”].

Returns i, j:

Indices of (x,y) in the image.

Return type:

tuple[ndarray[Any, dtype[Union[floating[Any], integer[Any]]]], ndarray[Any, dtype[Union[floating[Any], integer[Any]]]]]