用Python绘制一个仿黑洞图像
简介
黑洞图像大家都知道，毕竟前几年刚发布的时候曾火遍全网，甚至都做成表情包了。
问题在于，凭什么认为这就是黑洞的照片，而不是一个甜甜圈啥的给整模糊了得到的呢？有什么理论依据吗？
单位制
利用einsteinpy模块中的Shadow类，可以实现黑洞吸积盘的发射强度，换言之，用理论模拟一下黑洞的照片应该是什么样子的。
22年5月份发布的是人马座A*的照片，这个"黑洞"的质量为8.26×1036kg，距离地球26673光年。
根据源码推测，einstenipy中采用几何单位制，光速和万有引力常数设为1， 即c=G=1。现设几何单位制下的时间、距离以及质量单位为S,M,KG，国际制单位为s,m,kg，若以时间单位秒为基本单位，即1s=1S，由于几何单位制下光速Sc=1M/S，则可得到距离单位的换算关系
同理，根据万有引力常数可得到质量单位的换算关系
所以质量和距离可以分别表示为
观测绘图
尽管我们在地球上观测这个黑洞，但由于相机镜头等因素，实际上肯定不能用这么大的尺度来观测，所以稍微选一个差不多的距离就好
import astropy.units as u from einsteinpy.rays import Shadow mass = 1.98 * u.kg fov = 100 * u.km shadow = Shadow(mass=mass, fov=fov, n_rays=1000)
shadow就是其发射模型，其内部用于绘图的主要成员有
- fb1 x xx负半轴
- fb2 x xx正半轴
- intensity 发射强度
下面绘制一下不同距离处的发射强度
import matplotlib.pyplot as plt
plt.plot(shadow.fb1, shadow.intensity, "r")
plt.plot(shadow.fb2, shadow.intensity, "r")
plt.xlabel("Impact Paramter (b)")
plt.ylabel("Intensity (Emissivity)")
plt.title("Intensity Plot")
plt.grid()
plt.show()
结果如下图所示，果然中间凹了下去
einsteinpy内嵌了绘图函数，通过ShadowPlotter(shadow, is_line_plot=True)可以更加方便地绘制上述图像，如果将is_line_plot设为False，则可得到强度分布图
from einsteinpy.plotting import ShadowPlotter obj = ShadowPlotter(shadow=shadow, is_line_plot=False) obj.plot() obj.show() plt.show()
这样一看是不是觉得和前面的那个甜甜圈确有相似之处。
这个模型是根据Cosimo Bambi的文章编写的，论文地址在这：Cosimo Bambi, 10.1103/PhysRevD.87.107501。本来合计着把这篇文章的公式啥的粘过来解读一下，但考虑到贴在这也没人看，想想还是算了。
到此这篇关于用Python绘制一个仿黑洞图像的文章就介绍到这了,更多相关Python仿黑洞内容请搜索脚本之家以前的文章或继续浏览下面的相关文章希望大家以后多多支持脚本之家！
相关文章
本篇文章主要介绍了详解python OpenCV学习笔记之直方图均衡化，小编觉得挺不错的，现在分享给大家，也给大家做个参考。一起跟随小编过来看看吧2018-02-02
在本篇文章里小编给大家整理的是关于Python实现的简易FTP的实例内容，有兴趣的朋友们可以参考下。2020-05-05
这篇文章主要介绍了pandas计数 value_counts()的使用，文中通过示例代码介绍的非常详细，对大家的学习或者工作具有一定的参考学习价值，需要的朋友们下面随着小编来一起学习学习吧2019-06-06
编写Python脚本批量下载DesktopNexus壁纸的教程
这篇文章主要介绍了编写Python脚本批量下载DesktopNexus壁纸的教程,相较于普通的爬虫抓取,本文的下载壁纸教程还包括了设置所要下载的分辨率等功能的实现,需要的朋友可以参考下2015-05-05
在Python中有两种函数，一种是def定义的函数，另一种是lambda函数，也就是大家常说的匿名函数。这篇文章主要介绍了Python的Lambda函数用法,需要的朋友可以参考下2019-09-09
这篇文章主要为大家介绍了Python pydotplus安装及可视化图形创建教程示例详解,有需要的朋友可以借鉴参考下,希望能够有所帮助,祝大家多多进步,早日升职加薪2023-10-10
这篇文章主要介绍了Python NumPy矩阵对象详解及方法，文章围绕主题展开详细的内容戒杀，具有一定的参考价值，需要的朋友可以参考一下2022-09-09
日常工作、生活中我们经常会遇到一些复杂的事务关系,比如人物关系,那如何才能清楚直观的看清楚这些任务关系呢？所以小编给大家介绍了Python如何使用Networkx实现复杂的人物关系图,文中通过代码示例讲解的非常详细,需要的朋友可以参考下2023-11-11
最近在做一个网站的后端开发。因为初期只有我一个人做，所以技术选择上很自由。在 web 服务器上我选择了 Tornado。虽然曾经也读过它的源码，并做过一些小的 demo，但毕竟这是第一次在工作中使用，难免又发现了一些值得分享的东西2014-05-05
JSON（JavaScript Object Nonation）是利用键值对+嵌套来表示数据的一种格式，以其轻量、易解析的优点，这篇文章主要介绍了利用Python实现Shp格式向GeoJSON的转换,需要的朋友可以参考下2019-07-07
最新评论