Still trouble with rimexam

[karlglazebrook]

rimexam still does not work. Running it with the m51 fits image from iraf:

Current image /Users/karl/Dropbox/Templates/m51_iraf.fits
using model: <class 'astropy.modeling.functional_models.Gaussian1D'>
Name: Gaussian1D
N_inputs: 1
N_outputs: 1
Fittable parameters: ('amplitude', 'mean', 'stddev')
xc=440.9853	yc=408.6807
(443,411) mean=   441.04, fwhm=     2.75
using model: <class 'astropy.modeling.functional_models.Gaussian1D'>
Name: Gaussian1D
N_inputs: 1
N_outputs: 1
Fittable parameters: ('amplitude', 'mean', 'stddev')
xc=429.1555	yc=405.2307
(435,407) mean=   397.66, fwhm=   125.70
using model: <class 'astropy.modeling.functional_models.Gaussian1D'>
Name: Gaussian1D
N_inputs: 1
N_outputs: 1
Fittable parameters: ('amplitude', 'mean', 'stddev')
xc=440.9852	yc=408.6807
(438,408) mean=   441.04, fwhm=     2.75
using model: <class 'astropy.modeling.functional_models.Gaussian1D'>
Name: Gaussian1D
N_inputs: 1
N_outputs: 1
Fittable parameters: ('amplitude', 'mean', 'stddev')
xc=440.9852	yc=408.6807
(440,407) mean=   408.72, fwhm=     2.53
using model: <class 'astropy.modeling.functional_models.Gaussian1D'>
Name: Gaussian1D
N_inputs: 1
N_outputs: 1
Fittable parameters: ('amplitude', 'mean', 'stddev')
xc=384.6102	yc=253.0607
(386,253) mean=   253.48, fwhm=     6.81
using model: <class 'astropy.modeling.functional_models.Gaussian1D'>
Name: Gaussian1D
N_inputs: 1
N_outputs: 1
Fittable parameters: ('amplitude', 'mean', 'stddev')
xc=403.5343	yc=273.3214
(406,275) mean=   273.36, fwhm=     2.55
using model: <class 'astropy.modeling.functional_models.Gaussian1D'>
Name: Gaussian1D
N_inputs: 1
N_outputs: 1
Fittable parameters: ('amplitude', 'mean', 'stddev')
xc=437.9990	yc=25.4316
(442,23) mean=    24.80, fwhm=     5.25
using model: <class 'astropy.modeling.functional_models.Gaussian1D'>
Name: Gaussian1D
N_inputs: 1
N_outputs: 1
Fittable parameters: ('amplitude', 'mean', 'stddev')
xc=378.1463	yc=65.8309
(376,67) mean=    65.87, fwhm=     2.52
using model: <class 'astropy.modeling.functional_models.Gaussian1D'>
Name: Gaussian1D
N_inputs: 1
N_outputs: 1
Fittable parameters: ('amplitude', 'mean', 'stddev')
xc=440.9852	yc=408.6807
(440,405) mean=   408.72, fwhm=     2.53
using model: <class 'astropy.modeling.functional_models.Gaussian1D'>
Name: Gaussian1D
N_inputs: 1
N_outputs: 1
Fittable parameters: ('amplitude', 'mean', 'stddev')
xc=346.5188	yc=230.6472
(345,230) mean=   346.57, fwhm=     2.66
using model: <class 'astropy.modeling.functional_models.Gaussian1D'>
Name: Gaussian1D
N_inputs: 1
N_outputs: 1
Fittable parameters: ('amplitude', 'mean', 'stddev')
xc=464.5063	yc=61.1100
(460,62) mean=   464.52, fwhm=     2.58
using model: <class 'astropy.modeling.functional_models.Gaussian1D'>
Name: Gaussian1D
N_inputs: 1
N_outputs: 1
Fittable parameters: ('amplitude', 'mean', 'stddev')
xc=440.9852	yc=408.6807
(441,405) mean=   441.04, fwhm=     2.75
using model: <class 'astropy.modeling.functional_models.Gaussian1D'>
Name: Gaussian1D
N_inputs: 1
N_outputs: 1
Fittable parameters: ('amplitude', 'mean', 'stddev')
xc=440.9852	yc=408.6807
(441,405) mean=   408.72, fwhm=     2.53
using model: <class 'astropy.modeling.functional_models.Gaussian1D'>
Name: Gaussian1D
N_inputs: 1
N_outputs: 1
Fittable parameters: ('amplitude', 'mean', 'stddev')
xc=440.9852	yc=408.6807
(441,405) mean=   441.04, fwhm=     2.75
using model: <class 'astropy.modeling.functional_models.Gaussian1D'>
Name: Gaussian1D
N_inputs: 1
N_outputs: 1
Fittable parameters: ('amplitude', 'mean', 'stddev')
xc=440.9852	yc=408.6807
(441,405) mean=   408.72, fwhm=     2.53
xc=440.9849	yc=408.6810
Max. pix. flux =  3164.000
amp =  3509.763
fwhm =     2.851
Traceback (most recent call last):
  File "<stdin>", line 1, in <module>
  File "/opt/anaconda3/envs/imexam2/lib/python3.10/site-packages/imexam/connect.py", line 150, in imexam
    self._run_imexam()
  File "/opt/anaconda3/envs/imexam2/lib/python3.10/site-packages/imexam/connect.py", line 264, in _run_imexam
    self.exam.do_option(
  File "/opt/anaconda3/envs/imexam2/lib/python3.10/site-packages/imexam/imexamine.py", line 201, in do_option
    self.imexam_option_funcs[key][0](x, y, self._data)
  File "/opt/anaconda3/envs/imexam2/lib/python3.10/site-packages/imexam/imexamine.py", line 1256, in radial_profile
    ax.set_xlim(0, datasize, 0.5)
TypeError: _AxesBase.set_xlim() takes from 1 to 3 positional arguments but 4 were given

I tried a different FITS image and got the same result.

The offending line of code is:

ax.set_xlim(0, datasize, 0.5)
I think this should only take two args so I do not know what is going on. Removing the '0.5' does make it work.

[sosey]

hi @karlglazebrook, apologies, this project hasn't been supported for a while and I can't work on during my normal hours. I still need to get this repo transfered out of spacetelescope and into my personal repository so that I can work on it in my free time. In the meantime, space telescope is supporting JDAVIZ, and IMVIZ which supports radial profiles, also through photutils, and you might try it out? https://jdaviz.readthedocs.io/en/stable/imviz/plugins.html#plot-options

[karlglazebrook]

No worries. I thought it was 'official'! If not I will recalibrate my expectations :-)

Basically I am a dinosaur imexam user (it is basically the only thing I still use iraf for) and just want something that will interact with DS9 for the radial etc plots. I don't think the **VIZ tools will interact with DS9?

If this is now a private GitHub project maybe I can contribute...