![](https://img-blog.csdnimg.cn/20201014180756928.png?x-oss-process=image/resize,m_fixed,h_64,w_64)
python
文章平均质量分 58
我,贝贝
这个作者很懒,什么都没留下…
展开
-
卫星轨道资料洋面风的绘图
其中u和v是二维的numpy.ndarray格式的。但是卫星轨道资料并不是二维格点的。它在解码后是dataframe格式的。所以,要画图,需要把卫星轨道资料插值到格点上,然后再画图。这里用meteva库的cressman插值函数进行站点到格点的插值。我有几个困惑:1是洋面风的色标好丑;2是轨道之外虽然没有数据,但还是有点,可不可以让轨道之外没有点?原创 2022-12-07 15:09:44 · 450 阅读 · 1 评论 -
构造一个xarray.Dataset
按照xarray的官方文档来学习xarray和dataset、dataArray等数据的结构和使用。构造一个维度为lon、lat的数据得到的ds变量为原创 2022-12-06 16:27:20 · 933 阅读 · 0 评论 -
用Python读中法卫星CFOSAT轨道资料的时候遇到的时间转换问题
CFOSAT轨道文件是nc格式的,用xarray.open_dataset读取后,它的时间是b’20221115T083925:456’的格式的。仔细看一下,这个时间的写法就很奇怪,00时60分00秒为什么不写成01时00分00秒呢?这样更符合使用习惯啊。结果上面的转换出现了报错,发现报错出现在时间为b’20201004T006000:456’时。转换思路:Y年M月D日H时M分S分 = Y月M日00时00分00秒+ H小时+M分钟+S秒。也许就是这样的奇怪写法导致了字符串转换为时间戳的时候报错。原创 2022-11-15 09:31:28 · 279 阅读 · 1 评论 -
用meteva画一定经纬度范围内的格点实况产品的图
用meteva画格点实况产品原创 2022-09-27 20:02:32 · 888 阅读 · 0 评论 -
EC预报的读取
提示:文章写完后,目录可以自动生成,如何生成可参考右边的帮助文档文章目录前言一、选择xarray来解码EC文件二、解码所有类型的EC文件1.所有类型的EC文件示例2.grib1和grib2后缀名的读取3.bz2后缀名的读取总结前言EC预报有很多类型,如下表所示提示:以下是本篇文章正文内容,下面案例可供参考一、选择xarray来解码EC文件参考这个博文二、解码所有类型的EC文件1.所有类型的EC文件示例里面有混进一些中国模式。从后缀看有bz2、bin、grib1、grib2。.原创 2022-04-15 15:25:16 · 1830 阅读 · 1 评论 -
EC预报下载读取
ECMWF官方网站上的数据集页面用户指南不花钱的实时开放数据空间分辨率:0.4°时间分辨率:3h时间范围:近四天的数据格式:grib2不花钱的实时数据(0.4°)的地址(4天)数据说明1数据说明2MUSIC上的用了接口调用测试服务后,得到的结果如图所示。天擎上的数据来源: 借助中国科技网中美俄环球科教网络(Global Ring Network for Advanced Applications Development,GLORIAD)通过国际通信系统传输至国家气象信息中原创 2022-03-31 10:45:39 · 4012 阅读 · 3 评论 -
气象数据下载和读取
登陆国家气象业务内网,可以用手机验证码登陆,然后进气象大数据云平台(天擎)。不不不,结果证明,不是这个网址,而是本省的网址http://10.155.64.50:8088/cmadaas。然后在用户空间,申请一个业务用户。请信息中心给通过。2.申请国家海洋科学数据中心账号进入国家海洋科学数据中心,先用手机号注册个人账号,然后升级成单位账号。3.NDBC浮标资料下载进入NDBC浮标的主页,点击Web Data Guide4.ascat资料的下载ASCAT的下载网址,可以直接原创 2022-03-24 09:55:20 · 9599 阅读 · 3 评论 -
踩坑实录——关于我升级了conda以后被迫重装anaconda的故事
提示:文章写完后,目录可以自动生成,如何生成可参考右边的帮助文档文章目录前言一、pandas是什么?二、使用步骤1.引入库2.读入数据总结前言脑子一抽升级了Numpy到1.19,然后可能是兼容性问题吧,anaconda里面import任何库都找不到了。试过用 conda install --revision,没有办法回滚于是尝试重装重装一直失败,提示failed to create manus尝试各种办法都不成功尝试不安装anaconda官网的版本,而是安装清华的版本安装过程中出现无法原创 2022-03-24 09:01:29 · 559 阅读 · 0 评论 -
下载一段时间的ERA5的资料
主要参考了这篇文章网站注册后点击这里查看url和key我增加了定义时间段的部分和保存文件名的部分。以下是代码import datetimeimport timefrom datetime import timedeltaimport pandas as pd```def get_time_df(start_time_btc,end_time_btc): #转换为世界时 start_time_t_utc = datetime.datetime.strptime(start_t原创 2022-03-24 09:00:11 · 573 阅读 · 0 评论 -
用Python生成时间序列
有的时候需要一个时间序列,比如,需要下载某个时间段的数据,或者需要给某段数据加上时间列。举个例子,我想要下载2022年3月24日08时起一天的逐小时数据。那么,我可以用以下代码实现。import pandas as pdimport datetimea= pd.date_range(start =datetime.datetime.strptime('2022032408','%Y%m%d%H'),periods=24,freq='H')print(a,type(a))得到的结果...原创 2022-03-24 08:38:19 · 4358 阅读 · 0 评论 -
python的dataframe填补时间序列中缺少的部分
构造一个dataframeimport pandas as pdimport datetimea = pd.DataFrame({'datetime':['2021092900','2021092902','2021092903'],'value':[1,2,3]})可以看到作为一个时间序列,中间缺了一个时次。原创 2021-11-19 19:56:33 · 2964 阅读 · 1 评论 -
学习python的xarray库
对于初学者,推荐看大佬“摸鱼咯”的系列博文学习。一、读取文件:得到dataset#从文件中读取import xarray as xrimport cfgribds = xr.open_dataset('G:/NEW_RAIN_QC/INPUT/HRCLDAS/TAIR-2021060904.GRB2',engine='cfgrib')读取得到的ds的类型是xarray.Dataset二、读数据:从dataset到dataarray(一)层层查看其中的Dimensions,Coordin原创 2021-09-28 22:54:14 · 4237 阅读 · 0 评论 -
试用python库meteva(气象行业)
meteva的帮助文档看这里使用过程中发现的一些问题及其修改方式:(1)站点插值到格点时,插值函数报错:%matplotlib inline%load_ext autoreload%autoreload 2import meteva.base as mebimport numpy as npimport pandas as pd#读文件file = 'G:/multi_elements_and_qc/station_hour/output/huitu_PRS.csv'data = pd原创 2021-09-17 16:38:03 · 5694 阅读 · 4 评论 -
在windows下用Python读grib2格式的三维云实况产品文件
前面已经写过如何在windows下安装wgrib,见这个博文实际应用中,发现grib2文件的读法还是有挺多需要细写的,这里以国家气象信息中心的三维云实况产品(Z_NAFP_C_BABJ_20210830203447_P_3DCloudA_RT_CHN_0P05_HOR-CCP3-2021083020.GRB2)为例。方法一:用wgrib将grib2格式转换为csv格式读取查看文件信息,在命令提示符下或者anaconda下,输入命令wgrib2 -v0 Z_NAFP_C_BABJ_2021083020原创 2021-09-08 19:27:34 · 1765 阅读 · 0 评论 -
读FY3D的海面风场数据(HDF格式)
文件下载注册了中国气象数据网的账号(教育科研级别),然后登陆风云卫星遥感数据服务网,找到相应的数据,然后申请,然后得到下载文件的目录,用页面提示的使用*.bat的批量下载的方式进行下载(需要预先在电脑上安装wget)。HDFVIEW安装HDFVIEW,我装的是HDFView-3.1.0-win10vs14_64。FengYunToolkit在风云卫星遥感数据服务网的工具页面可以找到风云卫星数据处理包。文件读取...原创 2021-06-23 11:14:10 · 2378 阅读 · 1 评论 -
解码和读取GRB2格式的气象文件
解码和读取GRB2格式的气象文件使用pygrib新的改变功能快捷键合理的创建标题,有助于目录的生成如何改变文本的样式插入链接与图片如何插入一段漂亮的代码片生成一个适合你的列表创建一个表格设定内容居中、居左、居右SmartyPants创建一个自定义列表如何创建一个注脚注释也是必不可少的KaTeX数学公式新的甘特图功能,丰富你的文章UML 图表FLowchart流程图导出与导入导出导入使用pygrib这里使用的GRB2格式文件是国家气象信息中心的CMPAS降水融合资料、HRCLDAS陆面融合实况资料、CLD原创 2021-05-14 13:03:26 · 4733 阅读 · 0 评论 -
学习tkinter
学习使用tkinter学习参考的博文有:1.这位博主写得很全,还有例子,适合我这样的初学者2.tkinter路径选择功能的实现3.tkinter获取文本框的值4.这个页面的字符串转MD5工具的例子的tkinter的class的写法挺好5.如何通过循环生成多个复选框checkbutton其他问题的解决:1.设置messagebox的大小:window.option_add('*Dialog.msg.width', 50)2.messagebox的import问题:要直接import,这样写im原创 2021-03-20 09:43:18 · 118 阅读 · 0 评论 -
python的.py文件在anaconda环境中打包成.exe可执行文件
python的打包:我是在jupyter上编程的,先从jupyter下载代码到本地,存为.py格式。然后,打包方法参照这位作者的博文(非常感谢!)。其实,我是尝试了pyinstaller后,打包的exe运行会闪退,不会出现GUI界面,我也不知道怎么改,因此换用cx-freeze。具体步骤很简单:在anaconda上安装cx-freeze。pip install cx-freeze。然后,在anaconda下进行打包cxfreeze g:\a.py --base-name=Win32GUI,不知道为什么加原创 2021-03-19 16:44:48 · 1277 阅读 · 3 评论 -
python读取文件出现的问题
python读取文件中 的矩阵file = open(‘filename’) a =file.read() b =a.split(‘\n’) for i in range(len(b)): b[i] = b[i].split()#使用空格分开len(b[0]) import numpy as npb = np.array(b)#转成numpy形的b.shapet...原创 2020-02-27 12:33:44 · 205 阅读 · 0 评论 -
用python比较两个文本
# -*- coding: utf-8 -*-#!/usr/bin/pythonimport difflib def string_similar(s1, s2): return difflib.SequenceMatcher(None, s1, s2).quick_ratio()# 打开文件fo = open("aa.txt",encoding='UTF8')f2 =...原创 2020-02-26 21:17:57 · 836 阅读 · 1 评论