I have satellite image that I want to convert into numpy array (and then to Pandas). I already know how to do that, but the problem is that it does not preserve the coordinate data.

This is how I do it now:

import xarray
array = src.read()


I have also tried to follow this answer (keeping the coordinate system of raster files in the resulting raster file after operation with numpy) but I always get this error:

NameError: name 'gdal_array' is not defined

which does not allow me to open the image.

My end goal is to have Pandas table that contains the bands values together with the coordinate values.


2 Answers 2


You can use the open_rasterio and to_dataframe methods to accomplish that.

import rioxarray

rds = rioxarray.open_rasterio("file.tif")

See also: https://gis.stackexchange.com/a/358057/144357


I would read with gdal and dump that to a numpy array. ReadAsArray().

Coordinates and projection are obtained using the GetGeotransform() SetGeotransform() and the GetProjection() and SetProjection() applied to the data set.

Example here. https://pcjericks.github.io/py-gdalogr-cookbook/raster_layers.html

Your Answer

By clicking “Post Your Answer”, you agree to our terms of service and acknowledge that you have read and understand our privacy policy and code of conduct.

Not the answer you're looking for? Browse other questions tagged or ask your own question.