一个RS农民工

前言  考试考image的函数的keywords的中文解释，下午也用了用，虽然不大习惯，但是和MATLAB比一比读一读英文再实践实践还是可以理解的。语法  IDL在doc里面留下的语法是这样的:graphic = IMAGE(DATA[,X,Y] [,Keywords=value],[,Properties=value,][,Map Properties=value,][,Map Grid Properties=value]])graphic = Image(Filename[,Keywords=

前言   突然想起来电脑最近网卡不大好，怕回去写不出来，就把第四篇也一带写了，最后辐射定标的回学校再写吧，这个程序属于稍微困难点的程序，需要好好构思构思，而且我对读文件的操作也不大熟，尽量讲清楚！ 需求分析Created with Raphaël 2.2.0开始读取*.hdr文件line.indexof=’...‘wavelenth读取*.sli文件绘制波谱曲线结束yesno这个程序稍微复杂一些，先看看绘制出的效果图吧！二话不说开始干活！1.读取hdr文件，读出波长等各种信息  

前言   发现这个好像和我之前写的程序重了，之前写考核的时候密度分割部分直接抄的课上的程序，这次看懂了，挺简单的 。 需求分析Created with Raphaël 2.2.0开始b3,b4波段计算NDVI密度分割密度分割图结束这个代码有点太过简单了…不废话了，干活！1.预定义颜色索引矩阵，读入近红和红数据，计算NDVI dst = [[-1.0, 0.0, 255, 127, 36], $ [.0, 0.2, 255, 165, 79], $ [0.2,

前言   明天晚上坐火车回学校了，电脑可能就收起来了，只好白天写了，争取写博客的时候教学性更强一点让人更加容易看懂，毕竟这几篇文章可能对复习IDL考试有一些帮助，主要是有客户群体，哈哈！不过程序是应用导向的，写程序还是应该明白自己要干什么再写，要不然写起来很费时费劲，更别说考试的时候写了，背?头给你拧下来你也背不出来。 需求分析Created with Raphaël 2.2.0开始读入二进制文件创建颜色索引表并显示创建标签并显示结束  所以开始按照需求部分来按块写代码，里面的内容尽量也写详细一

前言   由于IDL要考试所以记录下5个程序和他们的语法过程，也权当是复习了，其实写过一遍之后一下就记住了，尽量注解详细一些，主要是我的老巢被人发现了，只能方便你我他了。 需求Created with Raphaël 2.2.0程序开始读入数据x,y绘制散点图计算误差根据误差绘图散点图程序结束  首先，定义好文件名，并且利用read_csv的function读入csv，读入的data是一个结构体，即struct，所以需要点号去访问，是个2个数组组合成的结构体看下图fn = 'C:\Users\

前言   前天在写IDL程序的时候遇到了NAN，虽然实际问题没能解决，但是大概已经有点眉目了。顺便记下遇到NAN，INF这些可能的操作。  NAN是not a number的缩写，在和浮点数比较的时候，它会出现警告，当时做的时候觉得反而是警告更好，NAN处理起来挺麻烦的。MATLAB处理这些异常时候挺熟练的，IDL度娘了半天没找到，只能自己读文档了，好在找到了文档。  inf就是无穷，也可以进行处理。创建这四种变量a = !VALUES.F_NAN ；NANIDL> b ...

前言  正巧IDL实验课考核的作业是利用4个Function和主Pro过程写一个遥感图像处理的代码，以前一直想着能不能计算NDVI和植被覆盖度(VFC)用IDL写出来，因为老是总强调处理遥感影像比较正统，一定程度上是这样说，因为IDL是Harris公司的，然后envi的底层是IDL，IDL的底层是cpp，写完了整个代码大概是140行，当然这个是仅仅利用2幅影像计算的，其实看大多数人的需求是长时间序列的，其实长时间序列只是进行了一个最大值滤波，之后再是几幅图像的计算。语法function name,..

前言  千呼万唤始出来，一直运行不成功原来是图像的问题，打不开第二波段的图像，人都给你整傻咯，IDL是真心的不好用而且还要和ENVI联动，就更蠢了，辐射定标ENVI自由的功能它不香吗，可是老师上课写了考试还考，只能含泪写完了它。流下了技术人debug一天的眼泪。 需求分析Created with Raphaël 2.2.0开始读取MTL文件gain 和 baisgain * DN + bais输出图像结束虽然流程看起来比较简单，但是中间的文本的操作真的是复杂到哭，真的是含泪写完了博客。1.打开tx