HAPI利用方法总结（HITRAN数据库，HITEMP数据库）

流逝的年华

之前我计较下温气鼓鼓体光谱特征是本人写的算法部分，查了很多多少文献，最初计较的成果也借算比力合意。厥后我才发明本来HITRAN民网上曾经供给了HAPI，只需求挪用HAPI就能够间接计较气鼓鼓体的光谱特征，而无需本人写具体的光谱特征计较办法。
HAPI (HITRAN Application Programming Interface) 是HITRAN online网站（https://hitran.org/）供给的python接心，能够完成长途会见HITRAN数据库，并能够逐止下载、过滤战处置份子战本子上的光谱数据的功用。

利用办法：

1. 间接从民网高低载hapi.py

2. 导进hapi模块

from hapi import *

复造代码

3. 创立文件夹，并指定该文件夹

首先，创立文件夹用于当地数据库的寄存。比方，文件夹名为data_hitemp。

然后利用以下语句选中该文件夹。

db_begin('data_hitemp')

复造代码

每次启动Python项目时，皆必需隐式指定命据库文件夹的称号。
4. 数据的获得

4.1 HITRAN数据库的获得：长途毗邻HITRAN网站下载
经由过程fetch()能够从HITRAN online主动下载HITRAN数据库的谱线数据。
fetch(TableName, M, I, numin, numax, ParameterGroups=[], Parameters=[])，此中TableName表示当地寄存的文件名；M表示HITRAN上份子的编号，如H2O的编号是1，CO2的编号是2，N2O的编号是4，等等；I表示HITRAN上同位素的编号；numin表示波数下限；numax表示波数上限；比方：

fetch('COC', 2, 1, 2000, 2100)

复造代码

上里代码表示将波数2000-2100（cm-1)的CO2的光谱数据寄存正在名为“COC”的文件中。此时能够正在data_hitemp文件夹中看到一个COC.data战一个COC.header文件。

4.2 HITEMP数据库的获得：间接下载HITEMP数据库
另外一种办法是间接下载好HITEMP数据库中的光谱数据，然后将其脚动寄存到data_hitemp文件夹中。从HITEMP数据库下载的数据格局是.par。

正在利用以下语句以后

db_begin('data_hitemp')

复造代码

能够看到，data_hitemp文件夹中的数据格局呈现了.header

5. 计较光谱特征

按照气鼓鼓体的温度，压力，光路少度，能够计较获得以下光谱特征：a)吸取系数（Absorption coefficient）；b)光谱吸取率（Absorption spectrum）；c)光谱透过率（Transmittance spectrum）；d)光谱辐射明度（Radiance spectrum）。HAPI供给的线形(line shape)包罗：Gaussian (Doppler)线形，Lorentzian线形，Voigt线形，Rautian线形，Speed-dependent Voigt线形，speed-dependent Rautian线形和Hartmann-Tran线形。
5.1 吸取系数
计较吸取系数可使用以下函数，别离代表差别的线形：

absorptionCoefficient_HT
absorptionCoefficient_Voigt
absorptionCoefficient_Lorentz
absorptionCoefficient_Doppler
absorptionCoefficient_SDVoigt

复造代码

上面以洛伦兹线形为例计较吸取系数：

nu,coef = absorptionCoefficient_Lorentz(SourceTables='02_2000-2125_HITEMP2010',HITRAN_units=False，Diluent={'air':1.0})
plot(nu,coef)

复造代码

此中，nu代表中间波数；coef代表吸取系数；Diluent表示气鼓鼓体混淆物的比例，Diluent={'self':A, 'air':B, 'CO2':C, ...}，且合意A+B+C+...=1。此中self用于计较自删宽半宽，air用于计较氛围删宽半宽。
数据库里只要氛围删宽半宽战自删宽半宽，所以其他组分除浓度中物资品种出有影响。
#若CO2露量为100%，则设置'self':1,另外一个参数设置air=0
#若CO2露量为1%，则设置'self'=0.01,另外一个设置为'air'=0.99
此时需求增长上面的代码

coef *= xco2 #xco2代表CO2的露量

复造代码

#若为杂氛围情况，设置'air'=1，便可
HITRAN_units: 假设为True，单位 cm^2/molec ；假设为False，单位cm-1 。为了后绝计较其他光谱特征如光谱吸取系数，光谱透过率，光谱辐射明度，必需将 HITRAN_units设为 False。
下左图是HITRAN_units=True，下左图是HITRAN_units=False.

完好的absorptionCoefficient_Lorentz参数以下：
absorptionCoefficient_Lorentz(SourceTables, Environment={'T':296.,'p':1.}, OmegaRange, OmegaStep, OmegaWing, IntensityThreshold, OmegaWingHW, GammaL='gamma_air', HITRAN_units, LineShift, File, Format, OmegaGrid, WavenumberRange, WavenumberStep, WavenumberWing, WavenumberWingHW, WavenumberGrid, Diluent={}, EnvDependences)
Environment={'T':296,'p':1}，默许温度T=296K，p=1atm，可按照实践温度战压力停止设置；
5.2 光谱吸取率，光谱透过率，光谱辐射明度
光谱吸取率，光谱透过率，光谱辐射明度的计较必需设置传输途径少度，默许的少度是1m。光谱吸取率战光谱透过率的默许参数是 Environment={'l': 100.0}，途径少度单位是cm。
光谱辐射明度的默许参数是Environment={'l': 100.0, 'T': 296.0}，除途径少度中，另有温度，那里的温度设置必需取之前正在absorptionCoefficient设置的温度连结分歧。

nu,absorp = absorptionSpectrum(nu,coef)
nu,trans = transmittanceSpectrum(nu,coef)
nu,radi = radianceSpectrum(nu,coef)

复造代码

需求留意的是，HAPI用户脚册中呈现了笔误，默许参数是Environment={'l': 100.0, 'T': 296.0}的该当为光谱辐射明度而非光谱透过率。

计较成果以下：
光谱吸取率：

光谱透过率：

光谱辐射明度：单位是W/cm^2/sr/cm^(-1)

期望能帮到需求的人，因为我打仗HAPI的工夫也很短，许多处所的了解能够没有到位，期望能战各人交流，增进我们配合前进。
更多闭于HAPI的利用办法能够下载用户利用脚册。
参考文章：
圆巍：HAPI——HITRAN数据库python接心 - 知乎 (zhihu.com)

免责声明：假如进犯了您的权益，请联络站少，我们会实时删除侵权内乱容，感谢协作！

HAPI利用方法总结（HITRAN数据库，HITEMP数据库）

浏览过的版块