ArcGIS教程：多值和值的提取至點(空間分析)

1、地理國情監(jiān)測云平臺北京數(shù)字空間科技有限公司地理國情監(jiān)測云平臺ArcGIS教程：多值和值的提取至點（空間分析）多值提取至點摘要在點要素類的指定位置提取一個或多個柵格像元值，并將值記錄到點要素類的屬性表中。用法可以將任意柵格組合（單波段或多波段）指定為輸入。將為每個輸入柵格提取一個像元值，并會將包含該提取值的新字段追加到輸入點要素類。默認(rèn)情況下，將根據(jù)輸入柵格的名稱來創(chuàng)建輸岀字段的名稱?；蛘吣部梢詾槊總€存儲柵格值的字段指定唯一名稱。當(dāng)輸入多波段柵格數(shù)據(jù)時，為所有波段添加的輸岀字段的名稱中將包含指示波段號的"blb2_, ,bn"前綴插值選項決定了從柵格中獲取值的方式。默認(rèn)選項

2、為使用被采樣像元的中心值。插值選項將使用雙線性插值法在像元中心內(nèi)插一個值。 Shapefile 格式的字段最大長度限制為10個字符。默認(rèn)情況下，追加到輸入shapefile 的輸出字段將被截斷并獲得唯一值。如果名稱很長或很相似，則可能很難區(qū)分各輸入柵格。在這種情況下，建議 您將要素轉(zhuǎn)換為文件地理數(shù)據(jù)庫。柵格中的NoData像元將被賦予空值。shapefile 不支持空值，而是將其賦予0（零）值。不支持將多點數(shù)據(jù)集作為輸入。語法bilinear_interpolate_values)ExtractMultiValuesToPoints (in_point_features, in_rasters

3、.參數(shù)說明數(shù)據(jù)類型n_point_features要添加?xùn)鸥裰档妮斎朦c要素。Feature Layern_rastersRaster, Output FieldName,.要基于輸入點要素的位置提取柵格值。您還可以為存儲柵格值的字段指定名稱。默認(rèn)情況下，將根據(jù) 輸入柵格數(shù)據(jù)集的名稱創(chuàng)建唯一的字段名稱。Value Tablebilinear_interpolate_values指定是否使用插值。Boolean可選） NONE不應(yīng)用任何插值法；將使用像元中心值。 BILINEAR 將使用雙線性插值法根據(jù)相鄰像元的有效 值計算像元值。除非所有相鄰像元都為NoData，否則會在插值時忽略NoData值

4、。代碼實例多值提取至點 （ExtractMultiValuesToPoints）示例 1（Python 窗口）將多個柵格的像元值提取到shapefile點要素類的屬性中。import arcpyfrom arcpy.sa import *from arcpy import envenv.workspace = "c:/sapyexamples/data"ExtractMultiValuesToPoints( "observers.shp" , "elevation" , "ELEV","costraste

5、r" , "COST', "flowdir" , "DIR" ,"NONE')多值提取至點 (ExtractMultiValuesToPoints)示例2(獨立腳本)使用插值法將多個柵格的像元值提取到shapefile點要素類的屬性中。# Name: ExtractMultiValuesToPoints_Ex_02.py# Description: Extracts the cells of multiple rasters as attributes in# an output point feature

6、class. This example takes a multiband IMG# and two GRID files as input.# Requirements: Spatial Analyst Extension# Import system modulesimport arcpyfrom arcpy import envfrom arcpy.sa import *# Set environment settingsenv.workspace = "C:/sapyexamples/data"# Set local variablesinPointFeatures

7、 = "poi.shp"inRasterList = "doqq.img" , "doqqval" , "redstd" , "focalstd","redmin" , "focalmin"# Check out the ArcGIS Spatial Analyst extension license arcpy.CheckOutExtension( "Spatial")# Execute ExtractValuesToPointsEx

8、tractMultiValuesToPoints(inPointFeatures, inRasterList,"BILINEAR")值提取至點摘要基于一組點要素提取柵格像元值，并將這些值記錄到輸岀要素類的屬性表。用法輸入點要素類中的所有字段均將包括在輸岀點要素類中。輸出要素類將添加一個名為RASTERVALC的新字段。對多波段柵格使用值提取至點時，RASTERVALU字段將包含輸入柵格的最后一個波段的值。要從多個柵格或某個多波段柵格數(shù)據(jù)集提取值，請使用多值提取至點工具。對于屬性表的 RASTERVALU字段，若賦值柵格中的像元值為 NoData，則該像元的字段值為 -999

9、9。插值選項決定了從柵格中獲取值的方式。默認(rèn)選項為使用被采樣像元的中心值。插值選項將使用雙線性插值法在像元中心內(nèi)插一個值。如果輸入柵格為浮點型，則生成的輸岀點數(shù)據(jù)集將僅包含輸入要素數(shù)據(jù)的屬性和像元值（由插值選項確定）。添加來自輸入柵格的屬性時，如果輸出點要素數(shù)據(jù)集是 shapefile ，則輸入柵格中不能存在名稱長度超過10個字符的字段。如果存在此類字段，則必須在運行該工具之前對其進(jìn)行重命名。不支持將多點數(shù)據(jù)集作為輸入。語法ExtractValuesToPoints (in_point_features, in_raster, out_point_features,interpolate_va

10、lues, add_attributes)參數(shù)說明數(shù)據(jù)類型in_point_features用于定位待提取柵格像元值的輸入點要素。Feature Layerin_raster待提取值的柵格數(shù)據(jù)集。其柵格數(shù)據(jù)可為整型或浮點型。Raster Layerout_point_features包含提取的柵格值的輸岀點要素數(shù)據(jù)集。Feature Classinterpolate_values (可選)指定是否使用插值。 NONE不應(yīng)用任何插值法；將使用像元中心值。 INTERPOLATE-將使用雙線性插值法根據(jù)相鄰像元的有效值計算像元值。除非所有相鄰像元都為NoData，否則會在插值時忽略 NoData值

11、。Booleanadd_attributes確定是否將柵格屬性寫入輸出點要素數(shù)據(jù)集。Boolean（可選）* VALUE_ONL* 僅將輸入柵格的值添加到點屬性。這是默認(rèn)設(shè)置。« ALL 輸入柵格的所有字段（“計數(shù)”除外）都將添加到點屬性。代碼實例值提取至點（ExtractValuesToPoints） 示例 1（Python 窗口）本示例根據(jù)點shapefile所定義的位置提取柵格的像元值，并創(chuàng)建這些值的輸岀點要素類。import arcpyfrom arcpy import envfrom arcpy.sa import *env.workspace = "C:/sap

12、yexamples/data"ExtractValuesToPoints("rec_sites.shp", "elevation" ,"C:/sapyexamples/output/outValPnts", "INTERPOLATE""VALUE_ONLY"值提取至點（ExtractValuesToPoints） 示例2（獨立腳本）本示例根據(jù)點shapefile所定義的位置提取柵格的像元值，并創(chuàng)建這些值的輸出點shapefile# Name: ExtractValuesToPoints

13、_Ex_02.py# Description: Extracts 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 variablesinPointFeatures = "rec_sites.shp"inRaster = "elevation"outPointFeatures = "C:/sapyexamples/output/extractvaluespts.shp"# Check out the ArcGIS Spatial Analyst extension licensearcpy.CheckOutExtension( "Spatial")# Execute ExtractValuesToPointsE