- 博客(23)
- 收藏
- 关注
RSS订阅
原创 ENVI/IDL混编 批量波段合成 简洁函数
环境：ENVI5.3+IDL8.5代码：此段代码主要用于sentinel-2波段批量合成以合成 band2,3,4,8,8a,11为例：1. 定义函数FUNCTION Layerstack, inpath, inname,outpath,outname;used for layer stacking outpath_name=FILEPATH(outname,ROOT_D...
2020-04-03 22:19:51
2171
1
原创 ENVI/IDL 编程：批量裁剪同一地区的多幅影像
问题描述：通常的批量裁剪方法是使用相同的矢量文件或者roi区域，分别对每景影像裁剪。可采用subset via rois等等方法，但由于多幅影像间（在坐标系相同，分辨率相同的情况下），并非完美配准，因此，裁剪到的影像会存在一行或者一列的差异，使得后续计算无法进行。因此，可采用另一种思路：利用裁剪区域的掩膜图mask（栅格），对多幅影像进行波段合成layer stacking，然后再掩膜掉无效区域。...
2019-02-14 10:09:33
5862
1
原创 Jupyterlab在windows下的安装
使用python 3.6为例，最近发现直接 pip install jupyterlab会报错 pywinty的问题，所以直接先装pywinty, 后装 jupyterlab，具体如下;pip install pywinty==1.1.0pip install jupyterlab==3.0.0无任何报错，即为成功。参考博客：https://blog.csdn.net/weixin_45661908/article/details/120713507...
2022-05-22 16:11:22
806
原创 Windows mmcv-full 1.5.0 详细安装过程
Openmmlab确实是个宝藏仓库，下面简要介绍windows下安装mmcv-full的过程。参考mmcv的官方文档（https://github.com/open-mmlab/mmcv）重要的是问题是 需要GPU cuda, pytorch 版本一致，与mmcv对应才行。步骤如下：1. 安装合适版本的cuda （此处不再介绍，自行百度）：本博主的cuda是10.12. 安装torch, torchvision版本：建议torch 1.6.0以上，本文测试了 torch1.7.0和1.8
2022-05-20 20:36:30
6382
2
原创 Google earth engine python本地环境搭载
配置的主要步骤如下：1. 安装python环境：建议安装发行版anaconda，采用其包管理器 conda，管理各种包，例如 6S包。当然也可以自行配置6S包，不过比较麻烦2.pip install earthengine-api3. 认证（需要科学上网）import ee# 注意需要设置本地代理服务器的端口，打开科学上网软件后，需要找到本地代理的端口号，如10809set http_proxy=http://127.0.0.1:10809set https_proxy=http..
2020-10-16 16:07:40
673
原创 ENVI/IDL 批量裁剪同一地区的多幅影像-第五篇
上一篇文章采用的是regrid的方法，适合将不同分辨率、不同坐标系的影像转到同一坐标系，而对于同一分辨率和相同坐标系的影像而言，会进行采样！这不同于裁剪。对此，下面提供了一种 利用参考影像的地理范围进行裁剪的方法：function y_geoclip, raster_base, raster, respath; Author: Yinxia Cao;Date: 2020.8.19; Define the spatial range of the subRectsubRect = [0,
2020-08-21 18:03:18
1456
原创 ENVI/IDL 批量裁剪同一地区的多幅影像-第四篇
之前采用的函数是envi_layer_stacking_doit，最近发现了一个十分高效的方法 RegridRaster task代码如下：function y_layer_stack,Raster_based,Raster_toresize,respath ;将raster_toresize 影像进行重采样，行列号以及投影与raster_based保持一致 ; Retrieve the grid parameters of the based raster PixelSize = R...
2020-07-01 15:17:27
1131
2
原创 ENVI/IDL 混合编程：批量裁剪影像-第三篇
今天发现之前处理多波段裁剪的代码有点问题，现在更新一下：批量对多波段影像进行掩膜处理的代码如下：;README;used for batch subset(loop pver all images);usage: function_subset(['img1.dat', 'img2.dat'], 'msk.dat', 'c:/data');requirement: ENVI/IDL 5.3 or higher;reference: the reference book of envi/id.
2020-06-19 21:36:56
1007
原创 ENVI/IDL混合编程-批量影像拼接
;README;used for batch mosaicing with sentinel-2 images (band by band);usage: function_mosaic(['/20200131_1/B1','/20200131_2/B1'],'/20200131/B1_mosaic');requirement: ENVI/IDL 5.3 or higher;referen...
2020-04-08 23:38:26
2203
1
原创 U盘格式化
问题描述：两年前，学校发了一个U盘（128G！），作为毕业礼物，甚是高兴，一直没舍得用。最近宅在家里，准备下载几个数据集来看看，随便一下就是多少个G，于是，我搬出了珍藏已久的U盘，怎知该U盘还未格式化。。。试了windows自带的格式化功能，但是失败。。。试了u盘启动盘制作软件，失败。。。解决方法：量产工具（顾名思义，从头开始制作U盘）具体步骤：1. 使用chipgenius查出U盘的...
2020-02-16 21:35:31
483
原创 利用AWS下载SpaceNet数据集
问题描述：下载SpaceNet的数据集，通常需要使用AWS服务器下载解决方法：具体可参考官方文档https://docs.amazonaws.cn/cli/latest/userguide/cli-services-s3-commands.html比较详细的博文https://blog.csdn.net/DynastyRumble/article/details/7664912...
2020-02-16 21:20:19
1590
16
原创 ENVI统计与matlab的区别
同样的数据，利用envi统计和matlab是一样的结果，但是为什么还是会有人觉得不一样呢？？？最近我也遇到了，其实是没有仔细看envi统计的区间，如下图所示：envi统计结果显示 308这个位置有694338个数值，难道意思就是指 等于308的值的个数为这么多吗？显然是错误的，从最左列可以看出，ENVI统计时，是按照4为间距进行统计，说明这其实是直方图统计...
2019-03-09 17:23:33
591
原创 ENVI/IDL 编程：批量使用伪不变特征法进行相对辐射校正--第二篇
重新改写了之前一篇的代码（https://blog.csdn.net/qq_41300185/article/details/87374427）因为之前的算法只能处理输入为：待纠正影像A-单波段，参考影像B-单波段那么，考虑其他可能的情况：1. A-单波段，B-单波段：参考之前的代码（https://blog.csdn.net/qq_41300185/article/details/8...
2019-03-04 16:23:44
833
1
原创 ENVI/IDL编程：批量重命名文件夹
今天，学习到一个批量重命名函数(移动文件的函数，实则是rename files)，如下：FILE_MOVE, oldpath, newpath需要注意的一点是：可能存在命名冲突的情况，在使用的时候，需进一步测试...
2019-03-02 22:57:45
1425
原创 ENVI/IDL 混合编程：批量裁剪影像-第二篇
今天，重新跑了一下，之前写的代码，发现，输入影像只能是单波段的，哈哈哈哈哈哈哈哈，之前忽略了这个问题，但问题不大。主要的波段合成函数：envi_layer_stacking_doit （由于版本限制，如果使用的是最新版ENVI 5.5，可以直接使用高级封装函数，但无奈我的版本是envi5.3 , IDL 8.5， 虽然envi改进了波段合成，使用了另外的函数，如下：重点：envi_l...
2019-03-02 22:00:00
2399
原创 Feature fusion methods in remote sensing
Today, I learned some feature fusion methods as follows (refer to the literature of qinrongjun 3D CD reviews):1. Post-refinement2. Direct feature fusion: considering all channels of information si...
2019-02-26 09:21:26
276
原创 ENVI/IDL编程：从硬盘上删除文件
很危险的函数 FILE_DELETE, filename，/RECYCLE保险的我加了一个recycle,避免彻底删除，将其放在trash中，仅对window有用，unix没有
2019-02-23 19:40:58
1792
转载 ENVI bandmath数据类型
主要的易错点：envi波段运算默认的输出数据的类型 是与 输入相同，意味着只有当你的空间是不变空间（矩阵论）时，结果才会正确。否则，就要实现转换数据类型，使用也很简单，比如，对于影像a, 将其转为16bits，则为fix(a)。(这个fix也是很奇葩了，不应该是int16吗，感觉更加直接）参看文章http://blog.sina.com.cn/s/blog_b27d56590102xa8u...
2019-02-22 22:14:32
3804
1
原创 Matlab编程：笔记一
终于，今天完成了函数整理工作，并完成了测试，主要的易错点有：1.函数路径：不得不吐槽一下，matlab的自定义函数 真的是很随意，随意到你可以进行改动，还没办法使用简单的方法进行固定，进行加密或者编成mex函数。。。很麻烦2.判断矩阵相同：isequal(a,b)：可能一开始会自然地使用a==b或者 a~=b，但很明显是错的（因为这其实是一个逻辑判断语言，最后的值是0（不相等）或者1（相等...
2019-02-20 21:59:30
177
原创 直方图匹配的原理
今天，重读了《遥感原理与应用》，发现里面写的影像辐射增强方法十分详细，且通俗简单。说到直方图匹配，其原理与直方图均衡化十分相似，如下：直方图均衡化的定义：将随机分布的直方图修改为均匀分布的直方图，其实质是对图像进行非线性拉伸，重新分配图像像元值，使一定灰度范围内的像元的数量大致相等。实现方法：1.需要已知图像均衡化后的灰度级m2.计算原影像的累计直方图，将累计直方图平均分成m份...
2019-02-19 22:16:03
3280
翻译 太阳、地球、卫星角度
参考： https://en.wikipedia.org/wiki/Azimuth涉及的角：高度角elevaion angle、天顶角zenith angle、方位角azimuth关系：elevation angle + zenith angle=90 °
2019-02-19 08:51:12
1063
原创 ENVI/IDL编程：批量使用伪不变目标法进行相对辐射校正
原理：参考博客http://blog.sina.com.cn/s/blog_764b1e9d0100pwby.html可直接利用的工具：http://blog.sina.com.cn/s/blog_764b1e9d0102y3b6.htmlIDL矩阵运算参考：http://blog.sciencenet.cn/home.php?mod=space&uid=346157&do...
2019-02-16 19:42:34
1602
3
原创 ENVI/IDL编程：相对辐射校正-直方图匹配
参考博文http://blog.sina.com.cn/s/blog_764b1e9d0102vqws.html 提供的envimosaicraster方法，编写批处理代码效果：测试结果与该博文的完全相同主要的方法：ENVIMosaicRaster数据格式：ENVI（ image与 image.hdr)PRO ZY3histogrammatch ;relative norm...
2019-02-14 16:08:59
4779
空空如也
空空如也
TA创建的收藏夹 TA关注的收藏夹
TA关注的人