New fitting module introduction - simple 1-temperature LTE case¶

RADIS has its own fitting feature, as shown in 1 temperature fit example, where you have to manually create the spectrum model, input the experimental spectrum and other ground-truths into numerous RADIS native functions, and adjust the fitting pipeline yourself.

Now with the new fitting module released, all you have to do is to prepare a .spec file containing your experimental spectrum, fill some JSON forms describing the ground-truth conditions just like how you fill a medical checkup paper, call the function fit_spectrum() and let it do all the work! If you are not satisfied with the result, you can simply adjust the parameters in your JSON such as slit and path_length, recall the function until the results are good.

Instruction:

Step 1: prepare a .spec file. Create a .spec file containing your experimental spectrum. You can do it with RADIS by saving a Spectrum object with store(). If your current data is not a Spectrum object, you can convert it to a Spectrum object from Python arrays or from text files, and then save it as .spec file as above.
Step 2: fill the JSON forms. There are 4 JSON forms you need to fill: experimental_conditions with ground-truth data about your experimental environment, fit_parameters with the parameters you need to fit (such as Tgas, mole fraction, etc.), bounding_ranges with fitting ranges for parameters you listed in fit_parameters, and fit_properties for some fitting pipeline references.
Step 3: call fit_spectrum() with the experimental spectrum and 4 JSON forms, then see the result.

This example features fitting an experimental spectrum with Tgas, using new fitting modules.

======================= COMMENCE FITTING PROCESS =======================

Successfully retrieved the experimental data in 0.0006549358367919922s.

Acquired spectral quantity 'radiance' from the spectrum.
NaN values successfully purged.Number of data points left: 16666 points.
Successfully refined the experimental data in 0.0011031627655029297s.
"tol" parameter spotted but "method" is not "lbfgsb"!
Using /home/docs/.radisdb/hitran/downloads__can_be_deleted/NH3


Data is fetched from http://hitran.org

