IDL使用shp文件裁剪tiff影像

IDL可以使用ENVI库中的函数来读取和处理遥感数据。如果你的shp文件是一个多边形，你可以使用ENVI的“ROI Mask”功能来创建一个掩膜，然后使用“Extract ROI”功能来提取掩模内的像素值。以下是一个简单的IDL代码：

; 定义shp文件名和tif文件名
shp_file = 'path/to/your/shapefile.shp'
tif_file = 'path/to/your/tiff/image.tif'

; 打开shp文件并获取多边形坐标
shp = SHP_OPEN(shp_file)
poly = SHP_READPOLY(shp, 0)

; 打开tif文件
envi_open_file, tif_file, r_fid=r_fid, $
                   r_var=r_var, $
                   quiet=1, $
                   /no_tiling

; 创建掩膜
mask = ENVIMaskRaster(r_var, poly, $
                      MASKNAME='MyMask', $
                      /INVERT_MASK)

; 提取掩膜内的像素值
roi = ENVIRasterSubset(r_var, $
                       SUB_RECT=mask.BBOX, $
                       MASK=mask, $
                       /NO_COPY)
                       
; 关闭文件
envi_file_query, r_fid, /close

这段代码首先打开shp文件，读取多边形坐标。然后，它打开tif文件，并使用ENVI库中的“ENVIMaskRaster”函数创建一个掩膜。掩膜是一个与原始图像大小相同的二进制数组，其中0表示像素在多边形外，1表示像素在多边形内。然后，它使用“ENVIRasterSubset”函数提取掩模内的像素值。最后，它关闭tif文件。

请注意，这仅适用于单波段tiff图像。如果你有多个波段，你需要分别处理每个波段。此外，如果你的shp文件包含多个多边形，你需要将代码修改为循环读取每个多边形并创建相应的掩模和ROI。