ARCGIS10柵格計(jì)算PYTHON代碼

1、使用python語言進(jìn)行柵格計(jì)算1重采樣通過更改像元大小和重采樣方法來更改柵格數(shù)據(jù)集用法：1）可更改像元大小，但柵格數(shù)據(jù)集的范圍將保持不變。2）本工具僅能輸出方形像元大小。3）可將輸出保存為 BIL、BIP、BMP、BSQ、DAT、GIF、GRID、IMG、JPEG、JPEG 2000、 PNG、TIFF格式或任意地理數(shù)據(jù)庫(kù)柵格數(shù)據(jù)集。4）有四個(gè)用于重采樣技術(shù)參數(shù)的選項(xiàng)：“最鄰近”選項(xiàng)，用于執(zhí)行最鄰近分配法，是速度最快的插值方法。此選項(xiàng)主要 用于離散數(shù)據(jù)（如土地使用分類），因?yàn)樗粫?huì)更改像元的值。最大空間誤差將 是像元大小的一半?！氨姅?shù)”選項(xiàng)，用于執(zhí)行眾數(shù)算法，可根據(jù)過濾器窗口內(nèi)的最常用值確定

2、像元的 新值。與最鄰近法一樣，此選項(xiàng)主要用于離散數(shù)據(jù)；但與“最鄰近”選項(xiàng)相比，“眾數(shù)”選項(xiàng)通?？缮筛交慕Y(jié)果?！半p線性”選項(xiàng)，用于執(zhí)行雙線性插值法，可根據(jù)四個(gè)最鄰近輸入像元中心的加 權(quán)平均距離確定像元的新值。此選項(xiàng)用于連續(xù)數(shù)據(jù)，并會(huì)生成平滑的數(shù)據(jù)?！叭巍边x項(xiàng)，用于執(zhí)行三次卷積插值法，可通過擬合穿過16個(gè)最鄰近輸入像 元中心的平滑曲線確定像元的新值。此選項(xiàng)適用于連續(xù)數(shù)據(jù)，盡管所生成的輸出 柵格可能會(huì)包含輸入柵格范圍以外的值。與通過運(yùn)行最鄰近重采樣算法獲得的柵 格相比，輸出柵格的幾何變形程度較小?！叭巍边x項(xiàng)的缺點(diǎn)是需要更多的處理 時(shí)間。在某些情況下，此選項(xiàng)會(huì)使輸出像元值位于輸入像元值范圍之

3、外。如果無 法接受此結(jié)果，請(qǐng)轉(zhuǎn)而使用“雙線性”選項(xiàng)?！半p線性”或“三次”選項(xiàng)不得用于分類數(shù)據(jù)，因?yàn)橄裨悼赡鼙桓摹?）輸出柵格數(shù)據(jù)集的左下角與輸入柵格數(shù)據(jù)集的左下犯具有相同的地圖空間坐 標(biāo)位置。語法：Resample_management （in_raster, out_raster, cell_size,resampling type）實(shí)例：ttResampletry:import arcpyarcpy. env. workspace = rD:/work”#Resample TIFF image to a higher resolutionarcpy.Resample_managemen

4、tCaa.tif, resample3. tif”, 4. 0e-006， “CUBIC）# NEAREST 最鄰近 | BILINEAR 雙線性 | CUBIC 三次 | MAJORITY 眾數(shù)except:print Resample example failed.print arcpy. GetMessages()2加逐個(gè)像元地將兩個(gè)柵格的值相加（求和）或與一定的數(shù)值相加。語法：Pl us 3d (in raster or constant!, in raster or constant2, out raster)實(shí)例：Name: jia. pyDescription: Adds the

5、 values of two rasters on a cell-by-cell basis.Requirements: 3D Analyst ExtensionAuthor: ESRIImport system modulesimport arcpyfrom arcpy import envSet environment settings env. workspace = D:/workSet local variables inRaster = bb.tif” inConstant = 100Check out the ArcGI 3D Analyst extension license

6、arcpy. CheckOutExtension(，Z3D，/)Execute Plusarcpy. Plus 3d (inRaster, inConstant, D:/work/output3. tif)結(jié)果：將柵格數(shù)據(jù)bb中每個(gè)像元的值加上100得到新的柵格數(shù)據(jù)3減逐個(gè)像元地從第一個(gè)輸入柵格的值中減去第二個(gè)輸入柵格的值或減去一個(gè)確定 的數(shù)值。語法：Minus_3d (in_raster_or_constant1, in_raster_or_constant2, out_raster)實(shí)例：Name: jian. pyDescription: Adds the values of two r

7、asters on a cell-by-cell basis.Requirements: 3D Analyst ExtensionAuthor: ESRIImport system modules import arcpy from arcpy import envSet environment settings env. workspace = D:/workSet local variables inRasterl = output. tif inRaster2 = output4. tif”Check out the ArcGI 3D Analyst extension licensea

