Python ArcGIS Raster用法及代碼示例

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

用法:

class arcgis.raster.Raster(path, is_multidimensional=False, extent=None, cmap=None, opacity=None, engine=None, gis=None)

Raster 對象是引用柵格的變量。可用於查詢柵格數據集的屬性。

用法：arcgis.raster.Raster(path, is_multidimensional=False,  engine=None, gis=None)

Raster 類可以與 arcpy 引擎或圖像服務器引擎一起使用。默認情況下，如果路徑是本地路徑，則 Raster 類使用 arcpy 引擎，否則將使用 image_server 引擎。

Parameter	Description
path	必需的字符串。輸入柵格。 例子： 路徑 = r”/path/to/raster” 路徑 = “https://myserver/arcgis/rest/services/ImageServiceName/ImageServer” 路徑 = “/fileShares/file_share_name/path/to/raster” 路徑 = “/cloudStores/cloud_store_name/path/to/raster” 路徑 = “https://sentinel-cogs.s3.us-west-2.amazonaws.com/sentinel-s2-l2a-cogs/43/M/BP/2021/6/S2A_43MBP_20210622_0_L2A/B08.tif” 使用數據存儲柵格或非影像服務 URL 時，應在活動 GIS 連接中啟用RasterRendering 服務
is_multidimensional	可選的布爾值。確定是否將輸入柵格視為多維。 如果輸入是多維的並且應作為多維處理，則指定 True，其中對數據集中的每個切片進行處理。如果輸入不是多維的，或者它是多維的並且不應作為多維處理，則指定 False。 默認為假
extent	可選字典。如果無法自動推斷輸入柵格的範圍，請傳入表示柵格範圍的字典，以便在 MapView 小部件上查看時使用。 例子： {“xmin”：-74.22655， “ymin”: 40.712216, “xmax”：-74.12544， “ymax”: 40.773941, “spatialReference”： { “wkid”: 4326 } }
cmap	可選的字符串。在 MapView 小部件中顯示 1 波段柵格時，要應用於柵格的 matplotlib 顏色圖。有關兼容值的列表，請參閱 arcgis.mapping.symbol.display_colormaps() 。
opacity	可選號碼。在 MapView 小部件中顯示柵格時，應用什麽不透明度。 0 完全透明，1 完全不透明。默認值：1
engine	可選字符串。要使用的後端引擎。可能的選項： “arcpy” : Use the arcpy engine for processing. “image_server” : Use the Image Server engine for processing.
gis	可選的地理信息係統。光柵對象的 GIS 。

例子：

# Example Usage

map = gis.map()

# Overlay an image service on the 'MapView' widget
service_url = gis.content.search("my_image_service", item_type="Imagery Layer")[0].url
raster = Raster(path=service_url, gis=gis)
map.add_layer(raster)

# Overlay .tif file present in user's registered fileShare datastore
# (Requires RasterRendering service to be enabled in the active GIS)
raster = Raster("/fileShares/data/Amberg.tif", gis=gis)
map.add_layer(raster)

# Overlay a publicly accesible Cloud-Optimized GeoTIFF
# (Requires RasterRendering service to be enabled in the active GIS)
raster = Raster("https://sentinel-cogs.s3.us-west-2.amazonaws.com/sentinel-s2-l2a-cogs/43/M/BP/2021/6/S2A_43MBP_20210622_0_L2A/B08.tif",
                gis=gis)
map.add_layer(raster)

# Overlay a local .tif file
raster = Raster(r"./data/Amberg.tif")
map.add_layer(raster)

# Overlay a 1-channel .gdb file with the "Orange Red" colormap at 85% opacity
raster = Raster("./data/madison_wi.gdb/Impervious_Surfaces",
                cmap = "OrRd",
                opacity = 0.85)
map.add_layer(raster)

# Overlay a local .jpg file by manually specifying its extent
raster = Raster("./data/newark_nj_1922.jpg",
                extent = {"xmin":-74.22655,
                          "ymin":40.712216,
                          "xmax":-74.12544,
                          "ymax":40.773941,
                          "spatialReference":{"wkid":4326}})
map.add_layer(raster)