Strane coordinate dopo la conversione dell'array numpy in raster

Jan 03 2021

Ho l'array NumPy e le sue coordinate del riquadro di delimitazione. Ho provato a convertirlo in raster usando rasterio, sulla base di questa risposta , e l'ho salvato come raster, ma quando uso rasterio.show le coordinate sono molto sbagliate.

Questo è lo script che ho usato:

bbox_coords_wgs84=[-101.7359960059834, 20.21904081937658, -100.5717967351885, 20.8312118894487]

#variables for the projection:
minx=bbox_coords_wgs84[0]
maxy=bbox_coords_wgs84[3]
pixel_size= 10

#according to the post on GIS SO:

import rasterio
from rasterio.transform import from_origin

transform=from_origin(minx,maxy,pixel_size,pixel_size)
crs_img='EPSG:4326'


with rasterio.open('test1.tif', 
                    'w',
                    driver='GTiff',
                    height=ndvi.shape[0],
                    width=ndvi.shape[1],
                    count=1,
                    dtype=ndvi.dtype,
                    crs=crs_img,
                    nodata=None, # change if data has nodata value
                    transform=transform) as dst:
        dst.write(ndvi, 1)
 
#display the results:

from matplotlib import pyplot
from rasterio.plot import show

src = rasterio.open('test1.tif')
show(src)

Come puoi vedere, i numeri non sono assolutamente le coordinate corrette.

Il mio obiettivo finale: essere in grado di riproiettare correttamente l'array NumPy in WGS84.

* Questo post si riferisce anche a questo post

Risposte

2 gene Jan 03 2021 at 21:44

Non si segnala che questa è una suite di riproiettare un array NumPy con trasformazione affine dove si utilizzarasterio.transform.from_bounds

Dal modulo rasterio.transform

rasterio.transform.from_bounds (ovest, sud, est, nord, larghezza, altezza)
Restituisce una trasformazione Affine dati limiti, larghezza e altezza.
Restituisce una trasformazione Affine per un raster georeferenziato dati i suoi limiti ovest, sud, est, nord e la sua larghezza e altezza in numero di pixel.

E

rasterio.transform.from_origin (west, north, xsize, ysize)
Restituisce una trasformazione Affine data in alto a sinistra e le dimensioni dei pixel.
Restituisce una trasformazione Affine per un raster georeferenziato date le coordinate del suo angolo superiore sinistro ovest, nord e le dimensioni dei pixel xsize, ysize.

Non è la stessa cosa ei risultati sono diversi

rasterio.transform.from_bounds( -101.7359960059834,20.21904081937658,-100.5717967351885,20.8312118894487,1103,2039)
Affine(0.0010554843796871222, 0.0, -101.7359960059834,
   0.0, -0.0003002310299519955, 20.8312118894487)

rasterio.transform.from_origin(-101.7359960059834,20.8312118894487,10,10)
Affine(10.0, 0.0, -101.7359960059834,
   0.0, -10.0, 20.8312118894487)

Nuovo

I quattro angoli del raster dal limite (larghezza = 1103, altezza = 2039)

fig,ax = plt.subplots()
ax.plot(0,0,'ro')
ax.plot(1103,0,'bo')
ax.plot(0,2039,'go')
ax.plot(1103,2039,'co')
plt.show()

La trasformazione

 trans = rasterio.transform.from_bounds(-101.7359960059834,20.21904081937658-100.5717967351885,20.8312118894487,1103,2039)

 
trans*(0,0)
(-101.7359960059834, 20.8312118894487)
trans*(1103,0) 
(-100.5717967351885, 20.8312118894487)
trans*(0,2039) 
(-101.7359960059834, 20.21904081937658)
trans*(1103,2039)
(-100.5717967351885, 20.21904081937658)

fig,ax = plt.subplots()
ax.plot(*(trans*(0,0)),'ro')
ax.plot(*(trans*(1103,0)),'bo')
ax.plot(*(trans*(0,2039)),'go')
ax.plot(*(trans*(1103,2039)),'co')
plt.show()