Skip to content

Commit

Permalink
Add files via upload
Browse files Browse the repository at this point in the history
  • Loading branch information
@piexpiex
piexpiex authored Dec 21, 2020
1 parent 4551d2b commit 1b94c5a
Show file tree
Hide file tree
Showing 2 changed files with 40 additions and 19 deletions.
29 changes: 20 additions & 9 deletions CFC_configuration/python_scripts/sex_analisis.py
View file
Original file line number Diff line number Diff line change
Expand Up @@ -487,6 +487,7 @@
delta_2_find_sources=catalog['DEdeg']
q_mode=catalog['q_mode']
class_sdss=catalog['class']
mode=catalog['mode']
if color=='u'or color=='U':
pmag=catalog['umag']
e_pmag=catalog['e_umag']
Expand Down Expand Up @@ -560,17 +561,20 @@
exit()

if sdss_key==0:
lista=[mag_sex,magerr_sex,ellongation,ellipticity,FWHM,pmag,e_pmag,SPREAD_VALUE,NUMBER_XMATCH,class_sdss,q_mode]
lista=[mag_sex,magerr_sex,ellongation,ellipticity,FWHM,pmag,e_pmag,SPREAD_VALUE,NUMBER_XMATCH,class_sdss,q_mode,mode]
if sdss_key==1:
lista=[mag_sex,magerr_sex,ellongation,ellipticity,FWHM,pmag,e_pmag,SPREAD_VALUE,NUMBER_XMATCH]

#search with 5 digits
#alpha_find_sources=np.around(alpha_find_sources,5)
#delta_find_sources=np.around(delta_find_sources,5)
Numbers,Numbers_ok=find_sources(alpha_find_sources,delta_find_sources,alpha_2_find_sources,delta_2_find_sources)

for k in range(len(lista)):
lista[k]=lista[k][Numbers]
lista[k]=lista[k][np.where(Numbers_ok==1)]
#bestmatch
#Numbers,Numbers_ok=find_sources(alpha_find_sources,delta_find_sources,alpha_2_find_sources,delta_2_find_sources)

#for k in range(len(lista)):
# lista[k]=lista[k][Numbers]
# lista[k]=lista[k][np.where(Numbers_ok==1)]

mag_sex=lista[0]
magerr_sex=lista[1]
Expand All @@ -585,6 +589,7 @@
if sdss_key==0:
class_sdss=lista[9]
q_mode=lista[10]
mode=lista[11]

if len(pmag)<6:
print('Not enough objects for the calibration')
Expand Down Expand Up @@ -629,6 +634,7 @@
if sdss_key==0:
class_sdss=lista[9]
q_mode=lista[10]
mode=lista[11]
#Morphology selection

median_FWHM=np.median(FWHM)
Expand Down Expand Up @@ -659,6 +665,7 @@
if sdss_key==0:
class_sdss=lista[9]
q_mode=lista[10]
mode=lista[11]

for j in range(len(NUMBER_XMATCH)):
source_flag[int(NUMBER_XMATCH[j]-1)]=5
Expand All @@ -671,9 +678,9 @@
if sdss_key==0:
for k in range(len(lista)):
try:
lista[k]=lista[k][np.where((class_sdss==6) & (q_mode=='+'))]
lista[k]=lista[k][np.where((class_sdss==6) & (q_mode=='+') & (mode==1.0))]
except:
lista[k]=lista[k][np.where((class_sdss==6) & (q_mode==1.0))]
lista[k]=lista[k][np.where((class_sdss==6) & (q_mode==1.0) & (mode==1.0))]

mag_sex=lista[0]
magerr_sex=lista[1]
Expand Down Expand Up @@ -827,8 +834,10 @@
c21 = fits.Column(name='Elongation',array=np.around(final_objects[:,ELONGATION],2), format='E')
c22 = fits.Column(name='Ellipticity',array=np.around(final_objects[:,ELLIPTICITY],2), format='E')
c23 = fits.Column(name='FWHM', unit='arcsec',array=np.around(3600*final_objects[:,FWHM_WORLD],2), format='E')
c24 = fits.Column(name='FLAGS',array=final_objects[:,FLAGS], format='E')
c25 = fits.Column(name='FLAGS_WEIGHT',array=final_objects[:,FLAGS_WEIGHT], format='E')

