Python changebox 예제들

예제 #1

from smooth_edge_bump import smooth_edge_bump
            
name='M31'
hdu = pyfits.open('resized_source/{0}_resized.fits'.format(name))
data = hdu[0].data.T
hdu.close()

plt.matshow(np.log10(data.T),origin='lower')
plt.colorbar()
plt.show()

data_fix=data
#==============================================================================
# if changebox
#==============================================================================
data_fix=changebox(data_fix,[249,249], [20, 26],[291, 184])
data_fix=changebox(data_fix,[178,102], [9, 8],[199, 85])
#==============================================================================
# if need smooth
#==============================================================================
data_fix=smooth_edge_bump(data_fix, r_out=83, thre=6447)

data_fix -= -434.893 # np.median(data_fix[0:35,0:35])  # subtract a "background" value; or -np.median(data_rehalf_sized[0:10,0:10])

plt.matshow(np.log10(data_fix.T),origin='lower')
plt.colorbar()
plt.show()

pyfits.PrimaryHDU(data_fix.T).writeto('fix/{0}_fix.fits'.format(name),clobber=True)

from shutil import copyfile

예제 #2

name = 'M51'
hdu = pyfits.open('resized_source/{0}_resized.fits'.format(name))
data = hdu[0].data.T
hdu.close()

plt.matshow(np.log10(data.T), origin='lower')
plt.colorbar()
plt.show()

data_fix = data
#==============================================================================
# if changebox
#==============================================================================
data_fix = changebox(data_fix,
                     center=[335, 306],
                     half_size=[5, 5],
                     shift_center=[350, 289])
data_fix = changebox(data_fix,
                     center=[261, 138],
                     half_size=[9, 8],
                     shift_center=[256, 125])
data_fix = changebox(data_fix,
                     center=[166, 95],
                     half_size=[9, 8],
                     shift_center=[150, 101])
data_fix = changebox(data_fix,
                     center=[353, 38],
                     half_size=[10, 10],
                     shift_center=[360, 54])

#==============================================================================

예제 #3

파일: cosmetology_NGC2397.py 프로젝트: dartoon/my_code

sys.path.insert(0, '../../evil_code/')
from changebox import changebox
from smooth_edge_bump import smooth_edge_bump

name = 'NGC2397'
hdu = pyfits.open('{0}_resized.fits'.format(name))
data = hdu[0].data.T
hdu.close()

plt.matshow(np.log10(data.T), origin='lower')
plt.colorbar()
plt.show()

data_fix = data
data_fix = changebox(data_fix,
                     center=[95, 346],
                     half_size=[7 / 2, 7 / 2],
                     shift_center=[75, 336])
data_fix = changebox(data_fix,
                     center=[112, 272],
                     half_size=[9 / 2, 9 / 2],
                     shift_center=[119, 279])
data_fix = changebox(data_fix,
                     center=[159, 16],
                     half_size=[10 / 2, 8 / 2],
                     shift_center=[181, 22])

data_fix = smooth_edge_bump(data_fix, r_out=27, thre=2.45)
data_fix = smooth_edge_bump(data_fix, r_out=59, thre=1.8)
data_fix = smooth_edge_bump(data_fix, r_out=94, thre=1.12)
data_fix = smooth_edge_bump(data_fix, r_out=118, thre=0.3)
data_fix = smooth_edge_bump(data_fix, r_out=126, thre=0.2)

예제 #4

파일: cosmetology_PGC55493.py 프로젝트: dartoon/my_code

name = 'PGC55493'
hdu = pyfits.open('{0}_resized.fits'.format(name))
data = hdu[0].data.T
hdu.close()

plt.matshow(np.log10(data.T), origin='lower')
plt.colorbar()
plt.show()

data_fix = data
#==============================================================================
# if changebox
#==============================================================================
data_fix = changebox(data_fix,
                     center=[94, 317],
                     half_size=[8 / 2, 7 / 2],
                     shift_center=[83, 310])
