Python绘制地磁场
文章目录
- 简介
- 磁场绘制
简介
pyIGRF
为国际参考磁场对Python的封装,可通过经纬高度以及时间来计算地磁场强度,使用方法简单粗暴,如下
import pyIGRF
pyIGRF.igrf_value(lat, lon, alt, date)
参数含义为
lat
纬度lon
经度alt
海拔date
日期,输入四位十进制整数表示年份
其返回值有七个,若写成
D, I, H, X, Y, Z, F = pyIGRF.igrf_value(lat, lon, alt, date)
各返回值含义为
D
偏角,以东向为正I
倾角,以竖直向下为正H
水平强度X
北向分量Y
东向分量Z
垂直分量,以下方为正F
总磁场强度
角度采取角度制, 磁场强度单位为nT
。
通过函数pyIGRF.igrf_variation
可以获取磁场的变化情况,其输入与igrf_value
相同,输出量的物理量也相同,但代表的是年变化量。
磁场绘制
接下来绘制一下某片区域的磁场分布,令经度范围在126附近,纬度范围在46附近,海拔由0到3000米。
地球建模采取WGS84,其长半轴为6378137km,扁率 e = 1 298.257223563 e=\frac{1}{298.257223563} e=298.2572235631,则短半轴为b = 6356752.314m。
import numpy as np
import matplotlib.pyplot as plt
xs, ys, zs = np.indices([5,5,5])
xs = (xs/500 + 125).flatten()
ys = (ys/500 + 45).flatten()
zs = (zs*600).flatten()
vals = [pyIGRF.igrf_value(y, x, z, 2022) for x,y,z in zip(xs,ys,zs)]
vals = np.array(vals)
Hx, Hy, Hz, H = vals[:,3:].TLAT = 6356752.314*np.pi/180
xs = (xs-np.min(xs))*np.sin(np.deg2rad(ys))*LAT
ys = (ys-np.min(ys))*LATax = plt.subplot(projection='3d')
ax.quiver(xs, ys, zs, Hx/300, Hy/300, Hz/300)
plt.show()
最终效果如图所示,可见随着高度的增加,地磁场是逐渐减弱的。
接下来可以查看某地随着高度变化,其地磁场分量的变化
hs = np.arange(10000)
vals = [pyIGRF.igrf_value(46, 125, h, 2022) for h in hs]
vals = np.array(vals).T
labels = ["North", "East", "up"]
for i in range(3):plt.plot(hs, vals[i+3], label=labels[i])plt.legend()
plt.show()
结果为
当然,这一局部地磁场分布其实看不出什么,接下来绘制一下全球的磁场分布情况
Ls, Bs = np.indices([180,90])
Ls = (Ls*2 - 180).flatten()
Bs = (Bs*2 - 90).flatten()
vals = [pyIGRF.igrf_value(y, x, 100, 2022) for x,y in zip(Ls,Bs)]
vals = np.array(vals).T
H = vals[-1]r = 6356752.314/2
Bs, Ls = np.deg2rad(Bs), np.deg2rad(Ls)
xs = r*np.cos(Bs)*np.cos(Ls)
ys = r*np.cos(Bs)*np.sin(Ls)
zs = r*np.sin(Bs)ax = plt.subplot(projection='3d')
cax = ax.scatter(xs, ys, zs, c=H, cmap='plasma')
plt.show()
效果如图所示
Python绘制地磁场相关推荐
- 怎么用python找论文_如何利用Python绘制学术论文图表方法
论文中图表形式多样,常用的处理工具有excel.MATLAB以及Python等,excel自处理的方法有两个缺陷: 1.当数据较多时,容易出现excel"翻白眼"的现象: 2.需要 ...
- python绘制条形图用什么函数_Python绘制正余弦函数图像完整代码
通过python绘制正弦和余弦函数,从默认的设置开始,一步一步地调整改进,让它变得好看,变成我们初高中学习过的图象那样.通过这个过程来学习如何进行对图表的一些元素的进行调整. 01. 简单绘图 mat ...
- python绘制所有特征的密度图(density plot)
python绘制所有特征的密度图(density plot) # python绘制所有特征的密度图: n = len(input_feats) #密度图kde,查查分布 cols = 4 rows = ...
- python绘制影像组学训练集、测试集对应的ROC曲线以及瀑布图(rad-score 瀑布图)
python绘制影像组学训练集.测试集对应的ROC曲线以及瀑布图(rad-score 瀑布图) # 所有数据的瀑布图 tagets_all = df['label'] tagets_all = df. ...
- 什么是折线图?怎样用Python绘制?怎么用?终于有人讲明白了(附代码)
来源:大数据DT 本文约5400字,建议阅读10分钟 本文为你介绍数据分析时经常用到的折线图,可以用来呈现哪些数据关系?在数据分析过程中可以解决哪些问题?怎样用Python绘制折线图? 01 概述 折 ...
- 柱状图、堆叠柱状图、瀑布图有什么区别?怎样用Python绘制?(附代码)
来源:大数据DT(ID:hzdashuju) 作者:屈希峰,资深Python工程师,知乎多个专栏作者 本文约8000字,建议阅读20分钟 柱状图是当前应用最广泛的图表之一,你几乎每天都可以在电子产品上 ...
- 层次聚类python_用python绘制层次聚类图
层次聚类(Hierarchical clustering)代表着一类的聚类算法,这种类别的算法通过不断的合并或者分割内置聚类来构建最终聚类. 聚类的层次可以被表示成树(或者树形图(dendrogram ...
- 利用Python绘制 3D 体素色温图
简 介: 测试了Matplotlib中的 体素绘制函数 voxels() ,但是无论是在 AI Studio中执行,还是直接在本地执行,都会出了错误.也就是无法通过 gca(projection='3 ...
- python绘制直方图显示数字_Python实现绘制双柱状图并显示数值功能示例
本文实例讲述了Python实现绘制双柱状图并显示数值功能.分享给大家供大家参考,具体如下: # -*- coding:utf-8 -*- #! python3 import matplotlib.py ...
最新文章
- 关于IE8.0 Internet Explorer 已经为了帮助保护您的计算机而关闭此网页 解决方法
- Lesson 7 (3) 深入理解PyTorch与PyTorch库架构
- 深度优先搜索(DFS)算法
- 1.4- 定时任务总结之九句箴言
- 当c语言老师遇到网络崩溃,网络教学搞不定?工科教授来支招!
- jmeter参数值只读取了第一个_基础入门篇Jmeter(五)_函数
- m1芯片Mac安装jekyll+搭建GitHub pages个人博客站点
- Theano 中文文档 0.9 - 7.2.1 起手式 —— 代数
- oracle设置no null,Oracle 在not in中使用null的问题
- ios html异步加载图片,iOS 异步加载本地图片
- 新型 JhoneRAT 恶意软件攻击中东地区
- shell脚本中执行命令_如何在Shell脚本中执行命令?
- Kattis - whatdoesthefoxsay —— 字符串
- 苹果Mac测试及维护工具:​​​​​​​​Techtool Pro
- tensorflow数据预处理
- VSCode下载安装和修改插件下载位置(配置右键菜单)
- android的opengl教程,android opengl 教程
- java毕业设计培训学校教学管理平台源码+lw文档+mybatis+系统+mysql数据库+调试
- 开放云将使业务焕然一新
- 手把手教你iPhone 3G手机软件开发 转帖
热门文章
- 数学建模的四大基本类型
- 代码实践:基于LSTM网络的DEAP情感数据集情感分类
- linux中重置是啥意思,linux – netstat -s中“由于意外数据而重置连接”是什么意思...
- ECshop导入淘宝数据包乱码问题解决方法
- 18 《社会性动物》 -豆瓣评分9.1
- softmax是什么?
- 记:《洛克菲勒留给儿子的38封信》-- 33
- Linux运维学习:中级进阶(2)——Linux企业常用服务(SSH、DHCP、FTP、NFS、DNS)
- Linux下 mdadm 模拟RAID 0, RAID1
- IT求职宝典-微软十佳金牌讲师孔文达作序推荐!