ENVI二次开发时的波段运算表达式语句合理性检测方法与ENVI检测方法探讨

合集 - 地理遥感生态网地理数据集5(100)
1.利用SuperView数据基于ENVI深度学习进行建筑物提取2.ENVI发布深度学习产品：ENVI Deep Learning Module3.澳大利亚市政局利用ENVI深度学习进行树木调查4.ENVI软件系列产品体系结构5.ENVI新机器学习：ENVITask 使用说明6.ENVI新机器学习：随机森林遥感分类工具操作手册7.ENVI新机器学习：异常探测分类工具操作手册8.新机器学习遥感分类工具9.ENVI深度学习2.0新特性10.IDL封装百度翻译API实现自动翻译和语种识别11.使用IDL进行氚羽流监测（五）—绘制污染物扩散带等值线图12.使用DL进行氚羽流监测（四）——绘制预期氚衰变图13.使用IDL进行氚羽流监测（三）——绘制含水层相对流速14.使用IDL进行氚羽流监测（二）——绘制3D地形和地下水位15.使用IDL进行氚羽流监测（一）——对研究区和环境进行建模16.【IDL】IDL8.8 安装中文语言包方法17.【IDL】IDL 8.7.1机器学习框架18.【IDL】IDL 8.7.1 新特性19.【IDL】获取參边形的最小外接矩形20.【IDL】图形保存图片时文字变黑块的临时解决方法21.【IDL】控制台使用Alt+/代码补全快捷键的方法22.【IDL】IDL 8.6新功能23.【IDL】 程序发布exe清晰步骤24.【IDL） 64位IDL中调用32位功能的方法25.【IDL】 去掉黄色警告图标的方法26.【IDL】 重写HDF文件中数据集的方法27.【IDL】无法使用格式快捷键的解决办法（Ctrl+ Shift+F）28.【IDL】IDL 8.5 中的颜色选择对话框29.【IDL】IDL 8.5 中的 Python Bridge30.【IDL】 开发环境（工作台）的语言切换万法31.【IDL】IDL8.4 文件监视系统体验32.【IDL】 批处理遇到错误时处理方法33.【IDL】开发游戏"2048"34.【IDL】IDL 8.4新特性介绍35.【IDL】IDL中乱码问题的解决方法36.【IDL】 自动构建泰森多边形（Voronoi Polygon）37.【IDL】 IDL 8.3 的新特性与新功能38.【IDL】 IDL 8.3 新功館39.2013开发竞赛ENV/IDL组作品欣赏——“伪装目标识别系统”40.【IDL】多元线性回归计算方法41.【IDL】 关键字继承42.【IDL】调用6S.exe生成查找表源码43.【IDL】使用WIDGET_TIMER进行动画控制44.【IDL】路径相关函数和方法汇总45.【IDL】绘图组件打开鼠标拖拽文件的方法46.【IDL】遠取文本中三維数組的方法47.【IDL】 IDL 8.2.2（SP2）新特性48.【IDL】IDL程序发布方法总结49.【IDL】宏命令菜单（Macros）的用处50.【IDL】 使用SPAWN启动路径中带有空格的文件51.【IDL】 实现双重排序功能（更新降序）52.【IDL】 C#调用ENVL_FX_SEGMENTONLY_DOIT需要注意的地方53.【IDL】 IDL 8.2 新增视频对象IDLffVideoWrite用法示例54.【IDL】生成等差/等比数列的方法55.【IDL】 判断双精度变量是合EQ56.【IDL】 File_Search函数用法详解57.【IDL】坐标轴中文标题乱码问题58.【IDL】C#调用IDL编写界面（含事件）的问题59.《IDL程序设计—数据可视化与ENV二次开发》目录60.《IDL程序设计—数据可视化与ENVI二次开发》已由高等教育出版社正式出版61.IDL Viewer 应用模版介绍（版本号更新为1.2）62.【IDL】日期和时间函数63.【IDL】数组求交集、并集和差集的函数64.【IL】IDL 8.2 新特性65.【IDL】 气象可视化程序配置（idL_guide5）66.【转】几个IDL在线教程网址67.【IDL】医学中IDL的应用68.【IDL】气象中IDL的应用69.【IDL】SIZE函数使用方法与参数含义详解70.【IDL】 DICOM扩展模块71.【IDL】几何图形数学运算函数72.【IDL】 弹出界面事件响应示例程序73.【IDL】鼠标信息显示界面（同理可实现鹰眼图等）74.【IDL】 IDL与C#混合编程技术75.【IDL〕编写程序启动界面76.【IDL】获取n*n数组的对角线元素的方法77.【IDL】同时显示不同分辨率图像的方法78.【IDL】帮助科学家了解太阳系的形成79.【IDL】 Polyval函数（多项式的估值）80.【IDL】 数组下标转换一维坐标索引函数：array_indices_reverse81.【IDL】数据处理时遇到了NAN怎么办
82.ENVI二次开发时的波段运算表达式语句合理性检测方法与ENVI检测方法探讨
83.【IDL】编写高效率的IDL程序84.【IDL】实现ENVI默认显示效果（线性拉伸算法）85.【IDL】双视图关联示例86.【IDL】图像真彩色、假彩色和伪彩色显示87.【IDL】国外学习与开源网站汇总88.【IDL】 如何保存IDL绘制的图像或图形89.【IDL】学习书籍《Modern IDL》90.【IDL】 读写Excel的方法总结91.【IDL】 绘制颜色棒 （colorbar）92.【IDL】 IDL8.1新功能介绍93.【IDL】 讨论浮点运算精度问题 （IDL与Excef和Matlab相比）94.【IDL】 图像显示程序（直接图形法与对象图形法混合显示）95.【IDL】 IDL Advanced及其详细功能介绍96.【IDL】自定义鼠标光标样式（直接图形法与对象图形法）97.【IDL】多条件判断的另类写法（case代替i）98.2014年IDL培训班素材包及录屏共享99.IDL中开发游戏“2048”100.ENVI 6.0自助试用许可开放申请
收起

