笔记
文章平均质量分 71
什么都不会的张同学
生活繁忙,随缘回复
展开
-
NCL绘制WRF domain区域并添加气象站点
原创 2024-05-13 18:17:51 · 471 阅读 · 0 评论 -
Fortran读取netcdf文件/WRF中的文件读取
一直很好奇WRF到底如何通过netcdf库读取netcdf文件,正巧有个机会,试了下fortran读取nc文件,总结一下。原创 2023-12-10 21:10:25 · 1140 阅读 · 1 评论 -
NCL读取hdf5文件——以AMSR2海冰密集度数据为例
由于大部分计算都是在服务器进行。因此考虑增进NCL和Shell脚本水平,之前的博文介绍了HDF5文件及其读写,但是由于服务器本身系统版本较老,且属于内网原因,决定直接使用自带的NCL脚本对其进行撰写记录。原创 2023-04-03 21:19:39 · 1130 阅读 · 0 评论 -
个人免费静态博客网站搭建步骤记录——Git+hexo
静态网站是指全部由HTML(标准通用标记语言的子集)代码格式页面组成的网站,所有的内容包含在网页文件中。网页上也可以出现各种视觉动态效果,如GIF动画、FLASH动画、滚动字幕等,而网站主要是静态化的页面和代码组成,一般文件名均以htm、html、shtml等为后缀。简单来说,与动态网页相比,静态网站的网页全由固定的html写成,内容稳定固定,没有接口与时间限制,如果网站创建者不进行修改,静态网页的内容与网址一般不会发生改变。原创 2022-11-09 11:46:01 · 1219 阅读 · 0 评论 -
Python气象绘图笔记——常用气象绘图函数脚本封装与使用记录
由于工作需要,将对我常用的python绘图脚本进行封装,为了防止代码丢失、忘记使用流程等,写个博客记录下。要加载的包。原创 2022-11-02 16:00:20 · 5013 阅读 · 0 评论 -
R语言空间插值/R语言离散数据网格化/R语言空间点插值/R语言nc日均转月均、日期转换
在我们日常使用数据时,我们获取的数据并不是连续的网格,而是网格状的离散数据。这时,我们就需要利用拥有的离散网格数据,自定义将其插值。通常插值我们需要shp文件给予边界,不过,只要知道了所需区域的经纬度范围,我们可以自行构建网格将点数据插值在自己定义的经纬度范围中。原创 2022-10-24 11:13:54 · 3386 阅读 · 0 评论 -
python气象数据处理实例:构造自定义时间索引(时间戳)/r站点日均转月均/dataframe添加时间索引(datetimeindex)
pandas作为python里常用的时间序列分析库,在气象数据处理中有着极为广泛地应用,其resample功能可以轻松地将不同频率时间序列转换,这正好满足气象中需要求日均、月均、甚至季节平均的需要,但是,resample需要其序列(series)或数据框(dataframe)具有时间索引,可在一些气象站点数据中,时间索引需要我们自己构建添加,而这些时间点并不均匀分布——换句话说,在某些时间点,数据缺失,这时需要我们自己构建时间日期,并添加至dataframe成为时间索引,再使用resample采样。原创 2022-10-09 16:05:59 · 2241 阅读 · 1 评论 -
python绘制北极风场与位势高度场/python cartopy极地投影重叠问题解决
python作为胶水语言,近年来在气象数据处理与可视化中有着极为广泛地应用,诸多第三方库极为便利地满足了气象研究者处理数据与绘图的要求,处理数据的包:xarray、pandas,读取不同气象格式的包:netcdf4、h5py,用于模式后处理的wrf-python、python-CDO等等。同样,python在气象可视化方面也有着许多包,一般而言,python的气象绘图是用matplotlib包绘图,用Basemao或cartopy包绘制地图底图,结合绘制。原创 2022-10-06 11:22:01 · 4476 阅读 · 6 评论 -
python绘制极地投影/cartopy20.0+解决python极地投影问题/python nc可视化——以北极10m风场为例
在之前的博客。中,我曾经提到过python cartopy在极地投影中难以添加坐标标签的问题,当时解决方法是使用文本的方式添加,但这种方式的麻烦且并不统一适用。而在之后的可视化过程,我发现python在极地投影时仍会出现一些奇怪的问题:比如,等值线扭曲,出现不规则多边形,风场分布不均匀等等。然而,这些问题都在最新版的cartopy21中解决,因此被cartopy极地投影折磨了数月而不得不转用m_map的我得到了解脱,下面将以极地2018极地春季的10m风场为例,给大家看看效果。原创 2022-10-04 20:13:08 · 4291 阅读 · 9 评论 -
WRF后处理:使用python paramiko库将服务器中wrfout文件批量下载至本地(sftp.get)/python通过sftp远程批量传输文件/sftp OSError: [Errno 22
我们在服务器上跑完WRF时,不免会有许多的wrfout文件,这时候就需要进行后处理。当然我们可以直接在服务器上进行后处理,但有时由于一些原因,比如笔者的服务器上无法连接conda channel(ˉ▽ˉ;)…,便将其下载至本地再进行后处理。通过xftp软件一个个下载未免有些困难,正巧python也有支持ssh与sftp的库,我们就用python批量下载吧。原创 2022-09-27 17:49:44 · 1471 阅读 · 0 评论 -
maltab datenum函数与正则表达式巧用:逐日数据转为逐月数据、日序转月序
在处理数据时,我们经常需要将逐日数据转为逐月数据,这就需要我们知道日序对应的月序,比方说,我有2019年一整年的天数据共365个,我知道这个文件代表的是2018年的第几天?这方法当然可行,但却是基于C语言的一种操作,实际上,利用matlab的datenum函数,我们只用一行代码就能实现这个过程。然后先批量读个数据文件(hdf格式),为了避免过于复杂,我们先算一年的吧,这时,i=1,year(i)=2005。当然,这种方法有一个要求,即,逐日数据需要完整,文件命名按照日序排列。原创 2022-09-02 19:21:00 · 2972 阅读 · 0 评论 -
R语言批量处理nc数据/R语言MK趋势检验/NCL绘图(nc可视化)
回所了,开始处理数据,由于cartopy绘制极地投影加标签实在是太麻烦了,就用R把nc数据处理了,再用ncl画图。原创 2022-08-26 11:23:58 · 4959 阅读 · 1 评论 -
Anaconda3 spyder无法打开问题解决实录
总的来讲这次的问题到解决方法,让我感到最大的问题在于andaconda的安装包太老了,清华源的anaconda虽然快,但安装起来问题极多,如果不嫌慢还是去官网下吧,下好了换了channel就行了。原创 2022-08-18 19:34:13 · 2725 阅读 · 0 评论 -
气象数据数据处理实例——matlab字符串切割匹配与R语言日期匹配(数据拼接)
这是两年半前大三时所写的代码和任务,如今看来其实很幼稚,当时得出的结论是1、在学习像R、matlab这种语言,一定要查明实现这种功能是否有类似函数解决2、R语言的处理数据类型时,涉及大量数据转换细节,只有多用才能较为熟练使用虽然当时处理时叫苦连天,不过从现在的目光看,这种类型的气象数据处理已经是很基本的操作了,之后我的数据处理变成了nc数据/卫星数据/栅格数据,对于气象站点数据的处理已经基本不做了,之后的方向也应该注重于卫星数据和模式上。...原创 2022-08-03 15:41:33 · 1396 阅读 · 0 评论 -
WRF后处理/Python处理nc数据与可视化/极坐标网格绘制(Cartopy、netcdf4)——以北极雪水当量数据为例
python处理NC数据,并可视化原创 2022-06-27 18:44:31 · 6448 阅读 · 3 评论 -
GITHUB无法打开与下载失败解决方法总结
学SNICAR模式,去GITHUB下源码,学校网速感人加上IP限制,折腾近一个小时才完成,网上文章太多,自己总结一下。下载单个文件通常下载提供都是整个项目下载,但对于过大项目,下载其中几个就够了,网上的方法花里胡哨,其实有很简单的方法。点进我们要单个下载的项目,比如这里我想先下nicarAD_v3.m 和 snicarAD_v3_drv.m 源文件研究一下。我们点进相应文件,看到右上角的raw,右键,点击将链接另存为即可。GITHUB无法打开与下载失败由于DMS解析的IP无法访问造成,我们需原创 2022-03-05 14:49:34 · 15848 阅读 · 1 评论 -
气候统计实习——EMD与ESMD方法介绍与应用
引言原创 2021-10-30 17:43:13 · 5028 阅读 · 20 评论 -
气候统计分析笔记之气象数据的预白化(WS2001)
“气候变化”这几年在公众视野频频中出现,恐怕连小学生都知道“温室气体会造成全球变暖,破坏环境,影响人类生存。”无论气象学家们承不承认,“气候变化”已经成为了当今气象学科不可回避的一个话题。与关注伴随而来的,是无数的质疑:“气候变化是否真的受到人类活动的影响?影响有多强?人类是否应该做出举措?”最尖锐,也最需要回答的问题是:”气候变化是否只是场骗局?“实际上,目前国际上对于气候变化的问题依然争论不休,虽然已有多项证据表明了气候变化的真实存在性与紧迫性。可相关科研人员也很难对这个问题给出十分确切的回答。原创 2021-10-24 21:21:11 · 3561 阅读 · 3 评论 -
基于Arcpy的栅格数据批量处理+简单说明
最近写毕业论文,把气象数据和卫星影像插值的插值,拼接的拼接,临到要计算了,问题成堆的出现,把人累的半死。本文的目的是给出一些常用的批量处理栅格数据的arcoy代码,做出简单说明,并说一些踩过的坑。1、代码批量裁剪(clip)代码:import arcpyarcpy.env.workspace="工作路径"rasters=arcpy.ListRasters("*","tif")mask="相应文件"#此处只能输入相应shp文件#否侧可能有:The clip feature is outsi原创 2021-01-26 15:31:48 · 2812 阅读 · 2 评论 -
matlab 批量修改文件名常见错误
问题我在使用matlab对删格文件进行计算时,发现由于命名问题,matlab读取文件顺序并不是按照顺序运行,例如,我的文件命名是“1.tif’,‘2.tif’,……‘11.tif’,‘12,tif’。那么,matlab批量读取的顺序便是’1.tif’,‘10.tif’,‘11.tif’,‘12.tif’,‘2.tif’……‘9.tif‘。为了解决这个问题,我需要将长度较小的文件前面加一个’0‘,这样,读取顺序会变为:’01.tif’,‘02.tif’……'12.tif’。批量修改文件名的操作有很多种原创 2021-01-15 11:42:56 · 1436 阅读 · 1 评论 -
基于MRT的MODIS数据批量拼接与投影——以MOD16为例
原本是打算尝试用python的,但是MRT毕竟是专门用来批量处理MODIS数据的工具,加之本人的python并不算熟练,与MOD16本身的一些数据特点,最终还是选择了MRT,当然,也是我比较懒惰的缘故。看了很多篇博文教程,过于零碎,在这里总结一下,顺便解释一下MOD16让我纠结的点。MRT的安装建议照着这篇做:MRT下载及安装数据MOD16的hv2504和hv2506两个通道,2005-2014年数据,时空分辨率1km,8day。纠结之处当我兴冲冲地用ENVI打开2005001的两景影像时,原创 2020-07-23 17:54:45 · 5324 阅读 · 2 评论 -
ANUSPLIN气象插值的那些坑
ANUSPLIN 是提供实用的转换分析和对多变量数据采用薄 盘光滑样条插值进行插值的工具。它提供了完整的统计分析、数据诊断以及空间分布标准误。同样也支持多种数据输入和表面查询功能。ANUSPLIN可以同时处理几个表面,因此在进行气候要素数据内插的批处理中显得更为方便。ANUSPLIN支持输入数据的多种格式变换,提供复杂的统计分析、透明的数据诊断过程和栅格化的拟合曲面以及标准误差曲面输出,因此目前已被广泛使用。对于样条空间插值的原理,请读者自行搜索。本文主要对ANUSILIN的使用步骤进行简单总结。数据原创 2020-07-10 15:13:03 · 6917 阅读 · 33 评论