- 博客(52)
- 收藏
- 关注
RSS订阅原创 根据shp统计栅格数据并进行shp数据选取
这里可以根据需求来进行统计，比如平均值，总和等。这一步主要将数据链接到目标shp中。三.按照目标条件来选择shp数据。一.根据shp数据统计栅格数据。二.将统计的dbf链接到shp。四.将选择的shp数据导出。
2024-07-10 13:28:16
204
原创 基于Matlab和Python泰勒图的绘制
泰勒图：泰勒图1常用于评价模型的精度，常用的精度指标有相关系数，标准差以及均方根误差(RMSE)。一般而言，泰勒图中的散点代表模型，辐射线代表相关系数，横纵轴代表标准差，而虚线代表均方根误差。这里介绍两种复现的代码，第一种是原始的泰勒图的代码(Python)，包括原始的误差项目，第二种是泰勒图种只显示均方根误差和R方改进的泰勒图(Matlab)。这个必须是高版本才能运行，比如我使用的是2022版本的，并且需要把泰勒图的包安装到matlab的bin文件中，并且要将代码文件也和泰勒图的代码包放在一起.
2024-07-10 13:13:19
337
原创 基于GEE平台提取研究区多年参数的均值（生成CSV文件）
GEE（全称Google Earth Engine）是大谷歌下属的一个可以批量处理卫星影像数据的工具，属于Google Earth一系列的工具。相比于ENVI等传统的处理影像工具，GEE可以快速、批量处理数量“巨大”的影像。通过GEE可以快速计算比如NDVI、LST等，可以预测作物相关产量，监测旱情长势变化，监测全球森林变化等。对于某些研究，需要长时间序列的数据，若下载遥感影像，再进行处理就会很繁琐，很复杂。直接通过GEE平台进行遥感影像的处理，可以节省大量的时间和精力。本文主要是基于GEE平台提取多年参数
2023-05-04 20:07:09
1439
原创 基于ArcGIS对遥感地物分类shp结果进行编辑和修改
我们采用遥感的手段进行土地分类，其结果可能都不是特别准确，那么就要在分类的结果上进行地物修改以达到分类结果更加精准的目的。那么对于有分类结果的数据，也就是有ENVI版本的.dat文件，那么可以在ENVI中进行手动修改，若没有ENVI版本的.dat文件，只有分类好的shp文件，那么只能在arcgis中对分类结果进行修改。1、若有ENVI版本的.dat文件的分类结果2、若没有ENVI版本的.dat文件的分类结果（只有shp的分类结果）：那么在arcgis中对shp进行编辑。（1）按照属性值将其显示在地图上
2022-09-18 20:52:51
3167
2
原创 借助excel工具进行多元线性回归模型的建立及案例分析
大家都知道两个变量进行相关性分析，最简单的是建立两个变量的散点图，通过计算R的平方来判定两个变量的相关性。但是实际的研究中，一个因变量是受到很多因素影响的。比如植被覆盖度不仅受到温度的影响，还有降水量和积温等影响，所以在分析相关性时简单建立两个变量的散点图是无法很好地分析他们的影响机制的。那么此时需要建立多个变量的多元线性回归方程。下面以植被覆盖度的案例来进行分析。案例：年尺度的综合因素对植被覆盖度的影响分析注： 因变量为植被覆盖度自变量为温度、降水量和积温工具：excel（1） 将数据进行归一
2022-05-22 12:17:26
9209
1
原创 PS初学者的一些实用技巧
最近接触PS比较多一些，我就分享几个我觉得有意思的PS技巧：1.调整像素：PS对写论文的人最重要的是调整像素，那么Ctrl+Alt+i，可以直接对像素进行调整。如图：2.画线条：对于PS如何采用不同的颜色进行线条的绘制。选择PS中的钢笔。红框里可以进行颜色、线条的形状、像素大小的更改。但需要注意的是，由于填充是白色，如果你在某幅影像上画线，其填充的颜色会遮挡你的内容，那么教大家一个办法解决。如图，将填充变为0%。3.画箭头：如图所示，找到自定义工具，其箭头的形状也可以自由选择。4.如何复制
2021-12-31 22:26:08
917
原创 计算栅格数据像元的趋势值及检验值P
本文章主要展示如何计算出栅格数据上每个像素的变化率及检验值P。与上篇的栅格数据之间相关性的计算是类似的代码，以下代码供大家参考学习：本代码是在matlab中运行的其输入和输出的数据为tif格式的栅格数据[a,R]=geotiffread('F:\SMC_V3.0\std_anlisyt\annual\x\FP_FVC_2013.8_class.tif_clp.tif');%先导入某个图像的投影信息，为后续图像输出做准确info=geotiffinfo('F:\SMC_V3.0\std_anlisyt
2021-12-23 21:32:16
3287
3
原创 ArcGIS10.2安装的详细步骤、注意事项及许可变永久操作
arcgis10.2版本可以联合其他的编程软件进行使用，因此说很方便，我强烈推荐各位小伙伴安装arcgis10.2版本。在安装之前要将防火墙关闭。下面就说一下安装arcgis10.2的详细步骤和注意事项：1.运行ESRI2.安装arcgis License manager一直next,其中注意修改路径到其他盘，尽量不要放在C盘。安装完毕后，启动License Manager，选择停止服务。如图：3.安装ArcGIS for Desktop安装期间一路下一步即可，什么也不用管。安装完成后启
2021-12-22 15:56:28
10097
原创 不同栅格数据之间的相关系数计算（输出为tif影像）
两类栅格数据之间的相关性分析终于搞出来了程序。程序的适用如下：比如一个植被指数NDVI与降水量的相关性分析、比如蒸散量与温度的相关性分析等等…只要大家觉得两个参数之间有什么可以分析的，都可以拿以下代码去干，哈哈哈。我这里是分析8年的植被覆盖度与年尺度的地表温度的例子。1.相关系数具体的计算公式为下：相关系数的范围是[-1,1]。2.代码如下（matlab）：[a,R]=geotiffread('F:\SMC_V3.0\std_anlisyt\annual\x\FP_FVC_2013.8_class
2021-11-02 21:24:43
23960
68
原创 栅格数据做相关性分析前的预处理（批量定义投影、栅格投影、重采样）
最近在做两个栅格数据的相关性分析，这个最主要的是两个栅格数的行列号是一致的（也就是说坐标系、分辨率一致），那么标准的栅格数据的分辨率是以米的形式存在的，但是有一些栅格数据是以度为单位存在的，所以在做分析前，一定要先把其转换为米，否则在重采样的时候就会出错。（相关文章——栅格数据以度为单位变为以米为单位的简单做法）当然，在变其分辨率的单位前，一定要将其坐标系转换为一致。所以就有以下的步骤（如何将两个栅格数据的坐标系、分辨率的单位和分辨率一致）1.首先，要看两个栅格数据的坐标系是否一致，如果不一致，要先统一到
2021-10-26 18:31:03
7117
9
原创 批量栅格投影（arcpy）
栅格投影的意义：将分辨率单位为度的变为米定义投影的意义：将某个栅格的坐标系变为想要的坐标系批量栅格投影的代码如下：# -*- coding: UTF-8 -*-#需要自行更改的如下：#inws：输入路径#outws：输出路径#Coordinate_System：目标投影坐标系import globimport osimport arcpy# 输入路径 应该注意，中文路径，会导致读不出文件inws = r"------"# 参考文件路径 使用栅格数据集（从其导入方形像元大
2021-10-26 18:25:18
3574
7
原创 批量定义投影(arcpy)
今天突然发现一些影像拖到ArcGIS中和shp文件不重合，实际上是投影不一样的问题，这些问题很常见。那么如何将大量的影像进行定义投影呢？那就用arcpy!哈哈哈！代码如下：# -*- coding: UTF-8 -*-import arcpyimport globimport os'''批量定义投影： 对输入文件夹下的数据，根据坐标系参考文件，进行批量定义投影 需要修改： inws：输入路径（必选） mask：参考范围的shp文件（必
2021-10-20 19:49:08
2909
1
原创 高分六号WFV预处理详细步骤（辐射校正、大气校正、几何校正）
高分六号WFV传感器获得的影像的预处理分为两个步骤：辐射校正和大气校正，辐射校正为了消除传感器的自身误差，大气校正为了消除大气的折射误差。一、辐射校正按照红框中信息进行更改。辐射校正前地物光谱反射率：校正后的地物光谱曲线的变化：二、大气校正ENVI中大气校正一般采用Flaash模型，但是这个模型需要设置许多的参数十分的麻烦，所以采用ENVI商店（app store,这个也需要下载的）里的Flaash Easy To Use 这个插件，可以快速地进行大气校正。大气校正后的地物反射率这
2021-09-18 18:04:20
9267
11
原创 在ArcGIS中提取所需地方区域的shp
最近在摸鱼的过程中，发现了如何从中国的一些省、市、县的行政边界提取你想要的区域的shp文件。因为大家也知道shp文件用来对影像的裁剪等等，作用极大，因此将这个小技巧分享给大家。提取步骤：1.将中国的省（市、县）的shp拖进ArcGIS中去。2.打开shp文件的属性表，选找到自己想要的区域，在最左边双击，如下图打对勾的地方3.那么你想要的区域就会在shp中显示。4.右击shp文件，将其数据导出。如下图所示5.将其保存一个地址即可。6.看一下结果图。这样就成功的将想要的shp文件进行
2021-09-05 16:32:15
34643
4
原创 免费的气象数据（GLDAS、ERA5）下载步骤及数据处理的详细过程
需要气象数据的小伙伴们看过来呦！首先声明，能够免费下载的气象数据一般分辨率比较低哟！一、GLDAS数据1.下载地址：（GLADS下载地址）GLADS数据分有3h、一天、一个月的产品，大家可以根据自己的需求进行下载。2.下载步骤：选择完数据产品、研究区范围、日期，将其生成下载的链接，网址会直接将链接放在文本里，用户可以进行下载。（1）首先，比较容易地下载步骤是将文本导入火狐浏览器中下载管理器中，具体可以参考——参考文章，此方法的缺点是一般不能链接到网址，下载总是出错。（2）采用wegt进行下载，下载
2021-09-02 22:40:53
20282
13
原创 阈值法批量生成二值化分类影像
这几天在弄一个阈值分类的事情。比如当NDVI大于0的时候是属于植被的，小于0为其他，以前做的时候就采用ENVI的决策树或者掩膜处理，之后就会生成一个二值化的影像。现在，想弄个程序，批量的得出分类之后的二值化影像。当然这个适用于一系列的影像的阈值分类。废话不说了，看程序。我是在IDL中实现的。pro yuzhi_class ;输出阈值分类过的二值化影像 compile_opt strictarr COMPILE_OPT idl2 envi, /restore_base_save_files
2021-08-02 15:12:04
1358
1
原创 土地转移矩阵的计算步骤
土地转移矩阵在遥感研究中非常的实用，可以计算几年内的土地的面积转移变化，废话不多说，看步骤！1.要找两年的土地利用数据，比如我这里是2013年和2019年MCD12Q1的500m的土地利用数据，然后要知道每个数值代表的是什么类，然后整理成像下面的表格一样。因为欧文这个研究区，只有这四类（15，16是一类），注意Excel保存为CSV格式！2.将每个土地利用数据建立文本字段“2013”和“2019”.将excel导入arcgis中，然后将每个土地利用数据都与这个excel相关联，2.将excel与各
2021-07-26 22:51:33
4579
原创 监督分类的具体步骤（依据高分辨率计算混淆矩阵）
好久没更新了，暑假的我貌似有些不务正业，晃晃悠悠的每天。最近河南暴雨、南京疫情，我的不务正业也随之停止，不过带球挺有意思。希望接下来的一切都能顺利！最近也在关注奥运会，昨天看了乒乓球赛，运动精神永远的神，又回想起来四年前代表学校第一次参加省运会，每赢一个球，队友的呐喊声，到现在都无法忘怀。作为女队队长，四年来大大小小的比赛，虽然没有怎么为母校带来很多荣誉，但是球队的凝聚力和队友的帮助，这辈子铭记于心。 坐在办公室窗户前，外面乌云密布，说下雨来着，也一直没有下，办公室还是熙熙攘攘的那几个人，很多次怀疑自己的能
2021-07-25 17:13:51
3801
原创 IDL多类影像批量计算注意的问题
这是学长在毕业之前为不上心的我解决的一个问题，在此祝学长们毕业快乐呀！真心感谢此为大神学长，希望学长以后一切可以顺利！1.问题： 我在使用IDL程序的时候，也没有报错，就是数值不太正确。一直不知道什么原因，后来发现，是因为两类的影像的像素大小不一样。2.解决方案：在以后处理影像的时候，要先拼接—投影—重采样—裁剪，按照这个顺序，就不会出错。顺序千万不能错，错了结果就不对。希望自己能不断进步吧 哈哈哈哈啊哈哈长路漫漫…唯有坚持…...
2021-06-30 20:43:19
300
原创 高程数据的下载以及运用高程数据进行坡度和坡向的计算
一、高程数据的下载1.下载参考链接——各种DEM数下载方法DEM数据质量对比——参考文章2.我直接在地理空间数据云里下载的，这个相对方便一些，下载可以根据自己所需的研究区域进行下载。二、DEM数据的预处理我采用的是NASA的DEM-V2数据（分辨率为30m），可以进行拼接和裁剪.这个数据下载下来就是TIF格式，可以直接用来做任何的处理。批量拼接——批量裁剪——...
2021-06-06 17:35:23
2728
原创 ArcGIS中cell size 是0.0000几变为以米为单位的解决方法
最近想将GLASS数据（分辨率为500m）的数据采样到30m，发现一直不能成功。那么故事就开始了…1.GLASS数据在arcgis中像元大小是用经纬度表示的，如下图所示。2.一般在arcgis中，这个像素的大小是以米为单位表示的。如下图所示的Landsat影像是30m的。那么为啥会出现这种情况？其实MODIS的数据在arcgis中都会出现这种情况。因为MODIS数据的投影坐标和landsat的不一样。要想将其经纬度显示变成标准的像元大小显示，就要重投影。如下图是MODIS数据的投影信息。下图是L
2021-06-04 17:01:22
25900
7
原创 土地利用分类数据类型和下载
土地利用分类数据目前有：MCD12Q1(500m分辨率)，GLC30m(目前有2010 2015 2017 2020-数据下载），FROM-GLC10m。主要介绍FROM-GLC10m。一、介绍在2017年GLC30m数据基础上，又提高了分类精度，其从30m提升到10m的土地利用数据——称为FROM-GLC10 - (2017 V0.1.3 with minor improvements)二、数据的下载下载链接——数据下载2.区域选择规则根据自己的经纬度范围选择需要的数据块它的数据块，经度间距
2021-06-02 21:34:18
5732
1
原创 ArcGIS中数据框的页面从竖向变为横向
现在自己太容易健忘了，这个小小的操作步骤总是记不住，所以记录下来，供大家参考大家在制图过程中，由于一些影像比较的宽，而ArcGIS中默认的页面是纵向的，所以不能满足我们的需求，因此我们要将页面调成横向的。1.刚开始的页面是如图这样的2.然后点击左上角的file中的Page and print setup3.点击Landscape，页面即可变成横向。变成横向的页面如下所示。这个操作虽然简单，但比较重要，希望对大家有所帮助~长路漫漫…唯有坚持…...
2021-05-27 08:57:05
6768
原创 ArcGIS中shp文件的裁剪
一般，我们都是用arcgis来进行栅格数据的裁剪，也就是卫星影像数据。但是在arcgis作图的时候，也需要对shp文件进行裁剪。步骤如下：首先创建一个shp文件，在Catalog中，连接一个文件，在此文件中创建shp。shp的格式是polygon（多边形）然后编辑，可以画出来多边形。然后再开始裁剪。注意：一般都是用研究区域shp去裁剪多边形。这样就可以了！小技巧，供大家分享！长路漫漫…唯有坚持…...
2021-05-12 16:28:44
28217
1
原创 总体分类精度和kappa系数计算实例详细介绍！
为啥要算总体分类精度和kappa系数呢？想必大家都知道是为了精度评价，当我们没有实测数据的时候，那么总体分类精度和kappa系数就派上用场了！我们没有实测数据，依旧能够评价自己的方法和模型的优良性。博客写的有点啰嗦~啊哈哈哈哈哈哈哈哈哈一、分类精度和kappa系数的计算公式首先先看一下总体分类精度和kappa系数的计算公式1.1总体分类精度：1.2kappa系数：我们先知道总体分类精度和kappa系数的公式，下面看具体怎么操作二、生成随机点在真彩色（或者假彩色）影像在随机生成随机点。
2021-05-10 16:22:34
19350
3
原创 哨兵2号数据的下载（各种方法）
哨兵二号是一个高分辨产品，对于分析地物很重要，下面我介绍2种下载哨兵2号的方法1.采用python代码进行下载在用此代码下载之前，要进入欧空局哥白尼数据中心——欧空局哥白尼数据中心网址先注册，注册挺简单。然后，将下载的区域，用框框画出来然后由于国内用户仅能下载Level-1C数据，而且一般都是2A级产品。所以选择如下：选择时间日期搜索出来的影像加入购物车，然后点击购物车里面的一个下载按钮，会出现一个文件，将其保存接下的代码，需要这个product.meta4文件接下来是下载数据的
2021-04-25 20:20:22
18643
5
原创 将一个文件夹里的所有影像按照天放进一个文件夹里
1.代码应用场景：比如现在我的文件夹里面有8个影像，分别是两天的数据。我想把其一天的数据放在一个文件夹里，方便下一步的处理。2.源代码如下：# -*-coding:utf-8-*-import osimport shutilimport re'''脚本说明： 根据文件夹下的文件名，来新建文件夹， 并把所有该天的文件移动到新建的文件夹'''path = r'E:\yz\MODIS\2019-1.taihu'# 列出文档file_list = os.listdi
2021-04-14 09:22:05
319
2
原创 利用python将文件夹按需进行分类
此次代码的作用是，由于MODIS的数据是8天尺度的，那么需要对应的CLDAS的数据也要是8天尺度的，那么就应该将CLADS的数据按照8天的形式进行分组，比如一年的数据有365个，那么8天一个文件夹，需要分成46组。这只是一个案例，大家可以根据自己的情况再分成不同的组数。看代码吧举个列子，比如要将两年的数据，每一年的数据按照2张影像为一组进行分类。我这里一年的数据有四张影像。注意：import shutil 不要采用python2.7，我用的是python3.7。from os import l
2021-04-13 09:48:21
788
3
原创 利用python批量从img影像中提取所需的波段
利用python批量从img影像中提取所需的波段现在我知道landsat数据下载，得到的波段是一张张单个的影像，但是MODIS影像下载出来的波段是合在一起的，如果处理的数据量很大，那么手动将其波段信息提取出来很麻烦，所以接下来的代码，可以将img影像中所需要的波段提取出来，废话不多说，看代码。1.注意事项：注：这个代码要用高版本的python，本人这个是3.7的，而且如果没有gdal numpy osgeo库需要自行配置，如果没有这几个库，代码无法运行成功。from osgeo import gda
2021-04-12 15:21:55
2813
原创 下载MODIS数据的巧妙办法——构造下载链接
大家在下载MODIS数据的时候，大多时候是一次性下载一年的，甚至是好多年，但是总有一些数据在下载过程中会错误，仅仅有10几个，20几张影像没有下载，再一一在NASA网址寻找十分麻烦，所以构造下载链接就可以解决这个问题，废话不多说，看步骤1.首先，要知道自己需要那几张影像，知道日期，那么我们在文本里输入自己想要的日期，如下图先写成这样的形式，方便我们接下来构造。然后最后写出这种的链接，直接在网址输入，便可以直接下载。比如如下为我构造出来，下载MYD09的链接。https://ladsweb.moda
2021-04-10 16:27:09
668
原创 Mathtype的安装、使用以及问题解决
首先，Mathtype是什么呢？是word中编辑公式的一个插件，大家可能会说，word中已经有编辑公式的，何必再加上一个这样的插件呢？word中带的编辑公式的，都是编辑一些比较基础的公式，对于科研狗肯定不行了，复杂公式一大堆，用Mathtype插件可以对编辑公式减少很多麻烦！废话不多说了，接下来就告诉大家如何安装和使用。1.首先是Mathtype的安装包，安装包下载的话...
2021-03-30 09:19:05
12114
31
空空如也
空空如也
TA创建的收藏夹 TA关注的收藏夹
TA关注的人