def source_finder(self, image=None, thresh=None, noise=None, lsmname=None, **kw): #TODO look for other source finders and how they operate thresh = thresh or self.pos_smooth image = image or self.imagename ext = utils.fits_ext(image) tpos = tempfile.NamedTemporaryFile(suffix="."+ext, dir=".") tpos.flush() # data smoothing mask, noise = utils.thresh_mask( image, tpos.name, thresh=thresh, noise=self.noise, sigma=True, smooth=True) lsmname = lsmname or self.poslsm # source extraction utils.sources_extraction( image=tpos.name, output=lsmname, sourcefinder_name=self.sourcefinder_name, blank_limit=self.noise/100.0, prefix=self.prefix, **kw)
def source_finder(self, image=None, thresh=None, prefix=None, noise=None, output=None, savemask=None, **kw): #kw.update(kwards) tpos = None naxis = self.header["NAXIS1"] boundary = numpy.array([self.locstep, self.cfstep]) #trim_box = (boundary.max(), naxis - boundary.max(), # boundary.max(), naxis - boundary.max()) trim_box = None # data smoothing if self.smoothing: ext = utils.fits_ext(image) tpos = tempfile.NamedTemporaryFile(suffix="."+ext, dir=".") tpos.flush() mask, noise = utils.thresh_mask(image, tpos.name, thresh=thresh, noise=self.noise, sigma=True, smooth=True, prefix=prefix, savemask=savemask) # using the masked image for forming islands kw["detection_image"] = tpos.name kw["blank_limit"] = self.noise/1.0e5 # source extraction utils.sources_extraction( image=image, output=output, sourcefinder_name=self.sourcefinder_name, trim_box=trim_box, prefix=self.prefix, **kw) if tpos: tpos.close()
def source_finder(self, image=None, thresh=None, prefix=None, noise=None, output=None, savemask=None, **kw): #kw.update(kwards) tpos = None naxis = self.header["NAXIS1"] boundary = numpy.array([self.locstep, self.cfstep]) #trim_box = (boundary.max(), naxis - boundary.max(), # boundary.max(), naxis - boundary.max()) trim_box = None # data smoothing if self.smoothing: ext = utils.fits_ext(image) tpos = tempfile.NamedTemporaryFile(suffix="." + ext, dir=".") tpos.flush() mask, noise = utils.thresh_mask(image, tpos.name, thresh=thresh, noise=self.noise, sigma=True, smooth=True, prefix=prefix, savemask=savemask) # using the masked image for forming islands kw["detection_image"] = tpos.name kw["blank_limit"] = self.noise / 1.0e5 # source extraction utils.sources_extraction(image=image, output=output, sourcefinder_name=self.sourcefinder_name, trim_box=trim_box, prefix=self.prefix, **kw) if tpos: tpos.close()