qq_41828818的博客

泰勒图：泰勒图1常用于评价模型的精度，常用的精度指标有相关系数，标准差以及均方根误差(RMSE)。一般而言，泰勒图中的散点代表模型，辐射线代表相关系数，横纵轴代表标准差，而虚线代表均方根误差。这里介绍两种复现的代码，第一种是原始的泰勒图的代码(Python)，包括原始的误差项目，第二种是泰勒图种只显示均方根误差和R方改进的泰勒图(Matlab)。这个必须是高版本才能运行，比如我使用的是2022版本的，并且需要把泰勒图的包安装到matlab的bin文件中，并且要将代码文件也和泰勒图的代码包放在一起.

代码功能主要实现的是横坐标是四类，图例为3类的柱状图。

最近在做两个栅格数据的相关性分析，这个最主要的是两个栅格数的行列号是一致的（也就是说坐标系、分辨率一致），那么标准的栅格数据的分辨率是以米的形式存在的，但是有一些栅格数据是以度为单位存在的，所以在做分析前，一定要先把其转换为米，否则在重采样的时候就会出错。（相关文章——栅格数据以度为单位变为以米为单位的简单做法）当然，在变其分辨率的单位前，一定要将其坐标系转换为一致。所以就有以下的步骤（如何将两个栅格数据的坐标系、分辨率的单位和分辨率一致）1.首先，要看两个栅格数据的坐标系是否一致，如果不一致，要先统一到

栅格投影的意义：将分辨率单位为度的变为米定义投影的意义：将某个栅格的坐标系变为想要的坐标系批量栅格投影的代码如下：# -*- coding: UTF-8 -*-#需要自行更改的如下：#inws：输入路径#outws：输出路径#Coordinate_System：目标投影坐标系import globimport osimport arcpy# 输入路径 应该注意，中文路径，会导致读不出文件inws = r"------"# 参考文件路径 使用栅格数据集（从其导入方形像元大

今天突然发现一些影像拖到ArcGIS中和shp文件不重合，实际上是投影不一样的问题，这些问题很常见。那么如何将大量的影像进行定义投影呢？那就用arcpy!哈哈哈！代码如下：# -*- coding: UTF-8 -*-import arcpyimport globimport os'''批量定义投影： 对输入文件夹下的数据，根据坐标系参考文件，进行批量定义投影 需要修改： inws：输入路径（必选） mask：参考范围的shp文件（必

大家在下载MODIS数据的时候，大多时候是一次性下载一年的，甚至是好多年，但是总有一些数据在下载过程中会错误，仅仅有10几个，20几张影像没有下载，再一一在NASA网址寻找十分麻烦，所以构造下载链接就可以解决这个问题，废话不多说，看步骤1.首先，要知道自己需要那几张影像，知道日期，那么我们在文本里输入自己想要的日期，如下图先写成这样的形式，方便我们接下来构造。然后最后写出这种的链接，直接在网址输入，便可以直接下载。比如如下为我构造出来，下载MYD09的链接。https://ladsweb.moda

最近在搞植被覆盖度，所以又和学长学到一个FVC（植被覆盖度）数据下载的方法。虽然，学长已经记录过了，但是…我记性不好，有些细节，还是自己来记录吧…哈哈哈哈啊哈哈哈…在此感谢学长…1.首先，GLASS数据的下载网址——美国马里兰大学2.选择第一个产品3.打开就是年份的文件4.由于此数据是MODIS数据为基础的，想要下载自己的区域，就要选择轨道号什么的，那么就麻烦了。那么就有了大神的代码了…此为…学长的代码博客——批量下载GLASS代码博客5.那么比如我选择的区域以及时间，跑出来的代码是6.将这

首先值得注意的是，对于python中的环境配置必须要有pip。pip是一个帮助python进行下载安装包的一个重要工具，所以python想要配置各种环境，必须有pip.1.下载pip教程——pip的下载链接2.有了pip就好办多了，不管配置什么环境都进入这个网站——python的各种环境的安装包3.在网站上下载好的包，比如下载pandas包点击下载，然后将其放在python下的scripts文件下。4.打开电脑中的开始——找到cmd命令提示符，一定要管理员身份运行，否则可能会出错。如果想找哪个盘

此次代码的作用是，由于MODIS的数据是8天尺度的，那么需要对应的CLDAS的数据也要是8天尺度的，那么就应该将CLADS的数据按照8天的形式进行分组，比如一年的数据有365个，那么8天一个文件夹，需要分成46组。这只是一个案例，大家可以根据自己的情况再分成不同的组数。看代码吧举个列子，比如要将两年的数据，每一年的数据按照2张影像为一组进行分类。我这里一年的数据有四张影像。注意：import shutil 不要采用python2.7，我用的是python3.7。from os import l

1.代码应用场景：比如现在我的文件夹里面有8个影像，分别是两天的数据。我想把其一天的数据放在一个文件夹里，方便下一步的处理。2.源代码如下：# -*-coding:utf-8-*-import osimport shutilimport re'''脚本说明： 根据文件夹下的文件名，来新建文件夹， 并把所有该天的文件移动到新建的文件夹'''path = r'E:\yz\MODIS\2019-1.taihu'# 列出文档file_list = os.listdi