iraf-community / iraf

IRAF - Image Reduction and Analysis Facility
https://iraf-community.github.io
Other
152 stars 52 forks source link

iraf disptrans EXTVER error #100

Open basriaz opened 4 years ago

basriaz commented 4 years ago

Is it possible to convert the wavelength scale to velocity scale for a MEF fits cube with dimensions [x,y,lambda] ?

When running disptrans,

disptrans diff_spec.fits test.fits angs

I get the error

Warning: Image header parameter not found (EXTVER)

I tried hedit diff_spec.fits EXTVER 3 add+

But got the same error when running disptrans again.

olebole commented 4 years ago

This seems to be an ancient problem; see iraf.net#136539 from 2006:

My question concerned disptrans. […] Input spectra (test0082): Output spectra (vdim0082): vtest0082 Output dispersion units (km/s 2.12125 micron): Warning: Image header parameter not found (EXTVER)As you see it crashes telling me that the image parameter EXTVER was not found. I have no clue what that image parameter would mean or what value it should have.

The message about EXTVER is actually bogus due to the way iraf error recovery works and is masking the real error (usually a missing header keyword or parameter problem). To see the real error just add a dummy keyword EXTVER=1 to the image and run the task again, you'll then need to address the error reported which will be the real problem. Note that EXTVER usually refers to multi-extension FITS files so the problem may just be that you need to add an extension number when calling the task (e.g. "disptrans test0082[1]"). Similar problems have also been reported when using a 2-D image, the solution there is to either SCOPY out a particular band/aperture, or use an image section. Hope this helps, I'll post a summary to a forum and you should followup there. Sorry again about the account problems, I see you have a login now.Cheers, -Mike

Or from 1996 (iraf.net#101761):

Date: Mon, 14 Oct 96 16:13:21 MST From: valdes (Frank Valdes) To: hamann@c[…]u Subject: Re: disptrans

Hello Fred,Your problem is that your data does not identify the units of the spectra; i.e. it does not know that the data is in angstroms. To fix this:

cl> wcsedit t076i222f unit angstroms 1 inter- update+

This will then define this data has being in angstroms. DISPTRANS will then know that you want to transform the angstroms to angstroms in vacuum.Let me know if you still have problems.Cheers,

Does this help you in some way?

basriaz commented 4 years ago

I ran

cl> wcsedit diff_H2-2_12.fits unit micron 3 inter- update+

cl> disptrans diff_H2-2_12.fits[3] test.fits km/s 2.12 micron

There's no error but there's no output file created.

What am I missing? Why did it not write the output?

olebole commented 4 years ago

Can you attach the relevant header of the file?

basriaz commented 4 years ago

header.txt

olebole commented 4 years ago

Thank you! I can't promise a schedule, but I will look into this.

basriaz commented 4 years ago

It was easy to do this unit conversion of a cube in python using the SpectralCube package.

On Fri, Mar 6, 2020 at 11:37 AM Ole Streicher notifications@github.com wrote:

Thank you! I can't promise a schedule, but I will look into this.

— You are receiving this because you authored the thread. Reply to this email directly, view it on GitHub https://github.com/iraf-community/iraf/issues/100?email_source=notifications&email_token=AFCZYHMUZLCZ6X7X7X6NNWDRGDG6FA5CNFSM4K5VM6K2YY3PNVWWK3TUL52HS4DFVREXG43VMVBW63LNMVXHJKTDN5WW2ZLOORPWSZGOEOA4VBA#issuecomment-595708548, or unsubscribe https://github.com/notifications/unsubscribe-auth/AFCZYHOKUUHQQHUT5IHBFSLRGDG6FANCNFSM4K5VM6KQ .