技术标签： IDL 经纬度 ENVI ENVI/IDL 波段值
以前，这一步，我都是手动导出的，详细请看另一篇博文：http://blog.chinaunix.net/uid-28490468-id-3462577.html
最近要处理较多的遥感数据，所以要把这一步通过代码自动化。
看了看help文档，发现ENVI支持一种叫IDL的语言，完全没有听说过，可能是我孤陋寡闻了，呵呵，有同感的举手。
有好玩的东西，那么就开始学呗。
打开ENVI+IDLX.X，编写IDL的IDE会自动和ENVI一起打开，这IDE很面熟，像eclipse。
IDL有自己的help文档，通过IDE可以打开。
另外网上也有些资料：
简洁明了： http://www2.geog.ucl.ac.uk/~mdisney/teaching/unix/idl/idl.html
一个国内的论坛，要注册，审核时间要几天（最好快点，呵呵）：http://bbs.esrichina-bj.cn/ESRI/forumdisplay.php?fid=28
各种资料胡乱看了一通，稍微总结几点：
1）IDL的数组是列优先的array[cols, rows]
2）数组是0开始的，和matlab不一样
3）procedure和function都叫routine，其中文件命名
File naming—a program file must be named the same as the main routine.
4）Primary routine must be the last in the file，这个关系到自动编译
5）对routine的参数传递有两种，argument和keyword
具体大家看help文档吧。
但是光有IDL这些基本的知识，根本无法处理像AVIRIS这种超光谱图片，我甚至尝试着去手动解析hdr头文件。
不过在愚蠢之余，试着google大法了一把，结果发现原来ENVI用IDL另外开发了一些API，特意用来处理硕大的遥感图片，
太激动了，呵呵。http://www.360doc.com/content/10/1101/20/472115_65790354.shtml
这些API在IDL的help文档里没有，要在ENVI的help文档里有，真是坑爹。
怎么使用这些API来完成我的任务呢？我是参考：http://www.360doc.com/content/10/0405/21/472115_21741788.shtml
用第一题的做法就能完成我的任务了，不过读取波段值可以用另一个函数，更方便。
具体API
下面是我的实验代码：（要把图片和对应的hdr头文件放在一个文件目录下）
PRO file_test
data_dir = "E:\Data\OilLeakData\data\f100524t01p00r11rdn_b\"
filename = data_dir + "f100524t01p00r11rdn_b_sc01_ort_img"
envi_open_file, filename, r_fid=fid, /no_realize, /no_interactive_query
if (fid eq -1) then return ;;check fid
envi_file_query, fid, ns=ns, nl=nl, nb=nb
pos = lindgen(nb);
xs = 225
xe = 344
ys = 6183
ye = 6193
iproj = envi_get_projection(fid=fid)
oproj = envi_proj_create(/geographic)
print, 'content'
for y = ys, ye do begin
for x = xs, xe do begin
data = envi_get_slice(fid=fid, line=y, pos=pos, $
xs=x, xe=x, /bip) ;; get bands' value
envi_convert_file_coordinates, fid, x, y, xmap, ymap, /to_map
envi_convert_projection_coordinates, xmap, ymap, iproj, oxmap, oymap, oproj ;;Lon = oxmap, Lat = oymap
help, x, xmap, oxmap, data
print, x, y, xmap, ymap, oxmap, oymap, data
endfor
endfor
envi_file_mng, id=fid, /remove ;;Don't forget to close file
END文章浏览阅读936次，点赞22次，收藏26次。React核心基础
文章浏览阅读2k次。linux系统查看磁盘空间的命令是【df -hl】，该命令可以查看磁盘剩余空间大小。如果要查看每个根路径的分区大小，可以使用【df -h】命令。df命令以磁盘分区为单位查看文件系统。本文操作环境：red hat enterprise linux 6.1系统、thinkpad t480电脑。(学习视频分享：linux视频教程)Linux 查看磁盘空间可以使用 df 和 du 命令。df命令df 以磁..._df -hl
文章浏览阅读923次。uses ComObj;var ExcelApp: OleVariant;implementationprocedure TForm1.Button1Click(Sender: TObject);const // SheetType xlChart = -4109; xlWorksheet = -4167; // WBATemplate xlWBATWorksheet = -4167_range[char(96 + acolumn) + inttostr(65536)].end[xlup]
文章浏览阅读2.3k次。上图为任务代码，在任务具体执行的方法中使用，一定要写在方法内使用SpringContextUtil.getBean()方法实例化Spring service类下边是ruoyi-quartz模块中util/SpringContextUtil.java(已改写)import org.springframework.beans.BeansException;import org.springframework.context.ApplicationContext;import org.s..._ruoyi-quartz无法引入ruoyi-admin的service
文章浏览阅读2w次，点赞10次，收藏77次。yum，全称“Yellow dog Updater, Modified”，是一个专门为了解决包的依赖关系而存在的软件包管理器。可以这么说，yum 是改进型的 RPM 软件管理器，它很好的解决了 RPM 所面临的软件包依赖问题。yum 在服务器端存有所有的 RPM 包，并将各个包之间的依赖关系记录在文件中，当管理员使用 yum 安装 RPM 包时，yum 会先从服务器端下载包的依赖性文件，通过分析此文件从服务器端一次性下载所有相关的 RPM 包并进行安装。_centos7配置yum源
文章浏览阅读828次，点赞21次，收藏8次。今天学长向大家分享一个毕业设计项目毕业设计 基于深度学习的抽烟行为检测算法实现(源码分享)毕业设计 深度学习的抽烟行为检测算法实现通过目前应用比较广泛的 Web 开发平台，将模型训练完成的算法模型部署，部署于 Web 平台。并且利用目前流行的前后端技术在该平台进行整合实现运营车辆驾驶员吸烟行为检测系统，方便用户使用。本系统是一种运营车辆驾驶员吸烟行为检测系统，为了降低误检率，对驾驶员视频中的吸烟烟雾和香烟目标分别进行检测，若同时检测到则判定该驾驶员存在吸烟行为。进行流程化处理，以满足用户的需要。
文章浏览阅读3.7k次，点赞3次，收藏14次。多个定时器同步触发启动是一种比较实用的功能，这里将对此做个示例说明。_stm32 定时器同步
文章浏览阅读348次。出处 : http://www.cnblogs.com/mythou/p/3187881.html本来想分析AppsCustomizePagedView类，不过今天突然接到一个临时任务。客户反馈说机器界面的图标很难点击启动程序，经常点击了没有反应，Boss说要优先解决这问题。没办法，只能看看是怎么回事。今天分析一下Launcher启动APP的过程。从用户点击到程序启动的流程，下面针对WorkSpa..._回调bubbletextview
文章浏览阅读6.2k次。Ubuntu 12 最快的两个源 个人感觉 163与cn99最快 ubuntu下包过慢 1、首先备份Ubuntu 12.04源列表 sudo cp /etc/apt/sources.list /etc/apt/sources.list.backup （备份下当前的源列表，有备无患嘛） 2、修改更新源 sudo gedit /etc/apt/sources.list （打开Ubuntu 12_un.12.cc
文章浏览阅读5.8k次，点赞6次，收藏86次。1.思路（1）动态添加路由肯定用的是addRouter，在哪用？（2）vuex当中获取到菜单，怎样展示到界面2.不管其他先试一下addRouter找到router/index.js文件，内容如下，这是我自己先配置的登录路由现在先不管请求到的菜单是什么样，先写一个固定的菜单通过addRouter添加添加以前注意：addRoutes()添加的是数组在export defult router的上一行图中17行写下以下代码var addRoute=[ { path:"/", name:"_vue动态路由权限
文章浏览阅读8.9k次。 关键词： JSTL 之变量赋值标签 /* * Author Yachun Miao * Created 11-Dec-06 */关于JSP核心库的set标签赋值变量，有两种方式: 1.日期" />2. 有种需求要把ApplicationResources_zh_CN.prope
文章浏览阅读3.1k次，点赞3次，收藏2次。1.1ZY5621概述ZY5621是VGA音频到HDMI转换器芯片，它符合HDMI1.4 DV1.0规范。ZY5621也是一款先进的高速转换器，集成了MCU和VGA EDID芯片。它还包含VGA输入指示和仅音频到HDMI功能。进一步降低系统制造成本，简化系统板上的布线。ZY5621方案设计简单，且可以完美还原输入端口的信号，此方案设计广泛应用于投影仪、教育多媒体、视频会议、视频展台、工业级主板显示、手持便携设备、转换盒、转换线材等产品设计上面。1.2 ZY5621 特性内置MCU嵌入式VGA_vga转hdmi带音频转换器,转接头拆解