mpdaf write subcube

emi1.write('mycube.fits') failed

In [96]: emi1.info() 160 X 10 X 20 cube (no name) .data(160,10,20) (erg/s/cm**2/Angstrom) fscale=1e-20, .var(160,10,20) center:(-30:00:00.45,01:20:00.438) size in arcsec:(2.000,4.000) step in arcsec:(0.200,0.200) rot:0.0 wavelength: min:9001.25 max:9200.00 step:1.25 Angstrom

In [97]: emi1.w emi1.wave emi1.wcs emi1.write

In [97]: emi1.write('mycube.fits')

IOError Traceback (most recent call last)

/data/MUSE/mpdaf/tutorial/data/ in ()

/usr/local/lib/python2.6/dist-packages/mpdaf-1.0.1-py2.6.egg/mpdaf/obj/cube.pyc in write(self, filename) 386 # for card in wcs_cards:

387 #                tbhdu.header.update(card.key, card.value, card.comment)

--> 388 cd = self.wcs.get_cd() 389 tbhdu.header.update('CTYPE1', wcs_cards['CTYPE1'].value, wcs_cards['CTYPE1'].comment) 390 tbhdu.header.update('CUNIT1', wcs_cards['CUNIT1'].value, wcs_cards['CUNIT1'].comment)

/usr/local/lib/python2.6/dist-packages/mpdaf-1.0.1-py2.6.egg/mpdaf/obj/coords.pyc in get_cd(self) 440 return cd 441 except: --> 442 raise IOError, 'No standard WCS' 443 444 def set_naxis1(self,n):

IOError: No standard WCS

was created with emi1 = obj[k1+1:k2+1,:,:] k1,k2 = sp1.wave.pixel([9000,9200], nearest=True)

In [109]: obj.info() 3601 X 10 X 20 cube (no name) .data(3601,10,20) (erg/s/cm**2/Angstrom) fscale=1e-20, .var(3601,10,20) center:(-30:00:00.45,01:20:00.438) size in arcsec:(2.000,4.000) step in arcsec:(0.200,0.200) rot:0.0 wavelength: min:4800.00 max:9300.00 step:1.25 Angstrom

In [111]: obj=cube[:,55-5:55+5,31-10:31+10] cube=Cube('Central_DATACUBE_FINAL_11to20_2012-05-16.fits')