【ArcGIS遇上Python】长时间序列(30年)每两组栅格数据对应做减法运算求物候参数
生活随笔
收集整理的這篇文章主要介紹了
【ArcGIS遇上Python】长时间序列(30年)每两组栅格数据对应做减法运算求物候参数
小編覺(jué)得挺不錯(cuò)的,現(xiàn)在分享給大家,幫大家做個(gè)參考.
Python語(yǔ)言在ArcGIS軟件中發(fā)揮著乾坤大挪移的作用,然而很多人并不知道它的厲害之處。今天我們利用Python語(yǔ)言來(lái)批處理遙感物候數(shù)據(jù),獲取物候期的參數(shù)。
說(shuō)明:
(1)Up_Path:始期數(shù)據(jù)路徑變量
(2)Down_Path:末期數(shù)據(jù)路徑變量
(3)outPath:相減輸出結(jié)果路徑變量
完整的Python代碼如下:
import arcpy arcpy.CheckOutExtension("spatial") arcpy.gp.overwriteOutput=1Up_Path="G:\\Phenology of 30 Years\\GIMMS 3g\\13Up\\2PointToRaster\\" Down_Path="G:\\Phenology of 30 Years\\GIMMS 3g\\14Down\\2PointToRaster\\" OutPath="G:\\Phenology of 30 Years\\GIMMS 3g\\15Length\\1Length\\" n=1983for i in range(0,30):OutName=OutPath+str(n)a=Up_Path+"p2r_"+str(n)b=Down_Path+"p2r_"+str(n)arcpy.gp.Minus_sa(b, a, OutName) #柵格相減,其實(shí)是字符串對(duì)應(yīng)的數(shù)據(jù)相減,但是不能直接用“-”來(lái)操作字符串print(str(n)+" has done") 創(chuàng)作挑戰(zhàn)賽新人創(chuàng)作獎(jiǎng)勵(lì)來(lái)咯,堅(jiān)持創(chuàng)作打卡瓜分現(xiàn)金大獎(jiǎng)總結(jié)
以上是生活随笔為你收集整理的【ArcGIS遇上Python】长时间序列(30年)每两组栅格数据对应做减法运算求物候参数的全部?jī)?nèi)容,希望文章能夠幫你解決所遇到的問(wèn)題。
- 上一篇: Java Eclipse开发环境搭建及注
- 下一篇: 【ArcGIS遇上Python】Pyth