# Calculate the S/N and coordinates
rdxqa = ReductionAssessment('SNRG', drpf, clobber=clobber, analysis_path=analysis_path)
# Peform the Voronoi binning to S/N>~10
binned_spectra = SpatiallyBinnedSpectra('VOR10', drpf, rdxqa, clobber=clobber,
analysis_path=analysis_path)
# Fit the stellar kinematics
stellar_continuum = StellarContinuumModel('GAU-MILESHC', binned_spectra, clobber=clobber,
guess_vel=vel, guess_sig=100.,
analysis_path=analysis_path)
# Get the emission-line moments
emission_line_moments = EmissionLineMoments('EMOMF', binned_spectra, clobber=clobber,
stellar_continuum=stellar_continuum,
redshift=nsa_redshift,
analysis_path=analysis_path)
# Get an estimate of the redshift of each bin using the first moment
# of the H-alpha emission line:
el_init_redshift = numpy.full(binned_spectra.nbins, nsa_redshift, dtype=float)
# HARDCODED FOR A SPECIFIC EMISSION-LINE MOMENT DATABASE
# TODO: Should
# - pass the EmissionLineMoments object to EmissionLineModel
# - include the channel used for this in EmissionLineModelDef
halpha_channel = 7
halpha_mom1_masked = emission_line_moments['ELMMNTS'].data['MASK'][:,halpha_channel] > 0
# - Use the 1st moment of the H-alpha line
el_init_redshift[ emission_line_moments['ELMMNTS'].data['BINID_INDEX'] ] \
= emission_line_moments['ELMMNTS'].data['MOM1'][:,halpha_channel] \
/ astropy.constants.c.to('km/s').value
binned_spectra = SpatiallyBinnedSpectra('Aperture', # Key for binning method
cube, # DRP data to bin
rdxqa, # Cube coordinates and S/N assessments
method_list=binning_method, # Methods to select from
output_path=output_path)
# The rest of this is just a single execution of the remaining
# analysis steps in
# $MANGADAP_DIR/python/mangadap/survey/manga_dap.py , with some
# simplifications
stellar_continuum = StellarContinuumModel('GAU-MILESHC', binned_spectra, guess_vel=vel,
guess_sig=100., output_path=output_path)
emission_line_moments = EmissionLineMoments('EMOMM', binned_spectra,
stellar_continuum=stellar_continuum,
redshift=nsa_redshift, output_path=output_path)
emission_line_model = EmissionLineModel('EFITM', binned_spectra,
stellar_continuum=stellar_continuum,
redshift=nsa_redshift, dispersion=100.0,
output_path=output_path)
spectral_indices = SpectralIndices('INDXEN', binned_spectra, redshift=nsa_redshift,
stellar_continuum=stellar_continuum,
emission_line_model=emission_line_model,
output_path=output_path)
construct_maps_file(cube, rdxqa=rdxqa, binned_spectra=binned_spectra,
stellar_continuum=stellar_continuum,
emission_line_moments=emission_line_moments,