从褶皱的多光谱合成图像计算NDSSI

Question

我不熟悉ArcPy和Python。我一直在尝试使用python脚本从文件夹中的100个复合光栅图像计算NDSSI

这个脚本应该从一个工作区获取光栅列表，计算NDSSI，并将结果以.tif格式存储在另一个文件夹中。然而，它不起作用

# Calculates NDSSI from multispectral composite images from a folder

import arcpy, string

from arcpy import env
from arcpy.sa import*

arcpy.CheckOutExtension("spatial")

env.workspace = r'E:\Landsat All Images\Processing\LT5\Composite'

outws = r'E:\Landsat All Images\Processing\LT5\NDSSI'

rasters = arcpy.ListRasters("*.tif")

for raster in rasters:


    Blue = raster + "\Band_1"
    NIR = raster + "\Band_4"

    Num = arcpy.sa.Float(Raster(Blue) - Raster(NIR))
    Denom = arcpy.sa.Float(Raster(Blue) + Raster(NIR))
    NIR_eq = arcpy.sa.Divide(Num, Denom)

    NIR_eq.save(outws)
print "Processing complete"