Rimozione del bordo esterno delle immagini satellitari Landsat utilizzando il file stac
Raccontare storie con le immagini satellitari è semplice. I paesaggi affascinanti fanno la maggior parte del lavoro. Tuttavia, visualizzarli richiede del lavoro come selezionando e ridimensionando i canali RGB. In questo articolo andremo oltre. Vedremo come sbarazzarci di quel brutto riquadro di delimitazione. Nello specifico, provvederemo a:
- Ritaglia e ruota le scene Landsat utilizzando Essere file
- Discutere su come possiamo farlo ma anche mantenere la geolocalizzazione dei pixel
Discuteremo le parti chiave del codice Python e potrai trovare il progetto completo su GitHub.
Iniziamo scaricando una scena Landsat. Puoi farlo usando il file EarthExplorer portale. In alternativa, se desideri utilizzare Python, l’articolo seguente ti guida attraverso il processo:
Alla fine, dovresti avere una cartella come Figura 1. Questi sono tutti i file disponibili per a Prodotto scientifico Landsat livello 2. Lavoreremo con i file evidenziati. Queste sono le 3 bande di luce visibile e la SR_stac file.
Questa scena particolare è stata scattata sopra Città del Capo, in Sud Africa. Per vedere questo visualizziamo le bande di luce visibile utilizzando il get_rgb funzione. Questo prende il nome/ID del file come parametro. Quindi caricherà le fasce (righe 8–10), le impilerà (riga 13), ridimensionarli (riga 14) e agganciarli (riga 17).
import tifffile as tiff
import numpy as np
data_file = "./data/"def get_rgb(ID):
# Load Blue (B2), Green (B3) and Red (B4) bands
R = tiff.imread(data_file +'{}/{}_SR_B4.TIF'.format(ID, ID))
G = tiff.imread(data_file +'{}/{}_SR_B3.TIF'.format(ID…
Fonte: towardsdatascience.com