t = fits.BinTableHDU.from_columns([c1,c2,c3,c4,c5,c6,c7,c8,c9,c10,c13,c14,c15,c16,c17,c18,c19,c20,c21,c22,c23],name='catalog')
t = fits.BinTableHDU.from_columns([c1,c2,c3,c4,c5,c6,c7,c8,c9,c10,c13,c14,c15,c16,c17,c18,c19,c20,c21,c22,c23,c24,c25],name='catalog')
t.writeto('catalogs_folder/CFC_catalogs/'+fichero[0:len(fichero)-5]+'_catalog.fits',overwrite=True)
votable1=Table.read('catalogs_folder/CFC_catalogs/'+fichero[0:len(fichero)-5]+'_catalog.fits')
votable1.write('catalogs_folder/CFC_catalogs/'+fichero[0:len(fichero)-5]+'_catalog.xml',table_id='table_id',format='votable',overwrite=True)
Expand Down Expand Up @@ -886,8 +895,10 @@
c20 = fits.Column(name='Ellipticity',array=np.around(final_objects[:,ELLIPTICITY],2), format='E')
c21 = fits.Column(name='FWHM', unit='arcsec',array=np.around(3600*total_objects[:,FWHM_WORLD],2), format='E')
c22 = fits.Column(name='source_type',array=source_flag, format='E')
c23 = fits.Column(name='FLAGS',array=total_objects[:,FLAGS], format='E')
c24 = fits.Column(name='FLAGS_WEIGHT',array=final_objects[:,FLAGS_WEIGHT], format='E')

t = fits.BinTableHDU.from_columns([c1,c2,c3,c4,c5,c6,c7,c8,c9,c10,c13,c14,c15,c16,c17,c18,c19,c20,c21,c22],name='catalog')
t = fits.BinTableHDU.from_columns([c1,c2,c3,c4,c5,c6,c7,c8,c9,c10,c13,c14,c15,c16,c17,c18,c19,c20,c21,c22,c23,c24],name='catalog')
t.writeto('catalogs_folder/CFC_sources/'+fichero[0:len(fichero)-5]+'_sources.fits',overwrite=True)
votable2=Table.read('catalogs_folder/CFC_sources/'+fichero[0:len(fichero)-5]+'_sources.fits')
votable2.write('catalogs_folder/CFC_sources/'+fichero[0:len(fichero)-5]+'_sources.xml',table_id='table_id',format='votable',overwrite=True)
Expand Down
30 changes: 20 additions & 10 deletions CFC_configuration/python_scripts/sex_analisis_saved.py
View file
Original file line number Diff line number Diff line change
Expand Up @@ -483,6 +483,7 @@
delta_2_find_sources=catalog['DEdeg']
q_mode=catalog['q_mode']
class_sdss=catalog['class']
mode=catalog['mode']
if color=='u'or color=='U':
pmag=catalog['umag']
e_pmag=catalog['e_umag']
Expand Down Expand Up @@ -556,17 +557,20 @@
exit()

if sdss_key==0:
lista=[mag_sex,magerr_sex,ellongation,ellipticity,FWHM,pmag,e_pmag,SPREAD_VALUE,NUMBER_XMATCH,class_sdss,q_mode]
lista=[mag_sex,magerr_sex,ellongation,ellipticity,FWHM,pmag,e_pmag,SPREAD_VALUE,NUMBER_XMATCH,class_sdss,q_mode,mode]
if sdss_key==1:
lista=[mag_sex,magerr_sex,ellongation,ellipticity,FWHM,pmag,e_pmag,SPREAD_VALUE,NUMBER_XMATCH]

#search with 5 digits
#alpha_find_sources=np.around(alpha_find_sources,5)
#delta_find_sources=np.around(delta_find_sources,5)
Numbers,Numbers_ok=find_sources(alpha_find_sources,delta_find_sources,alpha_2_find_sources,delta_2_find_sources)

for k in range(len(lista)):
lista[k]=lista[k][Numbers]
lista[k]=lista[k][np.where(Numbers_ok==1)]
#bestmatch
#Numbers,Numbers_ok=find_sources(alpha_find_sources,delta_find_sources,alpha_2_find_sources,delta_2_find_sources)

#for k in range(len(lista)):
# lista[k]=lista[k][Numbers]
# lista[k]=lista[k][np.where(Numbers_ok==1)]

mag_sex=lista[0]
magerr_sex=lista[1]
Expand All @@ -580,6 +584,7 @@
if sdss_key==0:
class_sdss=lista[9]
q_mode=lista[10]
mode=lista[11]