ENVI的波段运算功能很强大，若输入表达式中有错误则无法add。这个检测是怎么实现的呢？如果在IDL下进行ENVI二次开发的时候，如何实现这个功能。

首先启动ENVI+IDL，进行波段运算，输入b1*a，然后点击Add。



IDL的控制台上会输出如下信息：

% Variable is undefined: A.

ENVI>

由此，联想到ENVI是否进行了一个测试调用，根据调用的返回值来判断该表达式是否合法。在IDL下进行测试：

ENVI> str1 ='b1*a'

ENVI> varname = Var_Extract(str1)

ENVI> isWrite = Execute(varName+"=indgen(5,5)")

ENVI> print,isWrite

      1

ENVI> isWrite = Execute(varName+"=" + Str1)

% Variable is undefined: A.

ENVI> print,isWrite

      0

该错误提示与ENVI的波段运算错误提示是对应的。应该就是用的这个Execute返回值来判断表达式是否正确。

其中调用到的函数来自YangHS，如下：

Function VarFind,Str

 Compile_opt strictarr

 

  Len = StrLen(Str)

  RetValue = 'b'

  RefArr = Strtrim(Indgen(10),2)

 if Len eq 1 then Return,RetValue

 for i=1,Len-1 do begin

   Char = StrMid(Str,i,1)

   if (Where(RefArr eq Char))[0] eq -1 then break

   RetValue = RetValue+Char

 endfor

 ;

 Return,RetValue

End

Function Var_Extract,ExpStr

 Compile_opt strictarr

 

  ExpStr = STRLOWCASE(ExpStr)

  StrArray = StrSplit(ExpStr,'b',/extract)

  StrArray = 'b'+StrArray

 ;

 for i=0,N_Elements(StrArray)-1 do begin

   StrArray[i] = VarFind(StrArray[i])

 endfor

 ;

  LenArr = StrLen(StrArray)

  index = Where(LenArr ne 1,count)

 if count gt 0 then RetArr=StrArray[index] else RetArr=''

 ;

  uniqIndex = Uniq(RetArr)

  RetArr = RetArr[uniqIndex]

 Return,RetArr

