射电天文偏振线的绘制

#读取数据I是流量强度，PA是偏振角，PI是偏振强度
I_data =  pf.getdata('9000.I.fits')
PA_hdr = pf.getheader('9000.PA.fits')
PI_hdr = pf.getheader('9000.PI.fits')
I_hdr = pf.getheader('9000.I.fits')#print(PA_data.shape, PI_data.shape, I_data.shape)Nx = 4   #这两个值为取值，越小对偏振的取值次数就会越频繁，反之亦然
Ny = 4
scale = 0.0007  #需要对PI进行比例的缩放，具体缩放程度看PI和I的数值差距
sigma = 9.89*10**-5 #根据自己数据的大小设置sigma，并且设置不绘制小于3倍sigma的值偏振信号
M0, N0 =PA_data.shape
#print(M0,N0)PA_w = WCS(PA_hdr)
PI_w = WCS(PI_hdr)
I_w = WCS(I_hdr)for i in range(0,M0,Nx):for j in range(0,N0,Ny):#没有偏振角或者偏振小于sigmaif np.isnan(PA_data[i,j]) or PI_data[i,j]< sigma: continue    pa0 = np.radians(PA_data[i, j]) + np.pi/2.pi0 = np.sqrt(PI_data[i, j]* scale) / 2.l, b = PA_w.wcs_pix2world([[j+1,i+1,0]], 0)[0][0:2]y1 = b+pi0*np.cos(pa0)x1 = l+pi0*np.sin(pa0)/np.cos(np.radians(y1))y2 = b-pi0*np.cos(pa0)x2 = l-pi0*np.sin(pa0)/np.cos(np.radians(y2))px1, py1 = I_w.wcs_world2pix([[x1, y1, 0, 0]], 0)[0][0:2]px2, py2 = I_w.wcs_world2pix([[x2, y2, 0, 0]], 0)[0][0:2]plt.plot([px1,px2],[py1,py2],'green') '''
#画偏振的箭头指向，拿的别人的代码，还没有尝试过，先这里吧
ang0 = 252.
d_ang = 59.2208573466
ang = np.radians(ang0 + d_ang)
x1 = 68.77
y1 = 2.82
y2 = y1 + 0.5*np.cos(ang)
x2 = x1 + 0.5*np.sin(ang)/np.cos(np.radians(y2))
px1, py1 = I_data.convproj.topixel((x1,y1))
px2, py2 = I_data.convproj.topixel((x2,y2))
plt.arrow(px1, py1, px2-px1, py2-py1, width=3., fc='c')
'''#plt.show()

射电天文偏振线的绘制相关推荐

射电天文谱线接收机和终端系统
射电天文中用于接收来自射电天体的原子和分子谱线信号,并测定谱线的频率.轮廓.线宽和亮温度等参数的一种特殊的接收设备.谱线接收机不同于连续谱接收机,它具有很高的频率分辨率,可以测出射电源亮温度随频率变化 ...
java在圆上画刻度线_Java – 绘制标尺(带有90度角刻度线的线)
我建议你 >实现一个标尺绘图方法,从左到右绘制一个简单的水平标尺 >使用Math.atan2计算出所需的角度. >在调用ruler-drawing-method之前应用带有平移和旋转 ...
Canvas中任意波浪线的绘制
最近接触一个墨水屏的项目,有个需求就是会进行一些标记绘制,类似直线.自由曲线.波浪线等,所以就涉及到波浪线如何绘制. 本篇文章并不会告诉你在 Android 中如何利用 Canvas 去绘制线条. 水 ...
openpyxl 绘制饼形图_好享学丨快速上手Pythonmatplotlib 箱线图绘制，学术人必备
好享学是高下制图推出的关于数据可视化经验分享栏目,我们将定期与您分享各界优秀人士的制图经验,一同学习. 01. 引言箱线图(Boxplot) 是一种用作显示一组数据分散情况资料的统计图表,本期推文 ...
python画波浪线_PPT绘制波浪线的四种方法
在开始学习教程前,先欣赏两幅画.画面中的波浪线使用的恰到好处,给整个画面增添张力与活力.那么用PPT能绘制出柔美的波浪线吗?答案当然是肯定的!一共四种方法,本文图文详解使用PPT绘制波浪线的四种方法. ...
autocad.net 画多段线_AutoCAD2016绘制多段线图文教程
AutoCAD2016如何绘制多段线?AutoCAD作为一款功能强大且专业的计算机辅助设计软件,所拥有的用户非常之多,可应用的范围也很广泛,AutoCAD中多段线,又被称为多义线,是CAD绘图过程中常 ...
python画k线_python画k线,python绘制动态k线及均线
内容导航: Q1:k线图上怎么画线一般你看K线软件都有画图功能,以同花顺为例,在上方菜单栏找到画图,点击,会在菜单栏下方一行或者屏幕最右边一列出现相应的画图工具. Q2:求教使用python绘制K线 ...
r语言ggplot2 多线图绘制图例_R语言绘制箱线图示例
R语言绘制箱线图箱线图(又称盒须图)是常见的统计图形之一,通过绘制连续型变量的五数总括,即最小值.下四分位数(第25百分位数).中位数(第50百分位数).上四分位数(第75百分位数)以及最大值,描述 ...
python绘制箱线图-python 中箱线图绘制
一.箱线图图绘制参数详解 plt.boxplot(x, notch=None, sym=None, vert=None, whis=None, positions=None, widths=None, ...
Arcgis-基于tif栅格图和线矢量绘制道路形变剖面线
基于Arcgis,根据以下步骤进行绘制道路形变剖面线图. 提取线两个端点. 数据管理工具 → 要素 → 要素折点转点,输入线矢量,选择点类型(两个端点). 计算线的长度,确定分割的点数,添加字段计算即 ...

射电天文偏振线的绘制

射电天文偏振线的绘制相关推荐

最新文章

热门文章