1、地理國情監(jiān)測(cè)云平臺(tái)ArcGIS教程：柵格像元的提取一、按屬性提取摘要基于邏輯查詢提取柵格像元。插圖用法· 如果 Where 子句的求值結(jié)果是 True，則將為該像元位置返回初始輸入值。如果其求值結(jié)果是 False，則將為像元位置指定 NoData。· 對(duì)于輸出柵格，將舍棄輸入柵格中的所有額外項(xiàng)(值和計(jì)數(shù)除外)。· 如果在 Where 子句中指定了除輸入柵格的值以外的某一項(xiàng)，則將為此像元位置返回初始輸入值。· 當(dāng)輸入為多波段柵格時(shí)，將輸出一個(gè)新的多波段柵格。對(duì)輸入多波段柵格中的每一單個(gè)波段都會(huì)進(jìn)行相應(yīng)分析。默認(rèn)輸出格式是 Esri Grid 堆棧。請(qǐng)注意，

2、Esri Grid 堆棧的名稱不能以數(shù)字開頭、不能使用空格，并且長(zhǎng)度不能超過 9個(gè)字符。· 如果輸入內(nèi)容為根據(jù)多波段柵格(超過三個(gè)波段)創(chuàng)建的圖層，提取操作將只考慮圖層加載(符號(hào)化)的波段。因此，輸出的多波段柵格只會(huì)具有三個(gè)波段，分別對(duì)應(yīng)于輸入圖層顯示中使用的多波段柵格。· 如果輸入柵格為整型，那么輸出柵格也為整型。如果輸入柵格為浮點(diǎn)型，則輸出柵格也為浮點(diǎn)型。語法ExtractByAttributes (in_raster, where_clause)參數(shù)說明數(shù)據(jù)類型in_raster提取像元的輸入柵格。Raster Layerwhere_clause用于選擇柵格像元子集的

3、邏輯表達(dá)式。表達(dá)式遵循 SQL 表達(dá)式的一般格式。SQL Expression返回值名稱說明數(shù)據(jù)類型out_raster包含提取像元值的輸出柵格。Raster代碼實(shí)例按屬性提取 (ExtractByAttributes) 示例 1(Python 窗口)本示例將基于邏輯查詢從柵格中提取像元，其中高程大于 1,000 米。import arcpyfrom arcpy import envfrom arcpy.sa import *env.workspace = "C:/sapyexampl

4、es/data"attExtract = ExtractByAttributes("elevation", "VALUE > 1000")attExtract.save("c:/sapyexamples/output/attextract")  按屬性提取 (ExtractByAttributes) 示例 2(獨(dú)立腳本)本示例將基于邏輯查詢從柵格中提取像元，其中高程大于 1,000 米。# Name: ExtractByAttributes_Ex_02.py# Description: Extrac

5、ts the cells of a raster based on a logical query.# Requirements: Spatial Analyst Extension# Import system modulesimport arcpyfrom arcpy import envfrom arcpy.sa import *# Set environment settingsenv.workspace = "C:/sapyexamples/data"# Set local varia

6、blesinRaster = "elevation"inSQLClause = "VALUE > 1000"# Check out the ArcGIS Spatial Analyst extension licensearcpy.CheckOutExtension("Spatial")# Execute ExtractByAttributesattExtract = ExtractByAttributes(inRaster, inSQLClause)# Save the outputattExtract.s

7、ave("c:/sapyexamples/output/attextract02")二、按圓形區(qū)域提取摘要基于圓提取柵格像元。用法· 可通過像元的中心來確定該像元是位于圓的內(nèi)部還是位于圓的外部。如果中心位于圓弧的內(nèi)部，則即使部分像元落在圓外，也會(huì)將此像元視為完全處于圓內(nèi)。· 未選擇的像元位置被賦予 NoData 值。· 當(dāng)輸入為多波段柵格時(shí)，將輸出一個(gè)新的多波段柵格。對(duì)輸入多波段柵格中的每一單個(gè)波段都會(huì)進(jìn)行相應(yīng)分析。默認(rèn)輸出格式是 Esri Grid 堆棧。請(qǐng)注意，Esri Grid 堆棧的名稱不能以數(shù)字開頭、不能使用空格，并且長(zhǎng)度不能超過 9

