创建全色锐化的栅格数据集 (数据管理)
用法
-
此工具只能应用于完全重叠的区域。
-
可将输出保存为 BIL、BIP、BMP、BSQ、DAT、GIF、GRID、IMG、JPEG、JPEG 2000、PNG、TIFF 格式或任意地理数据库栅格数据集。
-
IHS 使用亮度、色调和饱和度颜色空间进行数据融合。Brovey 和 ESRI 使用基于光谱建模的算法进行数据融合。
-
“简单均值”方法使用原始的红色、绿色、蓝色值与全色像素值之间的平均像素值。对每种颜色（红色、绿色和蓝色）都执行此操作。
NewRedValue = ((OldRedValue + PanValue) / 2)
NewGreenValue = ((OldGreenValue + PanValue) / 2)
NewBlueValue = ((OldBlueValue + PanValue) / 2)
-
这四个权重值可用于调整全色锐化算法。所有权重相加后总和应等于 1。
-
在三波段栅格数据集中执行全色锐化将生成具有三个波段的栅格数据集。
在四波段栅格数据集中执行全色锐化将生成具有四个波段（三个已全色锐化的波段以及红外波段）的栅格数据集。
语法
| 参数 | 说明 | 数据类型 |
in_raster |
要全色锐化的输入栅格数据集。 这通常是低分辨率的多波段栅格数据集。 | Raster Layer |
red_channel （可选） |
要用红色显示的输入栅格波段。 | Long |
green_channel （可选） |
要用绿色显示的输入栅格波段。 | Long |
blue_channel （可选） |
要用蓝色显示的输入栅格波段。 | Long |
infrared_channel （可选） |
输入栅格的红外波段。并不是每个栅格数据集中都存在红外波段。 | Long |
out_raster_dataset |
输出栅格数据集。 以文件格式存储栅格数据集时，需要指定文件扩展名，具体如下：
以地理数据库形式存储栅格数据集时，不应向栅格数据集的名称添加文件扩展名。 将栅格数据集存储到 JPEG 文件、JPEG 2000 文件、TIFF 文件或地理数据库时，可以指定压缩类型和压缩质量。 | Raster Dataset |
in_panchromatic_image |
全色栅格数据集。这通常是高分辨率的灰度影像。 | Raster Dataset |
pansharpening_type |
全色锐化方法：
 旧版本：较旧的脚本和模型可能仍使用旧关键字。在 ArcGIS 10.0 版本中，此关键字被替换为“Simple_Mean”。这两个关键字都可继续使用，但为了清楚起见，最好将其更新为新关键字。 | String |
red_weight （可选） |
红色波段的权重值。 所有权重值（红、绿、蓝和红外）的总和应等于 1。 | Double |
green_weight （可选） |
绿色波段的权重值。 所有权重值（红、绿、蓝和红外）的总和应等于 1。 | Double |
blue_weight （可选） |
蓝色波段的权重值。 所有权重值（红、绿、蓝和红外）的总和应等于 1。 | Double |
infrared_weight （可选） |
红外波段的权重值。仅在存在红外波段并且使用红外选项时，此参数才有效。 所有权重值（红、绿、蓝和红外）的总和应等于 1。 | Double |
代码示例
这是 CreatePansharpenedRasterDataset 工具的 Python 示例。
import arcpy
arcpy.CreatePansharpenedRasterDataset_management("c:/data/rgba.img","3","2","1",
"4","output\\rgba_pan.tif",
"pan4.img","ESRI","0.16",
"0.16","0.16","0.52")这是 CreatePansharpenedRasterDataset 工具的 Python 脚本示例。
##====================================
##Create Pan-sharpened Raster Dataset
##Usage: CreatePansharpenedRasterDataset in_raster red_channel green_channel blue_channel {infrared_channel}
## out_raster_dataset in_panchromatic_image ESRI | IHS | Brovey |
## Simple_Mean {red_weight} {green_weight} {blue_weight} {infrared_weight}
try:
import arcpy
arcpy.env.workspace = r"\\MyMachine\PrjWorkspace\RasGP"
##3 Band RGB Pansharpen with Brovey algorithm
arcpy.CreatePansharpenedRasterDataset_management("rgb.img","3","2","1","1","output\\rgb_pan.img","pan.img","Brovey")
##4 Band Infrared Image Pansharpen with Customized ESRI algorithm
arcpy.CreatePansharpenedRasterDataset_management("rgba.img","3","2","1","4","output\\rgba_pan.tif","pan4.img","ESRI","0.16","0.16","0.16","0.52")
except:
print "Create Pan-sharpened Raster Dataset example failed."
print arcpy.GetMessages()