I have downloaded multiple products from EO-1 satellite and I have created digital rasters out of them. Each of them contains 10 or 242 bands (those are hyperspectral). I have also created a shapefile for each of them with multiple points. Now I would like to create a table/array with Python containing values of a pixel from each band with coordinates from a shapefile. I would like to have all bands in separate columns and all points in separate rows. I have managed to create such data frame in R, and I tried to do so with GDAL and numpy but I am far less familiar with Python.

  • Hey, welcome to GIS SE. Can you provide your efforts (i.e a code sample of what you have tried using Python)? Furthermore, check this answer I wrote sometime ago. It might be a good starting point for retrieving pixel values from xy coordinates. – Marcelo Villa Jul 3 '19 at 16:48
  • I think, it has already been answered here and here. – pan Jul 3 '19 at 17:02

Your Answer

By clicking “Post Your Answer”, you agree to our terms of service, privacy policy and cookie policy

Browse other questions tagged or ask your own question.