Beispiel #1
0
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)
Beispiel #2
0
    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/*'