- 博客(25)
- 资源 (1)
- 收藏
- 关注
RSS订阅原创 GEDI学习笔记 二：通过Subsetter.py代码对GEDI L1 L2级别光斑数据进行处理：属性筛选、区域裁剪和矢量格式转换（最后附上代码链接）
GEDI（Global Ecosystem Dynamics Investigation）卫星是美国国家航空航天局（NASA）开展的一项地球观测任务。GEDI卫星的主要目标是获取全球各地森林和陆地生态系统的三维结构和动态信息，以帮助科学家更好地了解地球生态系统的变化和功能。
2023-06-18 10:10:33
713
原创 GEDI学习笔记 一： 在EarthData上批量、快速下载GEDI数据方法（速度几M每秒）
GEDI（Global Ecosystem Dynamics Investigation）卫星是美国国家航空航天局（NASA）开展的一项地球观测任务。GEDI卫星的主要目标是获取全球各地森林和陆地生态系统的三维结构和动态信息，以帮助科学家更好地了解地球生态系统的变化和功能。在EarthData上快速下载GEDI数据方法
2023-05-28 09:47:33
3134
4
原创 Sentinel-2学习笔记 一：利用ENVI软件对Sentinel数据进行波段组合（本实验基于ENVI5.1）
一、利用ENVI软件对Sentinel数据进行波段组合（本实验基于ENVI5.1）前提：已经通过SNAP软件对Sentinel数据进行大气校正和重采样。Sentinel数据如下：启动ENVI后打开LayerStacking工具：2.导入数据3. 依次按照波段顺序导入后，投影信息选择WGS84，重采样方法选双线性：选择输出文件夹4 等待一段时间5. 结果...
2021-09-25 20:09:25
1058
1
原创 Sen2Cor+SNAP处理Sentinel数据
一、下载安装Sen2Cor1.在官网上下载Sen2Cor后，解压，并将其放到C:\Users2.打开cmd命令提示符（Window10系统，电脑左下角右键点击如下）3. 输入cmd4. 输入命令：cd+空格+Sen2Cor路径进行安装。安装之后，输入：L2A_Process --help，出现以下窗口则显示成功：二、Sen2Cor命令大气校正1.将需要大气校正的文件夹进行复制（路径不要太长，尽量不要有中文路径）2. 将路径复制下来例如：G:\Text\20210923\S2B_M
2021-09-25 19:50:31
4605
15
原创 Sentinel-2数据下载方法
一、Sentinel-2数据下载方法目前Sentinel-2数据下载包括以下5种方式，以从欧空局官网下载为例：（1）欧空局官网（欧空局官网）（2）USGS（3）Remote pixel（4）地理空间数据云（5）遥感集市1. 进入网站后先注册（最好使用科学上网，速度会快点，这里不翻墙也可以）2.选择传感器类型、参数、时间、区域等3.参数选择好之后，点击搜索，即可出现结果。注意，如果是L2A级别的影像，则是已经处理过的，不再需要进行大气校正。如果是L1C级别的影像，下载后则需要自己进
2021-09-25 19:25:46
10768
2
原创 01 Pytorch学习-创建自己的第一个神经网络模型
1.导入Pytorch框架import torchfrom torch.utils.data import Datasetfrom torchvision import datasetsfrom torchvision.transforms import ToTensorimport matplotlib.pyplot as plt;2.下载训练数据以及训练数据标training_data = datasets.FashionMNIST( root="data", train
2021-09-13 16:25:46
897
原创 如何将Excel中画好的图导入Word
在写文章时遇到过这样一个问题，在Excel中画好的图不知道怎样复制到Word中，直接复制粘贴过去的话，不是图片格式，导致后期更改不方便，博主通过摸索找到了一种很好的解决这个问题的方法：以如下一张图为例：如果直接复制粘贴到Word中，如下图所示，缩小会改变图像，使图像变形，不利于后期的更改，且不能保存为单独的图像按照如下方法则可以避免这个问题在Excel中选择图片后，点击 开始->复制为图片，选择“如打印效果”，点击确定。在Word中直接进行粘贴就可以了，接下来可以
2020-08-05 16:35:56
8727
原创 ENVI计算植被指数时大部分值变为0
给大家分享一个博主之前遇到的坑，通过已有多光谱数据计算NDVI植被指数，本来看着挺简单的，但有个小地方值得注意一下，不然一下子就入坑了（计算出来有很多无效值）。首先介绍一下NDVI植被指数计算公式NDVI = (nir-red)/(nir+red)主要是用ENVI的math band工具进行计算,主要用到近红外和红外两个波段。根据NDVI计算公式，在math band工具中输入如下公式进行计算b1选择nirb2选择red通过这样计算得出结果，发现结果很不理想，很多地
2020-06-08 17:50:05
5232
4
原创 投稿SCI，在最后一步校正时Word里面图片显示不出来解决办法
有些时候我们投稿SCI，需要最后一步在选校正了才能点击投稿，这时可能会发生下面这种情况：出现这样情况的原因是因为图片的格式是有问题的，需要对图片进行修改，具体问题在哪里，接下来给大家进行分析：以我自己的二值图像为例，格式为Tif，大部分图像都是对的，有个别几张图就是显示不出来，分别用Matlab进行读取后，发现两种图虽然显示没问题，但像素值是由去别的，如下图所示：不能正确显示的二值图像像素值能正确显示的二值图像像数值两张图片对比就很容易发现问题所在了，将二值图像的值全部换成255..
2020-05-25 17:15:28
2711
原创 MODIS数据批量下载方法，简单实用，速度在不使用科学上网的前提下也能达到1M左右
首先总结一下目前网上的几种下载方式以及方法链接，有兴趣的朋友可以了解一下。通过linux、shell、python批量下载https://my.oschina.net/u/2424163/blog/1799992.通过官网给出的wget命令批量下载https://ladsweb.modaps.eosdis.nasa.gov/tools-and-services/data-download...
2020-05-06 16:47:44
1562
1
原创 Matlab 3-D-Skel2Graph
Matlab Skel2Graph 3D()[A,node,link] = Skel2Graph(skel,THR);skel: 输入的3D 二值图像THR: 分支长度最小的阈值，用于过滤骨架化后的毛刺node有如下属性 ：idx :该节点立体像素值links:该节点的连接列表conn：该节点终点列表comX,comY,comZ：该节点所有立体像素的中心ep：如果该节点为终点节...
2020-04-25 19:25:40
535
1
原创 EARTH DATA网站注册走过的坑
1、需要下载遥感影像，在Earth data中下载时发现需要注册，注册就注册嘛，要求还挺多的。2、好不容易用户名和密码符合要求了，其他该填的也填了，最后人机验证也弄了，但就是一直提示出错，试了好多次都一样，反正信息都是对的，就是注册失败。3、停下来仔细思考一下，想着毕竟是国外的网站，需不需要连外网试试，果然，试一次成功了。（果然啥都得多试试，哈哈）4、所以总结起来就是，在Earth dat...
2019-10-18 12:52:45
8608
6
原创 Word里面如何在两个字正中间加入点
写论文时在作者栏有时需要在两个名字中间加入点，自己琢磨了一下，可以通过以下流程实现：1）在Word中点击“插入”——“符号”——“其他符号”2)在符号列表中下拉，找到所需要的点状符号（大致位置在后面，注意是正中间点）。3）结果...
2019-10-11 21:55:47
5585
原创 ArcGIS图像处理信息缺失
Arcgis中利用Raster to Polyline工具将栅格转为矢量时，发现有部分线段缺失。1、原始栅格图2、转换后的图用其它栅格图像转挺正常的，但用这个图像转就出错，实验了很久，发现原因是因为图像像素值的问题。其他二值栅格图像值包含0，1，而这张栅格图像值包含0，255，表面上看都属于二值图像，实则后者数据灰度图像，值范围0-255。在ArcGIS中的解决方法如下：在栅格图像上...
2019-09-30 11:41:56
2797
原创 Matlab二值图像保存避免信息缺失技巧
现象所谓二值图像，就是像素只包含0或者1的图像，由于论文需求，此处需要在MATLAB中将二值图像进行保存，二值图像如下所示：本来应该是个挺简单的操作，直接用imwrite（）函数保存就行了，于是我的操作如下：1.imwrite(Image_building,‘impervious.jpg’)但再次在MATLAB中导入图片的时候，问题就来了，二值图像值变得很奇怪：值在0-255之间，不再...
2019-09-05 14:54:09
1637
2
原创 Spectral（SPY）学习笔记(4)
Fisher线性判别数据读取显示：采用线性特征降维利用降维后的高光谱进行分类：目标检测（Target Detectors）估计背景统计信息。要将像素声明为异常，需要指定阈值RX分数。此处我们选择所有RX分数相对于背景的概率小于0.001的图像像素为异常像素:不需要为异常像素指定阈值，也可以简单地查看原始RX分数的图像，其中越亮的像素被认为“越异常”:对于样本图像，由于使用了线性...
2019-08-18 10:12:10
928
2
原创 Spectral（SPY）学习笔记(3)-高光谱降维（Dimensionality Reduction）
实现高光谱数据降维由于所谓的“维数诅咒”，处理数百个波段的高光谱图像可能会增加计算量，影响分类精度。为了缓解这些问题，通常需要降低数据的维数。主成分变换高光谱图像中的许多波段往往具有很强的相关性。主成分变换表示原始图像波段到一组新的、不相关的特征的线性变换。这些新特征对应于图像协方差矩阵的特征向量，其中相关特征值表示特征向量方向上的方差。在相对较少的主成分中(与原始频带数相比)可以捕获到非常...
2019-08-17 20:02:20
8168
1
原创 Spectral（SPY）学习笔记(2)-Supervised Classification
高光谱数据实现监督分类读取数据并显示-Unsupervised Classification
高光谱实现k-means聚类读取ENVI格式数据并显示## 开始K-means聚类参数说明：img：输入的高光谱影像5：第一个5表示k值5：第二个5表示迭代次数m:分类结果图像c:表示每个类别的波谱曲线显示分类结果...
2019-08-17 15:45:26
649
翻译 Python pyecart和matplotlib的区别
最近制图用到了这两个python库，发现制图功能都很强大，现总结下两者的区别以方便以后选用。1）pyechart更适用于做一些高大上的图，比如地图、地域分布、热点图、热力图等；而matplotlib更加适合于做一些常用的图，如折线图、柱状图、散点图等。2）pyechart主要可以做交互式的图，前者风格偏向数据工程demo展示；matplotlib主要可以做静态图，风格偏向学术。...
2019-07-19 15:43:19
2578
原创 ArcGIS对栅格数据进行掩膜提取后值变化以及不显示的情况
ArcGIS对栅格数据进行掩膜提取后，DN值变化以及图像不显示1、本人遇到的问题现有栅格数据以及矢量数据（红色框部分）本来想用矢量数据对栅格数据进行掩膜提取，得到研究区栅格数据，过程如下，在ArcGIS中AecToolbox>Spatical Analyst Tools>Extraction>Extract by Mask,过程很简单，运行也正常，但结果却出现错误，数据...
2019-07-03 21:11:01
25520
2
原创 ArcGIS中令栅格数据研究区中水体的DN值为零
1 数据准备1)研究区栅格数据2）研究区水体区域栅格数据合并起来是这样的：我们的目标是将图一研究区中的水体DN值变成0，而不是简单的剔除。2 数据处理利用ArcGIS中的Data Management Tools > Raster > Raster Dataset > Mosaic To New Raster工具Input Rasters中输入需要合并的栅格。O...
2019-06-30 10:58:55
1971
转载 常见的标准化方法
1 什么是数据标准化（Normalization）将数据按比例缩放，使之落入一个小的特定区间。在某些比较和评价的指标处理中经常会用到，去除数据的单位限制，将其转化为无量纲的纯数值，便于不同单位或量级的指标能够进行比较和加权。2 常用的标准化方法2.1：规范化方法也叫离差标准化，是对原始数据的线性变换，使结果落到[0,1]区间，转换函数如下:其中max为样本数据的最大值，min为样本数据...
2019-06-27 09:05:02
2302
原创 ArcGIS如何方便修改图例
ArcGIS修改图例小技巧加载图例后，右键单击选择Convert to Graphics，然后在右键单击，选择ungroup，相当于把图例的每一部分单独分离了，接下来就可以根据自己的额需求分别对每一部分进行编辑了。...
2019-06-25 09:20:43
21261
2
原创 U盘制作Ubuntu启动盘后容量突然变得只有几M
U盘制作Ubuntu启动盘后容量突然变得只有几M解决办法：进入计算机-管理-磁盘管理，找到U盘所在的磁盘，发现大部分空间都处于“未分配”的状态。右键单击选择重新新建卷，卷大小为未分配内存大小即可...
2019-06-24 09:38:05
2829
1
空空如也
TA创建的收藏夹 TA关注的收藏夹
TA关注的人