审查委员会-审查请求绘制使用Python

审查委员会这些项目中有一个是Python社区一直引以为傲的吗事实上的现在代码审查的标准。审查委员会有一个有趣的an亚洲金博宝非常有用的REST API您可以使用它来检索关于代码审查、注释、差异和存储在数据库中的每一个单独的信息。亚洲金博宝因此,我创建了一个小的Python脚本,名为rbstats它检索用户完成的审查请求的信息,然后使用matplotlib。我将展示在审查委员会制度Apache基金会的成员。

要使用该工具,只需使用Review Boars系统的API URL调用,即:

python rb_stats.py——max-results 80 https://reviews.apache.org/api

然后你会得到一个这样的图形点击放大):

其中“最热”点根据用户为另一个axis用户创建的代码检查的数量进行加权。还可以根据用户(例如用户)绘制统计数据benjaminhindman使用秋天colormap和400 max结果:

python rb_stats.py——max-results 400——来自用户benjaminhindman——colormap autumn https://reviews.apache.org/api

点击放大图片:

用Python实现黎曼函数的可视化

mpmpath这是黎曼ζ函数求值程序,我使用Matplotlib (the源代码)。

黎曼函数是一个解析函数它定义在复平面上有一个复变量表示为年代”。黎曼在数学中非常重要因为它和质数亚洲金博宝有很深的关系;函数为:

\泽塔(s) = \ sum_ {n = 1} ^ \ infty \压裂{1}{n ^ s} = \压裂{1}{1}^ s + \压裂{1}{2 ^ (s} + \压裂{1}{3 ^ s} + \ cdots

(年代)> 1

所以,让\泽塔(s) = v在哪里s = a + b\imathv = u + k\imath

第一个情节使用了三连音(x, y, z)绘制三维空间的坐标,其中各分量为:

  • x =再保险(v)(或x = u,之前定义);
  • y = Im (v)(或y = k,之前定义);
  • z = Im (s)(或z = b,之前定义);
  • 动画中使用的时间组件被调用\θ它由\θ=再保险(s)或者简单地\θ=一个

为了绘制动画,我使用了以下间隔:

  • 再保险(s):从\[0.01,10.0 \左右),每0.01步计算一次,如图右上角的图所示;亚洲金博宝
  • Im (s):从\[0.1,200.0 \左右),以每0.1步计算-如亚洲金博宝z轴。

此图使用固定的间隔(不自动缩放)(x, y, z)坐标。在哪里(s) = 1/2(\θ= 1/2)为黎曼函数的非平凡零存在时。

现在看到相同的图形,但这次使用自动缩放(自动调整大小)x, y坐标):

注意(x, y)自动伸缩功能。

现在从另一个使用2D空间的绘图,其中每个组件是:

  • x = Im (s)(或x = b,之前定义);
  • y = Im (v)(蓝色)和再保险(v)(绿色);
  • 动画中使用的时间组件被调用\θ它由\θ=再保险(s)或者简单地\θ=一个

为了绘制动画,我使用了以下间隔:

  • 再保险(s):从\[0.01,10.0 \左右),每0.01步计算一次——见顶部图标题;亚洲金博宝
  • Im (s):从\[0.1,50.0 \左右),以每0.1步计算-如亚洲金博宝x轴。

此图使用固定的间隔(不自动缩放)(x, y)坐标。在哪里(s) = 1/2(\θ= 1/2)为黎曼函数的非平凡零存在时。黎曼函数的前10个非平凡零被显示为一个红点,当两个级数,Im (v)再保险(v)在红线的红点上互相交叉((s) = 1/2)是函数的0的位置,注意实部和虚部是如何变成的再保险(s)增加。

现在看到相同的图形,但这次使用自动缩放(自动调整大小)y协调):

如果你对黎曼Zeta函数的更多形象化感兴趣,你会喜欢那篇很棒的论文j . Arias-de-Reyna被称为“黎曼泽塔函数的x射线”。

我一直很喜欢可视化对数学函数理解的影响。安斯科姆的四方是可视化的重要性的一个清晰的例子。

源代码用于创建的情节可在这里:

2D图形的源代码

3D图的源代码

希望你喜欢这篇文章!来制作我所使用的情节和视频matplotlib,mpmathMEncoder

- Christian S. Perone

使用Python使用近实时MODIS Aqua+Terra卫星图像探测实时海地地震数据亚洲金博宝

不久前,我发表了一篇关于使用Python和NASA MODIS公开的数据处理MODIS近实时卫星图像的文章,如果你对这些图像感兴趣,可以看看这里的帖子

现在我已经创建了一个名为pyearthquake自动检索任何MODIS子集图像从美国宇航局快速反应系统在一个简单的方法,只需使用子集,卫星和分辨率参数。这个包有一个模块可以从美国地质调查局获取实时地震数据美国地质勘探局网站并自动解析它,以及从USGS检索shakemaps和其他产品的功能。

检索、处理和绘制USGS数据

