麻辣GIS微信平台

更多 GIS 干货

微信关注不错过

IDL调用ENVI函数进行栅格转矢量操作

最近使用IDL用到了栅格数据和矢量数据相互转换的方法,在网上搜索了一下,找到这个方法,值得尝试一下:

COMPILE_OPT idl2
;ENVI调用初始化
ENVI,/restore_base_save_files
ENVI_BATCH_INIT
;打开图像文件  ;
ENVI_OPEN_FILE, file, r_fid=fid
IF (fid EQ -1) THEN BEGIN
ENVI_BATCH_EXIT
RETURN
ENDIF
;
ENVI_FILE_QUERY, fid, dims=dims,nb = nb
;对第一个波段进行计算
pos = [0]
;将灰度值为0的转换为vector
values = [0]
;
l_name = 'zeroValue'
evffile = FILE_DIRNAME(file)+'\img2vec.evf'
;
; 栅格转换为矢量
;
ENVI_DOIT, 'rtv_doit', $
fid=fid, pos=pos, dims=dims, $
IN_MEMORY = LINDGEN(N_ELEMENTS(values)), $
values=values, l_name=l_name, $
out_names=evffile
;evf转换为shp文件
shapefile = FILE_DIRNAME(file)+'\img2vec.shp'
EVF_ID = ENVI_EVF_OPEN(evffile)

ENVI_EVF_TO_SHAPEFILE, EVF_ID, shapefile
ENVI_EVF_CLOSE, EVF_ID

; 退出ENVI
envi_batch_exit

效果如图:

麻辣GIS-Sailor
作者:
GIS爱好者,学GIS,更爱玩GIS。
微博关注
手机阅读
赞赏支持
手机阅读
微信捐助麻辣GIS
微信打赏
支付宝捐助麻辣GIS
支付宝打赏

声明

1.本文所分享的所有需要用户下载使用的内容(包括但不限于软件、数据、图片)来自于网络或者麻辣GIS粉丝自行分享,版权归该下载资源的合法拥有者所有,如有侵权请第一时间联系本站删除。

2.下载内容仅限个人学习使用,请切勿用作商用等其他用途,否则后果自负。

3.如果本文对您有所帮助,请您对我们团队进行 打赏捐助让我们在传播3S的路上可以走得更远,不胜感激。PS:如果你是学生党,请优先把经费用于购买学习资料 以及 与小哥哥/小姐姐约会上:-)

如您有疑问,可在文末留言,或到麻辣GIS QQ群中提问。

QQ群(一):291616564 加入QQ群

QQ群(二):166408035 加入QQ群

QQ群(三):627853279 加入QQ群

QQ群(四):436386604 加入QQ群

QQ群(五):606176554 加入QQ群

QQ群(六):946178380 加入QQ群

QQ群(七):861271808 加入QQ群

新疆互助群:910717627 加入QQ群

微信公众号:malagis,扫描右边二维码直接关注。

已有 3 条评论

  1. 麻辣GIS-aaaaa
    1#
    aaaaa  · 2014-05-19 15:50

    楼主这转出来的矢量文件是多边形文件吧,能不能直接转出线文件来?

  2. 麻辣GIS-高同学
    2#
    高同学  · 2015-08-20 14:39

    博主您好,我用您的方法用二值图像提取的矢量边界用ENVI重新打开,边界在原图像右下方位移两个像素,不知是什么原因?望您给与回复

  3. 麻辣GIS-刘同学
    3#
    刘同学  · 2016-04-29 12:41

    这个做完能回到idl的模式吗

发表评论