计算全色锐化权重 (数据管理)

描述

为新的或自定义的传感器数据计算一组最佳的全色锐化权重。

使用方法

• 此工具将计算一组最佳的全色锐化权重，这些值可用于需要全色锐化权重的其他工具。

• 如果将栅格产品用作输入栅格，则应用栅格产品模板中的波段顺序。

语法

ComputePansharpenWeights(in_raster, in_panchromatic_image, {band_indexes})

参数	说明	数据类型
in_raster	具有全色波段的多光谱栅格。	Mosaic Dataset; Mosaic Layer; Raster Dataset; Raster Layer
in_panchromatic_image	与多光谱栅格相关联的全色波段。	Raster Layer
band_indexes （可选）	全色锐化权重的波段顺序。 如果将栅格产品用作 in_raster 参数，则应用栅格产品模板中的波段顺序。	String

派生输出

名称	说明	数据类型
out_string	输出全色锐化权重。	字符串

代码示例

ComputePansharpenWeights 示例 1（Python 窗口）

这是 ComputePansharpenWeights 工具的 Python 示例。

import arcpy
arcpy.ComputePansharpenWeights_management(
    "c:/data/rgb.tif", "c:/data/image.tif", "3 2 1 4")

ComputePansharpenWeights 示例 2（独立脚本）

这是 ComputePansharpenWeights 工具的 Python 脚本示例。

#Run Compute Pan Sharpen Weights tool using the bands 4,3,2,1 

import arcpy

InMSraster = "C:\\Landsat7\\L71046029_02920050705_MTL.txt\Multispectral" 
InPANraster = "C:\\Landsat7\\L71046029_02920050705_MTL.txt\Panchromatic"
band_index = "3 2 1 5"

arcpy.ComputePansharpenWeights_management(InMSraster, InPANraster, band_index)

ComputePansharpenWeights 示例 3（工作流脚本）

这是在其他工具中使用 ComputePansharpenWeights 输出的 Python 脚本示例。

#Compute the pansharpening weights and use the results in the 
#create pansharpening tool.

try:
    import arcpy
    
    InRGBraster = "C:\\temp\\rgb.img"
    InPanraster = "C:\\temp\\pan.tif"
    
    #Compute Pan Sharpen Weights  
    out_pan_weight = arcpy.ComputePansharpenWeights_management(
        InRGBraster, InPanraster, "3 2 1 4")
    
    #Get results 
    pansharpen_weights = out_pan_weight.getOutput(0)
    
    #Split the results string for weights of each band
    pansplit = pansharpen_weights.split(";")
    
    #Run the Create pan sharpened raster dataset tool. 
    arcpy.CreatePansharpenedRasterDataset_management(
        InRGBraster, "3", "2", "1", "4", "C:\\temp\\pansharpened_raster.tif",
        InPanraster, "Gram-Schmidt", pansplit[0].split(" ")[1],  
        pansplit[1].split(" ")[1], pansplit[2].split(" ")[1],
        pansplit[3].split(" ")[1])
    
except arcpy.ExecuteError:
    print(arcpy.GetMessages())
except Exception as err:
    print(err[0])

环境

此工具不使用任何地理处理环境。

许可信息

• Basic: 是
• Standard: 是
• Advanced: 是