Basics

All iSpec's functions described in this manual, and some additional ones that cannot be accessed through the visual interface, can be called from a Python script. In the iSpec directory there is a file named 'example.py' which implements an extensive list of examples:

  • read_write_spectrum()
  • convert_air_to_vacuum()
  • plot()
  • cut_spectrum_from_range()
  • cut_spectrum_from_segments()
  • determine_radial_velocity_with_mask()
  • determine_radial_velocity_with_template()
  • correct_radial_velocity()
  • determine_tellurics_shift_with_mask()
  • determine_tellurics_shift_with_template()
  • degrade_resolution()
  • smooth_spectrum()
  • resample_spectrum()
  • coadd_spectra()
  • merge_spectra()
  • normalize_spectrum_using_continuum_regions()
  • normalize_spectrum_in_segments()
  • normalize_whole_spectrum()
  • normalize_whole_spectrum_ignoring_prefixed_strong_lines()
  • normalize_whole_spectrum_with_template()
  • filter_cosmic_rays()
  • find_continuum_regions()
  • find_continuum_regions_in_segments()
  • find_linemasks()
  • fit_lines_determine_ew_and_crossmatch_with_atomic_data(use_ares=False)
  • fit_lines_determine_ew_and_crossmatch_with_atomic_data(use_ares=True)
  • fit_lines_already_crossmatched_with_atomic_data_and_determine_ew(use_ares=False)
  • fit_lines_already_crossmatched_with_atomic_data_and_determine_ew(use_ares=True)
  • calculate_barycentric_velocity()
  • estimate_snr_from_flux()
  • estimate_snr_from_err()
  • estimate_errors_from_snr()
  • clean_spectrum()
  • clean_telluric_regions()
  • adjust_line_masks()
  • create_segments_around_linemasks()
  • synthesize_spectrum(code="spectrum")
  • synthesize_spectrum(code="turbospectrum")
  • synthesize_spectrum(code="sme")
  • synthesize_spectrum(code="moog")
  • synthesize_spectrum(code="synthe")
  • add_noise_to_spectrum()
  • generate_new_random_realizations_from_spectrum()
  • precompute_synthetic_grid(code="spectrum")
  • precompute_synthetic_grid(code="turbospectrum")
  • precompute_synthetic_grid(code="sme")
  • precompute_synthetic_grid(code="moog")
  • precompute_synthetic_grid(code="synthe")
  • determine_astrophysical_parameters_using_synth_spectra(code="spectrum")
  • determine_astrophysical_parameters_using_synth_spectra(code="turbospectrum")
  • determine_astrophysical_parameters_using_synth_spectra(code="sme")
  • determine_astrophysical_parameters_using_synth_spectra(code="moog")
  • determine_astrophysical_parameters_using_synth_spectra(code="synthe")
  • determine_astrophysical_parameters_using_synth_spectra_and_precomputed_grid(code="spectrum")
  • determine_astrophysical_parameters_using_synth_spectra_and_precomputed_grid(code="turbospectrum")
  • determine_astrophysical_parameters_using_synth_spectra_and_precomputed_grid(code="sme")
  • determine_astrophysical_parameters_using_synth_spectra_and_precomputed_grid(code="moog")
  • determine_astrophysical_parameters_using_synth_spectra_and_precomputed_grid(code="synthe")
  • determine_abundances_using_synth_spectra(code="spectrum")
  • determine_abundances_using_synth_spectra(code="turbospectrum")
  • determine_abundances_using_synth_spectra(code="sme")
  • determine_abundances_using_synth_spectra(code="moog")
  • determine_abundances_using_synth_spectra(code="synthe")
  • determine_abundances_line_by_line_using_synth_spectra(code="spectrum")
  • determine_abundances_line_by_line_using_synth_spectra(code="turbospectrum")
  • determine_abundances_line_by_line_using_synth_spectra(code="sme")
  • determine_abundances_line_by_line_using_synth_spectra(code="moog")
  • determine_abundances_line_by_line_using_synth_spectra(code="synthe")
  • determine_loggf_line_by_line_using_synth_spectra(code="spectrum")
  • determine_loggf_line_by_line_using_synth_spectra(code="turbospectrum")
  • determine_loggf_line_by_line_using_synth_spectra(code="sme")
  • determine_loggf_line_by_line_using_synth_spectra(code="moog")
  • determine_loggf_line_by_line_using_synth_spectra(code="synthe")
  • determine_astrophysical_parameters_from_ew(code="moog", use_lines_already_crossmatched_with_atomic_data=True)
  • determine_astrophysical_parameters_from_ew(code="moog", use_lines_already_crossmatched_with_atomic_data=False)
  • determine_astrophysical_parameters_from_ew(code="width", use_lines_already_crossmatched_with_atomic_data=True)
  • determine_astrophysical_parameters_from_ew(code="width", use_lines_already_crossmatched_with_atomic_data=False)
  • determine_astrophysical_parameters_from_ew(code="spectrum", use_lines_already_crossmatched_with_atomic_data=True)
  • determine_astrophysical_parameters_from_ew(code="spectrum", use_lines_already_crossmatched_with_atomic_data=False)
  • determine_astrophysical_parameters_from_ew(code="turbospectrum", use_lines_already_crossmatched_with_atomic_data=True)
  • determine_astrophysical_parameters_from_ew(code="turbospectrum", use_lines_already_crossmatched_with_atomic_data=False)
  • determine_abundances_from_ew(code="spectrum", use_lines_already_crossmatched_with_atomic_data=True)
  • determine_abundances_from_ew(code="turbospectrum", use_lines_already_crossmatched_with_atomic_data=True)
  • determine_abundances_from_ew(code="spectrum", use_lines_already_crossmatched_with_atomic_data=True)
  • determine_abundances_from_ew(code="turbospectrum", use_lines_already_crossmatched_with_atomic_data=True)
  • determine_abundances_from_ew(code="moog", use_lines_already_crossmatched_with_atomic_data=True)
  • determine_abundances_from_ew(code="moog", use_lines_already_crossmatched_with_atomic_data=False)
  • determine_abundances_from_ew(code="width", use_lines_already_crossmatched_with_atomic_data=True)
  • determine_abundances_from_ew(code="width", use_lines_already_crossmatched_with_atomic_data=False)
  • determine_abundances_from_ew(code="moog", use_lines_already_crossmatched_with_atomic_data=True)
  • determine_abundances_from_ew(code="width", use_lines_already_crossmatched_with_atomic_data=True)
  • calculate_theoretical_ew_and_depth()
  • paralelize_code()
  • estimate_vmic_from_empirical_relation()
  • estimate_vmac_from_empirical_relation()
  • generate_and_plot_YY_isochrone()
  • interpolate_atmosphere(code="spectrum")
  • interpolate_atmosphere(code="turbospectrum")
  • interpolate_atmosphere(code="moog")

The easiest way to use the script is to duplicate it, and erase all the functions except the ones that fit the user's needs. The variable ispec_dir (inside the script) should point to the right directory where iSpec is installed, for instance:

#--- iSpec directory -------------------------------------------------------------
ispec_dir = '/home/virtual/iSpec/'
sys.path.insert(0, os.path.abspath(ispec_dir))
import ispec

The script can be executed by writing on a terminal "python example.py".