99视频精品,国产色在线视频,国产福利片在线观看

在ArcGIS中使用Python裁剪柵格數(shù)據(jù)的流程主要包括以下步驟：

設(shè)置工作空間：使用arcpy.env.workspace來指定數(shù)據(jù)的工作空間。
定義柵格數(shù)據(jù)和矢量數(shù)據(jù)路徑：使用變量來存儲柵格數(shù)據(jù)和矢量數(shù)據(jù)的路徑。
設(shè)置輸出柵格數(shù)據(jù)路徑：同樣使用變量來指定裁剪后的柵格數(shù)據(jù)輸出路徑。
使用arcpy.Clip_management函數(shù)裁剪柵格數(shù)據(jù)：該函數(shù)可以裁剪指定的柵格數(shù)據(jù)，輸入?yún)?shù)包括待裁剪柵格數(shù)據(jù)路徑、裁剪區(qū)域的矢量數(shù)據(jù)路徑、輸出柵格數(shù)據(jù)路徑以及其他裁剪選項(xiàng)（如裁剪幾何類型）。

使用的包：在ArcGIS中，用于操作GIS數(shù)據(jù)和執(zhí)行地理處理任務(wù)的主要包是arcpy。

主要函數(shù)：在裁剪柵格數(shù)據(jù)的過程中，主要用到的函數(shù)：

arcpy.Clip_management()：裁剪柵格數(shù)據(jù)的函數(shù)，它有多個參數(shù)，包括待裁剪柵格數(shù)據(jù)路徑、輸出柵格數(shù)據(jù)路徑、裁剪區(qū)域的矢量數(shù)據(jù)路徑等。

arcpy.management.Clip(in_raster,?rectangle,?out_raster,?{in_template_dataset},?{nodata_value},?{clipping_geometry},?{maintain_clipping_extent})

示例

import?arcpy

#?Set?workspace
arcpy.env.workspace?=?r"D:\Software\DataScience\ArcgisPro_Workspace\數(shù)據(jù)"

#?Define?paths
raster_path?=?"SoilType.tif"
shp_file?=?"珠三角_Project.shp"
out_raster_path?=?"output\SoilType_clip.tif"

#?Clip?raster?using?shapefile
arcpy.Clip_management(raster_path,?"#",?out_raster_path,?shp_file,?"#",?"ClippingGeometry")

當(dāng)然了，我們使用Python更多的是為了更加簡便的操作重復(fù)性的工作。如果是要批量裁剪，可以在以上代碼的基礎(chǔ)上補(bǔ)充循環(huán)處理的代碼：

獲取所有待裁剪的柵格數(shù)據(jù)文件列表：使用 arcpy.ListRasters() 函數(shù)來獲取所有待裁剪的柵格數(shù)據(jù)文件名列表。
遍歷柵格數(shù)據(jù)文件列表并進(jìn)行裁剪操作：使用循環(huán)遍歷所有柵格數(shù)據(jù)文件，每次都使用 arcpy.Clip_management() 函數(shù)來裁剪柵格數(shù)據(jù)。

import?arcpy
import?os

#?設(shè)置工作空間
arcpy.env.workspace?=?r"D:\Software\DataScience\ArcgisPro_Workspace\數(shù)據(jù)"

#?定義矢量數(shù)據(jù)路徑
shp_file?=?"珠三角_Project.shp"

#?設(shè)置輸出文件夾路徑
output_folder?=?r"D:\Software\DataScience\ArcgisPro_Workspace\output"

#?獲取所有tif文件列表
tif_files?=?arcpy.ListRasters("*",?"TIF")

#?遍歷所有tif文件并進(jìn)行裁剪操作
for?raster_path?in?tif_files:
????#?拼接輸出柵格數(shù)據(jù)路徑
????filename?=?os.path.splitext(os.path.basename(raster_path))[0]
????out_raster_path?=?os.path.join(output_folder,?f"{filename}_clip.tif")
????
????#?裁剪柵格數(shù)據(jù)
????arcpy.Clip_management(raster_path,?"#",?out_raster_path,?shp_file,?"#",?"ClippingGeometry")