if len(pmag)<6:
print('Not enough objects for the calibration')
Expand Down Expand Up @@ -624,6 +629,7 @@
if sdss_key==0:
class_sdss=lista[9]
q_mode=lista[10]
mode=lista[11]
#Morphology selection

median_FWHM=np.median(FWHM)
Expand Down Expand Up @@ -654,7 +660,7 @@
if sdss_key==0:
class_sdss=lista[9]
q_mode=lista[10]

mode=lista[11]
for j in range(len(NUMBER_XMATCH)):
source_flag[int(NUMBER_XMATCH[j]-1)]=5

Expand All @@ -666,9 +672,9 @@
if sdss_key==0:
for k in range(len(lista)):
try:
lista[k]=lista[k][np.where((class_sdss==6) & (q_mode=='+'))]
lista[k]=lista[k][np.where((class_sdss==6) & (q_mode=='+') & (mode==1.0))]
except:
lista[k]=lista[k][np.where((class_sdss==6) & (q_mode==1.0))]
lista[k]=lista[k][np.where((class_sdss==6) & (q_mode==1.0) & (mode==1.0))]

mag_sex=lista[0]
magerr_sex=lista[1]
Expand Down Expand Up @@ -823,8 +829,10 @@
c21 = fits.Column(name='Elongation',array=np.around(final_objects[:,ELONGATION],2), format='E')
c22 = fits.Column(name='Ellipticity',array=np.around(final_objects[:,ELLIPTICITY],2), format='E')
c23 = fits.Column(name='FWHM', unit='arcsec',array=np.around(3600*final_objects[:,FWHM_WORLD],2), format='E')
c24 = fits.Column(name='FLAGS',array=final_objects[:,FLAGS], format='E')
c25 = fits.Column(name='FLAGS_WEIGHT',array=final_objects[:,FLAGS_WEIGHT], format='E')

t = fits.BinTableHDU.from_columns([c1,c2,c3,c4,c5,c6,c7,c8,c9,c10,c13,c14,c15,c16,c17,c18,c19,c20,c21,c22,c23],name='catalog')
t = fits.BinTableHDU.from_columns([c1,c2,c3,c4,c5,c6,c7,c8,c9,c10,c13,c14,c15,c16,c17,c18,c19,c20,c21,c22,c23,c24,c25],name='catalog')
t.writeto('catalogs_folder/CFC_catalogs/'+fichero[0:len(fichero)-5]+'_catalog.fits',overwrite=True)
votable1=Table.read('catalogs_folder/CFC_catalogs/'+fichero[0:len(fichero)-5]+'_catalog.fits')
votable1.write('catalogs_folder/CFC_catalogs/'+fichero[0:len(fichero)-5]+'_catalog.xml',table_id='table_id',format='votable',overwrite=True)
Expand Down Expand Up @@ -883,8 +891,10 @@
c20 = fits.Column(name='Ellipticity',array=np.around(final_objects[:,ELLIPTICITY],2), format='E')
c21 = fits.Column(name='FWHM', unit='arcsec',array=np.around(3600*total_objects[:,FWHM_WORLD],2), format='E')
c22 = fits.Column(name='source_type',array=source_flag, format='E')
c23 = fits.Column(name='FLAGS',array=total_objects[:,FLAGS], format='E')
c24 = fits.Column(name='FLAGS_WEIGHT',array=final_objects[:,FLAGS_WEIGHT], format='E')

t = fits.BinTableHDU.from_columns([c1,c2,c3,c4,c5,c6,c7,c8,c9,c10,c13,c14,c15,c16,c17,c18,c19,c20,c21,c22],name='catalog')
t = fits.BinTableHDU.from_columns([c1,c2,c3,c4,c5,c6,c7,c8,c9,c10,c13,c14,c15,c16,c17,c18,c19,c20,c21,c22,c23,c24],name='catalog')
t.writeto('catalogs_folder/CFC_sources/'+fichero[0:len(fichero)-5]+'_sources.fits',overwrite=True)
votable2=Table.read('catalogs_folder/CFC_sources/'+fichero[0:len(fichero)-5]+'_sources.fits')
votable2.write('catalogs_folder/CFC_sources/'+fichero[0:len(fichero)-5]+'_sources.xml',table_id='table_id',format='votable',overwrite=True)
Expand Down

0 comments on commit 1b94c5a

Please sign in to comment.