IRAF can't deal with .fits files


by PLuz
Tags: astronomy, deal, files, fits, iraf
PLuz
PLuz is offline
#1
Oct5-13, 06:46 AM
P: 56
I've been (trying) to learn astronomy, more specifically image reduction, using IRAF. I have, though, been having problems with the fits files from the Sloan Digital Sky Survey site. For example from this archive. I can't do anything with them. Implot, imstat all give errors. Using python, numdisplay can't work with this files either, although ds9 works fine.

I assume that with IRAF everything is fine since it can work with .fits from other surveys. My insistence with the above archive is the availability of the files (I haven't been able to find other source of unprocessed files...)

Can somebody help me? Is something wrong with the files or it might be something with IRAF?
Phys.Org News Partner Astronomy news on Phys.org
Unique pair of supermassive black holes in an ordinary galaxy discovered
Red stars and big bulges: How black holes shape galaxies
Image: X-raying the cosmos
phyzguy
phyzguy is offline
#2
Oct5-13, 09:05 AM
P: 2,070
Quote Quote by PLuz View Post
I've been (trying) to learn astronomy, more specifically image reduction, using IRAF. I have, though, been having problems with the fits files from the Sloan Digital Sky Survey site. For example from this archive. I can't do anything with them. Implot, imstat all give errors. Using python, numdisplay can't work with this files either, although ds9 works fine.

I assume that with IRAF everything is fine since it can work with .fits from other surveys. My insistence with the above archive is the availability of the files (I haven't been able to find other source of unprocessed files...)

Can somebody help me? Is something wrong with the files or it might be something with IRAF?
I've never used IRAF, but Marthematica can open FITS files and allow you to examine them or plot the data inside. Python also has a module called Pyfits (available for download here). Here is some simple code I use to examine the FITS headers and data. One of these two should allow you to determine if the problem is with the files or with IRAF, although if ds9 works then the files are probably OK.

import pyfits

file = pyfits.open("your_file.fits")
header = file[0].header
data = file[0].data

print "Header items are:",header.keys()
for name in header.keys():
    print name, " is ",header[name]
print "\n"
print "Data size is:", data.shape
print "First data element is ",data[0,0]
PLuz
PLuz is offline
#3
Oct5-13, 09:24 AM
P: 56
How happy would I be if instead of IRAF I could just use Mathematica, matlab or even python, I hate that there isn't a viable alternative...

I just tried your suggestion and right after I enter the first command (after the import) pyraf tells me:

The following header keyword is invalid or follows an unrecognized non-standard convention:
SDSS                           / THIS HEADER CONTAINS SPECIAL SDSS KEYWORDS     
The following header keyword is invalid or follows an unrecognized non-standard convention:
UNSIGNED                       / UNSIGNED DATA - NO BSCALE OR BZERO KEYWORDS
So it seems the problem is in the file. My question is then changed to: How do we make IRAF deal with this errors?

phyzguy
phyzguy is offline
#4
Oct5-13, 12:40 PM
P: 2,070

IRAF can't deal with .fits files


Quote Quote by PLuz View Post
How happy would I be if instead of IRAF I could just use Mathematica, matlab or even python, I hate that there isn't a viable alternative...

I just tried your suggestion and right after I enter the first command (after the import) pyraf tells me:

The following header keyword is invalid or follows an unrecognized non-standard convention:
SDSS                           / THIS HEADER CONTAINS SPECIAL SDSS KEYWORDS     
The following header keyword is invalid or follows an unrecognized non-standard convention:
UNSIGNED                       / UNSIGNED DATA - NO BSCALE OR BZERO KEYWORDS
So it seems the problem is in the file. My question is then changed to: How do we make IRAF deal with this errors?
What are you trying to do with the data that you can't do with ds9 or Python? Have you tried looking at this site, which has tools for vieing the SDSS data?
Chronos
Chronos is offline
#5
Oct6-13, 01:56 AM
Sci Advisor
PF Gold
Chronos's Avatar
P: 9,183
You may find this useful http://www.twilightlandscapes.com/IR...Fintro_01.html
PLuz
PLuz is offline
#6
Oct6-13, 08:28 AM
P: 56
To phyzguy:
There's nothing specific that I want to do with the images that I can't do with python, the thing is that I don't have time to program all the tasks on python. Sure there are some repositories but they don't have all the tasks and time is an issue right now.

To Chronos:
That site is really a big help. Thank you for this reference it has proven to be very helpful.
irafnet
irafnet is offline
#7
Oct9-13, 01:47 PM
P: 1
The SDSS FITS files are in MEF (multi-extension FITS) format, and so you need to specify the image extension to the iraf task, e.g.

cl> imstat sdss.fits[1]

where the "[1]" specifies image extension one (the usual image, others include BINTABLES that cannot be used with image tasks). The FITSUTIL external package has some tasks for dealing with MEF that may help, esp the FXHEADER task that will list the extensions. If you just want to deal with the image you can copy it out of the file using

cl> imcopy sdss.fits[1] new.fits
matt.o
matt.o is offline
#8
Oct20-13, 03:09 AM
P: 391
If you want to combine python and iraf, you may find pyraf useful:

http://www.stsci.edu/institute/software_hardware/pyraf


Register to reply

Related Discussions
IRAF for the beginner Astrophysics 1
photometric data from SDSS fits files Astrophysics 1
Why do we use IRAF? General Astronomy 5
conversion of healpix fits file to normal fits file Cosmology 3