Rotazione di raster con Rasterio

Rotazione raster con Rasterio

Utilizzando Python per ruotare le immagini satellitari mantenendo l’accuratezza della geolocalizzazione

I dati raster sono simili ai dati delle normali immagini. Tuttavia, ogni pixel è collegato a una posizione sulla superficie terrestre. Questo complica le cose. Se vogliamo ruotare i dati, dobbiamo anche considerare il sistema di riferimento delle coordinate (CRS) sottostante. Distorcere un raster senza adeguare la sua geolocalizzazione porterà ad un’analisi spaziale inaccurata.

L’aggiustamento della geolocalizzazione non è semplice. Fortunatamente, Rasterio può aiutare. È una popolare libreria Python per l’analisi dei dati geospaziali. Utilizzeremo il pacchetto per:

• ruotare un raster
• e riproiettare l’immagine nel sistema di riferimento delle coordinate corretto.

Durante il processo, discuteremo il codice Python e potrai trovare il progetto completo su GitHub.

Questo articolo assume una conoscenza di base dei dati raster e del lavoro con il loro CRS. Dai un’occhiata all’articolo qui sotto se vuoi ripassare. Approfondisce maggiormente la riproiezione dei dati raster.

• Come imparerei Data Science con ChatGPT (se potessi ricominciare)
• Replicazione dei dati CDC tecniche, compromessi, intuizioni
• Impatto di ChatGPT sui copywriter difficoltà e speranze per un futuro integrato con l’IA

Come tracciare le coordinate sulle immagini satellitari Landsat con Python

Mappare le posizioni dei pixel alle coordinate geografiche utilizzando i metadati Landsat e Rasterio

towardsdatascience.com

Scaricare le scene Landsat

Per i nostri dati raster, lavoreremo con immagini satellitari. Nello specifico, scene Landsat. Puoi scaricarne una utilizzando il portale EarthExplorer. In alternativa, se vuoi utilizzare Python, l’articolo qui sotto ti guiderà nel processo:

Scaricare immagini satellitari Landsat con Python

Ottimizzare i download delle scene Landsat con il pacchetto Python landsatxplore

towardsdatascience.com

Alla fine, dovresti avere una cartella contenente tutti i file disponibili per un prodotto scientifico di livello 2 Landsat. Lavoreremo con la banda del rosso della luce visibile. Per le scene Landsat 8 o 9, questa è indicata dalla banda B4.

Apertura di un file raster

Utilizziamo il codice qui sotto per aprire e visualizzare questa banda. ID fornisce l’ID della scena Landsat per questa particolare scena…