BEGIN DOWNLOAD: NH3_1
bytes written to /home/docs/.radisdb/hitran/downloads__can_be_deleted/NH3/NH3_1.data
bytes written to /home/docs/.radisdb/hitran/downloads__can_be_deleted/NH3/NH3_1.data
bytes written to /home/docs/.radisdb/hitran/downloads__can_be_deleted/NH3/NH3_1.data
bytes written to /home/docs/.radisdb/hitran/downloads__can_be_deleted/NH3/NH3_1.data
bytes written to /home/docs/.radisdb/hitran/downloads__can_be_deleted/NH3/NH3_1.data
bytes written to /home/docs/.radisdb/hitran/downloads__can_be_deleted/NH3/NH3_1.data
bytes written to /home/docs/.radisdb/hitran/downloads__can_be_deleted/NH3/NH3_1.data
bytes written to /home/docs/.radisdb/hitran/downloads__can_be_deleted/NH3/NH3_1.data
bytes written to /home/docs/.radisdb/hitran/downloads__can_be_deleted/NH3/NH3_1.data
bytes written to /home/docs/.radisdb/hitran/downloads__can_be_deleted/NH3/NH3_1.data
bytes written to /home/docs/.radisdb/hitran/downloads__can_be_deleted/NH3/NH3_1.data
bytes written to /home/docs/.radisdb/hitran/downloads__can_be_deleted/NH3/NH3_1.data
bytes written to /home/docs/.radisdb/hitran/downloads__can_be_deleted/NH3/NH3_1.data
bytes written to /home/docs/.radisdb/hitran/downloads__can_be_deleted/NH3/NH3_1.data
bytes written to /home/docs/.radisdb/hitran/downloads__can_be_deleted/NH3/NH3_1.data
bytes written to /home/docs/.radisdb/hitran/downloads__can_be_deleted/NH3/NH3_1.data
bytes written to /home/docs/.radisdb/hitran/downloads__can_be_deleted/NH3/NH3_1.data
bytes written to /home/docs/.radisdb/hitran/downloads__can_be_deleted/NH3/NH3_1.data
bytes written to /home/docs/.radisdb/hitran/downloads__can_be_deleted/NH3/NH3_1.data
bytes written to /home/docs/.radisdb/hitran/downloads__can_be_deleted/NH3/NH3_1.data
bytes written to /home/docs/.radisdb/hitran/downloads__can_be_deleted/NH3/NH3_1.data
bytes written to /home/docs/.radisdb/hitran/downloads__can_be_deleted/NH3/NH3_1.data
bytes written to /home/docs/.radisdb/hitran/downloads__can_be_deleted/NH3/NH3_1.data
bytes written to /home/docs/.radisdb/hitran/downloads__can_be_deleted/NH3/NH3_1.data
bytes written to /home/docs/.radisdb/hitran/downloads__can_be_deleted/NH3/NH3_1.data
bytes written to /home/docs/.radisdb/hitran/downloads__can_be_deleted/NH3/NH3_1.data
bytes written to /home/docs/.radisdb/hitran/downloads__can_be_deleted/NH3/NH3_1.data
bytes written to /home/docs/.radisdb/hitran/downloads__can_be_deleted/NH3/NH3_1.data
bytes written to /home/docs/.radisdb/hitran/downloads__can_be_deleted/NH3/NH3_1.data
bytes written to /home/docs/.radisdb/hitran/downloads__can_be_deleted/NH3/NH3_1.data
bytes written to /home/docs/.radisdb/hitran/downloads__can_be_deleted/NH3/NH3_1.data
bytes written to /home/docs/.radisdb/hitran/downloads__can_be_deleted/NH3/NH3_1.data
bytes written to /home/docs/.radisdb/hitran/downloads__can_be_deleted/NH3/NH3_1.data
bytes written to /home/docs/.radisdb/hitran/downloads__can_be_deleted/NH3/NH3_1.data
bytes written to /home/docs/.radisdb/hitran/downloads__can_be_deleted/NH3/NH3_1.data
bytes written to /home/docs/.radisdb/hitran/downloads__can_be_deleted/NH3/NH3_1.data
bytes written to /home/docs/.radisdb/hitran/downloads__can_be_deleted/NH3/NH3_1.data
bytes written to /home/docs/.radisdb/hitran/downloads__can_be_deleted/NH3/NH3_1.data
bytes written to /home/docs/.radisdb/hitran/downloads__can_be_deleted/NH3/NH3_1.data
bytes written to /home/docs/.radisdb/hitran/downloads__can_be_deleted/NH3/NH3_1.data
bytes written to /home/docs/.radisdb/hitran/downloads__can_be_deleted/NH3/NH3_1.data
bytes written to /home/docs/.radisdb/hitran/downloads__can_be_deleted/NH3/NH3_1.data
bytes written to /home/docs/.radisdb/hitran/downloads__can_be_deleted/NH3/NH3_1.data
bytes written to /home/docs/.radisdb/hitran/downloads__can_be_deleted/NH3/NH3_1.data
bytes written to /home/docs/.radisdb/hitran/downloads__can_be_deleted/NH3/NH3_1.data
bytes written to /home/docs/.radisdb/hitran/downloads__can_be_deleted/NH3/NH3_1.data
bytes written to /home/docs/.radisdb/hitran/downloads__can_be_deleted/NH3/NH3_1.data
bytes written to /home/docs/.radisdb/hitran/downloads__can_be_deleted/NH3/NH3_1.data
bytes written to /home/docs/.radisdb/hitran/downloads__can_be_deleted/NH3/NH3_1.data
bytes written to /home/docs/.radisdb/hitran/downloads__can_be_deleted/NH3/NH3_1.data
bytes written to /home/docs/.radisdb/hitran/downloads__can_be_deleted/NH3/NH3_1.data
bytes written to /home/docs/.radisdb/hitran/downloads__can_be_deleted/NH3/NH3_1.data
bytes written to /home/docs/.radisdb/hitran/downloads__can_be_deleted/NH3/NH3_1.data
bytes written to /home/docs/.radisdb/hitran/downloads__can_be_deleted/NH3/NH3_1.data
bytes written to /home/docs/.radisdb/hitran/downloads__can_be_deleted/NH3/NH3_1.data
bytes written to /home/docs/.radisdb/hitran/downloads__can_be_deleted/NH3/NH3_1.data
bytes written to /home/docs/.radisdb/hitran/downloads__can_be_deleted/NH3/NH3_1.data
bytes written to /home/docs/.radisdb/hitran/downloads__can_be_deleted/NH3/NH3_1.data
bytes written to /home/docs/.radisdb/hitran/downloads__can_be_deleted/NH3/NH3_1.data
bytes written to /home/docs/.radisdb/hitran/downloads__can_be_deleted/NH3/NH3_1.data
bytes written to /home/docs/.radisdb/hitran/downloads__can_be_deleted/NH3/NH3_1.data
bytes written to /home/docs/.radisdb/hitran/downloads__can_be_deleted/NH3/NH3_1.data
bytes written to /home/docs/.radisdb/hitran/downloads__can_be_deleted/NH3/NH3_1.data
bytes written to /home/docs/.radisdb/hitran/downloads__can_be_deleted/NH3/NH3_1.data
bytes written to /home/docs/.radisdb/hitran/downloads__can_be_deleted/NH3/NH3_1.data
bytes written to /home/docs/.radisdb/hitran/downloads__can_be_deleted/NH3/NH3_1.data
bytes written to /home/docs/.radisdb/hitran/downloads__can_be_deleted/NH3/NH3_1.data
bytes written to /home/docs/.radisdb/hitran/downloads__can_be_deleted/NH3/NH3_1.data
bytes written to /home/docs/.radisdb/hitran/downloads__can_be_deleted/NH3/NH3_1.data
bytes written to /home/docs/.radisdb/hitran/downloads__can_be_deleted/NH3/NH3_1.data
bytes written to /home/docs/.radisdb/hitran/downloads__can_be_deleted/NH3/NH3_1.data
bytes written to /home/docs/.radisdb/hitran/downloads__can_be_deleted/NH3/NH3_1.data
bytes written to /home/docs/.radisdb/hitran/downloads__can_be_deleted/NH3/NH3_1.data
bytes written to /home/docs/.radisdb/hitran/downloads__can_be_deleted/NH3/NH3_1.data
bytes written to /home/docs/.radisdb/hitran/downloads__can_be_deleted/NH3/NH3_1.data
bytes written to /home/docs/.radisdb/hitran/downloads__can_be_deleted/NH3/NH3_1.data
bytes written to /home/docs/.radisdb/hitran/downloads__can_be_deleted/NH3/NH3_1.data
bytes written to /home/docs/.radisdb/hitran/downloads__can_be_deleted/NH3/NH3_1.data
bytes written to /home/docs/.radisdb/hitran/downloads__can_be_deleted/NH3/NH3_1.data
bytes written to /home/docs/.radisdb/hitran/downloads__can_be_deleted/NH3/NH3_1.data
bytes written to /home/docs/.radisdb/hitran/downloads__can_be_deleted/NH3/NH3_1.data
bytes written to /home/docs/.radisdb/hitran/downloads__can_be_deleted/NH3/NH3_1.data
bytes written to /home/docs/.radisdb/hitran/downloads__can_be_deleted/NH3/NH3_1.data
bytes written to /home/docs/.radisdb/hitran/downloads__can_be_deleted/NH3/NH3_1.data
bytes written to /home/docs/.radisdb/hitran/downloads__can_be_deleted/NH3/NH3_1.data
bytes written to /home/docs/.radisdb/hitran/downloads__can_be_deleted/NH3/NH3_1.data
bytes written to /home/docs/.radisdb/hitran/downloads__can_be_deleted/NH3/NH3_1.data
bytes written to /home/docs/.radisdb/hitran/downloads__can_be_deleted/NH3/NH3_1.data
bytes written to /home/docs/.radisdb/hitran/downloads__can_be_deleted/NH3/NH3_1.data
bytes written to /home/docs/.radisdb/hitran/downloads__can_be_deleted/NH3/NH3_1.data
bytes written to /home/docs/.radisdb/hitran/downloads__can_be_deleted/NH3/NH3_1.data
bytes written to /home/docs/.radisdb/hitran/downloads__can_be_deleted/NH3/NH3_1.data
bytes written to /home/docs/.radisdb/hitran/downloads__can_be_deleted/NH3/NH3_1.data
bytes written to /home/docs/.radisdb/hitran/downloads__can_be_deleted/NH3/NH3_1.data
bytes written to /home/docs/.radisdb/hitran/downloads__can_be_deleted/NH3/NH3_1.data
bytes written to /home/docs/.radisdb/hitran/downloads__can_be_deleted/NH3/NH3_1.data
bytes written to /home/docs/.radisdb/hitran/downloads__can_be_deleted/NH3/NH3_1.data
bytes written to /home/docs/.radisdb/hitran/downloads__can_be_deleted/NH3/NH3_1.data
bytes written to /home/docs/.radisdb/hitran/downloads__can_be_deleted/NH3/NH3_1.data
bytes written to /home/docs/.radisdb/hitran/downloads__can_be_deleted/NH3/NH3_1.data
bytes written to /home/docs/.radisdb/hitran/downloads__can_be_deleted/NH3/NH3_1.data
bytes written to /home/docs/.radisdb/hitran/downloads__can_be_deleted/NH3/NH3_1.data
bytes written to /home/docs/.radisdb/hitran/downloads__can_be_deleted/NH3/NH3_1.data
bytes written to /home/docs/.radisdb/hitran/downloads__can_be_deleted/NH3/NH3_1.data
bytes written to /home/docs/.radisdb/hitran/downloads__can_be_deleted/NH3/NH3_1.data
bytes written to /home/docs/.radisdb/hitran/downloads__can_be_deleted/NH3/NH3_1.data
bytes written to /home/docs/.radisdb/hitran/downloads__can_be_deleted/NH3/NH3_1.data
bytes written to /home/docs/.radisdb/hitran/downloads__can_be_deleted/NH3/NH3_1.data
bytes written to /home/docs/.radisdb/hitran/downloads__can_be_deleted/NH3/NH3_1.data
bytes written to /home/docs/.radisdb/hitran/downloads__can_be_deleted/NH3/NH3_1.data
bytes written to /home/docs/.radisdb/hitran/downloads__can_be_deleted/NH3/NH3_1.data
bytes written to /home/docs/.radisdb/hitran/downloads__can_be_deleted/NH3/NH3_1.data
bytes written to /home/docs/.radisdb/hitran/downloads__can_be_deleted/NH3/NH3_1.data
bytes written to /home/docs/.radisdb/hitran/downloads__can_be_deleted/NH3/NH3_1.data
bytes written to /home/docs/.radisdb/hitran/downloads__can_be_deleted/NH3/NH3_1.data
bytes written to /home/docs/.radisdb/hitran/downloads__can_be_deleted/NH3/NH3_1.data
bytes written to /home/docs/.radisdb/hitran/downloads__can_be_deleted/NH3/NH3_1.data
bytes written to /home/docs/.radisdb/hitran/downloads__can_be_deleted/NH3/NH3_1.data
bytes written to /home/docs/.radisdb/hitran/downloads__can_be_deleted/NH3/NH3_1.data
bytes written to /home/docs/.radisdb/hitran/downloads__can_be_deleted/NH3/NH3_1.data
bytes written to /home/docs/.radisdb/hitran/downloads__can_be_deleted/NH3/NH3_1.data
bytes written to /home/docs/.radisdb/hitran/downloads__can_be_deleted/NH3/NH3_1.data
bytes written to /home/docs/.radisdb/hitran/downloads__can_be_deleted/NH3/NH3_1.data
bytes written to /home/docs/.radisdb/hitran/downloads__can_be_deleted/NH3/NH3_1.data
bytes written to /home/docs/.radisdb/hitran/downloads__can_be_deleted/NH3/NH3_1.data
bytes written to /home/docs/.radisdb/hitran/downloads__can_be_deleted/NH3/NH3_1.data
bytes written to /home/docs/.radisdb/hitran/downloads__can_be_deleted/NH3/NH3_1.data
bytes written to /home/docs/.radisdb/hitran/downloads__can_be_deleted/NH3/NH3_1.data
bytes written to /home/docs/.radisdb/hitran/downloads__can_be_deleted/NH3/NH3_1.data
bytes written to /home/docs/.radisdb/hitran/downloads__can_be_deleted/NH3/NH3_1.data
bytes written to /home/docs/.radisdb/hitran/downloads__can_be_deleted/NH3/NH3_1.data
bytes written to /home/docs/.radisdb/hitran/downloads__can_be_deleted/NH3/NH3_1.data
bytes written to /home/docs/.radisdb/hitran/downloads__can_be_deleted/NH3/NH3_1.data
bytes written to /home/docs/.radisdb/hitran/downloads__can_be_deleted/NH3/NH3_1.data
bytes written to /home/docs/.radisdb/hitran/downloads__can_be_deleted/NH3/NH3_1.data
bytes written to /home/docs/.radisdb/hitran/downloads__can_be_deleted/NH3/NH3_1.data
bytes written to /home/docs/.radisdb/hitran/downloads__can_be_deleted/NH3/NH3_1.data
bytes written to /home/docs/.radisdb/hitran/downloads__can_be_deleted/NH3/NH3_1.data
bytes written to /home/docs/.radisdb/hitran/downloads__can_be_deleted/NH3/NH3_1.data
bytes written to /home/docs/.radisdb/hitran/downloads__can_be_deleted/NH3/NH3_1.data
bytes written to /home/docs/.radisdb/hitran/downloads__can_be_deleted/NH3/NH3_1.data
bytes written to /home/docs/.radisdb/hitran/downloads__can_be_deleted/NH3/NH3_1.data
bytes written to /home/docs/.radisdb/hitran/downloads__can_be_deleted/NH3/NH3_1.data
bytes written to /home/docs/.radisdb/hitran/downloads__can_be_deleted/NH3/NH3_1.data
bytes written to /home/docs/.radisdb/hitran/downloads__can_be_deleted/NH3/NH3_1.data
bytes written to /home/docs/.radisdb/hitran/downloads__can_be_deleted/NH3/NH3_1.data
bytes written to /home/docs/.radisdb/hitran/downloads__can_be_deleted/NH3/NH3_1.data
bytes written to /home/docs/.radisdb/hitran/downloads__can_be_deleted/NH3/NH3_1.data
bytes written to /home/docs/.radisdb/hitran/downloads__can_be_deleted/NH3/NH3_1.data
bytes written to /home/docs/.radisdb/hitran/downloads__can_be_deleted/NH3/NH3_1.data
bytes written to /home/docs/.radisdb/hitran/downloads__can_be_deleted/NH3/NH3_1.data
bytes written to /home/docs/.radisdb/hitran/downloads__can_be_deleted/NH3/NH3_1.data
bytes written to /home/docs/.radisdb/hitran/downloads__can_be_deleted/NH3/NH3_1.data
bytes written to /home/docs/.radisdb/hitran/downloads__can_be_deleted/NH3/NH3_1.data
bytes written to /home/docs/.radisdb/hitran/downloads__can_be_deleted/NH3/NH3_1.data
bytes written to /home/docs/.radisdb/hitran/downloads__can_be_deleted/NH3/NH3_1.data
bytes written to /home/docs/.radisdb/hitran/downloads__can_be_deleted/NH3/NH3_1.data
bytes written to /home/docs/.radisdb/hitran/downloads__can_be_deleted/NH3/NH3_1.data
bytes written to /home/docs/.radisdb/hitran/downloads__can_be_deleted/NH3/NH3_1.data
bytes written to /home/docs/.radisdb/hitran/downloads__can_be_deleted/NH3/NH3_1.data
bytes written to /home/docs/.radisdb/hitran/downloads__can_be_deleted/NH3/NH3_1.data
bytes written to /home/docs/.radisdb/hitran/downloads__can_be_deleted/NH3/NH3_1.data
bytes written to /home/docs/.radisdb/hitran/downloads__can_be_deleted/NH3/NH3_1.data
bytes written to /home/docs/.radisdb/hitran/downloads__can_be_deleted/NH3/NH3_1.data
bytes written to /home/docs/.radisdb/hitran/downloads__can_be_deleted/NH3/NH3_1.data
bytes written to /home/docs/.radisdb/hitran/downloads__can_be_deleted/NH3/NH3_1.data
bytes written to /home/docs/.radisdb/hitran/downloads__can_be_deleted/NH3/NH3_1.data
bytes written to /home/docs/.radisdb/hitran/downloads__can_be_deleted/NH3/NH3_1.data
bytes written to /home/docs/.radisdb/hitran/downloads__can_be_deleted/NH3/NH3_1.data
bytes written to /home/docs/.radisdb/hitran/downloads__can_be_deleted/NH3/NH3_1.data
bytes written to /home/docs/.radisdb/hitran/downloads__can_be_deleted/NH3/NH3_1.data
bytes written to /home/docs/.radisdb/hitran/downloads__can_be_deleted/NH3/NH3_1.data
bytes written to /home/docs/.radisdb/hitran/downloads__can_be_deleted/NH3/NH3_1.data
bytes written to /home/docs/.radisdb/hitran/downloads__can_be_deleted/NH3/NH3_1.data
bytes written to /home/docs/.radisdb/hitran/downloads__can_be_deleted/NH3/NH3_1.data
bytes written to /home/docs/.radisdb/hitran/downloads__can_be_deleted/NH3/NH3_1.data
bytes written to /home/docs/.radisdb/hitran/downloads__can_be_deleted/NH3/NH3_1.data
bytes written to /home/docs/.radisdb/hitran/downloads__can_be_deleted/NH3/NH3_1.data
bytes written to /home/docs/.radisdb/hitran/downloads__can_be_deleted/NH3/NH3_1.data
bytes written to /home/docs/.radisdb/hitran/downloads__can_be_deleted/NH3/NH3_1.data
bytes written to /home/docs/.radisdb/hitran/downloads__can_be_deleted/NH3/NH3_1.data
bytes written to /home/docs/.radisdb/hitran/downloads__can_be_deleted/NH3/NH3_1.data
bytes written to /home/docs/.radisdb/hitran/downloads__can_be_deleted/NH3/NH3_1.data
bytes written to /home/docs/.radisdb/hitran/downloads__can_be_deleted/NH3/NH3_1.data
bytes written to /home/docs/.radisdb/hitran/downloads__can_be_deleted/NH3/NH3_1.data
bytes written to /home/docs/.radisdb/hitran/downloads__can_be_deleted/NH3/NH3_1.data
bytes written to /home/docs/.radisdb/hitran/downloads__can_be_deleted/NH3/NH3_1.data
bytes written to /home/docs/.radisdb/hitran/downloads__can_be_deleted/NH3/NH3_1.data
Header written to /home/docs/.radisdb/hitran/downloads__can_be_deleted/NH3/NH3_1.header
END DOWNLOAD
                     Lines parsed: 76353
