从星系光谱计算积分线通量和等效宽度的软件包

pyspeclines的Python项目详细描述


#pyspeclines

一些有用的脚本,用于使用3种不同的方法从星系光谱计算积分线通量和等效宽度:

-数值(梯形)积分
-使用"标准"Levenberg-Marquardt算法进行高斯拟合
-使用MCMC算法进行高斯拟合安装包

yspeclines——升级
````

每个光谱都应该包含在一个单独的fits文件中,该文件**必须**包含以下列:

-``wl`、``ang``
-``flux``中的波长数组、``erg s^-1cm^-2ang^-1``
-``err``、``erg s^-1cm^-2ang^-1`

fits header

`` pyspeclines``还使用以下(**optional**)适合标题关键字:

-``redshift``,用于将光谱从观测帧红移到静止帧
-``ra`,对象在``deg``,在``icrs``帧中的赤经,用于校正银河吸收
-`` dec``,物体在``deg``中的偏角,在``icrs``框架中,用于校正银河吸收

r&davis(1998)由schlafly&finkbeiner(2011)重新校准,fitzpatrick(1999)的"r_v=3.1"消光曲线。

json配置文件

json文件允许您选择和配置要测量的发射线。在[pyspeclines/files]文件夹(https://github.com/jacopo-chevallard/pyspeclines/tree/master/pyspeclines/files)中提供了一些示例json文件。json文件包含一个"key:values"字典,其中"key"标记行(或一组行),而"values"包含多个条目。

686.0],
"wl-u-range":[4680.0,4681.0],
"continuum-left":[4672.0,4679.0],
"continuum-right":[4692.0,4700.0]
}
}
`````
-``wl-central``用作高斯拟合中直线中心的起点
-``wl-range``是r线的数值积分范围为左连续体的数值积分范围为右连续体的数值积分范围。"D是``[Continuum左[0],Continuum右[1]]````````


``行双:
```JSON
{
"SiI6716-SiI6731":{
"WL中心":[6716.0,6731.0],
"WL范围":[6710.0,6723.0],
"排除":[6710.0,6721.0,6721.0,6726.0,6740.0],
"Continuum连续统"。左":[6695.0,6705.0],
"continuum-right":[6740.0,6750.0]
}
}
````
-wrt在上面的例子中,``key``由两个标签组成,用下划线`````
-``exclude``允许定义区域(``[exclude[0],exclude[1]],[exclude[2],exclude[3]```)不包括在连续体配件中

-多个运动部件
``json
{
"oiii5007n-oiii5007b":{
"wl-central":[5007.0,5007.0],
"width":[100.0,400.0],
"wl-range":[5000.0,5014.0],
"exclude":[4995.0,5020.0],
"连续体左":[4990.0,5000.0],
"连续体右":[5020.0,5034.0]
}
}
`````
-``宽度``允许定义多个运动学组件,在本例中为"窄"(标记为"oiii5007n")和"宽"(标记为"oiii5007b")组件,其起始宽度**必须**设置为不同的值。

-具有多个运动学组件的多行/>"NII6548"是指"WLU中央":[6548.05,6563.0,6563.0,6563.0,6584.0],
"宽度":[100.0,100.0,400.0,100.0],
"排除":[6542.0,6580.0],
"WL范围":[6541.0,6575.0],
"左连续体":[6515.06542.0],
"C中心区":[6548.05,6563.0,6563.0,6563.0,6584.0],
"宽度":[100.0,100.0,100.0,400.0,100.0],
"排除":[6542.0,6580.0是的连续右":[6595.06610.0]
}
}
````
-在这种情况下,我们希望对不同的行(``nii6548``,``halphan``和``nii6584``)使用相同的宽度,对``halphab``使用不同的宽度。因此,我们对"nii6548"、"halphan"和"nii6584"使用相同的"width"值,因为这将在高斯拟合期间"绑定"它们的宽度,而"halphab"的"width"将保持独立。





nning
```
pyspeclines--help
``


如果光谱是在观察到的帧中提供的,则必须在包含对象红移的fits头中提供一个"redshift"关键字。

-使用数值积分计算通量和ewsle my_spectrum.fits--json文件发射线ews_config.json
```

-使用高斯拟合(levenberg marquardt)计算通量和ews
```
pyspeclines--file my_spectrum.fits--json文件发射线ews_config.json--gaussian拟合
````

-计算fl使用高斯拟合(mcmc)的uxes和ews
```
pyspeclines--file my戋spectrum.fits--json file emission戋lines戋ews戋config.json--gaussian fit--use pymc--mcmcmc samples 5000
```

欢迎加入QQ群-->: 979659372 Python中文网_新手群

推荐PyPI第三方库


热门话题
java加入条件查询出现错误“无法定位属性”   播放和停止音乐   java在两个活动之间传递int   java使Spring SAML SP与8080以外的其他端口一起工作   JavaSpring3:注入默认Bean,除非存在另一个Bean   java图像未在简单JFrame中显示   用户界面可以帮助我用这个类似的代码构造,但添加了需求(正确的日期格式,必须是18岁以上)Java GUI   带有WHERE、AND、OR子句的java Android Sqlite搜索查询   下载时出错'http://java.sun.com/xml/ns/persistence/persistence_2_1.xsd'到C://<path>   java使用迭代器对象循环throuh arraylist并检查它是否高于某个值   作为流的java祖先列表   java不能使用注释运行简单的servlet   Mule中的javaxmlxsd验证和路由   JAVA:当整数大于128时比较不起作用   Java+Web界面   java我如何无法在Gridlayout中再次设置可见?   java是初始化对象数组的正确方法吗?   在Java中使用jsonql数据源的json Jasper填充报告