8、個(gè)字符。· 如果輸入內(nèi)容為根據(jù)多波段柵格(超過三個(gè)波段)創(chuàng)建的圖層，提取操作將只考慮圖層加載(符號(hào)化)的波段。因此，輸出的多波段柵格只會(huì)具有三個(gè)波段，分別對(duì)應(yīng)于輸入圖層顯示中使用的多波段柵格。· 如果輸入柵格為整型，那么輸出柵格也為整型。如果輸入柵格為浮點(diǎn)型，則輸出柵格也為浮點(diǎn)型。語法ExtractByCircle (in_raster, center_point, radius, extraction_area)參數(shù)說明數(shù)據(jù)類型in_raster提取像元的輸入柵格。Raster Layercenter_point指示用于定義提取區(qū)域的圓的中心坐標(biāo) (x,y)&#

9、160;的點(diǎn)類。該類的形式為：·         Point (x, y)將指定坐標(biāo)使用與輸入柵格相同的地圖單位。Pointradius用于定義提取區(qū)域的圓半徑。將以地圖單位指定半徑，并且與輸入柵格的單位相同。Doubleextraction_area(可選)標(biāo)識(shí)要提取輸入圓內(nèi)部還是輸入圓外部的像元。·         INSIDE  指定應(yīng)選擇輸入圓內(nèi)部的像元并將其寫入輸出柵

10、格的關(guān)鍵字。圓形區(qū)域外部的所有像元都將在輸出柵格中獲得 NoData 值。·         OUTSIDE  指定應(yīng)選擇輸入圓外部的像元并將其寫入輸出柵格的關(guān)鍵字。圓形區(qū)域內(nèi)部的所有像元都將在輸出柵格中獲得 NoData 值。String返回值名稱說明數(shù)據(jù)類型out_raster包含提取像元值的輸出柵格。Raster代碼實(shí)例按圓形區(qū)域提取 (ExtractByCircle) 示例 1(Python 窗口)本示例將在某一點(diǎn)位置方圓 50

11、0 米的半徑范圍內(nèi)提取像元。import arcpyfrom arcpy import envfrom arcpy.sa import *env.workspace = "C:/sapyexamples/data"outExtCircle = ExtractByCircle("elevation", arcpy.Point(482838.823, 222128.982),         

12、;                       500, "INSIDE")outExtCircle.save("c:/sapyexamples/output/extcircle")按圓形區(qū)域提取 (ExtractByCircle) 示例 2(獨(dú)立腳本)本示例將在某一點(diǎn)位置方圓 1,000 米的半徑范圍內(nèi)提取像元。#

13、Name: ExtractByCircle_Ex_02.py# Description: Extracts the cells of a raster based on a circle.# Requirements: Spatial Analyst Extension# Import system modulesimport arcpyfrom arcpy import envfrom arcpy.sa import *# Set environment settingsenv.workspace = "

14、;C:/sapyexamples/data"# Set local variablesinRaster = ("elevation")centerPoint = arcpy.Point(482838.823, 222128.982)circRadius = 1000extractType = "INSIDE"# Check out the ArcGIS Spatial Analyst extension licensearcpy.CheckOutExtension("Spatial")# Exe

15、cute ExtractByCircleoutExtCircle = ExtractByCircle(inRaster, centerPoint, circRadius,                               extractType)# Save the outputoutEx

