胡修修-CSDN博客

原创【MATLAB】柱状堆叠图给不同数据赋不同颜色bar

但一句bar只能赋一个颜色，所以用‘stacked’关键字无法调整颜色。效果如下：data7是浅绿色部分；data6是深绿色部分......柱状堆叠图如上，加个‘stacked’关键字即可。想要给不同数据指定颜色，必须分开写——多少个数据就写多少句。也就是人为堆叠，把要放在最上面的数据写在第一行，依次往下。

2024-06-04 09:59:24 363

原创 matlab画相关系数图（长时序--分时期分月份）

a图为例子，纵坐标是1-12月，横坐标是不同的因子。如左下第一格表示的就是x1与AO因子在十二月的相关性。另外，多子图的colorbar的设置需要在最后单独写几行代码，可根据自己需求更改。上面的代码大部分都是重复的，只要弄懂第一个子图怎么画，后面的改变量就行。（a/b/c是x1因子，d/e/f是x2因子）祝大家画出好看的图。

2024-04-19 10:54:38 566

原创 matlab中设置坐标轴范围和刻度 xlim/xtick/xticklable

设置要显示坐标刻度。%X轴的数据显示范围。

2024-03-12 08:47:46 1891

原创 python--gdal:tif图像坐标/投影坐标/经纬度转换（理清i和j的顺序）

上一段代码先将图像坐标（即tran的6个参数）→转化为了投影坐标（即px、py）再→转为经纬度（lon、lat）检查投影信息方法：print（dataset.GetProjection()）--若输出为空字符串则无投影。px是我们要计算的x轴方向上的坐标，而col代表列数（XSize），trans[1]是x方向的分辨率。gdal能读出里面的基本地理编码信息，根据这6个基本的信息，可以推导出其他需要的地理信息。换句话说，col是x方向上的间隔个数，而trans[1]是x方向上每个间隔的宽度。

2024-03-07 11:25:29 1579 1

原创 Sentinel-1相关信息（数据类型、成像模式、分辨率......）

哨兵-1（Sentinel-1）SAR影像，包含S1A和S1B。每颗星重访周期是12天，两颗星轨道相位差180°（也就是彼此位于地球对面的两个位置）。由于两颗星的所有设置一致，两颗星可以直接混用，即时间分辨率可达到6天。提供SAR极化方式：HH（水平-水平）、VV（垂直-垂直）、VH（垂直-水平）、HV（水平-垂直）。

2024-01-04 15:47:56 4398

原创 python二进制文件转nc（以PIOMAS海冰厚度数据为例）

首先要明白，二进制文件是4个字节按顺序依次存储所有数据，因此heff.H1979是没有记录对应的地理信息的。数据的地理信息会在另一个文件中，也就是下图中的grid.dat。这是第二个要下载的文件。（二进制文件存储的数据和.dat文件存储的地理信息是一一对应的，及二进制第一个数据所在的位置就是.dat中第一个经纬度）路径如图所示：girds放下载的grid.bat；sit放下载的二进制海冰厚度数据；output是输出结果文件的地方。以其中的海冰厚度数据为例进行转化。

2023-11-02 17:28:54 505 1

原创 matlab创建矩阵、理解三维矩阵

全0矩阵：a = zeros(2,3,4)全1矩阵：a = ones(2,3,4)matlab创建矩阵、理解三维矩阵

2023-10-27 15:41:19 1931

原创如何分离一个要素的shp矢量文件：利用ArcGIS分割工具

如何分离一个要素的shp矢量文件：利用ArcGIS分割工具。下面介绍如何用ArcGIS对含有多个分离区域的一整个面要素进行分割。找到裁剪工具，中文路径对应：分析工具>>提取>>裁剪。

2023-10-25 15:01:29 1489

原创 ftp传输工具--FileZilla安装教程及使用

FileZill安装教程、使用

2023-06-15 11:11:14 370

原创 python中关于模块/包/from...import...一些概念

（若sea_ice_motion文件夹内还有子文件夹作为包，该子文件夹也必须包含__init__.py文件。例如下图，在sea_ice_motion文件夹下，有个seaicemotion.py，如果要引用seaicemotion.py这个模块，首先要在sea_ice_motion。如果要使用某个模块，就要去对应的包（包可以理解为路径或文件夹）下引用对应的模块。（1）和当前.py文件同一目录下（即同一个包）的.py文件。（3）其他包的模块（即其他文件夹的.py文件），其中包含类、方法、变量和常量。

2023-05-30 20:07:24 501

原创 gdal读取tif常用函数

【代码】gdal读取tif常用函数。

2023-05-29 15:14:05 841

原创 python处理NC文件时，常用的一些函数

参数pathname，定义了文件路径匹配规则，可以是绝对路径，也可以是相对路径。返回所有匹配的文件路径列表。

2023-05-29 09:40:20 426

原创 ERA5风场速度提取（某区域某时间段），u/v合成风向计算，python绘图。

（维度就是帮助我们将数据分层分类，以便更好的进行索引，以前的文章中有详细解释，可以看看）气象上定义正北方向为0°（即风从北吹向南，是数学坐标系中的 -90°），顺时针转动角度增大。（我这里下载了02-20年的10m风场的u分量+v分量+风速大小）风向Dir=0°（或360°）, u=0, v0，正南风；风向Dir=270°, u>0, v=0，正西风。风向Dir=90°, u

