卫星影像密集匹配之GDAL预处理
参考:https://gdal.org/index.html
主要步骤:
- PanSharpen
- Translate
- Gdaladdo
一、PanSharpen
全色锐化(PanSharpen)是合并高分辨率全色和低分辨率多光谱图像以创建单个高分辨率彩色图像的过程。使用高分辨率灰度波段中的空间信息和多光谱波段中的颜色信息来创建高分辨率彩色图像,从本质上提高数据集中颜色信息的分辨率以匹配全色波段的分辨率。
具体为:低分辨率色带 + 高分辨率灰度带 = 高分辨率彩色图像
作用:高分辨率彩色图像相较于高分辨率全色影像具有更丰富的图像纹理信息,在进行匹配时能够更好的提取同名点。大大提高影像的重建效果。
1
2
3
4
set ms_dir="IMG_BJ3N1_202205090240444_MS_SEN_PWOI_000020136_1_3_F_1_R1C1.TIF"
set pan_dir="IMG_BJ3N1_202205090240444_PAN_SEN_PWOI_000020136_1_3_F_1_P_R1C1.TIF"
set sharpen_name="IMG_BJ3N1_202205090240444_PANSHARPENED_SEN_PWOI_000020136_1_3_F_1_P_R1C1.TIF"
gdal_pansharpen %pan_dir% %ms_dir% -b 3 -b 2 -b 1 -r cubic -threads ALL_CPUS -co "RPCTXT=YES" -co "TILED=YES" %sharpen_name%
二、Translate
The gdal_translate utility can be used to convert raster data between different formats, potentially performing some operations like subsettings, resampling, and rescaling pixels in the process.
简而言之即缩放图像到Uint8,转换rgb通道等。
作用:方便软件处理
1
2
set sharpen_name_byte="IMG_BJ3N1_202205090240444_PANSHARPENED_SEN_PWOI_000020136_1_3_F_1_P_R1C1.TIF"
gdal_translate -ot Byte -scale 0 1000 0 255 -b 1 -b 2 -b 3 -co "TILED=YES" %sharpen_name% %sharpen_name_byte%
三、Gdaladdo
作用:生成影像概览。
1
gdaladdo -ro %sharpen_name_byte%