擅长:python、mysql、java
<p>下面的脚本根据4波段NAIP图像计算NDVI,其中波段4=nIR,波段3=Red。您需要空间分析员扩展。</p>
<p>请记住,陆地卫星TM波段4=近红外波段3=红色,陆地卫星8波段5=近红外波段4=红色。<a href="http://landsat.usgs.gov/band_designations_landsat_satellites.php" rel="nofollow">USGS Reference</a></p>
<pre><code># Calculates NDVI from multispectral imagery
import arcpy, string
from arcpy import env
from arcpy.sa import*
arcpy.CheckOutExtension("spatial")
env.workspace = r'C:\Your\workspace'
input = r'C:\Your\raster.tif'
result = "outputName.tif"
# You may need to change the band combinations.
# This is for 4-band NAIP imagery or Landsat TM.
NIR = input + "\Band_4"
Red = input + "\Band_3"
NIR_out = "NIR.tif"
Red_out = "Red.tif"
arcpy.CopyRaster_management(NIR,NIR_out)
arcpy.CopyRaster_management(Red, Red_out)
Num = arcpy.sa.Float(Raster(NIR_out) - Raster(Red_out))
Denom = arcpy.sa.Float(Raster(NIR_out) + Raster(Red_out))
NIR_eq = arcpy.sa.Divide(Num, Denom)
NIR_eq.save(result)
</code></pre>