python-海洋数据处理
文章平均质量分 58
简朴-ocean
一年级在读
海洋科学专业
Python学习记录
物理海洋方向
大气波动动力学
WRF数值模拟
公众号:简谱学记
展开
-
cdo | 常用命令
整理一下平时经常会使用的cdo命令。原创 2024-05-31 23:03:30 · 728 阅读 · 0 评论 -
CMIP6降水单位转换
CMIP6 降水单位转换原创 2024-04-30 11:39:12 · 492 阅读 · 1 评论 -
cdo 修改 calendar 为标准的格式
cdo 修改 calendar 为标准格式原创 2024-04-24 10:45:11 · 275 阅读 · 0 评论 -
cdo + shell 批量插值nc文件
近期,开始处理CMIP6模式数据,由于不同的模式数据的分辨率之间各不相同,为了方面后续分析,这里对于数据进行插值处理。原创 2024-04-23 16:54:41 · 430 阅读 · 0 评论 -
Python | 月平均气候态 | SST
Python | 月平均气候态 | SST原创 2024-04-12 16:29:14 · 512 阅读 · 2 评论 -
GPM降水数据下载-Linux系统上使用wget
如何下载GPM降水数据原创 2023-12-16 18:33:10 · 407 阅读 · 0 评论 -
python -- 替换netcdf文件中的时间
以上,通过shell脚本与python脚本的结合,初步实现了缺失时刻的nc数据的填充。对于海表面温度(SST)、海面风速等海洋要素特征来说,短时间内变化不那么剧烈,可以采用上述方法;当然,这里展示的只是通过复制前一时刻的数据,也可以通过插值的方法,这就需要额外编写相关的python脚本,然后加到对应的。原创 2023-09-21 10:01:01 · 260 阅读 · 0 评论 -
python - 快速搜索指定格式文件的方法
python - 快速搜索指定格式文件的方法,使用pathlib库进行实现原创 2023-09-19 17:06:35 · 457 阅读 · 0 评论 -
python -- 实现路径的匹配,剔除掉指定路径,并保存路径
【代码】python -- 实现路径的匹配,剔除掉指定路径,并保存路径。原创 2023-09-01 11:23:57 · 287 阅读 · 0 评论 -
关于多grib文件合并以及netcdf和grib的方法
介绍一些好用的工具处理netcdf和grib原创 2023-07-04 21:52:40 · 1371 阅读 · 0 评论 -
python&深度学习 -- 训练的数据集中出现NaN值应该怎么办?Filling gaps in a numpy array
但是,如果训练集中存在 NaN 值,则计算梯度时出现除以 0 的情况,从而导致梯度无法计算。因此,在深度学习中,通常会在训练集中删除所有包含 NaN 值的样本,或者使用插值或其他技术来填充 NaN 值,以确保数据集不包含 NaN 值。函数来计算每个 NaN 值到离它最近的非 NaN 值的距离,并将其替换为该非 NaN 值。模型无法收敛:如果训练集中存在 NaN 值,则可能会干扰损失函数的计算过程,从而影响模型参数的更新,使模型无法收敛。基于上述问题,这里介绍一个处理Nan值的简单方法,实现NaN值的填充。原创 2023-05-11 20:35:26 · 1319 阅读 · 1 评论 -
python--读取TRMM-3B43月平均降水绘制气候态空间分布图(陆地区域做掩膜)
python--读取TRMM-3B43月平均降水绘制气候态空间分布图(陆地区域做掩膜)原创 2023-05-05 13:30:47 · 1814 阅读 · 0 评论 -
python -- 检验GPM数据是否缺失
python -- 检验GPM数据是否缺失原创 2023-04-28 10:33:13 · 852 阅读 · 0 评论 -
windows & Linux :python 脚本 下载日本葵花8号卫星数据L1级产品,自定义分辨率、时间
近期需要用到日本葵花8号卫星数据,用于相关研究,而通过官方提供的下载方法,难以针对性的下载所需要的数据类型,因此这里编写了一个针对葵花8号卫星L1级数据产品的脚本下载,主要实现两个功能:1、自定义时间下载2、选择任意的分辨率下载nn=08代表葵花8号,nn=09表示葵花9号年月日小时分钟Rbb = R21,表示的是整个区域,波段从 “01” 到 “16”;Rbb=R14,表示的是日本区域,波段为14。原创 2023-04-13 11:04:17 · 1566 阅读 · 0 评论 -
python -- 批量读取多个文件,并将每个文件中相同变量累加
python -- 批量读取多个文件,并将每个文件中相同变量累加原创 2023-03-13 19:35:17 · 2826 阅读 · 3 评论 -
python--海温、OLR数据分布做显著性检验,绘制空间分布并打点
【图+代码】python--海温、OLR数据分布做显著性检验,绘制空间分布并打点。原创 2022-10-22 23:02:42 · 4823 阅读 · 5 评论 -
python--使用convolve 对二维数据进行平滑
python--使用convolve 对二维数据进行平滑。原创 2022-10-13 17:38:52 · 4273 阅读 · 4 评论 -
python--测试使用不同的方式计算位涡平流项的差异
python-计算位涡平流项原创 2022-10-13 17:22:46 · 611 阅读 · 0 评论 -
linux下使用.csh脚本自动化下载GFS气象预报数据
linux下使用.csh脚本自动化下载GFS气象数据。原创 2022-10-13 17:08:12 · 2768 阅读 · 0 评论 -
python--对站点数据做EOF并做插值绘制填色图
使用python对站点数据进行EOF处理,并使用两种插值方法到规则网格上。原创 2022-09-23 21:53:19 · 2915 阅读 · 4 评论 -
python读取excel数据、转换度分为度,绘制填色图
【代码】python读取excel数据、转换度分为度,绘制填色图。原创 2022-09-19 15:38:58 · 698 阅读 · 0 评论 -
python--计算纬度/经度格式的网格点之间的实际距离
一般在处理海洋、气象数据时,往往会通过一些收支方程来对某一气象要素进行诊断分析,而大部分的诊断方程往往都会包含对于变量的偏导。原创 2022-09-13 23:08:55 · 1211 阅读 · 1 评论 -
python--循环绘制ERA5风场的空间分布图
通常,在处理气象海洋资料时,经常会绘制风场的空间分布图进行简单分析,而常常需要连续绘制多天,并将多张子图绘制到同一个图片中,因此这就需要用到循环绘图。原创 2022-09-02 23:14:57 · 4430 阅读 · 1 评论 -
python -- 计算有效波高、平均波高
现有波面数据如下图所示,采用频率为20Hz,需要根据波面数据计算得到平均波高、有效波高等信息。原创 2022-09-02 23:02:31 · 2359 阅读 · 0 评论 -
python--转换wrf输出的风场数据为网页可视化的json格式
一般网页可视化风场中的数据都是json格式,而如果我们希望将wrf模式模拟输出的风场数据在网页中进行展示,这就需要先将wrfoutput数据转换为网页可以识别的json格式。原创 2022-09-01 11:05:08 · 1410 阅读 · 0 评论 -
Neither Quantity object nor its magnitude supports indexing
python 使用metpy给变量赋单位时,出现报错Error: Neither Quantity object nor its magnitude supports indexing原创 2022-08-16 18:52:50 · 432 阅读 · 0 评论 -
python--实现Lorenz 63模式
python 实现洛伦兹模式原创 2022-08-09 23:18:32 · 1561 阅读 · 0 评论 -
python--由wrfouput的数据计算位势涡度,并插值到指定压力层
通过python,计算wrf模式输出数据中的位涡原创 2022-07-26 21:14:02 · 2174 阅读 · 1 评论 -
python:将二维array of float32类型的数据转为float类型
通过python将二维array of float32类型的数据转为float类型原创 2022-07-15 13:07:05 · 2177 阅读 · 0 评论 -
ERA5-hourly-total precipitation的单位换算问题:m转为mm
Conversion table for accumulated variables (total precipitation/fluxes)原创 2022-07-08 14:31:26 · 3922 阅读 · 7 评论 -
python绘图--由逐日风场数据计算月平均风场数据并绘制二维填色图
利用python处理风场数据,将逐日每天四次的日资料处理为月平均数据,并绘制空间填色图。原创 2022-06-08 11:26:11 · 2379 阅读 · 1 评论 -
python 将垂直速度从压力坐标系转为高度坐标系(pa/s转为m/s)
前言最近在计算一下涡度的垂直项时,发现一个问题,如果使用笛卡尔坐标系下的涡度计算公式,会发现对于我们从ERA5下载的垂向速度ω,它的单位是:pa/s,而对于水平的速度u、v来说,它们的单位是:m/s,这在计算速度的旋度时就会出现问题,由于单位不一致引起计算无法进行,对于后续数据处理造成干扰。本来我在网上搜寻解决方案时,发现大部分的解答都只提供了解决思路(由公式转换),往往比较冗杂,不够简单。所辛,我在metpy这个工具包中发现了包装好的函数,只需要提供相应的数据,即可成功转换!这里展示由压力坐标向高度原创 2022-05-13 21:10:27 · 4239 阅读 · 8 评论 -
python 基于metpy计算位涡平流项(水平)
python 计算位涡的平流项原创 2022-05-13 20:49:25 · 3049 阅读 · 5 评论 -
使用salem处理wrfout数据,进行切片、并插值到等压面
使用salem处理wrfout数据,进行切片、并插值到等压面最近,有学习到通过salem处理wrfout数据,非常的简单快捷,读取的变量也比较方面。也可以快速出图,下面简单对比一下xarray和salem读取wrfout文件的区别:import xarray as xrfrom salem.utils import get_demo_fileds = xr.open_dataset(get_demo_file('wrfout_d01.nc'))可以发现,通过xarray 读取wrfout文件原创 2022-04-26 22:07:09 · 2551 阅读 · 0 评论 -
python 对海洋、气象数据进行滤波--带通滤波处理
python对于海洋、气象数据进行带通滤波处理原创 2022-04-12 22:08:31 · 5381 阅读 · 17 评论 -
python 对陆地数据进行掩膜的两种方法
1、使用python自带的库进行掩膜2、使用cdo+geopandas+Salem进行掩膜原创 2022-03-18 20:28:36 · 9985 阅读 · 20 评论 -
python 读取netcdf4文件的全过程(基础教程,看了包会)
python 读取netcdf4(nc)文件的完整教学过程,再基础不过拿到一个nc文件,如何使用python进行读取呢?本文带你完整走一遍流程。前期准备:xarray 、matplotlib、cartopy库的安装:conda install xarray matplotlib cartopypython安装一个nc文件数据python 安装的教程就不详细介绍了,不会的可以参考这个链接:Windows系统下python 及相关库、环境配置的安装教程下面开始正式数据读取,第一步,打开基于你原创 2022-02-15 20:27:37 · 15363 阅读 · 15 评论 -
python 基于cartopy库绘制台风路径(包含代码详细解释)
使用python基于cartopy库绘制台风路径使用python 绘制西太平洋进入我国的台风路径,文件为.dat格式,内容如下所示:代码如下:import netCDF4 as nc4import matplotlib.pyplot as pltimport numpy as npimport datetimeimport osimport cartopy.crs as ccrspath='E://' #文件路径files= os.listdir(path) #得到文件夹下的所有原创 2022-02-11 21:54:33 · 4601 阅读 · 5 评论 -
python--根据任意非网格经纬度坐标,找到均匀网格点上最接近的经纬度坐标
需求:根据非规则经纬度坐标,查找均匀网格点上最接近的经纬度坐标,并提取该点上的变量。思路:1、首先选取一个较大的经纬度范围,将你的非规则经纬度坐标包含进去;2、计算均匀网格点上经纬度与非规则经纬度坐标的绝对值距离;3、选择绝对值距离最短的点的索引坐标;4、根据索引坐标提取最邻近的经纬度坐标5、选取最近邻经纬度坐标对应的变量数据:海表面高度异常数据:sla分辨率:0.25X0.25数据类型:2017年月平均数据(netcdf)代码实现:主要用到几个库:xarray、numpy、ma原创 2022-02-07 21:03:57 · 4624 阅读 · 1 评论 -
python 批量读取地转流速日资料绘制气候态年平均海表流场(填色为流速)
利用python对海洋流速数据进行处理,绘制海表流场原创 2022-01-27 18:09:14 · 3482 阅读 · 3 评论