zgyin的个人博客分享 http://blog.sciencenet.cn/u/zgyin

博文

信噪比图当中寻找云核的原因

已有 2413 次阅读 2013-1-11 20:22 |系统分类:科研笔记

为什么要从S/N图中寻找云核,这种方法是可行的,理由如下:
1、Li, D., Velusamy, T., Goldsmith, P.F. & Langer, W., 2007, “Massive Quiescent Cores in Orion part II – Core Mass Function”, The Astrophysical Journal, Vol. 655, pp. 351.这篇文章就是从S/N图中寻找云核的
2、http://docs.jach.hawaii.edu/star/sun255.htx/node18.html 中讲到如果假设噪音恒定,从S/N中寻找核要比从flux图中寻找核要好,如果没有S/N图可以采用,网页中给的方法计算S/N图。

我今天呈现了数据的3D图(注意z轴的刻度

x轴是赤纬(dec),y轴是赤经(ra),z轴是流量(flux)

x轴是赤纬(dec),y轴是赤经(ra),z轴是信噪比(s/n)
对图进行分析,可以发现流量高的地方信噪比大,流量图类似与信噪比图,显而易见,通过信噪比图更容易找核
数据可视化的程序
pro plo
ih = readfits("/root/data/ROPH_SWAS.37708-37914.extended.fits",EXTEN_NO=0)
im = readfits("/root/data/ROPH_SWAS.37708-37914.extended.fits",EXTEN_NO=1)
in = readfits("/root/data/ROPH_SWAS.37708-37914.extended.fits",EXTEN_NO=3)
iflux=fltarr(n_elements(im[*,1]),n_elements(im[1,*]))
istn=fltarr(n_elements(im[*,1]),n_elements(im[1,*]))
for i=0,n_elements(im[*,1])-1 do begin
       for j=0,n_elements(im[1,*])-1 do begin
             if finite(im[i,j]) ne 0 then begin
                 if im[i,j] ge 2.0 then begin
                    iflux[i,j]= ih[i,j]
                    istn[i,j]=in[i,j]
                  endif
              endif
          endfor
endfor
set_plot,'ps'
device,file='/root/data/3dimage/3dflux.ps'
shade_surf,iflux
device,/close
set_plot,'ps'
device,file='/root/data/3dimage/3dstn.ps'
shade_surf,istn
device,/close
end
同时通过下面的程序计算了流量图有效值的平均值和信噪比图的有效值的平均值
pro aver
ih = readfits("/root/data/ROPH_SWAS.37708-37914.extended.fits",Header1,EXTEN_NO=0)
im = readfits("/root/data/ROPH_SWAS.37708-37914.extended.fits",Header2,EXTEN_NO=1)
in = readfits("/root/data/ROPH_SWAS.37708-37914.extended.fits",Header2,EXTEN_NO=3)
a=0.0d
b=0.0d
n=0l
for i=0,n_elements(im[*,1])-1 do begin
       for j=0,n_elements(im[1,*])-1 do begin
             if finite(im[i,j]) ne 0 then begin
                 if im[i,j] ge 2.0 then begin
                    a=a+ih[i,j]
                    b=b+in[i,j]
                    n++     
                  endif
              endif
          endfor
endfor
           
averflux=a/(n*1.0)
averstn=b/(n*1.0)
ratio=averstn/averflux

print,'average of flux =',averflux
print,'average of stn =',averstn
print,'averstn/averflux=',ratio
end
就算结果:
average of flux =     0.029920901
average of stn =      0.21340572
averstn/averflux=       7.1323292
从信噪比图中寻找到的云核



https://wap.sciencenet.cn/blog-783377-652132.html

上一篇:从信噪比图当中寻找云核
下一篇:fedora 键盘失效而鼠标没问题的解决方案
收藏 IP: 202.101.72.*| 热度|

0

该博文允许注册用户评论 请点击登录 评论 (0 个评论)

数据加载中...

Archiver|手机版|科学网 ( 京ICP备07017567号-12 )

GMT+8, 2023-2-6 01:01

Powered by ScienceNet.cn

Copyright © 2007- 中国科学报社

返回顶部