Python ArcGIS ImageryLayer.compute_stats_and_histograms用法及代码示例

本文简要介绍 python 语言中 arcgis.raster.ImageryLayer.compute_stats_and_histograms 的用法。

用法:

compute_stats_and_histograms(geometry, mosaic_rule=None, rendering_rule=None, pixel_size=None, time=None, process_as_multidimensional=False)

返回：

字典

compute_stats_and_histograms 方法计算给定范围内 ImageryLayer 对象的统计数据和直方图。

Parameter	Description
geometry	必需的 Geometry ( Polygon 或 Envelope )。定义几何图形的几何图形，在该几何图形中计算统计数据和直方图。
mosaic_rule	可选字典。在定义应如何镶嵌单个图像时指定镶嵌规则。未指定镶嵌规则时，将使用影像图层的默认镶嵌规则(如根资源中所宣传的：defaultMosaicMethod、mosaicOperator、sortField、sortValue)。
rendering_rule	可选字典。指定应如何呈现所请求图像的呈现规则。
pixel_size	可选字符串或字典。正在使用的像素级别(或正在查看的分辨率)。如果未指定像素大小，则 pixel_size 将默认为数据集的基本分辨率。镶嵌数据集中指定像素大小的栅格将用于直方图计算。 用法: 字典结构：pixel_size={point} 点简单语法：pixel_size='<x>,<y>' 例子： pixel_size={“x”: 0.18, “y”: 0.18} pixel_size='0.18,0.18'
time	可选 datetime.date、datetime.datetime 或时间戳字符串。计算统计数据和直方图的时刻或时间范围。指定为 datetime.date、datetime.datetime 或自纪元以来的毫秒数的时间戳语法：time=<timeInstant> 指定为 [<startTime>, <endTime>] 列表的时间范围 对于时间范围，<startTime> 或 <endTime> 之一可以为 None。为开始时间或结束时间指定的“无”值将分别表示开始时间或结束时间为无穷大。语法：time=[<开始时间>, <结束时间>] ;指定为 datetime.date、datetime.datetime 或时间戳 注意： 该参数是在 10.8 中添加的。
process_as_multidimensional	可选的布尔值。指定是否将影像服务作为多维影像服务进行处理。 False - Statistics and histograms of pixel values from only the first slice is computed. This is the default. True - The image service is treated as a multidimensional raster, and statistics and histograms of pixel values from all selected slices are computed. 注意： 在 10.9 中为使用 ArcObjects11 或 ArcObjectsRasterRendering 作为服务提供者的图像服务添加。

示例 1：

# Usage Example 1: Compute the stats and histograms in the specified area of interest for a time instant.

aoi = {
        "spatialReference": {"wkid": 32610},
        "xmax": 725000,
        "xmin": 720000,
        "ymax": 4300000,
        "ymin": 4250000
      }

aoi_geometry = Geometry(aoi)

comp_stats_hist_01 = img_lyr.compute_stats_and_histograms(geometry=aoi,
                                                          rendering_rule={"rasterFunction":None},
                                                          time="1326650400000")

示例 2：

# Usage Example 2: Compute the stats and histograms in the specified area of interest for a time extent.

aoi = {
        "spatialReference": {"wkid": 32610},
        "xmax": 725000,
        "xmin": 720000,
        "ymax": 4300000,
        "ymin": 4250000
      }

aoi_geometry = Geometry(aoi)

# If the datetime object is not in the UTC timezone, the API will internally convert it to the UTC timezone.
start = datetime.datetime(2012,1,15,18,0,0, tzinfo=datetime.timezone.utc)
end = datetime.datetime(2012,1,15,21,0,0, tzinfo=datetime.timezone.utc)

comp_stats_hist_02 = img_lyr.compute_stats_and_histograms(geometry=aoi,
                                                          rendering_rule={"rasterFunction":None},
                                                          time=[start,end])