我终于（有点）让这个脚本工作了，除了现在我收到一条新的错误消息。. 在“提取目录并识别中央星系”之后似乎失败了。
代码：
def objmask(inimgs, inwhts, thresh1='20.0', thresh2='2.0', tfdel=True,
xceng=65., yceng=65., outdir='.', tmpdir='tmp'):
print np.size(inimgs)
# initial detection of main galaxy with SExtractor for re-centering purposes
if outdir!='.':
if not os.path.exists(outdir):
os.makedirs(outdir)
if not os.path.exists(tmpdir):
os.makedirs(tmpdir)
for c in range(np.size(inimgs)):
print 'Creating Aperture Run:', c
subprocess.call(['sex',inimgs[c],'-c','gccg.sex',
'-CATALOG_NAME','_tmp_seobj'+str(c)+'.cat',
'-PARAMETERS_NAME','gccg_ell.param',
'-FILTER_NAME','gccg.conv',
'-STARNNW_NAME','gccg.nnw',
# '-CHECKIMAGE_TYPE','APERTURES',
'-VERBOSE_TYPE','QUIET',
'-DETECT_THRESH',thresh1,
'-ANALYSIS_THRESH',thresh2,
'-WEIGHT_IMAGE',inwhts[c]],
)
# extract catalog and identify central galaxy
secat=asciitable.read('_tmp_seobj'+str(c)+'.cat',
names=['flux','ferr','xmin','ymin','xmax','ymax',
'xc','yc','cxx','cyy','cxy'])
robj = np.sqrt((secat['xc']-xceng)**2.0+(secat['yc']-yceng)**2.0)
rmin = (robj==np.min(robj))
wrmin = np.where(robj==np.min(robj))
xc_min,yc_min = secat['xc'][rmin],secat['yc'][rmin]
print 'extract catalog complete'
# shift images and masks to a common center
hdu=pf.open(inimgs[c])
img=hdu[0].data
xdel=xceng-xc_min
ydel=yceng-yc_min
img_sh=shift(img,[ydel,xdel])
hdu2=pf.open(inwhts[c])
mask=hdu2[0].data
mask_sh=shift(mask,[yceng-yc_min,xceng-xc_min])
xmin=np.delete(secat['xmin']+xdel,wrmin)
xmax=np.delete(secat['xmax']+xdel,wrmin)
xcen=np.delete(secat['xc']+xdel,wrmin)
ymin=np.delete(secat['ymin']+ydel,wrmin)
ymax=np.delete(secat['ymax']+ydel,wrmin)
ycen=np.delete(secat['yc']+ydel,wrmin)
print 'shift complete'
# mask detected objects and write mask file
#retgen.pbar ((float(c)*ncols+1.0)/18.0)
omask=tgen.semask(img.shape,xmin,xmax,ymin,ymax,xcen,ycen,
np.delete(secat['cxx'],wrmin),
np.delete(secat['cyy'],wrmin),
np.delete(secat['cxy'],wrmin),2.0)
hdimsh=pf.PrimaryHDU(img_sh)
hdmsh=pf.PrimaryHDU(omask)
with warnings.catch_warnings():
warnings.simplefilter('ignore')
hdimsh.writeto(outdir+'/sh_img_'+str(c)+'.fits',
clobber=True)
hdmsh.writeto(tmpdir+'/_omask.fits',clobber=True)
print 'mask file complete'
# write combined g+r+i image and mask
if c==2:
cimg/=3.0
csub/=3.0
cmask/=3.0
hdcim_out=pf.PrimaryHDU(cimg)
hdmod_out=pf.PrimaryHDU(csub)
hdcm_out=pf.PrimaryHDU(cmask)
with warnings.catch_warnings():
warnings.simplefilter('ignore')
hdcim_out.writeto(outdir+'/im_comb.fits',
clobber=True)
hdmod_out.writeto(tmpdir+'/_im-mod_comb.fits',clobber=True)
hdcm_out.writeto(tmpdir+'/_comb_mask.fits',clobber=True)
错误：
In [15]: fetch_swarp4.objmask(['sciPHOTOf105w7.fits'],['whtPHOTOf105w7.fits'],thresh1='20.0',thresh2='2.0',tfdel=True,xceng=65.,yceng=65.,outdir='.',tmpdir='tmp')
1
Creating Aperture Run: 0
---------------------------------------------------------------------------
ValueError Traceback (most recent call last)
/home/vidur/se_files/<ipython-input-15-92433178d837> in <module>()
----> 1 fetch_swarp4.objmask(['sciPHOTOf105w7.fits'],['whtPHOTOf105w7.fits'],thresh1='20.0',thresh2='2.0',tfdel=True,xceng=65.,yceng=65.,outdir='.',tmpdir='tmp')
/home/vidur/se_files/fetch_swarp4.pyc in objmask(inimgs, inwhts, thresh1, thresh2, tfdel, xceng, yceng, outdir, tmpdir)
33 'xc','yc','cxx','cyy','cxy'])
34 robj = np.sqrt((secat['xc']-xceng)**2.0+(secat['yc']-yceng)**2.0)
---> 35 rmin = (robj==np.min(robj))
36 wrmin = np.where(robj==np.min(robj))
37 xc_min,yc_min = secat['xc'][rmin],secat['yc'][rmin]
/usr/lib/python2.7/dist-packages/numpy/core/fromnumeric.pyc in amin(a, axis, out)
1893 except AttributeError:
1894 return _wrapit(a, 'min', axis, out)
-> 1895 return amin(axis, out)
1896
1897
ValueError: zero-size array to minimum.reduce without identity
坦率地说，我完全不知道发生了什么，我完全迷失了。它在 rmin 失败，但我不知道错误是什么意思，也不知道为什么会这样。
PS 文件类型是“适合”文件，用于天文学。它基本上是一个可以用来推断大量天文数据的数组。