python绘制彩色地震剖面断层解释_python绘制地震散点图
本项目是利用五年左右的世界地震数据,通过python的pandas库、matplotlib库、basemap库等进行数据可视化,绘制出地震散点图。主要代码如下所示
from __future__ import division
import pandas as pd
from pandas import Series,DataFrame
import numpy as np
from matplotlib.patches import Polygon
chi_provinces = ['北京','天津','上海','重庆',
'河北','山西','辽宁','吉林',
'黑龙江','江苏','浙江','安徽',
'福建','江西','山东','河南',
'湖北','湖南','广东','海南',
'四川','贵州','云南','陕西',
'甘肃','青海','台湾','内蒙古',
'广西','西藏','宁夏','新疆',
'香港','澳门'] #list of chinese provinces
def is_in_china(str):
if str[:2] in chi_provinces:
return True
else:
return False
def convert_data_2014(x):
try:
return float(x.strip())
except ValueError:
return x
except AttributeError:
return x
def format_lat_lon(x):
try:
return x/100
except(TypeError):
return np.nan
df = pd.read_excel(r'C:/Users/GGWS/Desktop/shuju/201601-12.xls')
df = df.append(pd.read_excel(r'C:/Users/GGWS/Desktop/shuju/201201-12.xls'),ignore_index = True)
df = df.append(pd.read_excel(r'C:/Users/GGWS/Desktop/shuju/shuju.xls'),ignore_index = True)
df = df.append(pd.read_excel(r'C:/Users/GGWS/Desktop/shuju/201501-12.xls'),ignore_index = True)
df_2014 = pd.read_excel(r'C:/Users/GGWS/Desktop/shuju/201401-12.xls') #have to introduce statics of 2014 independently because the format and the type of data of specific column in this data set are different from others
df['longitude'] = df['longitude'].apply(convert_data_2014)
df['latitude'] = df['latitude'].apply(convert_data_2014)
df_2014['longitude'] = df_2014['longitude'].apply(convert_data_2014)
df_2014['latitude'] = df_2014['latitude'].apply(convert_data_2014)
df = df.append(df_2014,ignore_index = True)
df = df[['latitude','longitude','magnitude','referenced place','time']] #only save four columns as valuable statics
df[['longitude','latitude']] = df[['longitude','latitude']].applymap(format_lat_lon) #use function "applymap" to convert the format of the longitude and latitude statics
df = df.dropna(axis=0,how='any') #drop all rows that have any NaN values
format_magnitude = lambda x: float(str(x).strip('ML'))
df['magnitude'] = df['magnitude'].apply(format_magnitude)
#df = df[df['referenced place'].apply(is_in_china)]
lon_mean = (df['longitude'].groupby(df['referenced place'])).mean()
lat_mean = (df['latitude'].groupby(df['referenced place'])).mean()
group_counts = (df['magnitude'].groupby(df['referenced place'])).count()
after_agg_data = pd.concat([lon_mean,lat_mean,group_counts], axis = 1 )
after_agg_data.rename(columns = {'magnitude':'counts'} , inplace = True)
#aggregate after grouping the data
after_sorted_data = after_agg_data.sort_values(by = 'counts',ascending = False)
new_index = np.arange(len(after_sorted_data.index))
after_sorted_data.index = new_index
paint_data = after_sorted_data[after_sorted_data['counts']>=after_sorted_data['counts'][80]]
import matplotlib.pyplot as plt
from mpl_toolkits.basemap import Basemap
plt.figure(figsize=(16,8))
m = Basemap()
m.readshapefile(r'C:/Users/GGWS/Desktop/jb/gadm36_CHN_1', 'states', drawbounds=True)
ax = plt.gca()
'''
for nshape,seg in enumerate (m.states):
poly = Polygon(seg,facecolor = 'r')
ax.add_patch(poly)
'''
m.drawcoastlines(linewidth=0.5)
m.drawcountries(linewidth=0.5)
m.shadedrelief()
for indexs in df.index:
lon2,lat2 = df.loc[indexs].values[1],df.loc[indexs].values[0]
x,y = m(lon2,lat2)
m.plot(x,y,'ro',markersize = 0.5) #获取经度值
'''
for indexs in after_sorted_data.index[:80]:
lon,lat = after_sorted_data.loc[indexs].values[0],after_sorted_data.loc[indexs].values[1]
x,y = m(lon,lat)
m.plot(x,y,'wo',markersize = 10*(after_sorted_data.loc[indexs].values[2]/after_sorted_data.loc[0].values[2]))
'''
plt.title("Worldwide Earthquake")
plt.show()
#indexs-len(df.index)+80
效果如下
以上就是本文的全部内容,希望对大家的学习有所帮助,也希望大家多多支持我们。
本文标题: python绘制地震散点图
本文地址: http://www.cppcns.com/jiaoben/python/262492.html
python绘制彩色地震剖面断层解释_python绘制地震散点图相关推荐
- python绘制彩色地震剖面断层解释_绘制地震地质解释剖面.PPT
绘制地震地质解释剖面 * 第七章 地震资料解释 浅震勘探的最终目的:解决具体的地质问题. 资料采集是基础,处理工作是关键,资料解释是地震工作的中枢. 第一节 反射波资料的解释 由于反射法数据处理最终得 ...
- python绘制彩色地震剖面_用matplotlib绘制地震摆动轨迹
我正在尝试使用matplotlib重新创建上面的打印样式. 原始数据存储在2D numpy数组中,其中快轴是时间. 画线很容易.我在努力让阴影区域更有效. 我目前的尝试是这样的:import nump ...
- python画五角星为什么144度_python绘制五角星
分享一个使用Python绘制的五角星小案例,欢迎留言探讨指正,大佬勿喷,抱拳 下面呈上小代码 代码: import turtle #导入turtle库包 turtle.fillcolor(" ...
- 用python画一个简单卡通人物_Python绘制可爱的卡通人物 | 【turtle使用】-Go语言中文社区...
微信公众号:AI算法与图像处理如果你觉得对你有帮助,欢迎关注.转发以及点赞哦-( ̄▽ ̄-)~ Turtle库 简介 什么是Turtle 首先,turtle库是一个点线面的简单图像库,能够完成一些比较简 ...
- python绘制彩色地震剖面_在地图上绘制饼图时“爆炸”楔形图(Python,matplotlib)...
所以我成功地在地图上绘制了饼图作为标记轴向散射,但我遇到了一些问题,一些楔子"爆炸"出了饼图.我似乎在我的代码中找不到原因,也无法在网上找到解释.这段代码基于示例here,一位同事 ...
- python绘制条形图用什么函数_Python绘制正余弦函数图像完整代码
通过python绘制正弦和余弦函数,从默认的设置开始,一步一步地调整改进,让它变得好看,变成我们初高中学习过的图象那样.通过这个过程来学习如何进行对图表的一些元素的进行调整. 01. 简单绘图 mat ...
- python slice函数画高维图_Python 绘制 3 维以上的高维图
实用技巧. 我们的大脑通常最多能感知三维空间,超过三维就很难想象了.尽管是三维,理解起来也很费劲,所以大多数情况下都使用二维平面.不过,我们仍然可以绘制出多维空间,今天就来用 Python 的 plo ...
- 用python画数学函数图像教程_Python 绘制你想要的数学函数图形
Python 非常热门,但除非工作需要没有刻意去了解更多,直到有个函数图要绘制,想起了它.结果发现,完全用不着明白什么是编程,就可以使用它完成很多数学函数图的绘制. 通过以下两个步骤,就可以进行数学函 ...
- python用cartopy包画地图_python绘制地图的利器Cartopy使用说明
python绘制地图一般使用Basemap绘图包,但该包配置相对较繁琐,自定义性不强,这里介绍一个绘制地图的利器Cartopy,个人认为该工具方便.快捷,附上一些自己写的程序. 准备工作,工欲善其事, ...
最新文章
- Aras学习笔记(1)学习Aras已半年有余,也积攒一些学习笔记,今天起会陆续分享出来,有兴趣的朋友一起交流...
- python_web框架
- HDU 4279 占坑先。。。
- 如何在程序中生成崩溃转储dump文件以及如何分析dump
- 模式识别机器学习术语
- 人不能两次踏进同一条河流
- pytorch教程龙曲良01-05
- 编写一个汇编语言程序,将字符串Hello World中的全部小写字母转换为大写字母,并存放回原地址处。
- php传中文给Java_完美解决PHP中文乱码(转) - - JavaEye技术网站
- 繁忙SQL Server上的查询计划
- Eclipse快捷键之搜索
- CC26xx(CortexM3) UARTs
- 一文看懂量子十问(上篇)
- DNS 工作原理是什么,域名劫持、域名欺骗、域名污染又是什么
- ACM/ICPC 大赛常见英语词汇
- 暴雪这次可真不要脸。。。
- 一个关于家具的横幅广告
- rust服务器假人文件,最新武林最新GM命令大全(精睿提取) - Welcome to XiongTianQi.CN...
- 2021-06-18杭电ACM-LCY算法进阶培训班-专题训练16
- antdvue upload组件的customRequest自定义上传事件一直uploading处理方法
热门文章
- 根据ip反查公司信息和域名方法(python实现)
- QMediaPlayer视频播放器
- PHP 数据库读写分离 2020
- 淘宝店铺的运营思路和宝贝描述的要素与量子分析,开店必须要学的
- 转:WINFORM加Pannel背景图闪烁终极解决办法
- 在哪里尝试爆火的AI绘画?16款AI绘画工具大盘点!
- 人的一生就是在不断的学习中成长
- 国内无线通信模组企业(TOP 80)
- 产业互联网与消费互联网
- PB修改源码时报错:C0176: Badly ordered TYPE and VARIABLE declarations. Is this modified exported source?