博客
关于我
强烈建议你试试无所不能的chatGPT,快点击我
ENVI\IDL 重采样 栅格单元大小设置
阅读量:5300 次
发布时间:2019-06-14

本文共 2341 字,大约阅读时间需要 7 分钟。

在ENVI\IDL 下图像重采用是用Resize_doit 函数。该函数能用来改变图像的大小和对图像重采样。

该函数有个参数RFACT。

RFACT

Use this keyword to specify a two-element array holding the rebin factors for x and y. The values of RFACT reflect the IDL convention for resizing data. A value of 1 does not change the size of the data. Values less than 1 cause the size to increase; values greater than 1 cause the size to decrease.

从以上可以看出,当Refact的Value大于1 到时候就是 放大图像,小与1就是缩小图像。

通过官方给出的实例:

PRO EXAMPLE_RESIZE_DOITcompile_opt IDL2 ; First restore all the base save files.envi, /restore_base_save_files ; Initialize ENVI and send all errors; and warnings to the file batch.txtenvi_batch_init, log_file='batch.txt'; Open the input fileenvi_open_file, 'can_tmr.img', r_fid=fidif (fid eq -1) then begin   envi_batch_exit   returnendif ; Set the POS keyword to process all; spectral data. Output the result; to disk.envi_file_query, fid, dims=dims, nb=nbpos = lindgen(nb)out_name = 'testimg' ; Perform the resize calculation.; Make the output image twice as; large in both X and Y. Use; bilinear interpolation.envi_doit, 'resize_doit', $   fid=fid, pos=pos, dims=dims, $   interp=1, rfact=[.5,.5], $   out_name=out_name, r_fid=r_fidEND

由于RFACT=[*,*]是浮点数表示,其重采样比率的计算方法是:采样后分辨率/采样前分辨率。如果原图像是20m,要采样成100m分辨率,RFACT=[100/20,100/20],即RFACT=[5,5],这个好理解。

但是个问题: 1. 如果在不知道原始图像的分辨率的基础上,要定量的设置重采样到30m的时候怎么设置倍数。

问题的解决办法就是在代码中加入获取原始图像的像元大小的代码!

proj=envi_get_projection(fid=fid,PIXEL_SIZE=ps,units=units)

其中ps 为原图像的像元大小。

则可以将帮助中的程序更改如下:

PRO EXAMPLE_RESIZE_DOIT, pixelSizeAfterResamplecompile_opt IDL2 ; First restore all the base save files.envi, /restore_base_save_files; Initialize ENVI and send all errors; and warnings to the file batch.txtenvi_batch_init, log_file='batch.txt' ; Open the input fileenvi_open_file, 'can_tmr.img', r_fid=fidif (fid eq -1) then begin   envi_batch_exit   returnendif; Set the POS keyword to process all; spectral data. Output the result; to disk.envi_file_query, fid, dims=dims, nb=nbpos = lindgen(nb)out_name = 'testimg' ; Perform the resize calculation. ;Make the output image twice as; large in both X and Y. Use; bilinear interpolation.envi_doit, 'resize_doit', $   fid=fid, pos=pos, dims=dims, $   interp=1, rfact=[pixelSizeAfterResample/ps,pixelSizeAfterResample/ps], $   out_name=out_name, r_fid=r_fidEND

 

 

 

 

转载于:https://www.cnblogs.com/myyouthlife/archive/2012/07/01/2571785.html

你可能感兴趣的文章
VC++ 模块与资源分离
查看>>
CWorkBooks、CWorkBook、CWorkSheets、CWorkSheet、CRange
查看>>
Ocx控件注册不成功?可能是tlb文件导致~
查看>>
进制转换
查看>>
JS鼠标事件大全
查看>>
css的动画和过渡
查看>>
Tomcat最大连接数问题
查看>>
发布功能完成
查看>>
Java 学习笔记 二 -- JDBC事务
查看>>
修改jar包中的XML配置文件
查看>>
【计算机网络】 网络体系结构分类: 客户机/服务器体系和P2P
查看>>
BZOJ 1096: [ZJOI2007]仓库建设 动态规划 + 斜率优化
查看>>
From:51JS,JS牛人的作品
查看>>
Navicat for MySQL中文破解版免费下载
查看>>
Fragment
查看>>
flask-session源码流程
查看>>
dede文章内容页增加视频文件
查看>>
ie浏览器兼容性快速处理小招
查看>>
[转] maven打包可运行的fat-jar的简单方法
查看>>
C# - 非中断(正常)模式下的调试
查看>>