PROCESSED

Data is fetched from http://hitran.org

BEGIN DOWNLOAD: NH3_2
bytes written to /home/docs/.radisdb/hitran/downloads__can_be_deleted/NH3/NH3_2.data
bytes written to /home/docs/.radisdb/hitran/downloads__can_be_deleted/NH3/NH3_2.data
bytes written to /home/docs/.radisdb/hitran/downloads__can_be_deleted/NH3/NH3_2.data
bytes written to /home/docs/.radisdb/hitran/downloads__can_be_deleted/NH3/NH3_2.data
bytes written to /home/docs/.radisdb/hitran/downloads__can_be_deleted/NH3/NH3_2.data
bytes written to /home/docs/.radisdb/hitran/downloads__can_be_deleted/NH3/NH3_2.data
bytes written to /home/docs/.radisdb/hitran/downloads__can_be_deleted/NH3/NH3_2.data
bytes written to /home/docs/.radisdb/hitran/downloads__can_be_deleted/NH3/NH3_2.data
bytes written to /home/docs/.radisdb/hitran/downloads__can_be_deleted/NH3/NH3_2.data
bytes written to /home/docs/.radisdb/hitran/downloads__can_be_deleted/NH3/NH3_2.data
bytes written to /home/docs/.radisdb/hitran/downloads__can_be_deleted/NH3/NH3_2.data
bytes written to /home/docs/.radisdb/hitran/downloads__can_be_deleted/NH3/NH3_2.data
bytes written to /home/docs/.radisdb/hitran/downloads__can_be_deleted/NH3/NH3_2.data
bytes written to /home/docs/.radisdb/hitran/downloads__can_be_deleted/NH3/NH3_2.data
bytes written to /home/docs/.radisdb/hitran/downloads__can_be_deleted/NH3/NH3_2.data
bytes written to /home/docs/.radisdb/hitran/downloads__can_be_deleted/NH3/NH3_2.data
bytes written to /home/docs/.radisdb/hitran/downloads__can_be_deleted/NH3/NH3_2.data
bytes written to /home/docs/.radisdb/hitran/downloads__can_be_deleted/NH3/NH3_2.data
bytes written to /home/docs/.radisdb/hitran/downloads__can_be_deleted/NH3/NH3_2.data
bytes written to /home/docs/.radisdb/hitran/downloads__can_be_deleted/NH3/NH3_2.data
bytes written to /home/docs/.radisdb/hitran/downloads__can_be_deleted/NH3/NH3_2.data
bytes written to /home/docs/.radisdb/hitran/downloads__can_be_deleted/NH3/NH3_2.data
bytes written to /home/docs/.radisdb/hitran/downloads__can_be_deleted/NH3/NH3_2.data
bytes written to /home/docs/.radisdb/hitran/downloads__can_be_deleted/NH3/NH3_2.data
bytes written to /home/docs/.radisdb/hitran/downloads__can_be_deleted/NH3/NH3_2.data
bytes written to /home/docs/.radisdb/hitran/downloads__can_be_deleted/NH3/NH3_2.data
bytes written to /home/docs/.radisdb/hitran/downloads__can_be_deleted/NH3/NH3_2.data
bytes written to /home/docs/.radisdb/hitran/downloads__can_be_deleted/NH3/NH3_2.data
bytes written to /home/docs/.radisdb/hitran/downloads__can_be_deleted/NH3/NH3_2.data
bytes written to /home/docs/.radisdb/hitran/downloads__can_be_deleted/NH3/NH3_2.data
bytes written to /home/docs/.radisdb/hitran/downloads__can_be_deleted/NH3/NH3_2.data
bytes written to /home/docs/.radisdb/hitran/downloads__can_be_deleted/NH3/NH3_2.data
bytes written to /home/docs/.radisdb/hitran/downloads__can_be_deleted/NH3/NH3_2.data
bytes written to /home/docs/.radisdb/hitran/downloads__can_be_deleted/NH3/NH3_2.data
Header written to /home/docs/.radisdb/hitran/downloads__can_be_deleted/NH3/NH3_2.header
END DOWNLOAD
                     Lines parsed: 13662