End

那么继续探讨下表达式变量的初值设置，编写一简单的波段运算函数test_bandmathVar。

function test_bandmathVar,inVar

 help,inVar,out = outD

 print,outD

 print,inVar

 return,inVar

end

编译后在波段运算中输入test_bandmathVar(b1)：

控制台输出如下：

INVAR          FLOAT    = Array[5, 5]

     36.0000     33.0000     53.0000     2.00000     119.000

     100.000     17.0000     20.0000     168.000     125.000

     97.0000     126.000     248.000     17.0000     2.00000

     224.000     76.0000     81.0000     237.000     38.0000

     217.000     140.000     196.000     19.0000     109.000

INVAR          FLOAT    = Array[5, 5]

     144.000     192.000     86.0000     30.0000     235.000

     65.0000     32.0000     131.000     210.000     201.000

     186.000     120.000     19.0000     80.0000     170.000

     147.000     139.000     170.000     112.000     81.0000

     141.000     105.000     175.000     247.000     34.0000

INVAR          FLOAT    = Array[5, 5]

     143.000     41.0000     121.000     125.000     51.0000

     45.0000     12.0000     114.000     139.000     197.000

     147.000     102.000     212.000     241.000     225.000

     141.000     152.000     50.0000     51.0000     34.0000

     176.000     223.000     151.000     131.000     247.000

INVAR          FLOAT    = Array[5, 5]

     61.0000     164.000     64.0000     246.000     82.0000

     255.000     29.0000     114.000     47.0000     113.000

     85.0000     142.000     94.0000     229.000     144.000

     250.000     30.0000     93.0000     190.000     100.000

     163.000     50.0000     139.000     253.000     144.000

由上面的输出看来ENVI下初步给定的是5*5的浮点型随机数组，随机变量构建代码如下：

ENVI> s = float(fix(randomu(aa,5,5)*255))

ENVI> print,s

     188.000     146.000     175.000     95.0000     116.000

     233.000     227.000     152.000     159.000     110.000

     58.0000     188.000     180.000     230.000     160.000

     152.000     254.000     31.0000     41.0000     18.0000

     127.000     198.000     11.0000     66.0000     195.000

ENVI> s = float(fix(randomu(aa,5,5)*255))

ENVI> print,s

     181.000     27.0000     4.00000     238.000     205.000

      34.0000     191.000     167.000     60.0000     161.000

     233.000     90.0000     112.000     81.0000     149.000

     210.000     236.000     66.0000     13.0000     188.000

     65.0000     63.0000     135.000     16.0000     129.000

ENVI> s = float(fix(randomu(aa,5,5)*255))

ENVI> print,s

     22.0000     198.000     172.000     66.0000    0.000000

     206.000     210.000     101.000     53.0000     86.0000

     65.0000     31.0000     117.000     167.000     50.0000

     65.0000     175.000     17.0000     192.000     41.0000

     36.0000     67.0000     97.0000     196.000     88.0000

ENVI>

 

 

posted @   地理遥感生态网平台  阅读(42)  评论(0)  编辑  收藏  举报
相关博文：
·  【IDL】实现ENVI默认显示效果（线性拉伸算法）
·  《IDL程序设计—数据可视化与ENV二次开发》目录
·  ENVI二次开发时的波段运算表达式语句合理性检测方法与ENVI检测方法探讨
·  【ENVI入门系列】21.波段运算与波谱运算工具
·  ENVI二次开发综述与代码分享
阅读排行：
· winform 绘制太阳，地球，月球 运作规律
· AI与.NET技术实操系列（五）：向量存储与相似性搜索在 .NET 中的实现
· 超详细：普通电脑也行Windows部署deepseek R1训练数据并当服务器共享给他人
· 【硬核科普】Trae如何「偷看」你的代码？零基础破解AI编程运行原理
· 上周热点回顾（3.3-3.9）
点击右上角即可分享
微信分享提示
AI IDE Trae