全国高分辨率土地利用数据服务 土地利用数据服务 土地覆盖数据服务 坡度数据服务 土壤侵蚀数据服务 全国各省市DEM数据服务 耕地资源空间分布数据服务 草地资源空间分布数据服务 林地资源空间分布数据服务 水域资源空间分布数据服务 建设用地空间分布数据服务 地形、地貌、土壤数据服务 分坡度耕地数据服务 全国大宗农作物种植范围空间分布数据服务
多种卫星遥感数据反演植被覆盖度数据服务 地表反照率数据服务 比辐射率数据服务 地表温度数据服务 地表蒸腾与蒸散数据服务 归一化植被指数数据服务 叶面积指数数据服务 净初级生产力数据服务 净生态系统生产力数据服务 生态系统总初级生产力数据服务 生态系统类型分布数据服务 土壤类型质地养分数据服务 生态系统空间分布数据服务 增强型植被指数数据服务
多年平均气温空间分布数据服务 多年平均降水量空间分布数据服务 湿润指数数据服务 大于0℃积温空间分布数据服务 光合有效辐射分量数据服务 显热/潜热信息数据服务 波文比信息数据服务 地表净辐射通量数据服务 光合有效辐射数据服务 温度带分区数据服务 山区小气候因子精细数据服务
全国夜间灯光指数数据服务 全国GDP公里格网数据服务 全国建筑物总面积公里格网数据服务 全国人口密度数据服务 全国县级医院分布数据服务 人口调查空间分布数据服务 收入统计空间分布数据服务 矿山面积统计及分布数据服务 载畜量及空间分布数据服务 农作物种植面积统计数据服务 农田分类面积统计数据服务 农作物长势遥感监测数据服务 医疗资源统计数据服务 教育资源统计数据服务 行政辖区信息数据服务
Landsat 8 高分二号 高分一号 SPOT-6卫星影像 法国Pleiades高分卫星 资源三号卫星 风云3号 中巴资源卫星 NOAA/AVHRR MODIS Landsat TM 环境小卫星 Landsat MSS 天绘一号卫星影像
基于一组坐标点提取栅格像元。
· 未选择的像元位置被赋予 NoData 值。
· 当输入为多波段栅格时，将输出一个新的多波段栅格。对输入多波段栅格中的每一单个波段都会进行相应分析。
默认输出格式是 Esri Grid 堆栈。请注意，Esri Grid 堆栈的名称不能以数字开头、不能使用空格，并且长度不能超过 9 个字符。
· 如果输入内容为根据多波段栅格（超过三个波段）创建的图层，提取操作将只考虑图层加载（符号化）的波段。因此，输出的多波段栅格只会具有三个波段，分别对应于输入图层显示中使用的多波段栅格。
· 如果输入栅格为整型，那么输出栅格也为整型。如果输入栅格为浮点型，则输出栅格也为浮点型。
ExtractByPoints (in_raster, points, {extraction_area})
|
参数 |
说明 |
数据类型 |
|
in_raster |
提取像元的输入栅格。 |
Raster Layer |
|
points [point,...] |
点类对象的 Python 列表用于指示要提取栅格值的位置。 点对象均在 x,y 坐标对列表中指定。对象形式为： · [point(x1,y1), point(x2,y2),...] 点所使用的地图单位与输入栅格相同。 |
Point |
|
extraction_area (可选) |
标识是基于指定点位置（内部）提取像元还是基于点位置外部（外部）提取像元。 · INSIDE — 指定将所选点落入的像元写入输出栅格的关键字。方框区域外部的所有像元都将在输出栅格中获得 NoData 值。 · OUTSIDE — 指定应选择输入点外部的像元并将其写入输出栅格的关键字。 |
String |
|
名称 |
说明 |
数据类型 |
|
out_raster |
包含提取像元值的输出栅格。 |
Raster |
本示例将基于指定点坐标提取栅格的像元。
import arcpy
from arcpy import env
from arcpy.sa import *
env.workspace = "C:/sapyexamples/data"
pointList = [arcpy.Point(743050, 4321275), arcpy.Point(743100, 4321200),
arcpy.Point(734500,4322000)]
outPointExtract = ExtractByPoints("soil", pointList,"INSIDE")
outPointExtract.save("c:/sapyexamples/output/pntextract")
本示例将基于指定点坐标提取栅格的像元。
# Name: ExtractByPoints_Ex_02.py
# Description: Extracts the cells of a raster based on a set of points.
# Requirements: Spatial Analyst Extension
# Import system modules
import arcpy
from arcpy import env
from arcpy.sa import *
# Set environment settings
env.workspace = "C:/sapyexamples/data"
# Set local variables
inRaster = "soil"
pointList = [arcpy.Point(743050, 4321275),
arcpy.Point(743100, 4321200),
arcpy.Point(743500, 4322000),
arcpy.Point(742900, 4321800)]
# Check out the ArcGIS Spatial Analyst extension license
arcpy.CheckOutExtension("Spatial")
# Execute ExtractByPoints
outPointExtract = ExtractByPoints("soil", pointList,"INSIDE")
# Save the output
outPointExtract.save("c:/sapyexamples/output/pntext")
