PyFITS:文件已经存在

时间:2013-07-18 01:14:22

标签: python fits pyfits

我真的很接近完成一个大代码,但它的最后一部分似乎失败了,我不知道为什么。我在这里要做的是采用图像阵列,将其与不同的图像阵列进行比较,并且在初始图像阵列等于1的任何地方,我想在第二个图像阵列中屏蔽该部分。但是,我收到一个奇怪的错误:

代码:

maskimg='omask'+str(inimgs)[5:16]+'.fits'
newmaskimg=pf.getdata(maskimg)
oimg=pf.getdata(inimgs)
for i in range (newmaskimg.shape[0]):
    for j in range (newmaskimg.shape[1]):
        if newmaskimg[i,j]==1:
            oimg[i,j]=0
pf.writeto('newestmask'+str(inimgs)[5:16]+'.fits',newmaskimg)

错误:

/home/vidur/se_files/fetch_swarp10.py in objmask(inimgs, inwhts, thresh1, thresh2, tfdel, xceng, yceng, outdir, tmpdir)
    122             if newmaskimg[i,j]==1:
    123                 oimg[i,j]=0
--> 124     pf.writeto('newestmask'+str(inimgs)[5:16]+'.fits',newmaskimg)
    125 
    126 

/usr/local/lib/python2.7/dist-packages/pyfits/convenience.pyc in writeto(filename, data, header, output_verify, clobber, checksum)
    396         hdu = PrimaryHDU(data, header=header)
    397     hdu.writeto(filename, clobber=clobber, output_verify=output_verify,
--> 398                 checksum=checksum)
    399 
    400 

/usr/local/lib/python2.7/dist-packages/pyfits/hdu/base.pyc in writeto(self, name, output_verify, clobber, checksum)
    348         hdulist = HDUList([self])
    349         hdulist.writeto(name, output_verify, clobber=clobber,
--> 350                         checksum=checksum)
    351 
    352     def _get_raw_data(self, shape, code, offset):

/usr/local/lib/python2.7/dist-packages/pyfits/hdu/hdulist.pyc in writeto(self, fileobj, output_verify, clobber, checksum)
    651                     os.remove(filename)
    652                 else:
--> 653                     raise IOError("File '%s' already exists." % filename)
    654         elif (hasattr(fileobj, 'len') and fileobj.len > 0):
    655             if clobber:

IOError: File 'newestmaskPHOTOf105w0.fits' already exists.

3 个答案:

答案 0 :(得分:2)

如果您不关心覆盖现有文件,pyfits.writeto会接受clobber参数自动覆盖现有文件(它仍会输出警告):

pyfits.writeto(..., clobber=True)

顺便说一下,让我非常强调你上面发布的代码非常使用Numpy的正确方法。代码中的循环可以写成一行,速度会快几个数量级。例如,许多可能性之一就是这样写:

oimg[newmaskimg == 1] = 0

答案 1 :(得分:0)

是的,添加clobber = True。我之前在我的代码中使用了它,它运行得很好。或者,只需转到sudo rm path/to/file并删除它们,以便再次运行它。

答案 2 :(得分:0)

我遇到了同样的问题,事实证明,使用参数clobber仍然有效,但将来的AstroPy版本将不支持该问题。

参数overwrite的作用相同,并且不会发出错误消息。