16、tCircle.save("c:/sapyexamples/output/extcircle02")三、按掩膜提取摘要提取所定義掩膜區(qū)域內(nèi)的柵格像元。插圖用法· 按掩膜提取工具與設(shè)置掩膜環(huán)境大體相同，但二者的區(qū)別在于輸入掩膜只適用于即時(shí)實(shí)例，而在環(huán)境中設(shè)置的掩膜可在發(fā)生更改或被禁用之前應(yīng)用于所有工具。· 當(dāng)輸入為多波段柵格時(shí)，將輸出一個(gè)新的多波段柵格。對(duì)輸入多波段柵格中的每一單個(gè)波段都會(huì)進(jìn)行相應(yīng)分析。默認(rèn)輸出格式是 Esri Grid 堆棧。請(qǐng)注意，Esri Grid 堆棧的名稱不能以數(shù)字開頭、不能使用空格，并且長(zhǎng)度不能超過 9個(gè)字符。· 如果

17、輸入內(nèi)容為根據(jù)多波段柵格(超過三個(gè)波段)創(chuàng)建的圖層，提取操作將只考慮圖層加載(符號(hào)化)的波段。因此，輸出的多波段柵格只會(huì)具有三個(gè)波段，分別對(duì)應(yīng)于輸入圖層顯示中使用的多波段柵格。· 如果輸入掩膜是柵格，則會(huì)將非 NoData 輸入像元位置的值復(fù)制到輸出柵格。可用來創(chuàng)建掩膜柵格的工具包括條件函數(shù)、條件測(cè)試，以及提取分析工具集中的其他工具。· 當(dāng)為輸入柵格掩膜指定多波段柵格時(shí)，將只在運(yùn)算中使用第一個(gè)波段。· 如果輸入柵格為整型，那么輸出柵格也為整型。如果輸入柵格為浮點(diǎn)型，則輸出柵格也為浮點(diǎn)型。語法ExtractByMask (in_raster, in_mask_dat

18、a)參數(shù)說明數(shù)據(jù)類型in_raster提取像元的輸入柵格。Raster Layerin_mask_data用于定義提取區(qū)域的輸入掩膜數(shù)據(jù)。它可以是柵格或要素?cái)?shù)據(jù)集。當(dāng)輸入掩膜數(shù)據(jù)為柵格時(shí)，將在輸出柵格中為掩膜數(shù)據(jù)中的 NoData 像元指定 NoData 值。Raster Layer | Feature Layer返回值名稱說明數(shù)據(jù)類型out_raster包含提取像元值的輸出柵格。Raster代碼實(shí)例按掩膜提取 (ExtractByMask) 示例 1(Python 窗口)本示例將從柵格內(nèi)提取由輸入面 shapefile 要素類所定義的掩膜內(nèi)的像元。im

19、port arcpyfrom arcpy import envfrom arcpy.sa import *env.workspace = "C:/sapyexamples/data"outExtractByMask = ExtractByMask("elevation", "mask.shp")outExtractByMask.save("C:/sapyexamples/output/maskextract")按掩膜提取 (E

20、xtractByMask) 示例 2(獨(dú)立腳本)本示例將從柵格內(nèi)提取由輸入面 shapefile 要素類所定義的掩膜內(nèi)的像元。# Name: ExtractByMask_Ex_02.py# Description: Extracts the cells of a raster that correspond with the areas#    defined by a mask.# Requirements: Spatial Analyst Extension# Import system modulesimport arcpyfrom ar

21、cpy import envfrom arcpy.sa import *# Set environment settingsenv.workspace = "C:/sapyexamples/data"# Set local variablesinRaster = "elevation"inMaskData = "mask.shp"# Check out the ArcGIS Spatial Analyst extension licensearcpy.Ch

22、eckOutExtension("Spatial")# Execute ExtractByMaskoutExtractByMask = ExtractByMask(inRaster, inMaskData)# Save the outputoutExtractByMask.save("C:/sapyexamples/output/extractmask")四、用點(diǎn)提取摘要基于一組坐標(biāo)點(diǎn)提取柵格像元。用法· 未選擇的像元位置被賦予 NoData 值。· 當(dāng)輸入為多波段柵格時(shí)，將輸出一個(gè)新的多波段柵格。對(duì)輸入多波段柵格中的每一單個(gè)波段都會(huì)

