Python滞后相关系数（Lagged correlation）代码分享，气象相关

最近研究南海海温时用到了这个，向大家进行一个代码分享

需要：catorpy，numpy，scipy

import matplotlib.pyplot as plt
import numpy as np
from scipy.stats import pearsonrdef LCC(data1, data2, n, path):
#   作滞后相关，n>0时，data2先于data1发生
#   如data2为海温，data1为降水，n为1——>LCC为海温关于当年降水和来年降水的相关a = n * 2 + 1b = len(data1)x = np.arange(-n, n + 1, 1)r = np.zeros((a))p = np.zeros((a))for i in range(a):if i < (n):r[n - i - 1], p[n - i - 1] = pearsonr(data1[:(b - i - 1)], data2[i + 1:])else :r[i], p[i] = pearsonr(data1[x[i]:], data2[:b - x[i]])
#   画图fig = plt.figure()ax = fig.add_subplot(111)ax.axhline(0,linewidth=3, color='k')   #零刻度线cs = ax.bar(x[n:], r[n:], align='center',edgecolor='k',linewidth=3,color = 'white')
#            从无滞后相关开始绘图，对齐方式为中心对齐count = nfor bar, height in zip(cs, r):         #开始显著性检验if p[count] < 0.01:bar.set(hatch='///')           #通过99%的显著性检验为斜线if p[count] < 0.005:bar.set(hatch='||')            #通过99.5%的显著性检验为竖线if p[count] < 0.001:bar.set(hatch='--')            #通过99.9%的显著性检验为横线count = count+1ax.bar_label(cs)                       #标注plt.title('cor_max='+ str(np.max(r))+'\n'+'cor_min='+str(np.min(r)))plt.show()plt.savefig(path)plt.close()

Python滞后相关系数（Lagged correlation）代码分享，气象相关相关推荐

python扫雷游戏实验分析_用python写扫雷游戏实例代码分享
扫雷是一个非常经典的WIN游戏,我们教给大家用python语言来写出这个游戏,以下是全部实例代码: #!/usr/bin/python #coding:utf-8 #python 写的扫雷游戏 imp ...
python画正方形的代码_python绘制正方形螺旋线的代码分享
python绘制正方形螺旋线的代码分享发布时间:2020-04-28 09:56:56 来源:亿速云阅读:219 作者:小新这篇文章主要为大家详细介绍了python绘制正方形螺旋线的代码分享,文 ...
圣诞节要来了，程序员女友用Python画圣诞树的浪漫代码作为最好的礼物
我从来没有过过一个真正意义上的圣诞节. 没有拥有过圣诞树, 没有和恋人一起看过雪,也没有收过精心准备的圣诞礼物. 我给自己安排很多与节日无关的事来度过那一天. 尽管如此,我还是很期待圣诞节. 哈哈哈哈 ...
电脑代码python_python发代码分享 -电脑资料
复制代码代码如下: import urllib.parse,os.path,time,sys,re,urllib.request from http.client import HTTPSConnec ...
基于python的音频播放器_基于python实现音乐播放器代码实例
基于python实现音乐播放器代码实例,一首,函数,按钮,布局,音乐基于python实现音乐播放器代码实例易采站长站,站长之家为您整理了基于python实现音乐播放器代码实例的相关内容. 核心播放 ...
摸鱼气象Python教案代码分享
摸鱼气象Python是大气科学领域内的一套有视频讲解,项目实战类的长期更新的教程. 其中涵盖了大量气象领域Python初学者需要用到的编程案例,如: 数据读取.数据处理.数值计算.数据可视化.数据分析 ...
python登录代码思路_用python登录Dr.com思路以及代码分享
用python登录Dr.com思路以及代码分享发布于 2014-08-28 22:31:52 | 192 次阅读 | 评论: 0 | 来源: 网友投递 Python编程语言Python 是一种面向对 ...
python turtle画画 30排以内_Python竟能画这么漂亮的花,帅呆了(代码分享)
阅读本文大概需要3分钟关于函数和模块讲了这么久,我一直想用一个好玩有趣的小例子来总结一下,同时也作为实战练习一下. 趣味编程其实是最好的学习途径,回想十几年前我刚毕业的时候,第一份工作就给手机上写a ...
python爬取网页书籍名称代码_python爬取亚马逊书籍信息代码分享
我有个需求就是抓取一些简单的书籍信息存储到mysql数据库,例如,封面图片,书名,类型,作者,简历,出版社,语种. 我比较之后,决定在亚马逊来实现我的需求. 我分析网站后发现,亚马逊有个高级搜索的功能 ...

Python滞后相关系数（Lagged correlation）代码分享，气象相关

Python滞后相关系数（Lagged correlation）代码分享，气象相关相关推荐

最新文章

热门文章