ARCGIS10柵格計算PYTHON代碼

ARCGIS10柵格計算PYTHON代碼使用python語言進行柵格計算1重采樣通過更改像元大小和重采樣方法來更改柵格數(shù)據(jù)集用法：1）可更改像元大小，但柵格數(shù)據(jù)集的范圍將保持不變。2）本工具僅能輸出方形像元大小。3）可將輸出保存為BIL、BIP、BMP、BSQ、DAT、GIF、GRID、IMG、JPEG、JPEG2000、PNG、TIFF格式或任意地理數(shù)據(jù)庫柵格數(shù)據(jù)集。4）有四個用于重采樣技術(shù)參數(shù)的選項：“最鄰近”選項，用于執(zhí)行最鄰近分配法，是速度最快的插值方法。此選項主要用于離散數(shù)據(jù)（如土地使用分類），因為它不會更改像元的值。最大空間誤差將是像元大小的一半。與最鄰近法一樣，此選項主要用于離散數(shù)據(jù)；但與“最鄰近”選項相比，“眾數(shù)”選項通?？缮筛交慕Y(jié)果。16個最鄰近輸入像元中心的平滑曲線確定像的值。與通過運行最鄰近重采樣算法獲得的柵格相比，輸出柵格的幾何變形程度較小。輸入像元值范圍之外。如果無法接受此結(jié)果，請轉(zhuǎn)而使用“雙線性”選項?！半p線性”或“三次”選項不得用于分類數(shù)據(jù)，因為像元值可能被更改。5）置。語法：Resample_management(in_raster,out_raster,{cell_size},{resampling_type})實例：#Resampletry:importarcpyarcpy.env.workspace=r\##ResampleTIFFimagetoahigherresolutionarcpy.Resample_management(\\#{NEAREST最鄰近|BILINEAR雙線性|CUBIC三次|MAJORITY眾數(shù)}except:print\printarcpy.GetMessages()2加逐個像元地將兩個柵格的值相加（求和）語法：Plus_3d(in_raster_or_constant1,in_raster_or_constant2,out_raster)實例：#Name:jia.py#Description:Addsthevaluesoftworastersonacell-by-cellbasis.#Requirements:3DAnalystExtension#Author:ESRI#Importsystemmodulesimportarcpyfromarcpyimportenv#Setenvironmentsettingsenv.workspace=\#SetlocalvariablesinRaster=\inConstant=100#CheckouttheArcGI3DAnalystextensionlicensearcpy.CheckOutExtension(\#ExecutePlusarcpy.Plus_3d(inRaster,inConstant,\結(jié)果：將柵格數(shù)據(jù)bb中每個像元的值加上100得到新的柵格數(shù)據(jù)3減語法：Minus_3d(in_raster_or_constant1,in_raster_or_constant2,out_raster)實例：#Name:jian.py#Description:Addsthevaluesoftworastersonacell-by-cellbasis.#Requirements:3DAnalystExtension#Author:ESRI#Importsystemmodulesimportarcpyfromarcpyimportenv#Setenvironmentsettingsenv.workspace=\#SetlocalvariablesinRaster1=\inRaster2=\#CheckouttheArcGI3DAnalystextensionlicensearcpy.CheckOutExtension(\#ExecutePlusarcpy.Minus_3d(inRaster1,inRaster2,\結(jié)果:對柵格數(shù)據(jù)output.tif進行與output4.tif逐個像元值的減法運算得出新的柵格數(shù)據(jù)。4乘語法：Times_3d(in_raster_or_constant1,in_raster_or_constant2,out_raster)實例：#Name:cheng.py#Description:Multipliesthevaluesoftworastersonacell-by-cellbasis.#Requirements:3DAnalystExtension#Author:ESRI#Importsystemmodulesimportarcpyfromarcpyimportenv#Setenvironmentsettingsenv.workspace=\#SetlocalvariablesinRaster=\inConstant=2#CheckouttheArcGIS3DAnalystextensionlicensearcpy.CheckOutExtension(\#ExecuteTimesarcpy.Times_3d(inRaster,inConstant,\結(jié)果：逐個將每個像元值乘以2輸出新的柵格。5除將兩個柵格的值逐個像元地相除或除以一個定值。當一個數(shù)除以零時，輸出結(jié)果為NoData。輸入的數(shù)據(jù)類型將決定輸出的數(shù)據(jù)類型：如果兩個輸入均為整數(shù)，則將執(zhí)行整數(shù)除法，并且輸出結(jié)果為整數(shù)。例如，如果3除以2，則輸出為1。如果兩個輸入中的任何一個為浮點型，則將執(zhí)行浮點型除法，并且輸出結(jié)果也為浮點型值。例如，如果3除以2.0，則輸出為1.5。語法：Divide_3d(in_raster_or_constant1,in_raster_or_constant2,out_raster)實例：#Name:chu.py#Description:Dividesthevaluesoftworastersonacell-by-cellbasis#Requirements:3DAnalystExtension#Author:ESRI#Importsystemmodulesimportarcpyfromarcpyimportenv#Setenvironmentsettingsenv.workspace=\#SetlocalvariablesinRaster1=\inRaster2=\#CheckouttheArcGIS3DAnalystextensionlicensearcpy.CheckOutExtension(\#ExecuteDividearcpy.Divide_3d(inRaster1,inRaster2,timesout2.tifoutput.tif6轉(zhuǎn)為整型通過截斷將柵格的每個單元值轉(zhuǎn)換為整型。1