23、進(jìn)行相應(yīng)分析。默認(rèn)輸出格式是 Esri Grid 堆棧。請(qǐng)注意，Esri Grid 堆棧的名稱不能以數(shù)字開頭、不能使用空格，并且長(zhǎng)度不能超過 9個(gè)字符。· 如果輸入內(nèi)容為根據(jù)多波段柵格(超過三個(gè)波段)創(chuàng)建的圖層，提取操作將只考慮圖層加載(符號(hào)化)的波段。因此，輸出的多波段柵格只會(huì)具有三個(gè)波段，分別對(duì)應(yīng)于輸入圖層顯示中使用的多波段柵格。· 如果輸入柵格為整型，那么輸出柵格也為整型。如果輸入柵格為浮點(diǎn)型，則輸出柵格也為浮點(diǎn)型。語法ExtractByPoints (in_raster, points, extraction_area)參數(shù)說明數(shù)據(jù)類型in_raster提取像元的輸

24、入柵格。Raster Layerpointspoint,.點(diǎn)類對(duì)象的 Python 列表用于指示要提取柵格值的位置。點(diǎn)對(duì)象均在 x,y 坐標(biāo)對(duì)列表中指定。對(duì)象形式為：·         point(x1,y1), point(x2,y2),.點(diǎn)所使用的地圖單位與輸入柵格相同。Pointextraction_area(可選)標(biāo)識(shí)是基于指定點(diǎn)位置（內(nèi)部）提取像元還是基于點(diǎn)位置外部（外部）提取像元。·    &#

25、160;    INSIDE  指定將所選點(diǎn)落入的像元寫入輸出柵格的關(guān)鍵字。方框區(qū)域外部的所有像元都將在輸出柵格中獲得 NoData 值。·         OUTSIDE  指定應(yīng)選擇輸入點(diǎn)外部的像元并將其寫入輸出柵格的關(guān)鍵字。String返回值名稱說明數(shù)據(jù)類型out_raster包含提取像元值的輸出柵格。Raster代碼實(shí)例用點(diǎn)提取 (ExtractByPoints) 示例 1(Python

26、窗口)本示例將基于指定點(diǎn)坐標(biāo)提取柵格的像元。import arcpyfrom arcpy import envfrom arcpy.sa import *env.workspace = "C:/sapyexamples/data"pointList = arcpy.Point(743050, 4321275), arcpy.Point(743100, 4321200),         

27、    arcpy.Point(734500,4322000)outPointExtract = ExtractByPoints("soil", pointList,"INSIDE")outPointExtract.save("c:/sapyexamples/output/pntextract")用點(diǎn)提取 (ExtractByPoints) 示例 2(獨(dú)立腳本)本示例將基于指定點(diǎn)坐標(biāo)提取柵格的像元。# Name: ExtractByPoints_Ex_02.py# Description: Extract

28、s the cells of a raster based on a set of points.# Requirements: Spatial Analyst Extension# Import system modulesimport arcpyfrom arcpy import envfrom arcpy.sa import *# Set environment settingsenv.workspace = "C:/sapyexamples/data"# Set local variab

29、lesinRaster = "soil"pointList = arcpy.Point(743050, 4321275),             arcpy.Point(743100, 4321200),             arcpy.Point(743500, 4322000), 

30、;            arcpy.Point(742900, 4321800)# Check out the ArcGIS Spatial Analyst extension licensearcpy.CheckOutExtension("Spatial")# Execute ExtractByPointsoutPointExtract = ExtractByPoints("soil", pointList,"INSIDE

31、")# Save the outputoutPointExtract.save("c:/sapyexamples/output/pntext")五、多邊形提取摘要基于多邊形提取柵格像元。用法· 若要基于要素類中的多邊形提取像元，而不是提供一系列 x,y 坐標(biāo)對(duì)，則可以使用按掩膜提取工具。· 可通過像元的中心來確定該像元是位于多邊形的內(nèi)部還是多邊形的外部。如果中心位于多邊形弧的內(nèi)部，則即使部分像元落在多邊形之外，也會(huì)將此像元視為完全處于多邊形之內(nèi)。· 多邊形最多可以有 1,000 個(gè)折點(diǎn)。多邊形折點(diǎn)必須按順時(shí)針順序輸入。如果要使用多個(gè)多

