- 博客(291)
- 资源 (36)
- 收藏
- 关注
RSS订阅原创 python 语言始终是中文
出现中文输出的原因可能与你的 系统区域设置 或 Python环境配置 有关。即使使用 calendar 模块，某些环境下 month_name 仍会返回本地化的月份名称（如中文）。最近在python绘图时发现，在使用pandas进行时间显示的格式时，明明设置的转换为英文，但是始终呈现的是中文的形式。方法二：通过 locale 设置强制使用英语环境。方法三：使用 pandas 的字符串接口（推荐）方法一：强制指定英文月份（不依赖系统设置）
2025-04-01 19:18:23
235
原创 Python | 计算散度
简单记录一下计算散度的方法，方便以后查找复制包括三种计算方式，。其中，numpy和metpy的方法进行了比较，结果比较一致。windspharm的方法里面包含了两种方法， 一致是直接调用函数实现，另一个是先计算梯度再相加。使用windspharm最简单，缺点就是需要使用全球的数据作为输入，而且在Linux上安装，我这里用的是模式数据就没去和metpy和numpy进行验证，只是记录作为一种方法。比较意外的是，同样的数据，numpy竟然比metpy还快一点。
2025-03-25 09:38:04
137
原创 笔记 | SVD | 奇异值分解
其原理是要求两个场的展开系数的协方差最大。比如对某海域的海表温度(SST)和海表气压(SLP)做SVD分析，是为了提取出这两个变量之间相互耦合(协方差矩阵)的最重要的空间模态和时间演变信息。其中的空间模态就是SVD得到的奇异向量(左右奇异向量分别对应SST和SLP)，它们有各自的时间系数表征了时间演变特征。EOF实际上是SVD的特殊情况，对于SVD针对的是两个物理时空场，而EOF是一个物理时空场，如m表示空间点，n表示时间点，那么。那么，同样的道理，根据原始数据和奇异向量也可以用来得到对应的时间序列，即。
2025-03-25 09:36:53
691
原创 Python | kelvin波的超前滞后合成图
记录一下kelvin波的超前滞后合成图主要通过kelvin波数据以及原始降水异常数据，选择同一个纬度点，比如说赤道上；然后选择一个超前滞后的时间，这里选择为10day最后得到的数组是一个leadtimexlon的数组，下面绘图就简单了。
2025-03-19 11:07:55
312
原创 Python | eps矢量图不支持透明度设置
今天绘图时遇到个问题，发现使用python绘图中使用了透明度参数alpha时，如果你要保存的图片格式为矢量图eps。初步试了一下，有点效果，但是对于多种透明度颜色时还是有较大限制的。而且，会导致保存出的图片占用较大空间。其他方法就简单了，先转换为其他格式，比如说pdf，再使用工具转换为eps格式。转换的方法的效果较好，能压缩图片的占用空间，画质还没有太大损害。或者，直接保存为png，再从ai里面另存为eps格式。相关问题也早有人提出，并给出了可能的解决方案.如果是png或者pdf格式是没有问题的。
2025-03-19 11:05:54
240
原创 当ssh出现：REMOTE HOST IDENTIFICATION HAS CHANGED!
这种情况一般是你重装了ssh或者认为修改了你的密钥，导致无法在其他机器上通过ssh远程连接此机器。
2024-12-25 15:56:58
611
原创 Linux 编译 convert_geotiff 时遇到的几个问题
这里对编译好的数据进行测试，新建一个文件夹，选择一个任意的tiff文件，调用相关命令。使脚本生效，再次编译，显示成功，这里有个警告但是不影响后续使用。这里我这里出现报错了，所以分开来debug，首先。所以，这里找到生成的Makefile文件，通常，会指定编译器和链接器的选项，你可以在其中找到。这个文件，这里通过locate命令进行查找，，在这里，原本的我的LIBS后面为。测试程序，是可以正常有提示信息的。然后，就发现编译好了，存在了。所以这里需要在环境配置脚本。注意，这里需要管理员权限。
2024-12-04 17:42:23
583
原创 Python | 同时绘制带有投影与非投影的子图
经常需要在一张图片上，同时绘制两个子图，一个是带有投影的空间分布图，另一个为不带投影的曲线图。比如说EOF的空间模态以及PC序列。在python中，如果在一张画布上同时绘制以上两个图，通常会使得两张图片的长宽比非常的不协调。虽然可以通过手动设置子图的大小来调整，但是如果画布的大小改变，则又需要进行细调，感觉非常笨重麻烦。今天在读文献的时候，又看到这种配置的图片，索性花时间来实现让其聪明的解决。
2024-12-03 10:57:11
712
1
原创 Python | 自定义colormap | 设置对称colormap
设置第一个颜色为白色设置中间的颜色为白色以上函数适用于设置的colorbar的颜色数量为偶数的设置情况，常用于以0为中心的正负对称colormap以及希望设置第一个颜色为白色的colormap的情况。
2024-12-03 10:53:55
376
原创 Python | 爬取微信公众号推文 | 仅供测试
之前一直想偷懒写个脚本，方便后续整理推文，昨天测试了半天，后续因为太频繁结果被封了ip。今天重新修改了代码，终于算是可以了。以下是抓取的结果：非常简单，只有两部分内容， 标题和link简单调试了一个脚本方便自己的后续推文整理，当然，爬取的内容可能不仅仅是原创推文，以及一些图片的文字。可以再进行相应调整。该教程为个人学习记录，本人公开的代码以及爬取的文本为本人个人公众号，相关代码仅供测试与学习！参考方法：https://github.com/wnma3mz/wechat_articles_spider?
2024-12-03 10:53:25
1786
2
原创 Python | 涡旋识别
听别人开题报告，介绍一个涡旋识别的方法，是基于开源的一个python的库。感觉挺有意思的，打算来试试效果以下windows的python版本为3.9，Linux的python版本为3.8。
2024-10-31 11:31:42
813
2
原创 Python | seaborn | heatmap
使用python基于seaborn绘制好看的热力图heatmap。。最无脑的操作就是将数据处理为的格式，这一点在xarray读取的nc中很方便。可以借助这个函数以下使用的测试数据是noaa的sst-anomaly的月平均数据，简单将其处理为二维数组，这里简单将其处理为-20°S-20°N，120°E-180°E区域平均内2000-2019年共20年的数据，同时为了满足为二维矩阵，将其按年份和月份分组，并计算每个月的平均值。然后调用进行绘图,以下是绘图结果。
2024-10-25 09:34:51
619
原创 Python | basemap空间绘图 | cartopy | geoviews
又处理了两周的cmip6数据，真的是daily的资料属实恶心人┗|｀O′|┛ 嗷~~。白天下数据，晚上运行插值脚本。第二天再下数据，再运行插值脚本。两周时间我觉得我的shell编程又进步明显。批量处理数据这一块不能说相当拿手，但是十分熟练了。峰回路转，抽空整理并复习了一下netcdf官网关于basemap绘制空间分布图的几种类型，为后续绘图做准备，发现basemap画图有他好用的道理，确实相比cartopy要方便不少。同时发现了两个有意思的库，geoviews以及holoviews。
2024-10-16 14:20:09
748
原创 Python | 由高程计算坡度和坡向
之前参加一个比赛，提供了中国的高程数据，可以基于该数据进一步计算坡度和坡向进行相关分析。对于坡度和坡向，这里分享一个找到的库，可以方便快捷的计算。RichDEM，官网地址如下相关代码和测试数据放到了GitHub中，感兴趣的可以试试。
2024-10-07 22:32:44
740
原创 Linux 系统上运行python 绘图显示找不到字体该怎么办
找到虚拟环境对应的matplotlib绘图库下面的字体所在的文件夹，手动将下载的缺失字体丢进去运行下面的代码即可。
2024-09-21 10:59:47
280
原创 python | x-y 网格切片
通常， 我们处理的毕竟完善的nc产品，一般呈现未的维度，且lon和lat都是规则的网格，我们可以方便的使用选择合适的区域进行切片。但是，部分nc产品比如卫星轨道或者模式输出的数据，可能并没有处理为规则的网格产品，其经纬度信息以x和y表示了分辨率对应的网格格点数量。如下面这种情况：那么，在这种分布下，如何进行区域的切片呢。在python中实现对于x-y网格的区域截取，相关数据和代码放到了GitHub上，感兴趣的朋友可以试试。
2024-09-20 10:58:38
337
原创 Python | 读取.dat 文件
使用matlab可以输出为.dat或者.mat形式的文件，之前介绍过读取.mat后缀文件，今天正好把.dat的读取也记录一下。复习了一下使用pandas读取.dat文件的相关函数，以及pandas的一些基础命令，绘图多y轴的方法。
2024-09-20 10:57:56
4028
原创 Python | 绘制核密度图
台风天，适合敲代码。前两天正好看到一个有意思的核密度图，使用的是seaborn绘制的。之前了解过这个包，但是一致没有去绘制相关的图，这次正好去学习一下相关的函数。绘制结果如下所示：主要两个有意思的地方在一张图上同时绘制两次填色图，一个有意思的实现在核密度填色图的两侧分别绘制数量占比的柱形图。
2024-09-20 10:56:13
828
原创 Python | GPCP | 趋势分析 | 气候态
最近，开始对于CMIP6的一些数据进行评估。Talor图是一个很好展现模式间误差的方式，这里简单记录一下在python中的实现方式。
2024-09-09 16:49:28
1155
原创 Python | 泰勒图
最近，开始对于CMIP6的一些数据进行评估。Talor图是一个很好展现模式间误差的方式，这里简单记录一下在python中的实现方式。主要为半图的画法效果大致下面这个样子这边在原本代码基础上，主要是设置了三个坐标label的颜色以及内部网格线的颜色的区别。
2024-09-03 10:37:18
1660
1
原创 识别热带赤道波动的艺术
天气尺度(时间) to 行星尺度(空间)，纬向传播的低纬扰动波动的周期从几天到几周线性波动理论的波解还包括热带扰动，东风波，季节内振荡KW：频率波数滤波器 (FWF) 方法比空间投影方法显示出更精细的结构和更慢的传播MRG：方法之间存在很大差异和复杂的垂直结构对于Kelvin 波来说， 作为结构最简单的波动，使用频率-波数滤波可以得到更精细的结构，有利于分析波动内部动力学以及其与背景气流的相互作用。
2024-08-29 12:00:16
889
原创 Python | Windows | 安装一个常用的海洋气象环境
经常会看到一些交流群里讨论关于python虚拟环境中某些库安装不成功的问题，比如说关于地学数据处理和绘图的salem、geopandas、关于wrf输出的wrf-python，以及一些仿照matlab绘图风格的basemap。这里，将我常用的一个安装方法分享给大家，便于在安装的虚拟环境出现问题时，可以比较便捷、快速的安装常用的库。包括常用到的以上提供了一个关于海洋-气象专业中常用的python库的便捷安装方法，比如salem、geopandas、basemap，可以非常方便的处理不同虚拟环境的问题。
2024-08-26 21:08:25
1108
1
原创 Python | 处理海洋2C 数据 | 非标准时间格式
最近，师弟在用Python读取某海洋2C数据时，突然冒出一个报错：“时间单位的参考日期无效，当前日期 00:00:00.0这让我回想起，似乎在很久很久以前，我处理SMAP和Argo数据时也遇到过类似的问题。为了不让未来的师弟师妹们再次被这种问题困扰，我决定在这里记录一下这个坑。希望未来再看到这条记录时，能少叹气、多微笑，毕竟已经知道怎么解决了嘛！
2024-08-26 21:07:28
1137
原创 Python | Linux | 解析Himawari-8/9 | Standard Data
之前一个相关的工作需要解析Himawari-8/9 Standard Data文件，因为他是二进制的，之前没有处理过，导致完全摸不着头脑。但是需要注意的是，一个时刻的文件包含的信息可能比较少，需要同时读取多个时刻的文件合成的nc文件才能获得有效的结果。说的比较含糊，总之是谁用谁知道。第二和第三个替换为你的netcdf中的include和lib的安装路径，意思就是告诉系统你的netcdf库安装的位置。下面找一个具体的数据进行测试，这部分的代码我通过python来实现，在python中调用刚刚编译的程序。
2024-08-26 21:06:37
1134
原创 对流耦合的热带波动 | MJO
这不太可能，因为虽然这种大尺度运动在理论上是可能的，但 J. Bjerknes（1938 年）和 E. Höiland （1939 年）已经证明，条件不稳定性有利于尽可能小尺度的积云对流。热带大气的垂直热力结构显然更有利于小尺度积云对流，而不利于热带气旋尺度的对流环流，为什么气旋会在这种条件不稳定的热带大气中形成？有人提出，气旋是通过一种次级不稳定性发展起来的，在这种次级不稳定性中，现有的积云对流在低层水平辐合区增强，而在低层辐散区减弱。第四种情况很可能是冷核东风波的原因，但几乎可以肯定不是飓风的原因。
2024-08-11 12:43:43
425
原创 如何解决 windows11系统 使用中电脑突然自动休眠的问题
注意：干净启动是通过一组最少的驱动程序和启动程序来启动 Windows，以便你可以确定后台程序是否干扰你的游戏或程序，并帮助您找出问题的原因。（请一定要先选择隐藏所有microsoft服务，再点全部禁用，不然可能会导致出现无法进入系统等不可预见的问题）如果以上两步过程中，故障消失，请尝试逐个移除您电脑上已经安装的杀毒软件／系统优化软件／系统监控软件直至问题消失。4.进入服务选项卡-点按左下角的隐藏所有microsoft服务，然后点按全部禁用。2.打开系统配置，选择常规选项卡-有选择的启动。
2024-08-11 12:42:37
2011
原创 Linux | 安装lb-toolkits 1.2.4库
最近又需要下载葵花的数据，之前分享过一次代码。今天发现之前的环境不小心被我删了，而运行相关的代码需要安装lb-toolkits这个库，今天正好记录了一下安装lb-toolkits的过程。这里安装的版本是1.2.4，别问为什么是这个版本，因为只安到了这个版本才成功。本来是直接通过命令直接安装的，但是发现一直报错。
2024-07-09 21:13:23
578
原创 如何将overleaf的latex格式转换为word
其中output.docx是你输出的word格式，main.tex 是你的overleaf格式。运行代码之后就可以看见生成的word了。但是，文件的参考链接是没有附带在文件夹里面的。
2024-07-09 11:38:32
7963
2
原创 Centos7 安装老版本的chrome
我这里是centOS 7。然后在安装最新版的google-chrome时，总是会报错显示存在依赖环境的问题，使得无法安装成功chrome。本来以为是yum的源的问题以及缺失相关依赖库的问题，尝试了换源发现还是存在这个依赖的环境；而且貌似是centOS7 根本不支持这个GLIBC_2.25。属于是搞了一晚上啥也没搞定，血压倒是上升了挺大。犯了很多国外的帖子，这个问题不少人也遇到过，但是基本上没有给出解决的办法。要么就是瞎写一通根本不行。
2024-07-04 10:45:05
3221
9
原创 Python | 计算位涡平流项
最近忙着复习、考试…都没怎么空敲代码，还得再准备一周考试。。。等考完试再慢慢更新了，今天先来浅更一个简单但是使用的在做动力机制分析时，我们常常需要借助收支方程来诊断不同过程的贡献，其中最常见的一项就包括水平平流项，如下所示，其中var表示某一个变量，V表示水平风场。−V⋅∇var−V⋅∇var以位涡的水平平流项为例，展开为−u∂pv∂xv∂pv∂y−u∂x∂pvv∂y∂pv位涡受背景气流的调控作用。
2024-06-30 19:59:41
1057
3
原创 如何在python中读取matlab输出的.mat文件
经常会使用matlab处理数据，比如之前在进行空间滤波发现matlab自带的空间滤波函数。但是，由于后续使用python绘图更熟练一点，所以涉及到在python中读取matlab输出的文件，一般以.mat文件为主。下面测试了几种读取方式，选了一个比较好用的方式。这里简单将读取的.mat1、array数组(如网格数据)2、str字符串列表(如网页链接)总体上，比较了对于纯数组array以及纯字符str输出的mat文件在python中的读取方式。对于array输出的mat文件，推荐使用。
2024-06-30 19:57:57
1019
natural-earth.zip
2021-09-22
HDFView-3.1.1.msi
2021-04-23
matlab-读取nc-GUI (2).zip
2020-06-06
气象数据处理方法、海洋数据分析方法
2020-06-06
空空如也
TA创建的收藏夹 TA关注的收藏夹
TA关注的人