下载 Landsat 场景
我们首先下载陆地卫星场景。您可以使用EarthExplorer门户来执行此操作。
数据下载后,您应该有一个下图所示的文件夹。这些是Landsat 2 级科学产品的所有可用文件。我们将处理突出显示的文件。这些是 3 个可见光波段和SR_stac文件。
加载图像和 stac 文件
这个特殊的场景是在南非开普敦上空拍摄的。为了看到这一点,我们使用get_rgb函数可视化可见光波段。这将文件名/ ID 作为参数。
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, ID))
B = tiff.imread(data_file +'{}/{}_SR_B2.TIF'.format(ID, ID))
# Stack and scale bands
RGB = np.dstack((R, G, B))
RGB = np.clip(RGB*0.0000275-0.2, 0, 1)
# Clip to enhance contrast
RGB = np.clip(RGB,0,0.2)/0.2
return RGB
我们使用此函数来获取图 1中下载的场景的 RGB 可视化(第 4-5 行)。您可以在图 2中看到生成的图像。如图所示,卫星图像通常在边界框内旋转。这是图像周围的黑色边框。
import matplotlib.pyplot as plt
# Get RGB image
ID = "LC08_L2SP_175083_20131218_20200912_02_T1"
RGB= get_rgb(ID)
# Display band
fig, ax = plt.subplot