PROCESSED
Added HITRAN-NH3 database in /home/docs/radis.json

Commence fitting process for LTE spectrum!

[[Fit Statistics]]
    # fitting method   = least_squares
    # function evals   = 16
    # data points      = 1
    # variables        = 1
    chi-square         = 8.9285e-14
    reduced chi-square = 8.9285e-14
    Akaike info crit   = -28.0469393
    Bayesian info crit = -30.0469393
[[Variables]]
    Tgas:  1000.10724 +/- 100.024360 (10.00%) (init = 700)

Successfully finished the fitting process in 4.655549049377441s.
/home/docs/checkouts/readthedocs.org/user_builds/radis/checkouts/latest/radis/misc/curve.py:268: UserWarning:

First spectrum has more resolution than 2nd. Reverse your spectra in interpolation/comparison for a better accuracy

/home/docs/checkouts/readthedocs.org/user_builds/radis/checkouts/latest/radis/misc/curve.py:268: UserWarning:

First spectrum has more resolution than 2nd. Reverse your spectra in interpolation/comparison for a better accuracy

/home/docs/checkouts/readthedocs.org/user_builds/radis/checkouts/latest/radis/misc/curve.py:241: UserWarning:

Presence of NaN in curve_divide!
Think about interpolation=2


======================== END OF FITTING PROCESS ========================


