GDAL修复Landsat ETM+影像条带
Landsat7 ETM+卫星影像由于卫星传感器故障,导致此后获取的影像出现了条带。如下图所示, 影像中均匀的布满条带。
使用GDAL修复影像条带的代码如下:
def gdal_repair(tif_name, out_name, bands):
"""
tif_name(string): 源影像名
out_name(string): 输出影像名
bands(integer): 影像波段数
"""
# 打开影像文件
tif = gdal.Open(tif_name)
# 根据文件类型获取对应的驱动程序
driver = gdal.GetDriverByName('GTiff')
# 根据指定文件的驱动程序,使用现有数据集创建新的可写数据集
# 所有支持创建新文件的驱动程序都支持该`CreateCopy()`方法, # 但仅`Create()`部分支持该方法
# CreateCopy的第一个参数为目标文件名,第二个参数为源数据集
# 第三个参数的值是`0`或`1`,值是`0`。即使无法将原始数据准确地转换为目标数据,程序仍将执行
new_img = driver.CreateCopy(out_name, tif, 0)