2023-04-12 20:38:25 8919 3

原创三维数组的理解和加减乘除运算python

二维数组可以理解为一个坐标系，通过（x,y）来索引。而三维数组其实就是将多个坐标系叠在一起，通过（z,x,y）来索引。如图：A和B分别是两个二维坐标，将它们叠在一起就变成了三维坐标。而z就是来表示每个二维坐标的位置的序号，也叫“通道”。（例：要找到a12，即先找到A（在第一通道，z=0），再找第一行第二列（x=0,y=1））

2023-03-31 10:29:27 1179

原创 nc文件的索引方法（loc关键字索引、位置索引）

一、glob.glob()获取文件路径二、索引NC文件时间

2023-03-30 15:59:24 585

原创问题解决：input should be an array, list, tuple or scalar

在一个画图程序中出现了input should be an array, list, tuple or scalar的报错。2.创建一个数组，将lons和lats放进数组。1.在18行和19行后面加上.data。

2023-03-30 15:41:05 224

原创 ArcGIS对矢量数据进行缓冲区分析和叠加分析具体操作

1、缓冲区分析技术缓冲区是地理空间，目标的一种影响范围或服务范围在尺度上的表现。它是一种因变量，由所研究的要素的形态而发生改变。从数学的角度来看，缓冲区是给定空间对象或集合后获得的它们的领域，而邻域的大小由邻域的半径或缓冲区建立条件来决定，因此对于一个给定的对象A，它的缓冲区可以定义为：（d一般是指欧式距离，也可以是其它的距离， r为邻域半径或缓冲区建立的条件）缓冲区建立的形态多种多样，这是根据缓冲区建立的条件来确定的，常用的对于点状要素有圆形，也有三角形、矩形和环形等；

2023-02-28 16:37:01 13061

原创 python读取二进制文件（简单小白版）以PIOMAS的海温数据读取为例

python读取二进制文件（简单小白版）以PIOMAS的海温数据和经纬度数据读取为例

2023-02-23 10:46:22 1418 4

原创海冰厚度数据来源及产品数据介绍（卫星数据、验证数据、模式数据）

海冰厚度数据来源及产品数据介绍（卫星数据、验证数据、模式数据）

2023-02-22 15:01:26 1360

原创 excel制作图表，设置坐标格式、双坐标、图例等。

excel制作图表（折线图、柱状图），设置横坐标、双左边、图例等

2023-02-17 19:22:09 6136

原创 OpenStreetMap与Google Map 的区别及其优缺点

OpenStreetMap与Google Map 的区别及其优缺点

2023-02-16 15:11:58 1110

原创利用ArcGIS绘制空气温度等值线（间隔1°C）和等值面的具体步骤

利用ArcGIS绘制空气温度等值线（间隔1°C）和等值面的具体步骤。

2023-02-16 15:10:39 1764

原创南北极地区常用地图投影及其特点

南北极地区常用的投影是等角投影（又称正形投影、相似投影），由于该投影在无穷小范围内使地图上的图形同椭球面上的原形保持相似，因而在极地得到广泛的应用。

2023-02-16 15:04:10 2730

原创点云数据的语义分割算法综述总结大全（传统方法+基于深度学习的方法）

近些年基于深度学习的点云语义分割研究开始被广泛关注，并取得了极大的进展。我围绕点云数据的语义分割，从传统方法和基于深度学习方法两方面，来对近些年点云语义分割的研究做出系统的归纳和总结。因此恬不知耻的取名未——点云数据的语义分割算法综述总结大全（传统方法+基于深度学习的方法）哈哈哈哈。

2023-01-19 20:23:39 11732 2

原创 ArcGIS提取栅格格网点经纬度，导出为CSV/TXT（面转栅格，再栅格转点）

ArcGIS提取栅格格网经纬度，导出为CSV/TXT（面转栅格，再栅格转点）

2023-01-19 20:07:08 10328 1

原创【python】wget 批量下载文件（多年多月多日时间序列文件）

【python】wget 批量下载文件（多年多月多日时间序列文件）

2023-01-05 08:57:29 672

原创 SAR的三种不同极化方式分别是什么（如何理解单极化、多极化、全极化）？

SAR的三种不同极化方式分别是什么（如何理解单极化、多极化、全极化）？1.什么是SAR极化方式？2.什么是SAR中单极化、双极化、全极化？3.如何选择SAR雷达卫星极化方式？

2022-12-04 14:18:44 10981 2

原创【气象】python计算平流并cartopy绘图（计算差分、弧度制转化、广播运算）

【气象数据处理】python计算平流并cartopy绘图（计算差分、弧度制转化、广播运算）

2022-11-28 14:18:07 2401

原创 python简单读取和索引.nc文件(气象小白入门版)

python简单读取和索引.nc文件(气象小白入门）这里的顺序指的是：一个是属性存放的顺序（时间是第一个，高度是第二个...），还有每个属性中的数据存放的顺序（时间属性中，1948排在1949前面）。print(file.air)运行结果如下图所示，可以看到属性存放的顺序是：时间-->高度-->纬度-->经度，索引的时候一定要按照这个顺序写约束。，有四个，每个维度的数字就代表有几个值，比如level:17就是有17个高度值，想要具体看看是什么值，可以print一下。，若下载的文件是温度数据，变量就是温度

2022-11-26 10:58:42 14630 6

Mluoo的博客