Implementing .shp boundaries for batch cropping TIF images with Python, processed with some errors

Causes

Because the .nc original image to tiff format is also processed by the code, there is some error in the data inside after cropping. No reason for this has been found

Ideas

The following code can be compiled with the GUI compiler that comes with ARCGIS, and can also be compiled with the python that comes with ARCGIS PRO, which is also on the Internet to find the relevant code for modification, the source code is now impossible to find, there are a lot of the Internet, and you should have viewed in the process of searching again.

Realization: python code

import arcpy
import os
 
arcpy.CheckOutExtension("spatial")
 
arcpy.env.workspace = "H:/LAB/LABDATA/ECMWF/3/TIF/si10" # 工作目录
 
rasters = arcpy.ListRasters("*", "TIF") # 获取所有的TIF影像
 
mask= "H:/LAB/LABDATA/landsat/shp/caohai2013.shp" # 获取掩膜边界
 
for raster in rasters:
 
    print(raster)
 
    out= "H:/LAB/LABDATA/ECMWF/3/TIF/caohai/si10/"+"clip_"+raster # 输出路径及文件名
 
    arcpy.Clip_management(raster,"#", out, mask, "0", "ClippingGeometry") # arcpy调用ARCMAP里的Clip工具
 
    print("clip_"+raster+"  has done")
 
print("All done")

The results are for reference only, if no errors are reported in the middle, when you see this information as shown, it should be correct!

This blog first appeared on CSDN