#data_fix=changebox(data_fix,center=[112,272], half_size=[9/2, 9/2], shift_center=[119, 279])
#data_fix=changebox(data_fix,center=[159,16], half_size=[10/2, 8/2], shift_center=[181,22])

#==============================================================================
# if need smooth
#==============================================================================
data_fix = smooth_edge_bump(data_fix, r_out=125, thre=0.15)

plt.matshow(np.log10(data_fix.T), origin='lower')
plt.colorbar()
plt.show()

data_fix -= 0.06  # np.median(data_fix[0:35,0:35])  # subtract a "background" value; or -np.median(data_rehalf_sized[0:10,0:10])
pyfits.PrimaryHDU(data_fix.T).writeto('fix/{0}_fix.fits'.format(name),

예제 #5

name = 'NGC6782'
hdu = pyfits.open('resized_source/{0}_resized.fits'.format(name))
data = hdu[0].data.T
hdu.close()

plt.matshow(np.log10(data.T), origin='lower')
plt.colorbar()
plt.show()

data_fix = data
#==============================================================================
# if changebox
#==============================================================================
data_fix = changebox(data_fix,
                     center=[76, 238],
                     half_size=[18, 18],
                     shift_center=[101, 279])
data_fix = changebox(data_fix,
                     center=[47, 276],
                     half_size=[18, 18],
                     shift_center=[78, 313])
data_fix = changebox(data_fix,
                     center=[260, 356],
                     half_size=[10, 9],
                     shift_center=[282, 328])
data_fix = changebox(data_fix,
                     center=[354, 249],
                     half_size=[9, 9],
                     shift_center=[363, 225])
data_fix = changebox(data_fix,
                     center=[178, 9],

예제 #6

name = 'NGC3370'
hdu = pyfits.open('resized_source/{0}_resized.fits'.format(name))
data = hdu[0].data.T
hdu.close()

plt.matshow(np.log10(data.T), origin='lower')
plt.colorbar()
plt.show()

data_fix = data
#==============================================================================
# if changebox
#==============================================================================
data_fix = changebox(data_fix,
                     center=[291, 96],
                     half_size=[5, 5],
                     shift_center=[278, 81])
data_fix = changebox(data_fix,
                     center=[348, 160],
                     half_size=[1, 1],
                     shift_center=[343, 154])
#==============================================================================
# if need smooth
#==============================================================================
data_fix = smooth_edge_bump(data_fix, r_out=59, thre=2.1)
data_fix = smooth_edge_bump(data_fix, r_out=75, thre=1.05)
data_fix = smooth_edge_bump(data_fix, r_out=98, thre=1.04)
data_fix = smooth_edge_bump(data_fix, r_out=124, thre=0.47)

data_fix -= 0.0965518  # np.median(data_fix[0:35,0:35])  # subtract a "background" value; or -np.median(data_rehalf_sized[0:10,0:10])

예제 #7

파일: cosmetology_NGC1309.py 프로젝트: dartoon/my_code

name = 'NGC1309'
hdu = pyfits.open('resized_source/{0}_resized.fits'.format(name))
data = hdu[0].data.T
hdu.close()

plt.matshow(np.log10(data.T), origin='lower')
plt.colorbar()
plt.show()

data_fix = data
#==============================================================================
# if changebox
#==============================================================================
data_fix = changebox(data_fix,
                     center=[153, 235],
                     half_size=[10, 9],
                     shift_center=[156, 76],
                     y_re=-1)
data_fix = changebox(data_fix,
                     center=[285, 86],
                     half_size=[11, 10],
                     shift_center=[272, 56])
data_fix = changebox(data_fix,
                     center=[283, 277],
                     half_size=[6, 6],
                     shift_center=[259, 281])
data_fix = changebox(data_fix,
                     center=[291, 33],
                     half_size=[8, 9],
                     shift_center=[257, 18])
#==============================================================================

예제 #8

파일: cosmetology_NGC1376.py 프로젝트: dartoon/my_code

name = 'NGC1376'
hdu = pyfits.open('resized_source/{0}_resized.fits'.format(name))
data = hdu[0].data.T
hdu.close()

plt.matshow(np.log10(data.T), origin='lower')
plt.colorbar()
plt.show()

data_fix = data
#==============================================================================
# if changebox
#==============================================================================
data_fix = changebox(data_fix,
                     center=[242, 255],
                     half_size=[6, 9],
                     shift_center=[258, 243])
data_fix = changebox(data_fix,
                     center=[42, 252],
                     half_size=[6, 7],
                     shift_center=[29, 235])
data_fix = changebox(data_fix,
                     center=[294, 129],
                     half_size=[5, 8],
                     shift_center=[293, 105])
data_fix = changebox(data_fix,
                     center=[281, 204],
                     half_size=[1, 1],
                     shift_center=[283, 200])
data_fix = changebox(data_fix,
                     center=[35, 275],

예제 #9

name = 'NGC5584'
hdu = pyfits.open('resized_source/{0}_resized.fits'.format(name))
data = hdu[0].data.T
hdu.close()

plt.matshow(np.log10(data.T), origin='lower')
plt.colorbar()
plt.show()

data_fix = data
#==============================================================================
# if changebox
#==============================================================================
data_fix = changebox(data_fix,
                     center=[282, 158],
                     half_size=[6, 9],
                     shift_center=[281, 137])
data_fix = changebox(data_fix,
                     center=[235, 8],
                     half_size=[4, 4],
                     shift_center=[238, 18])
data_fix = changebox(data_fix,
                     center=[170, 278],
                     half_size=[6, 5],
                     shift_center=[154, 277])
data_fix = changebox(data_fix,
                     center=[95, 187],
                     half_size=[3, 4],
                     shift_center=[103, 186])
data_fix = changebox(data_fix,
                     center=[159, 241],

예제 #10

파일: cosmetology_NGC278.py 프로젝트: dartoon/my_code

sys.path.insert(0, '../../evil_code/')
from changebox import changebox
from smooth_edge_bump import smooth_edge_bump

name = 'NGC278'
hdu = pyfits.open('{0}_resized.fits'.format(name))
data = hdu[0].data.T
hdu.close()

plt.matshow(np.log10(data.T), origin='lower')
plt.colorbar()
plt.show()

data_fix = data
data_fix = changebox(data_fix,
                     center=[102, 228],
                     half_size=[5 / 2, 5 / 2],
                     shift_center=[117, 235])
data_fix = changebox(data_fix,
                     center=[233, 201],
                     half_size=[9 / 2, 9 / 2],
                     shift_center=[226, 214])
#
data_fix = smooth_edge_bump(data_fix, r_out=85, thre=0.9)
data_fix = smooth_edge_bump(data_fix, r_out=115, thre=0.4)

plt.matshow(np.log10(data_fix.T), origin='lower')
plt.colorbar()
plt.show()

data_fix -= 0.0620392  # np.median(data_fix[0:35,0:35])  # subtract a "background" value; or -np.median(data_rehalf_sized[0:10,0:10])
pyfits.PrimaryHDU(data_fix.T).writeto('fix/{0}_fix.fits'.format(name),

예제 #11

파일: cosmetology_NGC1084.py 프로젝트: dartoon/my_code

sys.path.insert(0, '../../evil_code/')
from changebox import changebox
from smooth_edge_bump import smooth_edge_bump

name = 'NGC1084'
hdu = pyfits.open('{0}_resized.fits'.format(name))
data = hdu[0].data.T
hdu.close()

plt.matshow(np.log10(data.T), origin='lower')
plt.colorbar()
plt.show()

data_fix = data
data_fix = changebox(data_fix,
                     center=[72, 205],
                     half_size=[6 / 2, 6 / 2],
                     shift_center=[79, 198])
data_fix = changebox(data_fix,
                     center=[81, 167],
                     half_size=[6 / 2, 6 / 2],
                     shift_center=[59, 158])
data_fix = changebox(data_fix,
                     center=[335, 316],
                     half_size=[8 / 2, 7 / 2],
                     shift_center=[345, 306])
data_fix = changebox(data_fix,
                     center=[380, 255],
                     half_size=[6 / 2, 8 / 2],
                     shift_center=[382, 242])
data_fix = changebox(data_fix,
                     center=[99, 112],

예제 #12

import sys
sys.path.insert(0,'../../evil_code/')
from changebox import changebox
from smooth_edge_bump import smooth_edge_bump
            
name='NGC4639'
hdu = pyfits.open('{0}_resized.fits'.format(name))
data = hdu[0].data.T
hdu.close()

plt.matshow(np.log10(data.T),origin='lower')
plt.colorbar()
plt.show()

data_fix=data
data_fix=changebox(data_fix,center=[5,51], half_size=[10/2, 26/2], shift_center=[6, 89])
#data_fix=changebox(data_fix,center=[112,272], half_size=[9/2, 9/2], shift_center=[119, 279])
#data_fix=changebox(data_fix,center=[159,16], half_size=[10/2, 8/2], shift_center=[181,22])

data_fix=smooth_edge_bump(data_fix, r_out=143, thre=0.22)

plt.matshow(np.log10(data_fix.T),origin='lower')
plt.colorbar()
plt.show()

data_fix -= 0.0728352 # np.median(data_fix[0:35,0:35])  # subtract a "background" value; or -np.median(data_rehalf_sized[0:10,0:10])
pyfits.PrimaryHDU(data_fix.T).writeto('fix/{0}_fix.fits'.format(name),clobber=True)

from shutil import copyfile
copyfile('cosmetology.py', 'fix/cosmetology_{0}.py'.format(name))

예제 #13

파일: cosmetology_ESO498G5.py 프로젝트: dartoon/my_code

name = 'ESO498G5'
hdu = pyfits.open('{0}_resized.fits'.format(name))
data = hdu[0].data.T
hdu.close()

plt.matshow(np.log10(data.T), origin='lower')
plt.colorbar()
plt.show()

data_fix = data
#==============================================================================
# if changebox
#==============================================================================
data_fix = changebox(data_fix,
                     center=[99, 269],
                     half_size=[15 / 2, 17 / 2],
                     shift_center=[119, 291])
data_fix = changebox(data_fix,
                     center=[255, 298],
                     half_size=[22 / 2, 21 / 2],
                     shift_center=[276, 263])
data_fix = changebox(data_fix,
                     center=[197, 231],
                     half_size=[7 / 2, 6 / 2],
                     shift_center=[204, 225])
data_fix = changebox(data_fix,
                     center=[215, 219],
                     half_size=[7 / 2, 6 / 2],
                     shift_center=[218, 211])
data_fix = changebox(data_fix,
                     center=[243, 75],

예제 #14

파일: cosmetology_NGC6503.py 프로젝트: dartoon/my_code

name = 'NGC6503'
hdu = pyfits.open('resized_source/{0}_resized.fits'.format(name))
data = hdu[0].data.T
hdu.close()

plt.matshow(np.log10(data.T), origin='lower')
plt.colorbar()
plt.show()

data_fix = data
#==============================================================================
# if changebox
#==============================================================================
data_fix = changebox(data_fix,
                     center=[65, 250],
                     half_size=[7, 6],
                     shift_center=[83, 259])
data_fix = changebox(data_fix,
                     center=[106, 204],
                     half_size=[7, 6],
                     shift_center=[118, 209])
data_fix = changebox(data_fix,
                     center=[235, 225],
                     half_size=[3, 3],
                     shift_center=[244, 225])
data_fix = changebox(data_fix,
                     center=[181, 256],
                     half_size=[3, 3],
                     shift_center=[189, 258])
data_fix = changebox(data_fix,
                     center=[14, 95],