- 博客(18)
- 问答 (1)
- 收藏
- 关注
RSS订阅原创 Python绘制GPM二级数据降水剖面图
我是学到一点就迫不及待想分享给大家，如有问题咨询，我会尽可能回复，私信最多一天回。图中黑线为0摄氏度等温线，红色为地形数据，可依据需要进行保留或删除。GPM2级数据丰富和强大，大部分都进行了注释。BlueWhiteOrangeRed等颜色网站如链接所示。参考链接在上一篇有附上。
2023-04-02 20:36:32
1341
转载 python读取GPM二级降水等卫星数据及思路
对于代码结果如图所示，之前网上找了很久没找到，在气象家园看到了大佬写的。有些人可能没注册，无法访问。本代码进行了详细注释。卫星扫描数据与网格数据最大不同，是在数据处理方面。如果不想下载可私信我，我将文件发送给您。文章以转载形式表达对原文的敬意。
2023-03-25 18:44:03
1917
9
原创 Python 使用seaborn画箱型图
文中使用seaborn画两个箱型图，seaborn不能画两个y轴的箱型图表示两个物理量，可以使用matplotlib画两个箱型图，中间大段全是重复。
2023-03-05 16:41:03
891
原创 python检查下载文件缺漏
如图我所下载的modis文件使用官网自动下载，总是有些文件漏掉，根据文件名发现其特点为中间A后面数字为年份＋天数，通过代码自动读取，当有遗漏时，代码会自动将所有遗漏的爆出来，根据文件名时间进行排查。代码很简单，只是对于字符串的运用。由于各种原因，常常自己所要下载的大量文件无法一次性完整下载完毕，且文件多有遗漏，利用python可以非常简单的进行排查遗漏的文件。
2023-02-07 20:15:42
401
原创 python区域掩膜白话，显示目标区域
一个简单的使用，就是读取文件，然后设定一个目标区域（我代码中设定的是亚洲区域）,然后进行填色图绘制，shp2clip函数需要防止，填色好的区域，ax，以及大的shp地图文件，最后是目标区域，从大的shp文件中裁剪。region：只想要一个区域就只放一个区域名字，想要多个区域就放多个区域名字，非常简单好用。查看的结果如下，因为不同shp文件里的排列顺序不同，图中是第三个位置是国家名字，所以是record[2],要是第四个是国家名字，就把2改成3！这是选择shp文件的代码。
2022-11-20 12:37:54
941
2
原创 Python 读取不规则位数的txt文件
针对txt文件中的数据个数位数不统一情况，并且存储不是每行数据个数都一致的情况，利用多种情况读取空格进行判断，只要读取到时空格则判断与上一次为同一个数据，并且利用多重循环，首选循环每一个大数据组，在循环每一个小数据组的每一行，在循环每一行的每一个读取到数组中。
2022-11-18 13:08:04
1048
原创 Python 处理dat文件并画图
思路：首先每行读取数据按照经纬度，数据是由北向南，由西向东的。通过np.delete()删除每行年份和月份。然后将字符转化成数字读取转化成三维数据，最后进行数据处理后绘图就行。上面没有经纬度，只有对应经纬度的数据，因此先进行数据处理。
2022-10-17 17:58:30
3024
5
原创 Python选取nc文件特定经纬度，时间和气压
有时候下载数据会下载很多时间点很多气压场的数据，通过筛选特定时间点和特定气压场的数据进行数据处理。
2022-09-08 20:29:49
6008
原创 Python用Basemap画位势高度距平图和5880线
本文数据为NCEP/NCAR 再分析数据，感觉这个真是比ERA5好用，就是分辨率和内容类型比ERA5差。链接如下：https://psl.noaa.gov/data/gridded/data.ncep.reanalysis.html。本文比较初级，就是画了等值线还有一个画出目标区域的线，shp文件网上有很多，想要可以发邮箱在评论区，我把我收藏的都给可以奉上。文中我把一个大文件裁剪成30个小文件，毫无必要，这个可以再优化下。希望能给大家提供帮助。...
2022-06-09 21:09:13
2916
原创 C语言实现矛盾方程求解
矛盾方程组，非常好理解，其实为翻译上述图片，当你图片看懂了，原理就懂了：代码实现：#include<stdio.h> # define M 100int main(){ int i,j,n,w1,w2,k1,k2,t,r,k,w,e,o1,o2,e1;double x[M],y[M],a[M][M],b[M][M],c[M],d[M],h[M][M],max,f,s; //录入数据 printf("输入未知数的个数n："); scanf("%d",&n...
2021-12-17 21:32:08
1435
原创 C语言实现拉格朗日插值法和牛顿插值法
拉格朗日插值法的实现就是通过构造函数，通过已知几对数据然后估算输入对应的数据所得得值。通过例子可知道拉格朗日法,假设有三组数据分别为（x0,y0）,(x1,y1),(x2,y2)，与第一个y值相乘得数的分母是(x0-x1)*（x0-x2),分子是：（x-x1)*(x-x2),第二个数就是分母是（x1-x0)*（x1-x2),分子是（x-x0)*(x-x2),第三个数分子是：（x-x0)*(x-x1),分母是（x2-x0)*(x2-x1).代码实现：#include<stdio.h..
2021-12-17 21:22:30
8313
2
原创 C语言实现方程组LU分解法和列主元LU分解法
根据老师所讲LU分解法和列主元LU分解法本质就是高斯消去法和列主元高斯消去法，将L部分变为都变为0，则就化为高斯消去法，我们在方程求回代的过程中只用到了U和y，不需要用到L，那个L部分实际保存了每一次运算的倍数关系，详细可自己举数学例子验证便可知。以例子解释规律，我们以（1，1）表示第一行第一列数,则（3，1）表示第三行第一列，那么第一行所有数保持不变，第一列数全部除上第一行第一列。然后原来（2，2）-(2,1)*(1,2)即：5-0.7*7。 6-0.7*8，5-0.7*7而第二列为（6-0.
2021-12-17 20:39:44
3155
原创 C语言实现高斯消去法和列主元高斯消去法
本篇主要实现高斯消去法和列主元高斯消去法 高斯消去法和列主元高斯消去法都是为了解线性方程组的有效方法，但列主元高斯消去法是高斯消去法的一个优化版本，强烈建议后面许多地方用到解方程组时，都用列主元高斯消去法。 高斯消去法：我个人觉得，例子比数学公式更好让人理解本质。上述是线性代数的高斯消元法，而我们只需要用程序将这个代码实现出来。上述方程组，第二行减去第一行，第三行减去两倍第一行。那么可以推断，若第三行要消去第一个数，则让第一行每个数除上第一行第一...
2021-12-17 19:57:53
8423
2
原创 C语言编写Picard迭代和牛顿迭代法
为了完成学校计算方法的作业，花了两个多星期将所有实验报告编写出来，看到同桌找代码的辛苦，CSDN也没有这类代码，我只是一时兴起，想分享我写的代码与大家共飨，有错误还请多多指正，代码都不难，也是为了让CSDN帮我保存下，我辛辛苦苦做出来的，倘若能给别人带来一丝一毫的帮助，我都万分荣幸。 1：什么是Picard迭代 先弄一个满足f(x) = 0的函数将其转化为x=正在上传…重新上传取消...
2021-12-17 18:59:23
4127
4
空空如也
怎么用python打开grib文件，window系统
2021-08-14
TA创建的收藏夹 TA关注的收藏夹
TA关注的人