32、邊形，則第一個(gè)折點(diǎn)和最后一個(gè)折點(diǎn)必須相同，以使多邊形閉合。如果最后的點(diǎn)不相同，多邊形將自動(dòng)閉合。多邊形的弧之間可以相交，但不建議您使用過于復(fù)雜的多邊形。· 未選擇的像元位置被賦予 NoData 值。· 當(dāng)輸入為多波段柵格時(shí)，將輸出一個(gè)新的多波段柵格。對(duì)輸入多波段柵格中的每一單個(gè)波段都會(huì)進(jìn)行相應(yīng)分析。默認(rèn)輸出格式是 Esri Grid 堆棧。請(qǐng)注意，Esri Grid 堆棧的名稱不能以數(shù)字開頭、不能使用空格，并且長(zhǎng)度不能超過 9個(gè)字符。· 如果輸入內(nèi)容為根據(jù)多波段柵格(超過三個(gè)波段)創(chuàng)建的圖層，提取操作將只考慮圖層加載(符號(hào)化)的波段。因此，輸出的多波段柵格只會(huì)具有

33、三個(gè)波段，分別對(duì)應(yīng)于輸入圖層顯示中使用的多波段柵格。· 如果輸入柵格為整型，那么輸出柵格也為整型。如果輸入柵格為浮點(diǎn)型，則輸出柵格也為浮點(diǎn)型。語法ExtractByPolygon (in_raster, polygon, extraction_area)參數(shù)說明數(shù)據(jù)類型in_raster提取像元的輸入柵格。Raster Layerpolygonpoint,.用于定義要提取的輸入柵格區(qū)域的一個(gè)或多個(gè)多邊形。每個(gè)多邊形部分都是由點(diǎn)類定義的一系列折點(diǎn)。作為可選的，可使用多邊形類來定義各個(gè)多邊形部分。點(diǎn)將指定為 x,y 坐標(biāo)對(duì)。對(duì)象形式為：· 

34、0;       point(x1,y1), point(x2,y2),point(xn,yn),.point(x1,y1), point(x'1,y'1), point(x'2,y'2),point(x'n,y'n),.,point(x'1,y'1)請(qǐng)注意，最后一個(gè)坐標(biāo)應(yīng)與第一個(gè)坐標(biāo)相同，從而使多邊形閉合。Pointextraction_area(可選)標(biāo)識(shí)要提取輸入多邊形內(nèi)部還是輸入多邊形外部的像元。·    

35、;     INSIDE  指定應(yīng)選擇輸入多邊形內(nèi)部的像元并將其寫入輸出柵格的關(guān)鍵字。多邊形區(qū)域外部的所有像元都將在輸出柵格中獲得 NoData 值。·         OUTSIDE  指定應(yīng)選擇輸入多邊形外部的像元并將其寫入輸出柵格的關(guān)鍵字。多邊形區(qū)域內(nèi)部的所有像元都將在輸出柵格中獲得 NoData 值。String返回值名稱說明數(shù)據(jù)類型out_raster包含提

36、取像元值的輸出柵格。Raster代碼實(shí)例多邊形提取 (ExtractByPolygon) 示例 1(Python 窗口)本示例將基于定義的多邊形坐標(biāo)提取柵格的像元。import arcpyfrom arcpy import envfrom arcpy.sa import *polyPoints = arcpy.Point(743050, 4321275), arcpy.Point(743100, 4321200),       

37、0;     arcpy.Point(743500, 4322000),arcpy.Point(742900, 4321800)env.workspace = "C:/sapyexamples/data"extPolygonOut = ExtractByPolygon("soil", polyPoints, "INSIDE")extPolygonOut.save("c:/sapyexamples/output/extpoly")多邊形提