Residual history:

[2.2244774290346794e-05, 2.224477385113398e-05, 6.249578796174649e-06, 6.249580089204621e-06, 6.706460908961787e-07, 6.706472220243689e-07, 3.284683533263895e-07, 3.2846783001901297e-07, 7.091278321790844e-07, 3.0428299211950907e-07, 3.042832337621719e-07, 3.090570218882054e-07, 2.99575991127554e-07, 2.9957609158733825e-07, 2.9880649987254274e-07, 2.988064553529746e-07, 2.9880649987254274e-07]

Fitted values history:

[700.0]
[700.0000104308128]
[1075.6382249193928]
[1075.6382409476514]
[1007.358215806425]
[1007.3582308172322]
[998.6117336654348]
[998.6117485459092]
[1007.8655628712253]
[1000.9251909668825]
[1000.9252058818302]
[999.2889522100642]
[1000.5161312776779]
[1000.51614618653]
[1000.107238643362]
[1000.1072535461212]
[1000.107238643362]

Total fitting time:
655549049377441 s

import astropy.units as u

from radis import load_spec
from radis.test.utils import getTestFile
from radis.tools.new_fitting import fit_spectrum

# ------------------------------------ Step 1. Load experimental spectrum ------------------------------------ #


# Load an experimental spectrum. You can prepare yours, or fetch one of them in the radis/test/files directory.
my_spec = getTestFile("synth-NH3-1-500-2000cm-P10-mf0.01-p1.spec")
s_experimental = load_spec(my_spec)


