sigma1 = noise.noise_by_data(FB,'FALSE') print clumps1 cmd = '%s/findclumps in=%s out=%s outcat=%s method=FellWalker deconv=True config=%s RMS=%s '%(cupdir,MSK,clumps1,cat1,config,sigma1) os.system(cmd) #TEST 2 #### FINDBACK + MAKESNR print '2 FINDBACK + MAKESNR' FBsnr = 'FW/temp/'+file+'FBSNR.sdf' cmd = '%s/makesnr in=%s out=%s'%(kapdir,FB,FBsnr) os.system(cmd) #snrMSK = 'FW/temp/'+file+'SNR_MSK.sdf' #mult = '%s/mult in1=%s in2=%s out=%s'%(kapdir,FBsnr,mask,snrMSK) #os.system(mult) clumps2 = 'FW/'+region+'clumpsSNR.sdf' cat2 = 'FW/cat/'+region+'clump_catSNR.fit' sigma2 = noise.noise_by_snr(FBsnr,'FALSE') cmd = '%s/findclumps in=%s out=%s outcat=%s method=FellWalker deconv=True config=%s RMS=%s '%(cupdir,FBsnr,clumps2,cat2,config,sigma2) os.system(cmd) print clumps2 #TEST 3 #### MAKESNR + FINDBACK print '3 MAKESNR + FINDBACK' snr = 'FW/temp/'+file+'SNR.sdf' cmd = '%s/makesnr in=%s out=%s'%(kapdir,map,snr) #os.system(cmd) snrFB = 'FW/temp/'+file+'SNRFB.sdf' print 'Findback box size ',box cmd = '%s/findback in=%s out=%s box=%s sub=TRUE %s'%(cupdir,snr,snrFB,box,'> /dev/null') #os.system(cmd) #snrFBMSK = 'FW/temp/'+file+'SNRFB_MSK.sdf' #mult = '%s/mult in1=%s in2=%s out=%s'%(kapdir,snrFB,mask,snrFBMSK)
snr = 'FW/temp/'+file+'SNR.sdf'
cmd = '%s/makesnr in=%s out=%s minvar=0'%(kapdir,i,snr)
os.system(cmd)
#box = 120
#print box
#snrFB = 'FW/temp/'+file+'snrFB.sdf'
#cmd = '%s/findback in=%s out=%s box=%s sub=TRUE %s'%(cupdir,snr,snrFB,box,'> /dev/null')
#os.system(cmd)
#PART 2 #### FIND clumps
print '==============================================='
print 'Finding clumps in: '+snr
sigma = noise.noise_by_snr(snr,'FALSE')
clumps = 'FW/clumps/'+region+'clumps.sdf'
cat = 'FW/cat/'+region+'clump_cat.fit'
config = '^inparams.txt'
cmd = '%s/findclumps in=%s out=%s outcat=%s method=FellWalker deconv=True config=%s RMS=%s '%(cupdir,snr,clumps,cat,config,sigma)
os.system(cmd)
print 'CLUMP MAP:'+clumps
j = j + 1
#break
tempfolder = 'FW/temp/*'
