- 博客(37)
- 收藏
- 关注
RSS订阅原创 色标在matplotlib和plottable中
是这样的，我有一个数组[-4.4, -2.8, -2.6, -2.2, -1.1, 1.1, 1.2, 1.3, 3.6, 6.0, 6.4, 12.3]，它需要绘制散点图，点的颜色来代表数值大小；同时，也需要在plottable上作为一列显示，同样用颜色来代表数值的大小。还有两个额外的需求：一、不同的数值的颜色要有区分度；二、0值需要是白色，正值为暖色系，负值为冷色系，从视觉上明显区分正负值。
2024-01-09 14:27:36
913
原创 读HDF5格式的文件
HDF（Hierarchical Data Format）指一种为存储和处理大容量科学数据设计的文件格式及相应库文件。其中是层级结构的意思，也就是说，HDF的文件结构是分层的，像文件夹一样，上层文件夹套下层文件夹。
2023-09-22 15:53:26
345
原创 卫星轨道资料洋面风的绘图
其中u和v是二维的numpy.ndarray格式的。但是卫星轨道资料并不是二维格点的。它在解码后是dataframe格式的。所以，要画图，需要把卫星轨道资料插值到格点上，然后再画图。这里用meteva库的cressman插值函数进行站点到格点的插值。我有几个困惑：1是洋面风的色标好丑；2是轨道之外虽然没有数据，但还是有点，可不可以让轨道之外没有点？
2022-12-07 15:09:44
431
1
原创 构造一个xarray.Dataset
按照xarray的官方文档来学习xarray和dataset、dataArray等数据的结构和使用。构造一个维度为lon、lat的数据得到的ds变量为
2022-12-06 16:27:20
892
原创 用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
272
1
原创 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
1684
原创 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
3746
3
原创 气象数据下载和读取
登陆国家气象业务内网，可以用手机验证码登陆，然后进气象大数据云平台（天擎）。不不不，结果证明，不是这个网址，而是本省的网址http://10.155.64.50:8088/cmadaas。然后在用户空间，申请一个业务用户。请信息中心给通过。2.申请国家海洋科学数据中心账号进入国家海洋科学数据中心，先用手机号注册个人账号，然后升级成单位账号。3.NDBC浮标资料下载进入NDBC浮标的主页，点击Web Data Guide4.ascat资料的下载ASCAT的下载网址，可以直接
2022-03-24 09:55:20
9310
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
530
原创 下载一段时间的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
559
原创 用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
4313
原创 python的dataframe填补时间序列中缺少的部分
构造一个dataframeimport pandas as pdimport datetimea = pd.DataFrame({'datetime':['2021092900','2021092902','2021092903'],'value':[1,2,3]})可以看到作为一个时间序列，中间缺了一个时次。
2021-11-19 19:56:33
2902
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
4057
原创 试用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
5019
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
1615
原创 WPS：找文件的历史版本
如何找到历史版本点击“首页”——“最近”——选择文件——右键“历史版本”，然后就可以在右栏看到历史版本，点击打开自己想要的历史版本。如何让WPS自己保存点击“首页”——“设置”——“设置”——“文档云同步”选择开...
2021-09-08 10:59:58
4556
原创 在WSL中安装ifort
首先，在intel的官网下载适合linux系统的OneAPI，这是链接我采取了直接下载*.sh文件，然后在ubuntu内键入命令的方式进行安装。sudo bash l_HPCKit_p_2021.3.0.3230.sh 并且参考安装指南，采用非互动式（静默）的安装方式，键入的命令实际上为sudo bash l_HPCKit_p_2021.3.0.3230.sh -s -a --silent --eula accept安装结束后的环境变量设置参考这个博文。但是，在键入命令source ~/
2021-06-29 16:55:59
844
1
原创 读FY3D的海面风场数据（HDF格式）
文件下载注册了中国气象数据网的账号（教育科研级别），然后登陆风云卫星遥感数据服务网，找到相应的数据，然后申请，然后得到下载文件的目录，用页面提示的使用*.bat的批量下载的方式进行下载（需要预先在电脑上安装wget）。HDFVIEW安装HDFVIEW，我装的是HDFView-3.1.0-win10vs14_64。FengYunToolkit在风云卫星遥感数据服务网的工具页面可以找到风云卫星数据处理包。文件读取...
2021-06-23 11:14:10
2188
1
原创 解码和读取GRB2格式的气象文件
解码和读取GRB2格式的气象文件使用pygrib新的改变功能快捷键合理的创建标题，有助于目录的生成如何改变文本的样式插入链接与图片如何插入一段漂亮的代码片生成一个适合你的列表创建一个表格设定内容居中、居左、居右SmartyPants创建一个自定义列表如何创建一个注脚注释也是必不可少的KaTeX数学公式新的甘特图功能，丰富你的文章UML 图表FLowchart流程图导出与导入导出导入使用pygrib这里使用的GRB2格式文件是国家气象信息中心的CMPAS降水融合资料、HRCLDAS陆面融合实况资料、CLD
2021-05-14 13:03:26
4403
原创 今日困惑：fortran不能读带中文的文件了？
我的fortran是VS2017+intel fortran 2019，今天发现只要读带中文的文件，就会读取错误，为啥？
2021-04-17 19:00:42
632
2
原创 安装intel fortran
主要参考了这个博文安装好后又遇到了““无法查找或打开PDB文件””的问题，解决方法参照这个博文，然后按照这个博文的一个评论，按Ctrl+F5调试运行程序（不使用本地调试），然后成功了。...
2021-03-26 07:51:51
1079
原创 学习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
107
原创 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
1232
1
原创 窗口排列
1.动态切换WIN7系统下，点住WINDOWS，然后点击TAB，停下在哪个窗口，就打开哪个窗口2.窗口排列鼠标右键任务栏，“堆叠显示窗口”是上下分布显示窗口，“并排显示窗口”是左右排列显示窗口。只显示未被最小化到任务栏的窗口。...
2020-12-17 15:05:35
1097
原创 WSL及库的安装
安装WSL以及所需的库安装WIN10系统安装WSL安装WSL的相关库安装WIN10系统我的目的是用WIN10自带的LINUX系统，即WSL。同事跟我说最方便安装的方式是在microsoft store里面直接安装。但是我的笔记本上找不到microsoft store，可能因为我的笔记本上预装的系统是win10神州政府版。期间也试过自己不通过microsoft store来装ubuntu，但是失败了。所以，决定先重装WIN10系统。然后，在网上找了几个一键安装系统的安装包，雨林木风的，深度系统的，最
2020-12-05 23:00:01
456
原创 在c++中拼接字符串，strcat的用法
（1）strcat(str1,str2),将str2拼接到str1后面，返回str1(2)str1的声明方式必须是char str1[100]=""(3)str2的声明方式可以是char str2[10]=“123” 或者char* str2=“123”
2020-04-01 08:45:40
1700
1
原创 Microft Visual Studio 2010报错“无法解析的外部符号”
我的程序是VERAQC_TEST.f90是主程序program；而vera_qc.f90是module，所有的功能块，各种subroutine和function都在里面。从报错的内容来看，是module里面对function和subroutine的link出了问题。可是，从语法上看，没有任何问题。试着把报错的subroutine和function剪切到program文件中。就可以转通了。...
2020-03-09 20:39:29
259
原创 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
200
原创 用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
825
1
原创 已知文件url，批量下载文件
打开迅雷，新建任务，将url批量复制到地址框内（最多1000条）。开始下载。 如果下载速度长时间为0。重启迅雷，重新开始下载。
2020-02-22 23:46:25
2482
原创 第一篇：希望养成记录的好习惯
先简单记录一些代码激活anaconda 里的ncl环境conda activate NCL_py在python里引入ncl库import Ngl , Nio
2020-02-20 21:33:39
201
空空如也
空空如也
TA创建的收藏夹 TA关注的收藏夹
TA关注的人