IDL
江北20190411
。
展开
-
图像客观质量评价指标IDL实现
;用于计算原始影像与经过处理后影像的平均绝对误差 MAEfunction calMAE, image_or, image_deal nDim = size(image_or, /n_dimension) nDim2 = size(image_deal, /n_dimension) if nDim le 2 and nDim2 le 2 then begin valid_data_or = image_or[where(~finite(image_or, /nan))] val原创 2020-06-27 17:26:11 · 835 阅读 · 0 评论 -
cmd调用Python脚本
cmd调用首先win + R 打开cmd窗口,cd到脚本所在的文件夹cd /d E:\pytest\银行系统练习然后输入python banksystem.py banksystem为脚本名IDL中打开cmd调用使用 & 可以一行执行多个命令spawn, 'cd /d E:\pytest\银行系统练习 & python banksystem.py'...原创 2020-05-21 21:24:36 · 466 阅读 · 0 评论 -
MODIS产品质量控制文件使用方法
官方关于产品质量控制的说明(机翻)质量指标 在生产过程中生成的CoreMetadata.0全局属性QA 中的元数据对象以及质量控制(QC)SDS中给出,或者在数据产品的产品后科学和质量检查中给出。CoreMetadata.0中的 QA元数据对象全局属性是AutomaticQualityFlag和ScienceQualityFlag及其相应的说明。根据运行LST算法期间遇到的数据条件,根据规则设置AutomaticQualityFlag。此质量检查标志的设置是完全自动化的。设置它的规则是自由的;几乎所有数原创 2020-05-18 00:51:13 · 7214 阅读 · 14 评论 -
IDL自定义滑动窗口函数
许多地表参数反演算法中是以整景影像中的每个小窗口为单位进行计算。这时候就需要对整景影像划分为数个固定窗口大小影像。以下函数实现输入想要设定窗口的大小,影像列行数,返回一个分块影像的左上角和右下角坐标的查找表数组。有这样一个数组,无论是影像计算还是裁剪深度学习样本都可以适用。原创 2020-05-11 23:04:10 · 1368 阅读 · 0 评论 -
ASTER连续光谱向传感器通道光谱转换
数据说明 ASTER光谱库包含三个子库,分别是 Johns Hopkins University (JHU) 光谱库、Jet Propulsion Laboratory (JPL) 光谱库和 United States Geological Survey (USGS) 光谱库。ASTER 光谱库包含多达 2300 余种典型地物的光谱,包括矿物质、岩石、土壤、植被、水体、冰雪等各类自然地物以及人造地物等类别,其波谱范围从 0.4 μm 到 15.4 μm。问题描述 但由原创 2020-05-09 16:17:46 · 1710 阅读 · 13 评论 -
矩匹配算法影像去噪
2020-5-2-9:19常规的矩匹配算法能起到去除条带噪声的效果,但去噪后影像灰度信息损失过大。IDL关键代码: dim = size(data, /dimensions) row = dim[1] col = dim[0] re_mean = mean(mean(data, dimension=1)) re_stev = mean(stddev(data, dimension=1)) or_mean = mean(data, dimension=1) or_st原创 2020-05-03 21:59:05 · 2688 阅读 · 1 评论 -
像元二分模型计算植被覆盖度
前言关于改进的像元二分模型第一次出现是在李苗苗老师的博士论文中,但她后来发表的一篇期刊文章把这一部分单独摘了出来并做了应用案例,期刊链接在下方。密云水库上游植被覆盖度的遥感估算论文要点计算公式为:Fc =(NDVI -NDVIsoil) (NDVIveg -NDVIsoil)问题的关键在于NDVIveg 与NDVIsoil如何取值。 其中 ,NDVIsoil 为完全是裸土或无植...原创 2020-05-01 00:43:52 · 11650 阅读 · 0 评论 -
6S模型逐像元大气校正MERSI2
前言FY3D-MERSI2 6S大气校正IDL调用6S模型之前完成了6S模型源码中添加MERSI2的光谱响应函数和实现了IDL调用6S模型,然后进行了应用同一参数进行大气校正。这次实现利用IDL调用6S模型建立大气校正系数查找表然后对MERSI2进行逐像元大气校正。由于MERSI2自带的角度信息是逐像元的数据,而我建立的查找表步长比较大,因此在进行校正前,先对角度信息和MODIS AOD...原创 2020-04-28 18:15:34 · 1524 阅读 · 0 评论 -
IDL调用6S模型
进行逐像元的大气校正前需要针对研究区域不同参数建立查找表,在这之前还需要对各参数进行敏感性分析,选择性忽略一部分参数从而提高大气校正效率。 这里省略了针对MERSI2的各参数大气校正的敏感性分析,仅仅为了实现IDL调用6S模型并提取输出结果到txt文档这个过程。 igeom = 0 ;自定义几何条件 vaz = 74 ;卫星方位角 saz = 195 month = 5 ;月份 day = 17 ;日...原创 2020-04-26 18:45:38 · 2996 阅读 · 11 评论 -
IDL矢量/掩膜裁剪影像(ENVI API编程)
前言当我们想要用一个矢量文件去裁剪影像时,可以使用ENVIVectorMaskRaster函数。maskedRaster = ENVIVectorMaskRaster(raster, shp_file)**图中标记红圈处是我们裁剪出的影像,非常小的一块。但是我却发现用ENVI-IDL裁剪出的影像仍然保留着原来影像的行列数,并不像ENVI中Subset Data from ROIS工具一...原创 2020-04-22 21:13:02 · 5688 阅读 · 6 评论 -
IDL实现FY3D MERSI2 250M/1KM数据预处理
结合我自己的需求只处理了Band3.4.15.18.24.25,后续获取光谱响应函数后会加上6s大气校正部分。由于构建250M的GLT文件时间比较长,所以没有添加构建GLT部分,构建GLT的方法可以参考ENVI-IDL技术殿堂官方博客。代码末删除临时文件的方法比较暴力,带有[envitempfile]的文件都会删除,结合实际情况,谨慎使用!原创 2020-04-19 22:09:00 · 3762 阅读 · 6 评论 -
MCTK批处理MODIS L2 swath产品
前言ENVI+IDL版本:5.3MCTK版本:V2.1.10ps:从2.0.0版本开始,MCTK使用完全相同的多线程、严格的swath投影引擎作为EPOC和VCTK,自动执行领结校正和支持所有三种标准插值方法,所以早期keywords中的/bowtie不再有效。而且插值方法keywords不再是7.8.9而是对应着0.1.2。(更多有关版本特点和keywords的详细介绍可参见User’s...原创 2020-04-09 09:06:14 · 3047 阅读 · 11 评论 -
IDL环境栅格文件批量转带地理信息的tiff格式
ENVI+IDL5.3博主在IDL中进行栅格运算后,应用raster.Export, outname, 'tiff'.直接将影像输出,但是丢失了投影等信息。以下代码用于解决在IDL中运算生成的dat文件批量转换成带有地理信息的tiff文件格式。(前提是,dat文件的hdr头文件中必须带有地理信息)pro write_tiff_withmapinfo;*****************...原创 2020-04-04 10:08:44 · 3001 阅读 · 2 评论 -
IDL interpol函数插值波谱响应函数
interpol函数介绍一维数组插值result=interpol(v,x,xout,/lsquadratic,/quadratic,/spline,/NAN)v为输入的数组,x为v中各个元素对应的横坐标值,xout为插值结果数组中各个元素对应的横坐标值,是一个输出的参数(可参考下面代码中的注释帮助理解)。/lsquadratic,/quadratic,/spline是不同的插值算法。/NA...原创 2020-04-04 08:52:06 · 3264 阅读 · 8 评论 -
IDL where函数
where函数介绍函数WHERE()能返回数组中满足指定条件的元素下标调用格式:Result=Where(数组表达式[,count][,Complement=变量1][,/L64][,NCOMPLEMENT=变量2])其中,关键字count返回符合指定条件的元素个数;变量1为不满足条件的数组元素下标;变量2为不满足条件的数组元素个数。where函数举例a=[[1,2],[3,4]]i...原创 2020-04-03 21:12:31 · 2791 阅读 · 0 评论