Changeset 69 for trunk

Timestamp:

Sep 25, 2014 12:51:44 PM (10 years ago)

Author:

JohnLightfoot

Message:

clean central quarter only

File:

• trunk/cleanimage.py (modified) (2 diffs)

trunk/cleanimage.py

@@ -34 +34 @@
         residualimage = np.zeros(np.shape(dirtyimage), np.float)
 
-#        for iwn,wn in enumerate(wavenumber):
-#            print 'cleaning %s' % wn
-#            max = np.max(dirtyimage[:,:,iwn])
-#            components, residual = pclean.hogbom(dirty=dirtyimage[:,:,iwn],
-#              psf=dirtybeam[:,:,iwn],
-#              window=True, gain=0.1, thresh=max/100.0, niter=5000)
-#            cleanimage[:,:,iwn] = components
-#            residualimage[:,:,iwn] = residual
-
         # calculate clean image for each wn
         jobs = {}
+        # clean only central quarter
+        imsize = np.shape(dirtyimage)[:2]
+        window = [slice(imsize[0]/4, imsize[0]*3/4), slice(imsize[1]/4, imsize[1]*3/4)]
         for iwn,wn in enumerate(wavenumber):
             # submit jobs
-            print 'starting clean for', wn
             immax = np.max(dirtyimage[:,:,iwn])
-            indata = (dirtyimage[:,:,iwn], dirtybeam[:,:,iwn], True, 0.1,
+            indata = (dirtyimage[:,:,iwn], dirtybeam[:,:,iwn], window, 0.1,
               immax/100.0, 5000,)
             jobs[wn] = self.job_server.submit(pythonclean.hogbom,
@@ -59 +52 @@
 
         self.result['dirtyimage'] = dirtyimage
-        self.result['cleanimage'] = cleanimage
+        self.result['cleanimage'] = cleanimage + residualimage
         self.result['residualimage'] = residualimage
         self.result['spatial axis [arcsec]'] = spatial_axis