低调的大耳朵图图
学遥感的小朋友
展开
-
基于Matlab和Python泰勒图的绘制
泰勒图:泰勒图1常用于评价模型的精度,常用的精度指标有相关系数,标准差以及均方根误差(RMSE)。一般而言,泰勒图中的散点代表模型,辐射线代表相关系数,横纵轴代表标准差,而虚线代表均方根误差。这里介绍两种复现的代码,第一种是原始的泰勒图的代码(Python),包括原始的误差项目,第二种是泰勒图种只显示均方根误差和R方改进的泰勒图(Matlab)。这个必须是高版本才能运行,比如我使用的是2022版本的,并且需要把泰勒图的包安装到matlab的bin文件中,并且要将代码文件也和泰勒图的代码包放在一起.原创 2024-07-10 13:13:19 · 177 阅读 · 0 评论 -
基于Python代码画柱状图(附全部代码)
代码功能主要实现的是横坐标是四类,图例为3类的柱状图。原创 2023-11-13 21:31:09 · 422 阅读 · 0 评论 -
栅格数据做相关性分析前的预处理(批量定义投影、栅格投影、重采样)
最近在做两个栅格数据的相关性分析,这个最主要的是两个栅格数的行列号是一致的(也就是说坐标系、分辨率一致),那么标准的栅格数据的分辨率是以米的形式存在的,但是有一些栅格数据是以度为单位存在的,所以在做分析前,一定要先把其转换为米,否则在重采样的时候就会出错。(相关文章——栅格数据以度为单位变为以米为单位的简单做法)当然,在变其分辨率的单位前,一定要将其坐标系转换为一致。所以就有以下的步骤(如何将两个栅格数据的坐标系、分辨率的单位和分辨率一致)1.首先,要看两个栅格数据的坐标系是否一致,如果不一致,要先统一到原创 2021-10-26 18:31:03 · 6858 阅读 · 9 评论 -
批量栅格投影(arcpy)
栅格投影的意义:将分辨率单位为度的变为米定义投影的意义:将某个栅格的坐标系变为想要的坐标系批量栅格投影的代码如下:# -*- coding: UTF-8 -*-#需要自行更改的如下:#inws:输入路径#outws:输出路径#Coordinate_System:目标投影坐标系import globimport osimport arcpy# 输入路径 应该注意,中文路径,会导致读不出文件inws = r"------"# 参考文件路径 使用栅格数据集(从其导入方形像元大原创 2021-10-26 18:25:18 · 3523 阅读 · 7 评论 -
批量定义投影(arcpy)
今天突然发现一些影像拖到ArcGIS中和shp文件不重合,实际上是投影不一样的问题,这些问题很常见。那么如何将大量的影像进行定义投影呢?那就用arcpy!哈哈哈!代码如下:# -*- coding: UTF-8 -*-import arcpyimport globimport os'''批量定义投影: 对输入文件夹下的数据,根据坐标系参考文件,进行批量定义投影 需要修改: inws:输入路径(必选) mask:参考范围的shp文件(必原创 2021-10-20 19:49:08 · 2848 阅读 · 1 评论 -
下载MODIS数据的巧妙办法——构造下载链接
大家在下载MODIS数据的时候,大多时候是一次性下载一年的,甚至是好多年,但是总有一些数据在下载过程中会错误,仅仅有10几个,20几张影像没有下载,再一一在NASA网址寻找十分麻烦,所以构造下载链接就可以解决这个问题,废话不多说,看步骤1.首先,要知道自己需要那几张影像,知道日期,那么我们在文本里输入自己想要的日期,如下图先写成这样的形式,方便我们接下来构造。然后最后写出这种的链接,直接在网址输入,便可以直接下载。比如如下为我构造出来,下载MYD09的链接。https://ladsweb.moda原创 2021-04-10 16:27:09 · 648 阅读 · 0 评论 -
GLASS数据的批量下载(以植被覆盖度FVC为例)
最近在搞植被覆盖度,所以又和学长学到一个FVC(植被覆盖度)数据下载的方法。虽然,学长已经记录过了,但是…我记性不好,有些细节,还是自己来记录吧…哈哈哈哈啊哈哈哈…在此感谢学长…1.首先,GLASS数据的下载网址——美国马里兰大学2.选择第一个产品3.打开就是年份的文件4.由于此数据是MODIS数据为基础的,想要下载自己的区域,就要选择轨道号什么的,那么就麻烦了。那么就有了大神的代码了…此为…学长的代码博客——批量下载GLASS代码博客5.那么比如我选择的区域以及时间,跑出来的代码是6.将这原创 2021-03-26 20:44:28 · 6481 阅读 · 7 评论 -
Python中opencv pandas 等库的配置(超级简单的方法)!
首先值得注意的是,对于python中的环境配置必须要有pip。pip是一个帮助python进行下载安装包的一个重要工具,所以python想要配置各种环境,必须有pip.1.下载pip教程——pip的下载链接2.有了pip就好办多了,不管配置什么环境都进入这个网站——python的各种环境的安装包3.在网站上下载好的包,比如下载pandas包点击下载,然后将其放在python下的scripts文件下。4.打开电脑中的开始——找到cmd命令提示符,一定要管理员身份运行,否则可能会出错。如果想找哪个盘原创 2021-03-05 21:26:18 · 471 阅读 · 0 评论 -
利用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 · 769 阅读 · 3 评论 -
将一个文件夹里的所有影像按照天放进一个文件夹里
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 · 296 阅读 · 2 评论