8、rcpy.CheckOutExtension(3D)Execute Plusarcpy. Minus_3d(inRasterl, inRaster2, Di/work/outputG. tif)結(jié)果:對(duì)柵格數(shù)據(jù)output, tif進(jìn)行與output4. tif逐個(gè)像元值的減法運(yùn)算得出新 的柵格數(shù)據(jù)。4乘將兩個(gè)柵格的值逐個(gè)像元地相乘或乘以一個(gè)定值。語法：Times_3d (in_raster_or_constantl, in_raster_or_constant2, out_raster)實(shí)例：Name: cheng. pyDescription: Multiplies the values

9、of two rasters on a cell-by-cell basis.Requirements: 3D Analyst ExtensionAuthor: ESRIImport system modules import arcpy from arcpy import envSet environment settings env. workspace = D:/workSet local variables inRaster = output. tif inConstant = 2Check out the ArcGIS 3D Analyst extension license arc

10、py. CheckOutExtension(Spatial)Execute Timesarcpy. Times 3d(inRaster, inConstant, D:/work/timesout2. tif)結(jié)果：逐個(gè)將每個(gè)像元值乘以2輸出新的柵格。5除將兩個(gè)柵格的值逐個(gè)像元地相除或除以一個(gè)定值。當(dāng)一個(gè)數(shù)除以零時(shí)，輸出結(jié)果為NoData。輸入的數(shù)據(jù)類型將決定輸出的數(shù)據(jù)類型：如果兩個(gè)輸入均為整數(shù)，那么將執(zhí)行整數(shù)除法，并且輸出結(jié)果為整數(shù)。例如，如果 3除以2,那么輸出為k如果兩個(gè)輸入中的任何一個(gè)為浮點(diǎn)型，那么將執(zhí)行浮點(diǎn)型除法，并且輸出結(jié)果也為 浮點(diǎn)型值。例如，如果3除以2.0,那么輸出為1.5。語

11、法：Dividc_3d (inrasterorconstantl, in_raster_or_constant2, out_raster) 實(shí)例：Name: chu.pyDescription: Divides the values of two rasters on a cel 1-by-cel 1 basisRequirements: 3D Analyst ExtensionAuthor: ESRIImport system modules import arcpyfrom arcpy import envSet environment settings env. workspace =

12、D:/workSet local variablesinRasterl = z，timesout2. tifinRaster2 = output. tifCheck out the ArcGIS 3D Analyst extension license arcpy. CheckOutExtens ionExecute Dividearcpy. Divide_3d（inRasterl, inRaster2, DVwork/outdivide. tif）結(jié)果：將timesout2. tif逐個(gè)像元的與output. tif的像元值相除輸出新柵格。6轉(zhuǎn)為整型通過截?cái)鄬鸥竦拿總€(gè)單元值轉(zhuǎn)換為整型。1輸

13、入值可為正數(shù)，也可為負(fù)數(shù)。2如果相比于截?cái)喔Ｍ捎蒙崛敕绞?，那么?qǐng)?jiān)趫?zhí)行操作之前添加0.5個(gè)輸入 柵格。3整型柵格值的最大支持范圍是從-2, 147, 483, 648（由2的31次確定的最小值） 到2, 147,483, 647 （由2的31次減去1確定的最大值）。如果將轉(zhuǎn)為整型用在 浮點(diǎn)型柵格上，且該浮點(diǎn)型柵格具有所含值超出該范圍的像元，那么這些像元在輸 出柵格中的值將為NoDatao4將分類（離散）數(shù)據(jù)作為整型柵格存儲(chǔ)所占用的磁盤空間要明顯少于將同樣的 信息作為浮點(diǎn)型柵格存儲(chǔ)所占用的磁盤空間。如有可能，建議使用該工具將浮點(diǎn) 型柵格轉(zhuǎn)換為整型。語法：Int_3d (in_raster_o

14、r_constant, out_raster)實(shí)例：Name: int. pyDescription: Converts each cell value of a raster to an integer by truncationRequirements: 3D Analyst ExtensionAuthor: ESRTImport system modules import arcpy from arcpy import envSet environment settings env. workspace =，D:/work，/Set local variables inRaster =

15、timesoutl.tif”Check out the ArcGIS 3D Analyst extension licensearcpy. CheckOutExtension(3D)Execute Intarcpy. Tnt 3d(inRaster, D:/work/outint. tif)7轉(zhuǎn)為浮點(diǎn)型將每個(gè)柵格像元的值轉(zhuǎn)換為浮點(diǎn)型表達(dá)形式。1輸入值可為正數(shù)，也可為負(fù)數(shù)。2如果對(duì)已經(jīng)是浮點(diǎn)型的輸入執(zhí)行轉(zhuǎn)為浮點(diǎn)型操作，那么輸出值將與輸入值相同。語法：Float_3d (in raster_or_constant, out_raster)實(shí)例：Name: float, pyDescription: Converts each cell value of a raster into a floating-point representationRequirements: Spatial Analyst ExtensionAuthor: ESRIImport system modules import arcpyfrom arcpy import env from arcpy. sa import *Set environment settings