# ------------------------------------ Step 2. Fill ground-truths and data ------------------------------------ #


# Experimental conditions which will be used for spectrum modeling. Basically, these are known ground-truths.
experimental_conditions = {
    "molecule": "NH3",  # Molecule ID
    "isotope": "1",  # Isotope ID, can have multiple at once
    "wmin": 1000,  # Starting wavelength/wavenumber to be cropped out from the original experimental spectrum.
    "wmax": 1050,  # Ending wavelength/wavenumber for the cropping range.
    "wunit": "cm-1",  # Accompanying unit of those 2 wavelengths/wavenumbers above.
    "mole_fraction": 0.01,  # Species mole fraction, from 0 to 1.
    "pressure": 1e6
    * u.Pa,  # Total pressure of gas, in "bar" unit by default, but you can also use Astropy units.
    "path_length": 10
    * u.mm,  # Experimental path length, in "cm" unit by default, but you can also use Astropy units.
    "slit": "1 nm",  # Experimental slit, must be a blank space separating slit amount and unit.
    "offset": "-0.2 nm",  # Experimental offset, must be a blank space separating offset amount and unit.
    "databank": "hitran",  # Databank used for the spectrum calculation. Must be stated.
}

# List of parameters to be fitted, accompanied by their initial values.
fit_parameters = {
    "Tgas": 700,  # Gas temperature, in K.
}

