请选择 进入手机版 | 继续访问电脑版
爱气象,爱气象家园! 

气象家园

 找回密码
 立即注册

QQ登录

只需一步,快速开始

新浪微博登陆

只需一步, 快速开始

搜索
查看: 9849|回复: 8

MeteoInfoLab脚本示例:散点密度图

[复制链接]

新浪微博达人勋

发表于 2017-3-6 10:26:27 | 显示全部楼层 |阅读模式

登录后查看更多精彩内容~

您需要 登录 才可以下载或查看,没有帐号?立即注册 新浪微博登陆

x
本帖最后由 MeteoInfo 于 2017-3-6 11:34 编辑

当绘制大量散点的时候会有很多点压盖在一起,这时候可以通过散点密度色阶来区分散点分布的密度。散点密度是指点在某个位置出现的频次,某个位置通常是一个小区域,因此可以将散点分布的整个大区域分成小的矩形区域,然后计算每个矩形区域中点的个数。MeteoInfoLab中的griddata函数的功能是将散点格点化,其中method参数设为'inside_count'即可计算每个格点的散点数目,其返回值前三个参数是格点数据,第4个参数是散点数据的密度值(散点所在网格的密度即该散点的密度)。argsort函数获得数组排序的索引数组,可以用来对散点数据数组进行排序,以保证密度大的点后绘制(不会被密度小的点压盖)。

  1. #Create normal distribution point data
  2. x = random.randn(10000)
  3. y = x * 3 + random.randn(10000)

  4. #Calculate point density
  5. x_g = linspace(x.min(), x.max(), 100)
  6. y_g = linspace(y.min(), y.max(), 100)
  7. z1,x1,y1,z = griddata((x,y), None, xi=(x_g, y_g), method='inside_count')

  8. #Sort the points by density, so that the densest points are plotted last
  9. idx = argsort(z)
  10. x = x[idx]
  11. y = y[idx]
  12. z = z[idx]

  13. #Plot scatter point colored by density
  14. ls = scatter(x, y, s=4, c=z, edgecolor=None, cnum=20)
  15. colorbar(ls)
  16. title('Scatter density plot example')


scatter_density.png

评分

参与人数 1金钱 +20 贡献 +2 收起 理由
methouse + 20 + 2 很给力!

查看全部评分

密码修改失败请联系微信:mofangbao

新浪微博达人勋

发表于 2017-3-24 09:59:34 | 显示全部楼层
正在学习,谢谢楼主分享。楼主有处理过闪电数据吗,有脚本可以学习下的?
密码修改失败请联系微信:mofangbao

新浪微博达人勋

 楼主| 发表于 2017-3-24 10:34:02 | 显示全部楼层
qdsonglinyy 发表于 2017-3-24 09:59
正在学习,谢谢楼主分享。楼主有处理过闪电数据吗,有脚本可以学习下的?

参考这个脚本汇总贴:http://bbs.06climate.com/forum.p ... &extra=page%3D1
密码修改失败请联系微信:mofangbao

新浪微博达人勋

发表于 2017-3-28 11:33:36 | 显示全部楼层
多谢分享!
密码修改失败请联系微信:mofangbao
回复

使用道具 举报

新浪微博达人勋

发表于 2017-12-10 11:00:31 | 显示全部楼层
学习了,谢谢
密码修改失败请联系微信:mofangbao

新浪微博达人勋

发表于 2018-8-16 21:25:58 | 显示全部楼层
哇哦,楼主写得太好了{:5_213:}{:5_213:}{:5_213:}{:5_213:}
密码修改失败请联系微信:mofangbao

新浪微博达人勋

发表于 2018-12-20 15:09:07 | 显示全部楼层
看来很不错的样子
密码修改失败请联系微信:mofangbao

新浪微博达人勋

发表于 2020-8-5 16:48:37 | 显示全部楼层
感谢分享
{:eb502:}{:eb502:}{:eb502:}{:eb502:}{:eb502:}{:eb502:}{:eb502:}{:eb502:}
密码修改失败请联系微信:mofangbao
回复 支持 反对

使用道具 举报

新浪微博达人勋

发表于 2021-3-16 15:38:27 | 显示全部楼层
{:5_213:}太强了,正准备用
密码修改失败请联系微信:mofangbao
回复 支持 反对

使用道具 举报

您需要登录后才可以回帖 登录 | 立即注册 新浪微博登陆

本版积分规则

Copyright ©2011-2014 bbs.06climate.com All Rights Reserved.  Powered by Discuz! (京ICP-10201084)

本站信息均由会员发表,不代表气象家园立场,禁止在本站发表与国家法律相抵触言论

快速回复 返回顶部 返回列表