【IDL代码库】栅格转矢量内部接口

地质少侠

自ENVI 5.2版本开始提供了ENVIClassificationToShapefileTask，可以实现由分类图像转换为shapefile矢量文件的功能。

但是，在使用过程中发现，当栅格图斑数量过多时，此Task效率较低。

于是尝试用旧版本的RTV_doit来实现，结果发现不会输出dbf文件，这样矢量文件也就没有了属性信息。很不开心。跟美国技术支持反馈后，他们给了我一个内部接口，可以实现正常转换。但是当转换的某类像元个数为0时（虽然这种情况极少出现），此接口会报错。

参考代码如下：

1. PRO test_raster_to_vector_doit
   
2.   COMPILE_OPT idl2
   
3.   envi, /restore_base_save_files
   
4.   ENVI_BATCH_INIT
   
5. 
   
6.   file = 'D:\Temp\class.dat' ;输入分类图像
   
7.   out_shpfile = 'D:\Temp\vector.shp' ;输出矢量文件路径
   
8.   idx = [1,2,3,4,5] ;输出类别索引
   
9. 
   
10.   ENVI_OPEN_FILE, file, r_fid=fid
    
11. 
    
12.   ENVI_FILE_QUERY, fid, dims=dims, $
    
13.     class_names=class_names, sname=sname
    
14.   pos = [0]
    
15. 
    
16.   l_name = class_names[idx]
    
17.   evf_file = envi_get_tmp()
    
18.   envi_hourglass
    
19.   envi_check_save, /vector
    
20.   rtv_class_setup, fid=fid, pos=pos, dims=dims, $
    
21.     class_names=l_name, class_ptr=idx,          $
    
22.     out_name=evf_file, in_memory=0,             $
    
23.     output_flag=1, sname=sname
    
24. 
    
25.   ;输出为shapefile
    
26.   evf_id = ENVI_EVF_OPEN(evf_file)
    
27.   ENVI_EVF_TO_SHAPEFILE, evf_id, out_shpfile
    
28. 
    
29.   ;关闭evf文件，并删除
    
30.   ENVI_EVF_CLOSE, evf_id
    
31.   FILE_DELETE, evf_file, /quiet
    
32. END

复制代码

jimi21

谢谢分享！

潜龙

谢谢分享资料！