# List of bounding ranges applied for those fit parameters above.
# You can skip this step and let it use default bounding ranges, but this is not recommended.
# Bounding range must be at format [<lower bound>, <upper bound>].
bounding_ranges = {
    "Tgas": [
        500,
        2000,
    ],
}

# Fitting pipeline setups.
fit_properties = {
    "method": "least_squares",  # Preferred fitting method from the 17 confirmed methods of LMFIT stated in week 4 blog. By default, "leastsq".
    "fit_var": "radiance",  # Spectral quantity to be extracted for fitting process, such as "radiance", "absorbance", etc.
    "normalize": False,  # Either applying normalization on both spectra or not.
    "max_loop": 150,  # Max number of loops allowed. By default, 200.
    "tol": 1e-15,  # Fitting tolerance, only applicable for "lbfgsb" method.
}

"""

For the fitting method, you can try one among 17 different fitting methods and algorithms of LMFIT,
introduced in `LMFIT method list <https://lmfit.github.io/lmfit-py/fitting.html#choosing-different-fitting-methods>`.

You can see the benchmark result of these algorithms here:
`RADIS Newfitting Algorithm Benchmark <https://github.com/radis/radis-benchmark/blob/master/manual_benchmarks/plot_newfitting_comparison_algorithm.py>`.

"""


# ------------------------------------ Step 3. Run the fitting and retrieve results ------------------------------------ #


# Conduct the fitting process!
s_best, result, log = fit_spectrum(
    s_exp=s_experimental,  # Experimental spectrum.
    fit_params=fit_parameters,  # Fit parameters.
    bounds=bounding_ranges,  # Bounding ranges for those fit parameters.
    model=experimental_conditions,  # Experimental ground-truths conditions.
    pipeline=fit_properties,  # Fitting pipeline references.
    fit_kws={"gtol": 1e-12},
)


# Now investigate the result logs for additional information about what's going on during the fitting process

print("\nResidual history: \n")
print(log["residual"])

print("\nFitted values history: \n")
for fit_val in log["fit_vals"]:
    print(fit_val)

print("\nTotal fitting time: ")
print(log["time_fitting"], end=" s\n")

Total running time of the script: (0 minutes 7.062 seconds)

RADIS

Navigation

Quick search

New fitting module introduction - simple 1-temperature LTE case¶