Python 包围盒裁剪卫星场景

下载 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.subplo