让我们开始讨论公共可用的目录与实时地震数据从USGS网站(他们是CSV格式,还有其他XML和KMZ格式,如果你感兴趣,但我重点在CSV格式文件创建Python模块:

USGS站点注意:文件不包括在内(例如,过去一天的文件不包括过去一小时的事件)。

M1 +地震(过去时)M1 + PAST_HOUR

这是过去一小时地震数据的全球目录;

M1 +地震(过去的)M1 + PAST_DAY

这是包含过去一天地震数据的全球目录;

M1 +地震(过去7天)-M1 + PAST_7DAY

这是世界范围内的地震目录,包括过去7天的地震数据;

要从这些目录中获取和处理数据,你可以使用我创建的pyearthquake包,让我们先安装它:

easy_install pyearthquake

这个命令将使用easy_install从PyPI衍生库自动安装包,它也将自动解决需求,但如果您遇到一些问题,以下是显式需求:matplotlib> = 0.99.0,numpy1.3.0 > =版本,公益诉讼> = 1.1.6和基础图> = 0.99.4。basemap包是matplotlib的map工具包。

让我们现在,例如,获取过去一小时地震数据的目录,我们可以使用Python解释器提示符,如下所示:

>>> catalog = usgs.retrieve_catalog("M1+PAST_HOUR")在0x00BD0FD0>处的USGSCatalog实例

所有可用的目录是:M1+PAST_HOUR, M1+PAST_DAY和M1+PAST_7DAY。但如果您键入了错误的目录名称,它将显示所有可用的名称。

您现在可以使用“catalog”对象来查看从USGS站点检索的数据:

> > >列表(目录)[{“Src”:“词”、“地区”:“加州中部”,“朗”:“-118.2063”,“Datetime”:“2010年1月16日星期六,17:12:20 UTC”、“深度”:“6.00”,“版本”:“1”,“纬度”:“36.2061”,“Eqid”:“10530221”,“大小”:“1.9”,“望远镜”:“27”},{“Src”:“正义与发展”,“地区”:“南阿拉斯加”、“朗”:“-150.0829”,“Datetime”:“2010年1月16日星期六,16:47:05 UTC”、“深度”:“34.90”,“版本”:“1”,“纬度”:“61年。4521”、“Eqid:“10029267”、“大小”:“2.3”、“望远镜”:“27”}]行目录> > >:…打印行“级”,行“深度”,行“Datetime”,行“深度”,行(“区域”)……2010年1月16日星期六16:47:05协调世界时34.90阿拉斯加南部

如您所见,“catalog”是一个可交互的对象,您可以在其中获取地震的所有数据,如震级、纬度、经度、深度、日期、地区名称等……

让我们看看过去7天的地震数据:

>>> catalog = usgs.retrieve_catalog("M1+PAST_7DAY"

在过去的7天里,地球上发生了1.142次1级以上的地震,如果我们想要找出上次海地悲剧发生的事件,我们可以将震级按6.0进行筛选:

[size = 6] [event for event in catalog if float(event[" size "]
对于大小为6的行:…打印行“Eqid”,行“级”,行“深度”,行“Datetime”,行“深度”,行(“区域”)……2010年1月14日星期四14:03:40 UTC 32.40马里亚纳群岛以南2010年1月12日星期二21:53:10国际标准时间13:00海地地区2010年1月10日星期日00:27:39 UTC 29.30加利福尼亚北部近海

如我们所见,编号为“2010rja6”的地震,是美国地质勘测局用来确定1月12日发生在海地的7.0级地震,这场地震摧毁了该地区及不幸的海地人民。

美国地质勘探局还提供了Shakemaps,这是一种计算机自动生成的地图,显示了地震地区的潜在破坏情况。

美国地质调查局6级、4级、5级。retrieve_shakemap (haiti_eq“INSTUMENTAL_INTENSITY”)

函数的第二个参数是shakemap类型,你也可以从USGS获得其他产品,比如媒体地图,峰值地面加速度。这些地图在包中有:仪器强度,裸露,装饰,不确定度和PEAK_GROUND_ACC。

packas具有绘制地震地图的功能(使用Matplotlib),参见如何绘制过去7天的地震地图:

> > > usgs.plot_events(目录)
< mpl_toolkits.basemap。Basemap对象0x01E6BDB0>

现在在海地地区,

这些点的颜色与震级有关,我这里使用的是Matplotlib的JET(黄色/红色的点是震级更高的地震)彩色地图,见此链接使用更多的彩色地图。

从MODIS卫星获取最后的图像和绘制地震图

现在我们来谈谈modis"的模块pyearthquake模块。该模块可用于获取MODIS Aqua/Terra卫星上最新的或存档的图像(高达250m分辨率)。亚洲金博宝首先,你必须知道你想从MODIS网站检索的子集,子集是可用的这里是NASA的基地。为了绘制海地区域,我将使用名为“GreaterAntilles”的子集,它覆盖了海地区域。要获得今天的卫星图像,我们可以这样做:

>>> from pyearthquake import modis >>> import datetime >>> now = date .date .now() > >0 >1 bmap = modis。get_modis_subset(现在,…“GreaterAntilles”……satellite_name =亚洲金博宝“terra”……分辨率= " 250 ")

这个命令将获得今天从“terra”卫星上拍摄的“GreaterAntilles”子集的合成图像,分辨率为250米(事实上,这将需要一段时间,因为我使用的是更高的分辨率,最大值为250米)。亚洲金博宝结果如下:

下面是太子港附近的海地地区:

这是MODIS有趣的部分,这张照片是今天拍摄的。

现在,让我们把美国地质勘探局的地震信息和海地最近的图像结合起来:

导入modis catalog = usgs.retrieve_catalog("M1+PAST_7DAY") now = date .date . time.now() bmap = modis。get_modis_子集(现在,“GreaterAntilles”,satellite_name=“terra”,亚洲金博宝分辨率=“250m”,show=False)bmap plot_events(目录)

这是来自Terra卫星的今日MODIS图像中过去7天地震图的结果:亚洲金博宝

你也可以使用这个软件包来绘制世界上其他地区的地震,只要从MODIS快速反应系统站点中选择另一个子集就可以了。

- Christian S. Perone

这里显示的所有图像和数据都来自MODIS响应系统或美国地质调查局。

本福德定律在Twitter上的应用分析金宝搏188官方

本福德定律是我们无法解释的奇怪现象之一,当亚洲金博宝我们发现越来越多的符合这个规律的现象时,我们很惊讶。两人(西蒙·纽科姆- 1881和弗兰克本福德(1938年)同样地在翻阅一本对数表格书的时候;他们注意到这本书的开头几页比结尾几页脏。

目前还没有先天的标准告诉我们什么时候数据集会或不会遵守本福德定律。正因为如此,我做了一个分析金宝搏188官方Twitter公共时间线

获取公共金宝搏188官方时间轴的Twitter API对于这种分析毫无用处,因为在API文档中,他们说公共时间轴缓存60秒!60秒是永恒的,每小时有150个请求。所以,它没有帮助,但是有一个alpha测试API它有来自公共时间轴的漂亮而有用的数据流;亚洲金博宝Twitter流媒体API中有很多方法,其中最有趣的是“Fi金宝搏188官方rehose”,它返回所有公共状态,但这个方法只对intere$ting的人可用,我不是其中之一。但是,我们有" Spritzer ",它返回所有公共状态的一部分,因为它只是我们目前可用的,所以它一定是有用的,它是一个相当不错的数据流=)

所以,我有Spritzer实时数据流和处理每一个新的状态,它到达一个正则表达式来找到状态中的所有数字;如果status为“I have 5 dog and 3 cats”,则收集的数量应为“[5,3]”。所有这些累积的数字然后与本福德定律进行核对。我使用Matplotlib绘制两条曲线(本福德定律和Twitter状态使用数字分布),以经验观察它们之间的相关性。金宝搏188官方你可以在视频的右上角看到皮尔森的相关性也是两个发行版之间的。

这里是视频(我只看到创建后的视频,但曲线的颜色是反向见图例):

该视频代表了15分钟(3.160个捕获状态)的Twitter公共时间线。金宝搏188官方在视频的最后,你可以看到皮尔逊相关性是0。95。看来我们又找到了本福德的另一个儿子。

这个用来处理Twitter Spritzer数据流并实时绘制金宝搏188官方相关图的小工具完全是用Python编写的,我将进行清理并在我有时间的时候发布到这里。该工具已经生成了1823张使用ffmpeg合并的png图像。

我希望你喜欢

引用这篇文章为:Christian S. Perone,《本福德定律应用于推特的分析》金宝搏188官方亚洲金博宝未发现的地域11/08/2009,//www.cpetem.com/2009/08/an-analysis-of-benfords-law-applied-to-金宝搏188官方twitter/

更新:11/08用户“poobare”引用了一篇关于本福德定律的有趣论文Reddit,这是的链接

更多关于本福德定律的帖子

质数和本福德定律

Delicious.com,根据本福德定律检查用户数量

本福德定律符合Python和苹果股票价格

接受近实时MODIS图像和来自美国宇航局Aqua和Terra卫星的火灾数据亚洲金博宝

Modis快速反应系统网站:

开发MODIS快速响应系统是为了几乎实时地每天提供地球陆地的卫星图像。真彩色,像照片一样的图像和假彩色图像在收集几小时内就可获得,这使得系统成为一个有价值的资源(…)。
中分辨率成像光谱仪(MODIS)飞行在NASA的Aqua和Terra卫星上,这是NASA中心的亚洲金博宝国际地球观测系统的一部分。这两颗卫星都绕地球运行,每天都能看到地球的大部分地方。亚洲金博宝板载亚洲金博宝时,MODIS在早上看到地球,而阿卡MODIS在下午绕地球运行。

少数人知道的一件事是,每个人都能获得实时的中分辨率图像从两个美国太空总署的卫星(Aq亚洲金博宝ua和Terra)环绕地球每天从南极到北极,捕捉图像和数据最喜欢热异常的地球的一部分。亚洲金博宝

在这篇文章中,我将展示如何使用PythonMatplotlib+技术工具包

Aqua和Terra绕轨道亚洲金博宝运行

“Terr亚洲金博宝a”轨道

阅读更多