网站设计制作用软件,中英文企业网站制作,公司起名大全免费版,wordpress 文档预览在之前的分享中#xff0c;我们介绍了如何使用ArcPy将GEE下载的遥感影像转为单波段文件。基于前面创建的单波段文件#xff0c;我们可以一次性计算多种遥感指数#xff0c;例如NDVI、EVI、NDSI等。我这里直接在ArcGIS Pro中自带的Notebook进行的运行。如下图所示#xff0c…在之前的分享中我们介绍了如何使用ArcPy将GEE下载的遥感影像转为单波段文件。基于前面创建的单波段文件我们可以一次性计算多种遥感指数例如NDVI、EVI、NDSI等。我这里直接在ArcGIS Pro中自带的Notebook进行的运行。如下图所示我们可以直接在ArcGIS Pro菜单栏中分析选项卡下的Python新建Notebook也可以直接在Catalog中新建一个notebook文件。 新建.ipynb文档后将以下代码复制进文档可以将代码均放入一个单元格点击运行将会一次性计算所有的指数也可以将代码分别放入不同的单元格只计算需要用到的指数。需要注意的是这里需要导入arcpy.sa我没深入分析这个问题按说ArcGIS Pro自带的Notebook已经默认导入了ArcPy包但是我不导入arcpy.sa确实会报错。
# 导入sa模块用于Power
from arcpy.sa import *
#读取单波段数据并栅格化
inPath E:/输入的单波段影像路径/
outPath E:/输出结果指数路径/NIR arcpy.ia.Raster(inPathNIR.tif)
R arcpy.ia.Raster(inPathR.tif)
G arcpy.ia.Raster(inPathG.tif)
B arcpy.ia.Raster(inPathB.tif)# 计算NDVI
NDVI (NIR-R)/(NIRR)
NDVI.save(outPathNDVI.tif)
# 计算EVI
EVI 2.5*(NIR-R)/(NIRR*6-7.5*B1)
EVI.save(outPathEVI.tif)
# 计算IFe2O3
IFe2O3 R/NIR
IFe2O3.save(outPathIFe2O3.tif)
# 计算NDSI
NDSI (R-NIR)/(RNIR)
NDSI.save(outPathNDSI.tif)
# 计算corsi
CORSI (BG)/(RNIR)*(NIR-R)/(NIRR)
CORSI.save(outPathCORSI.tif)
# 计算BI
BI Power((Power(R,2)Power(NIR,2)),0.5)
BI.save(outPathBI.tif)
# 计算SI1
SI1 R*NIR/G
SI1.save(outPathSI1.tif)
#需要计算其他指数的话直接补充其他计算公式即可如下所示我将导入单波段的代码与栅格化的代码各放入了一个单元格中计算各个指数的代码也各在一个单元格中这样需要计算哪个指数时仅仅运行某一个指数即可。 更多分享见合集ArcGIS_Pro_Python
