Python ArcGIS ImageryLayer.compute_histograms用法及代碼示例

本文簡要介紹 python 語言中 arcgis.raster.ImageryLayer.compute_histograms 的用法。

用法:

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

返回：

一個字典

compute_histograms 方法返回給定範圍內為 ImageryLayer 計算的所有柵格波段的直方圖列表。

Parameter	Description
geometry	必需的 Geometry ( Polygon 或 Envelope )。定義在其中計算直方圖的幾何圖形的幾何圖形。
mosaic_rule	可選字符串。在定義應如何鑲嵌單個圖像時指定鑲嵌規則。未指定鑲嵌規則時，將使用影像圖層的默認鑲嵌規則(如根資源中所宣傳的：defaultMosaicMethod、mosaicOperator、sortField、sortValue)。有關詳細信息，請參閱Mosaic rule objects 幫助。
rendering_rule	可選規則。指定應如何處理請求的圖像的呈現規則。響應是更新的層信息，反映了由呈現規則定義的自定義處理。例如，如果renderingRule 包含一個attributeTable 函數，則響應將指示“hasRasterAttributeTable”: true；如果 renderingRule 包含改變波段數的函數，則響應將指示正確的 bandCount 值。
pixel_size	可選字符串或字典。正在使用的像素級別(或正在查看的分辨率)。如果未指定像素大小，則 pixel_size 將默認為數據集的基本分辨率。 pixel_size 參數的結構與ArcGIS REST API 返回的點對象的結構相同。除了字典結構之外，您還可以使用逗號分隔的字符串指定像素大小。 用法: 字典結構：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 - The histogram 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 histograms of pixel values from all selected slices are computed. 注意： 在 10.9 中為使用 ArcObjects11 或 ArcObjectsRasterRendering 作為服務提供者的圖像服務添加。

示例 1：

# Usage Example 1: Compute the 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_hist_01 = img_lyr.compute_histograms(geometry=aoi,
                                          rendering_rule={"rasterFunction":None},
                                          time="1326650400000")

示例 2：

# Usage Example 2: Compute the 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_hist_02 = img_lyr.compute_histograms(geometry=aoi,
                                          rendering_rule={"rasterFunction":None},
                                          time=[start, end])