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)