38、取 (ExtractByPolygon) 示例 2(獨(dú)立腳本)本示例將基于定義的多邊形坐標(biāo)提取柵格的像元。# Name: ExtractByPolgyon_Ex_02.py# Description: Extracts the cells of a raster based on a polygon.# Requirements: Spatial Analyst Extension# Import system modulesimport arcpyfrom arcpy import envfrom arcpy.sa import&#

39、160;*# Set environment settingsenv.workspace = "C:/sapyexamples/data"# Set local variablesinRaster = "soil"polyPoints = arcpy.Point(743050, 4321275), arcpy.Point(743100, 4321200),             arcpy.Po

40、int(743500, 4322000),arcpy.Point(742900, 4321800)# Check out the ArcGIS Spatial Analyst extension licensearcpy.CheckOutExtension("Spatial")# Execute ExtractByPolygonextPolygonOut = ExtractByPolygon(inRaster, polyPoints, "INSIDE")# Save the outputextPolygonOut.save(

41、"c:/sapyexamples/output/extpoly02")六、用矩形提取摘要基于矩形提取柵格像元。用法· 可通過像元的中心來確定該像元是位于矩形的內(nèi)部還是矩形的外部。如果中心位于矩形輪廓的內(nèi)部，則即使部分像元落在矩形之外，也會(huì)將此像元視為完全處于矩形內(nèi)部。· 未選擇的像元位置被賦予 NoData 值。· 當(dāng)輸入為多波段柵格時(shí)，將輸出一個(gè)新的多波段柵格。對(duì)輸入多波段柵格中的每一單個(gè)波段都會(huì)進(jìn)行相應(yīng)分析。默認(rèn)輸出格式是 Esri Grid 堆棧。請(qǐng)注意，Esri Grid 堆棧的名稱不能以數(shù)字開頭、不能使用空格，并且長(zhǎng)度不能超過 9個(gè)字

42、符。· 如果輸入內(nèi)容為根據(jù)多波段柵格(超過三個(gè)波段)創(chuàng)建的圖層，提取操作將只考慮圖層加載(符號(hào)化)的波段。因此，輸出的多波段柵格只會(huì)具有三個(gè)波段，分別對(duì)應(yīng)于輸入圖層顯示中使用的多波段柵格。· 如果輸入柵格為整型，那么輸出柵格也為整型。如果輸入柵格為浮點(diǎn)型，則輸出柵格也為浮點(diǎn)型。語法ExtractByRectangle (in_raster, rectangle, extraction_area)參數(shù)說明數(shù)據(jù)類型in_raster提取像元的輸入柵格。Raster Layerrectangleextent用于定義待提取區(qū)域的矩形。可以使用范圍對(duì)象來指定坐標(biāo)。對(duì)象形式為：

43、3;         Extent(XMin，YMin，XMax，YMax)其中 XMin 和 YMin 定義待提取區(qū)域左下方的坐標(biāo)，XMax 和 YMax 定義右上方的坐標(biāo)。將指定坐標(biāo)使用與輸入柵格相同的地圖單位。Extentextraction_area(可選)標(biāo)識(shí)要提取輸入矩形內(nèi)部還是輸入矩形外部的像元。·         INSI

44、DE  指定應(yīng)選擇輸入矩形內(nèi)部的像元并將其寫入輸出柵格的關(guān)鍵字。矩形區(qū)域外部的所有像元都將在輸出柵格中獲得 NoData 值。·         OUTSIDE  指定應(yīng)選擇輸入矩形外部的像元并將其寫入輸出柵格的關(guān)鍵字。矩形區(qū)域內(nèi)部的所有像元都將在輸出柵格中獲得 NoData 值。String返回值名稱說明數(shù)據(jù)類型out_raster包含提取像元值的輸出柵格。Raster代碼實(shí)例用矩形提取 (ExtractByRectangle) 示例 1(Python 窗口)本示例會(huì)將矩形范圍以外的像元提取到新柵格。import ar