tyut_qc-优快云博客

原创关于哨兵2-LIC级数据辐射定标几个参数的获取

Sentinel-2L1C级数据是经正射校正和亚像元级几何精校正后的大气表观反射率（TOA）产品，在我们需要使用辐射率（Radiance）数据时需要一些转换，所以该数据的辐射定标不在是简单的乘以一个Gain，加上一个Offset那么简单。而是利用公式关系反算辐射率数据。其详细计算公式参见(ENVI官博)http://blog.sina.com.cn/s/blog_764b1e9d0102y...

2019-03-19 15:50:24 8004 14

原创 IDL之创建隐藏式窗体（进度条）

WIDGET_CONTROL这个方法中有个map关键字，将map设为0即隐藏主窗体。下面以IDL中的Itools进度条为例，我们知道进度条需要依赖一个WIDGET_BASE顶级容器，但一般不希望该窗体显示。pro IDLProgressbar tlb=WIDGET_BASE(xsize=400,ysize=400) ;窗口居中显示 DEVICE,get_screen_size=...

2019-03-19 14:58:42 1304

原创 IDL之多波段运算(Map Variables to Multiband Images)

IDL波段运算，常规用法就是选择要运算的几个波段，经过运算后，输出为单波段结果。可我们经常也遇到这样的场景，如大气校正后，结果被放大10000倍，需要将每个波段都除以10000，再合成所有波段。ENVI中可以使用Map Variables to Multiband Images这个功能，可以选择几个波段或全部的文件作为运算波段，结果输出为对应多波段影像。 ENVI中有两个函数，即...

2018-11-01 10:22:54 5776 4

原创 DevExpress TreeList 节点查找(类Word导航栏标题查找逻辑)

前段时间写项目文档，发现Word导航栏标题的节点搜索还是蛮有意思的，想着自己也实现一下这样的逻辑。下图就是Word节点搜索时的图。仔细观察搜索时的行为，可以分析需要实现以下功能：1、节点递归查找2、颜色标注匹配节点3、若匹配节点其父节点尚未展开，则需颜色标注该父节点，表示其子项存在匹配项4、展开节点时，若节点并非是匹配节点，则取消颜色标注5、收回节...

2018-10-31 16:38:18 3256

原创 IDL矢量裁剪影像(解决复杂多边形、空心、岛等裁剪问题)

IDL矢量裁剪影像的代码已经有不少人写了，之前的项目中，同事也是参考了ENVI官博中的写法。后来，客户那边反馈裁剪的不准，才发现了问题，就是这些代码用简单的矢量裁剪影像没有问题，但客户那边的是河流的矢量(如下图)，中间有些挖空，在转换为ROI的时候，往往不能正确裁剪。矢量代码裁剪结果 ...

2018-08-28 18:00:01 3707 1

原创利用IDL输出影像快视图、拇指图

大家做影像管理方面的业务，一定会需要针对影像生成快视图、拇指图，代码方面用的最多的应该是GDAL库了，其实如果安装了ENVI，利用IDL生成的话也很方便。下面我就测试数据(GF2（高分二号）影像)生成了快视图和拇指图，主要函数就是 stretch_doit、congrid。其中，快视图最大尺寸1024，拇指图最大尺寸128，影像按照321波段真彩色进行拉伸显示，输出格式为jpg,大家可以...

2018-07-03 10:25:49 4374 2

原创 Arcgis Engine 实现单波段栅格图层色带渲染(包括色带反转)

最近有同事问我ArcMap中的色带反转怎么实现，刚好我在整理以前的项目，找到了一些源码备忘。先说下思路，简单来说就是得到栅格图层的IRasterRenderer,转换为IRasterStretchColorRampRenderer,在将这个接口转换为IRasterStretch,IRasterStretch中有个Invert的方法可以控制色带是否反转...

2018-06-25 13:43:21 6569 1

原创 HSV与RGB色彩模型互相转换

HSV与RGB是ArcGIS中用的最多的两种色彩模型，这两种色彩模型基本信息就不做过多介绍了。下面是互相转换的代码1、HSV转RGB Arcgis Engine中可以将 Symbol中的Color先转为 IHsvColor这个接口，在这个接口中获取H、S、V三个值，我这里方便计算，将S、V值除以了100。 IHsvColor pHsvColor = pSet...

2018-06-25 13:15:20 2198

原创 Arcgis Engine 自定义栅格分类渲染(排除不想要的值)

在做栅格分类渲染的时候，大家有可能想排除一些不想要的值，如背景值或异常值，在ArcMap中对应的功能就是Exclusion，如下图。那么，二次开发的时候怎么写呢？，其实，只要利用IRasterDataExclusion这个接口就可以实现，根据官方开发文档查了下分级的类RasterClassifyColorRampRendererClass所对应的接口就一目了然了。 ...

2018-06-11 17:22:51 2198

原创 IDL调用ENVI-FLAASH大气校正异常退出解决办法

最近把IDL编写的FLAASH大气校正模块嵌到了项目里，今天几类卫星的测试数据收到了，就着手测试了大气校正功能是否正常。测试了GF1(高分一号)、GF2(高分二号)、ZY3(资源三号-01a)数据都没有问题，可是到了BJ2（北京二号）大气校正界面弹出了几秒就莫名其妙的关闭了，百思不得其解，于是记录下排除错误和解决问题的方法。在ENVI中打开BJ2辐射定标数据，发现BJ2缺少...

2018-05-18 15:25:44 4860 6

原创 IDL调用FLAASH大气校正，不启动ENVI，显示自带进度条

自从ENVI官博放出了IDL调用FLAASH大气校正的代码，相信不少人已经在批处理业务或者系统建设中已经用上了FLAASH大气校正模块，当然那段代码针对不同的ENVI/IDL版本还是有点问题，（大家有问题可以一起讨论下）。emmm，回归正题。因为ENVI5.0后，发布了ENVI面向对象开发新方式，启动ENVI就是 e=envi(/headless)，/headless的意思就是不启动ENVI软...

2018-05-16 15:33:22 4418 25

原创关于ENVI-FLAASH大气校正报 .. bin\flaash\scriptfile.006系统找不到指定的路径的错误

完整的错误信息如下图：本人的软件环境为ENVI5.3.1、IDL8.5，之前在ENVI5.1、IDL8.3下同样遇到过这样的错误，只不过由于事情比较多，都是让遇到这个问题的同事把ENVI卸载掉，重新安装。最近由于项目需要，需要利用IDL调用FLAASH模块，再次遇到了这个问题。。。从这个问题出发，既然报找不到script.006这个文件，我就到这个目录下找...

2018-05-16 14:55:46 3080

tyut_qc的博客