Tarea 5

Luego de descargar las imágenes en los filtros F475W y F850LP del objeto VCC1316 (M87) se siguen los pasos de la primera tarea para generar el catálogo. De Sirianni et. al (2005) se obtiene la escala de WFC (0.05''/px) y los zeropoint (según la tabla 10) y se ejecuta Sextractor. Se corrige por apertura según la tabla 6 (2 pixeles de radio a escala de 0.05''/px corresponden a 0.1'') para longitudes de onda efectivas iguales a las longitudes de onda centrales. Se corrige por reddening según la tabla 14 y el valor E(B-V) de 0.023 de Côté et al. (2006), tabla 1.


In [23]:
from astropy.io import fits
import numpy as np

f475 = fits.open('hst_9401_02_acs_wfc_f475w_drz.fits')
f850 = fits.open('hst_9401_02_acs_wfc_f850lp_drz.fits')

f475[1].writeto('sci_f475w_m87.fits',clobber=True)
f475[2].writeto('invvar_f475w_m87.fits',clobber=True)

f850[1].writeto('sci_f850lp_m87.fits',clobber=True)
f850[2].writeto('invvar_f850lp_m87.fits',clobber=True)

f475w.close()
f850lp.close()

In [28]:
!sextractor sci_f475w_m87.fits -c f475w.sex
!sextractor sci_f850lp_m87.fits -c f850lp.sex


> 
----- SExtractor 2.19.5 started on 2016-11-24 at 18:32:27 with 1 thread

> Setting catalog parameters
> Reading detection filter
> Initializing catalog
> Looking for sci_f475w_m87.fits
----- Measuring from: sci_f475w_m87.fits [1/1]
      "Unnamed" / no ext. header / 4238x4213 / 32 bits (floats)
> Looking for invvar_f475w_m87.fits
----- Weighting from: invvar_f475w_m87.fits [1/1]
      "Unnamed" / no ext. header / 4238x4213 / 32 bits (floats)
Detection+Measurement image: > Setting up background maps
> Setting up background map at line:   64
> Setting up background map at line:  128
> Setting up background map at line:  192
> Setting up background map at line:  256
> Setting up background map at line:  320
> Setting up background map at line:  384
> Setting up background map at line:  448
> Setting up background map at line:  512
> Setting up background map at line:  576
> Setting up background map at line:  640
> Setting up background map at line:  704
> Setting up background map at line:  768
> Setting up background map at line:  832
> Setting up background map at line:  896
> Setting up background map at line:  960
> Setting up background map at line: 1024
> Setting up background map at line: 1088
> Setting up background map at line: 1152
> Setting up background map at line: 1216
> Setting up background map at line: 1280
> Setting up background map at line: 1344
> Setting up background map at line: 1408
> Setting up background map at line: 1472
> Setting up background map at line: 1536
> Setting up background map at line: 1600
> Setting up background map at line: 1664
> Setting up background map at line: 1728
> Setting up background map at line: 1792
> Setting up background map at line: 1856
> Setting up background map at line: 1920
> Setting up background map at line: 1984
> Setting up background map at line: 2048
> Setting up background map at line: 2112
> Setting up background map at line: 2176
> Setting up background map at line: 2240
> Setting up background map at line: 2304
> Setting up background map at line: 2368
> Setting up background map at line: 2432
> Setting up background map at line: 2496
> Setting up background map at line: 2560
> Setting up background map at line: 2624
> Setting up background map at line: 2688
> Setting up background map at line: 2752
> Setting up background map at line: 2816
> Setting up background map at line: 2880
> Setting up background map at line: 2944
> Setting up background map at line: 3008
> Setting up background map at line: 3072
> Setting up background map at line: 3136
> Setting up background map at line: 3200
> Setting up background map at line: 3264
> Setting up background map at line: 3328
> Setting up background map at line: 3392
> Setting up background map at line: 3456
> Setting up background map at line: 3520
> Setting up background map at line: 3584
> Setting up background map at line: 3648
> Setting up background map at line: 3712
> Setting up background map at line: 3776
> Setting up background map at line: 3840
> Setting up background map at line: 3904
> Setting up background map at line: 3968
> Setting up background map at line: 4032
> Setting up background map at line: 4096
> Setting up background map at line: 4160
> Filtering background map(s)
> Computing background d-map
> Computing background-noise d-map
(M+D) Background: 0.167588   RMS: 0.0131859  / Threshold: 0.0395578  
> Scanning image
> Line:   25  Objects:        0 detected /        0 sextracted
> Line:   50  Objects:        3 detected /        0 sextracted
> Line:   75  Objects:        6 detected /        0 sextracted
> Line:  100  Objects:       13 detected /        0 sextracted
> Line:  125  Objects:       28 detected /        0 sextracted
> Line:  150  Objects:       39 detected /        0 sextracted
> Line:  175  Objects:       48 detected /        0 sextracted
> Line:  200  Objects:       60 detected /        0 sextracted
> Line:  225  Objects:       69 detected /        0 sextracted
> Line:  250  Objects:       74 detected /        0 sextracted
> Line:  275  Objects:       82 detected /        0 sextracted
> Line:  300  Objects:       96 detected /        0 sextracted
> Line:  325  Objects:      106 detected /        0 sextracted
> Line:  350  Objects:      120 detected /        0 sextracted
> Line:  375  Objects:      137 detected /        0 sextracted
> Line:  400  Objects:      152 detected /        0 sextracted
> Line:  425  Objects:      165 detected /        0 sextracted
> Line:  450  Objects:      177 detected /        0 sextracted
> Line:  475  Objects:      191 detected /        0 sextracted
> Line:  500  Objects:      205 detected /        0 sextracted
> Line:  525  Objects:      215 detected /        0 sextracted
> Line:  550  Objects:      225 detected /        0 sextracted
> Line:  575  Objects:      241 detected /        0 sextracted
> Line:  600  Objects:      267 detected /        0 sextracted
> Line:  625  Objects:      279 detected /        0 sextracted
> Line:  650  Objects:      290 detected /        0 sextracted
> Line:  675  Objects:      303 detected /        0 sextracted
> Line:  700  Objects:      314 detected /        0 sextracted
> Line:  725  Objects:      325 detected /        0 sextracted
> Line:  750  Objects:      338 detected /        0 sextracted
> Line:  775  Objects:      347 detected /        0 sextracted
> Line:  800  Objects:      355 detected /        0 sextracted
> Line:  825  Objects:      365 detected /        0 sextracted
> Line:  850  Objects:      381 detected /        0 sextracted
> Line:  875  Objects:      390 detected /        0 sextracted
> Line:  900  Objects:      400 detected /        0 sextracted
> Line:  925  Objects:      415 detected /        0 sextracted
> Line:  950  Objects:      426 detected /        0 sextracted
> Line:  975  Objects:      443 detected /        0 sextracted
> Line: 1000  Objects:      457 detected /        0 sextracted
> Line: 1025  Objects:      473 detected /        0 sextracted
> Line: 1026  Objects:      473 detected /        0 sextracted
> Line: 1050  Objects:      496 detected /        4 sextracted
> Line: 1075  Objects:      508 detected /       10 sextracted
> Line: 1100  Objects:      523 detected /       19 sextracted
> Line: 1125  Objects:      533 detected /       30 sextracted
> Line: 1150  Objects:      540 detected /       38 sextracted
> Line: 1175  Objects:      563 detected /       50 sextracted
> Line: 1200  Objects:      581 detected /       61 sextracted
> Line: 1225  Objects:      598 detected /       69 sextracted
> Line: 1250  Objects:      620 detected /       74 sextracted
> Line: 1275  Objects:      639 detected /       81 sextracted
> Line: 1300  Objects:      661 detected /       88 sextracted
> Line: 1325  Objects:      683 detected /       95 sextracted
> Line: 1350  Objects:      699 detected /      106 sextracted
> Line: 1375  Objects:      717 detected /      112 sextracted
> Line: 1400  Objects:      743 detected /      124 sextracted
> Line: 1425  Objects:      759 detected /      135 sextracted
> Line: 1450  Objects:      788 detected /      147 sextracted
> Line: 1475  Objects:      816 detected /      163 sextracted
> Line: 1500  Objects:      839 detected /      174 sextracted
> Line: 1525  Objects:      857 detected /      185 sextracted
> Line: 1550  Objects:      878 detected /      195 sextracted
> Line: 1575  Objects:      897 detected /      214 sextracted
> Line: 1600  Objects:      921 detected /      236 sextracted
> Line: 1625  Objects:      945 detected /      248 sextracted
> Line: 1650  Objects:      964 detected /      257 sextracted
> Line: 1675  Objects:      995 detected /      270 sextracted
> Line: 1700  Objects:     1027 detected /      281 sextracted
> Line: 1725  Objects:     1062 detected /      290 sextracted
> Line: 1750  Objects:     1112 detected /      302 sextracted
> Line: 1775  Objects:     1159 detected /      311 sextracted
> Line: 1800  Objects:     1245 detected /      321 sextracted
> Line: 1825  Objects:     1382 detected /      331 sextracted
> Line: 1850  Objects:     1559 detected /      342 sextracted
> Line: 1875  Objects:     1705 detected /      352 sextracted
> Line: 1900  Objects:     1815 detected /      363 sextracted
> Line: 1925  Objects:     1895 detected /      373 sextracted
> Line: 1950  Objects:     1972 detected /      386 sextracted
> Line: 1968  Objects:     2031 detected /      400 sextracted
> Line: 1975  Objects:     2045 detected /      405 sextracted
> Line: 2000  Objects:     2126 detected /      420 sextracted
> Line: 2025  Objects:     2212 detected /      436 sextracted
> Line: 2050  Objects:     2289 detected /      450 sextracted
> Line: 2075  Objects:     2391 detected /      464 sextracted
> Line: 2100  Objects:     2471 detected /      476 sextracted
> Line: 2125  Objects:     2547 detected /      488 sextracted
> Line: 2150  Objects:     2659 detected /      496 sextracted
> Line: 2175  Objects:     2735 detected /      518 sextracted
> Line: 2200  Objects:     2794 detected /      537 sextracted
> Line: 2225  Objects:     2848 detected /      555 sextracted
> Line: 2250  Objects:     2895 detected /      574 sextracted
> Line: 2275  Objects:     2933 detected /      593 sextracted
> Line: 2300  Objects:     2986 detected /      611 sextracted
> Line: 2325  Objects:     3051 detected /      633 sextracted
> Line: 2350  Objects:     3154 detected /      654 sextracted
> Line: 2375  Objects:     3303 detected /      672 sextracted
> Line: 2400  Objects:     3420 detected /      692 sextracted
> Line: 2425  Objects:     3519 detected /      710 sextracted
> Line: 2450  Objects:     3586 detected /      735 sextracted
> Line: 2475  Objects:     3637 detected /      759 sextracted
> Line: 2500  Objects:     3680 detected /      784 sextracted
> Line: 2525  Objects:     3720 detected /      800 sextracted
> Line: 2527  Objects:     3723 detected /      800 sextracted
> Line: 2550  Objects:     3779 detected /      821 sextracted
> Line: 2575  Objects:     3812 detected /      835 sextracted
> Line: 2600  Objects:     3840 detected /      860 sextracted
> Line: 2625  Objects:     3864 detected /      880 sextracted
> Line: 2650  Objects:     3889 detected /      898 sextracted
> Line: 2675  Objects:     3906 detected /      923 sextracted
> Line: 2700  Objects:     3927 detected /      952 sextracted
> Line: 2725  Objects:     3953 detected /      977 sextracted
> Line: 2750  Objects:     3972 detected /     1006 sextracted
> Line: 2775  Objects:     3997 detected /     1040 sextracted
> Line: 2800  Objects:     4021 detected /     1093 sextracted
> Line: 2825  Objects:     4042 detected /     1167 sextracted
> Line: 2840  Objects:     4051 detected /     1200 sextracted
> Line: 2850  Objects:     4066 detected /     1238 sextracted
> Line: 2875  Objects:     4094 detected /     1302 sextracted
> Line: 2900  Objects:     4116 detected /     1357 sextracted
> Line: 2925  Objects:     4140 detected /     1398 sextracted
> Line: 2950  Objects:     4159 detected /     1446 sextracted
> Line: 2975  Objects:     4182 detected /     1480 sextracted
> Line: 3000  Objects:     4200 detected /     1517 sextracted
> Line: 3025  Objects:     4218 detected /     1546 sextracted
> Line: 3050  Objects:     4233 detected /     1574 sextracted
> Line: 3075  Objects:     4251 detected /     1593 sextracted
> Line: 3083  Objects:     4258 detected /     1600 sextracted
> Line: 3100  Objects:     4277 detected /     1617 sextracted
> Line: 3125  Objects:     4286 detected /     1646 sextracted
> Line: 3150  Objects:     4306 detected /     1679 sextracted
> Line: 3175  Objects:     4327 detected /     1702 sextracted
> Line: 3200  Objects:     4345 detected /     1730 sextracted
> Line: 3225  Objects:     4365 detected /     1753 sextracted
> Line: 3250  Objects:     4381 detected /     1776 sextracted
> Line: 3275  Objects:     4398 detected /     1797 sextracted
> Line: 3300  Objects:     4411 detected /     1823 sextracted
> Line: 3325  Objects:     4425 detected /     1860 sextracted
> Line: 3350  Objects:     4443 detected /     1908 sextracted
> Line: 3375  Objects:     4459 detected /     1963 sextracted
> Line: 3396  Objects:     4476 detected /     2000 sextracted
> Line: 3400  Objects:     4479 detected /     2013 sextracted
> Line: 3425  Objects:     4493 detected /     2091 sextracted
> Line: 3450  Objects:     4510 detected /     2133 sextracted
> Line: 3475  Objects:     4523 detected /     2173 sextracted
> Line: 3500  Objects:     4544 detected /     2212 sextracted
> Line: 3525  Objects:     4557 detected /     2247 sextracted
> Line: 3550  Objects:     4569 detected /     2280 sextracted
> Line: 3575  Objects:     4581 detected /     2321 sextracted
> Line: 3600  Objects:     4595 detected /     2348 sextracted
> Line: 3625  Objects:     4606 detected /     2371 sextracted
> Line: 3650  Objects:     4618 detected /     2392 sextracted
> Line: 3667  Objects:     4625 detected /     2400 sextracted
> Line: 3675  Objects:     4633 detected /     2407 sextracted
> Line: 3700  Objects:     4645 detected /     2426 sextracted
> Line: 3725  Objects:     4665 detected /     2449 sextracted
> Line: 3750  Objects:     4684 detected /     2473 sextracted
> Line: 3775  Objects:     4695 detected /     2498 sextracted
> Line: 3800  Objects:     4700 detected /     2516 sextracted
> Line: 3825  Objects:     4713 detected /     2533 sextracted
> Line: 3850  Objects:     4727 detected /     2559 sextracted
> Line: 3875  Objects:     4737 detected /     2589 sextracted
> Line: 3900  Objects:     4747 detected /     2615 sextracted
> Line: 3925  Objects:     4763 detected /     2630 sextracted
> Line: 3950  Objects:     4771 detected /     2656 sextracted
> Line: 3975  Objects:     4785 detected /     2678 sextracted
> Line: 4000  Objects:     4796 detected /     2690 sextracted
> Line: 4025  Objects:     4809 detected /     2710 sextracted
> Line: 4050  Objects:     4813 detected /     2725 sextracted
> Line: 4075  Objects:     4819 detected /     2745 sextracted
> Line: 4100  Objects:     4823 detected /     2762 sextracted
> Line: 4125  Objects:     4827 detected /     2775 sextracted
> Line: 4150  Objects:     4829 detected /     2797 sextracted
> Line: 4156  Objects:     4830 detected /     2800 sextracted
> Line: 4175  Objects:     4830 detected /     2814 sextracted
> Line: 4200  Objects:     4833 detected /     2834 sextracted
> Line: 4213  Objects:     4833 detected /     3200 sextracted
      Objects: detected 4833     / sextracted 3306            

> Closing files
> 
> All done (in 2.6 s: 1646.3 lines/s , 1291.8 detections/s)
> 
----- SExtractor 2.19.5 started on 2016-11-24 at 18:32:30 with 1 thread

> Setting catalog parameters
> Reading detection filter
> Initializing catalog
> Looking for sci_f850lp_m87.fits
----- Measuring from: sci_f850lp_m87.fits [1/1]
      "Unnamed" / no ext. header / 4238x4213 / 32 bits (floats)
> Looking for invvar_f850lp_m87.fits
----- Weighting from: invvar_f850lp_m87.fits [1/1]
      "Unnamed" / no ext. header / 4238x4213 / 32 bits (floats)
Detection+Measurement image: > Setting up background maps
> Setting up background map at line:   64
> Setting up background map at line:  128
> Setting up background map at line:  192
> Setting up background map at line:  256
> Setting up background map at line:  320
> Setting up background map at line:  384
> Setting up background map at line:  448
> Setting up background map at line:  512
> Setting up background map at line:  576
> Setting up background map at line:  640
> Setting up background map at line:  704
> Setting up background map at line:  768
> Setting up background map at line:  832
> Setting up background map at line:  896
> Setting up background map at line:  960
> Setting up background map at line: 1024
> Setting up background map at line: 1088
> Setting up background map at line: 1152
> Setting up background map at line: 1216
> Setting up background map at line: 1280
> Setting up background map at line: 1344
> Setting up background map at line: 1408
> Setting up background map at line: 1472
> Setting up background map at line: 1536
> Setting up background map at line: 1600
> Setting up background map at line: 1664
> Setting up background map at line: 1728
> Setting up background map at line: 1792
> Setting up background map at line: 1856
> Setting up background map at line: 1920
> Setting up background map at line: 1984
> Setting up background map at line: 2048
> Setting up background map at line: 2112
> Setting up background map at line: 2176
> Setting up background map at line: 2240
> Setting up background map at line: 2304
> Setting up background map at line: 2368
> Setting up background map at line: 2432
> Setting up background map at line: 2496
> Setting up background map at line: 2560
> Setting up background map at line: 2624
> Setting up background map at line: 2688
> Setting up background map at line: 2752
> Setting up background map at line: 2816
> Setting up background map at line: 2880
> Setting up background map at line: 2944
> Setting up background map at line: 3008
> Setting up background map at line: 3072
> Setting up background map at line: 3136
> Setting up background map at line: 3200
> Setting up background map at line: 3264
> Setting up background map at line: 3328
> Setting up background map at line: 3392
> Setting up background map at line: 3456
> Setting up background map at line: 3520
> Setting up background map at line: 3584
> Setting up background map at line: 3648
> Setting up background map at line: 3712
> Setting up background map at line: 3776
> Setting up background map at line: 3840
> Setting up background map at line: 3904
> Setting up background map at line: 3968
> Setting up background map at line: 4032
> Setting up background map at line: 4096
> Setting up background map at line: 4160
> Filtering background map(s)
> Computing background d-map
> Computing background-noise d-map
(M+D) Background: 0.175354   RMS: 0.0133454  / Threshold: 0.0400361  
> Scanning image
> Line:   25  Objects:        0 detected /        0 sextracted
> Line:   50  Objects:        2 detected /        0 sextracted
> Line:   75  Objects:        4 detected /        0 sextracted
> Line:  100  Objects:        8 detected /        0 sextracted
> Line:  125  Objects:       18 detected /        0 sextracted
> Line:  150  Objects:       27 detected /        0 sextracted
> Line:  175  Objects:       37 detected /        0 sextracted
> Line:  200  Objects:       49 detected /        0 sextracted
> Line:  225  Objects:       56 detected /        0 sextracted
> Line:  250  Objects:       60 detected /        0 sextracted
> Line:  275  Objects:       71 detected /        0 sextracted
> Line:  300  Objects:       81 detected /        0 sextracted
> Line:  325  Objects:       89 detected /        0 sextracted
> Line:  350  Objects:      101 detected /        0 sextracted
> Line:  375  Objects:      112 detected /        0 sextracted
> Line:  400  Objects:      122 detected /        0 sextracted
> Line:  425  Objects:      137 detected /        0 sextracted
> Line:  450  Objects:      145 detected /        0 sextracted
> Line:  475  Objects:      157 detected /        0 sextracted
> Line:  500  Objects:      166 detected /        0 sextracted
> Line:  525  Objects:      174 detected /        0 sextracted
> Line:  550  Objects:      181 detected /        0 sextracted
> Line:  575  Objects:      195 detected /        0 sextracted
> Line:  600  Objects:      214 detected /        0 sextracted
> Line:  625  Objects:      224 detected /        0 sextracted
> Line:  650  Objects:      233 detected /        0 sextracted
> Line:  675  Objects:      241 detected /        0 sextracted
> Line:  700  Objects:      253 detected /        0 sextracted
> Line:  725  Objects:      264 detected /        0 sextracted
> Line:  750  Objects:      272 detected /        0 sextracted
> Line:  775  Objects:      282 detected /        0 sextracted
> Line:  800  Objects:      292 detected /        0 sextracted
> Line:  825  Objects:      299 detected /        0 sextracted
> Line:  850  Objects:      311 detected /        0 sextracted
> Line:  875  Objects:      321 detected /        0 sextracted
> Line:  900  Objects:      328 detected /        0 sextracted
> Line:  925  Objects:      338 detected /        0 sextracted
> Line:  950  Objects:      348 detected /        0 sextracted
> Line:  975  Objects:      360 detected /        0 sextracted
> Line: 1000  Objects:      372 detected /        0 sextracted
> Line: 1022  Objects:      387 detected /        0 sextracted
> Line: 1025  Objects:      390 detected /        3 sextracted
> Line: 1050  Objects:      405 detected /        4 sextracted
> Line: 1075  Objects:      414 detected /        9 sextracted
> Line: 1100  Objects:      428 detected /       14 sextracted
> Line: 1125  Objects:      436 detected /       22 sextracted
> Line: 1150  Objects:      445 detected /       27 sextracted
> Line: 1175  Objects:      462 detected /       38 sextracted
> Line: 1200  Objects:      485 detected /       48 sextracted
> Line: 1225  Objects:      502 detected /       54 sextracted
> Line: 1250  Objects:      528 detected /       60 sextracted
> Line: 1275  Objects:      553 detected /       67 sextracted
> Line: 1300  Objects:      586 detected /       69 sextracted
> Line: 1325  Objects:      603 detected /       75 sextracted
> Line: 1350  Objects:      623 detected /       86 sextracted
> Line: 1375  Objects:      644 detected /       91 sextracted
> Line: 1400  Objects:      676 detected /       99 sextracted
> Line: 1425  Objects:      694 detected /      107 sextracted
> Line: 1450  Objects:      726 detected /      118 sextracted
> Line: 1475  Objects:      755 detected /      130 sextracted
> Line: 1500  Objects:      798 detected /      137 sextracted
> Line: 1525  Objects:      835 detected /      146 sextracted
> Line: 1550  Objects:      888 detected /      155 sextracted
> Line: 1575  Objects:      943 detected /      172 sextracted
> Line: 1600  Objects:      991 detected /      189 sextracted
> Line: 1625  Objects:     1050 detected /      195 sextracted
> Line: 1650  Objects:     1114 detected /      203 sextracted
> Line: 1675  Objects:     1178 detected /      212 sextracted
> Line: 1700  Objects:     1267 detected /      225 sextracted
> Line: 1725  Objects:     1386 detected /      236 sextracted
> Line: 1750  Objects:     1539 detected /      246 sextracted
> Line: 1775  Objects:     1714 detected /      254 sextracted
> Line: 1800  Objects:     1889 detected /      266 sextracted
> Line: 1825  Objects:     2035 detected /      274 sextracted
> Line: 1850  Objects:     2166 detected /      283 sextracted
> Line: 1875  Objects:     2282 detected /      292 sextracted
> Line: 1900  Objects:     2380 detected /      299 sextracted
> Line: 1925  Objects:     2497 detected /      306 sextracted
> Line: 1950  Objects:     2583 detected /      316 sextracted
> Line: 1975  Objects:     2682 detected /      329 sextracted
> Line: 2000  Objects:     2790 detected /      342 sextracted
> Line: 2025  Objects:     2889 detected /      357 sextracted
> Line: 2050  Objects:     3009 detected /      372 sextracted
> Line: 2075  Objects:     3100 detected /      381 sextracted
> Line: 2100  Objects:     3214 detected /      388 sextracted
> Line: 2125  Objects:     3350 detected /      399 sextracted
> Line: 2132  Objects:     3383 detected /      400 sextracted
> Line: 2150  Objects:     3486 detected /      411 sextracted
> Line: 2175  Objects:     3595 detected /      427 sextracted
> Line: 2200  Objects:     3703 detected /      451 sextracted
> Line: 2225  Objects:     3810 detected /      466 sextracted
> Line: 2250  Objects:     3883 detected /      488 sextracted
> Line: 2275  Objects:     3958 detected /      507 sextracted
> Line: 2300  Objects:     4040 detected /      528 sextracted
> Line: 2325  Objects:     4136 detected /      554 sextracted
> Line: 2350  Objects:     4263 detected /      576 sextracted
> Line: 2375  Objects:     4407 detected /      595 sextracted
> Line: 2400  Objects:     4601 detected /      617 sextracted
> Line: 2425  Objects:     4812 detected /      636 sextracted
> Line: 2450  Objects:     5026 detected /      657 sextracted
> Line: 2475  Objects:     5192 detected /      681 sextracted
> Line: 2500  Objects:     5332 detected /      708 sextracted
> Line: 2525  Objects:     5432 detected /      737 sextracted
> Line: 2550  Objects:     5539 detected /      763 sextracted
> Line: 2575  Objects:     5626 detected /      786 sextracted
> Line: 2588  Objects:     5659 detected /      800 sextracted
> Line: 2600  Objects:     5694 detected /      814 sextracted
> Line: 2625  Objects:     5774 detected /      847 sextracted
> Line: 2650  Objects:     5833 detected /      873 sextracted
> Line: 2675  Objects:     5875 detected /      909 sextracted
> Line: 2700  Objects:     5923 detected /      947 sextracted
> Line: 2725  Objects:     5963 detected /      991 sextracted
> Line: 2750  Objects:     6004 detected /     1059 sextracted
> Line: 2775  Objects:     6048 detected /     1135 sextracted
> Line: 2797  Objects:     6076 detected /     1200 sextracted
> Line: 2800  Objects:     6086 detected /     1213 sextracted
> Line: 2825  Objects:     6106 detected /     1279 sextracted
> Line: 2850  Objects:     6136 detected /     1329 sextracted
> Line: 2875  Objects:     6154 detected /     1372 sextracted
> Line: 2900  Objects:     6173 detected /     1403 sextracted
> Line: 2925  Objects:     6197 detected /     1444 sextracted
> Line: 2950  Objects:     6215 detected /     1480 sextracted
> Line: 2975  Objects:     6237 detected /     1529 sextracted
> Line: 3000  Objects:     6257 detected /     1565 sextracted
> Line: 3019  Objects:     6269 detected /     1600 sextracted
> Line: 3025  Objects:     6273 detected /     1608 sextracted
> Line: 3050  Objects:     6288 detected /     1657 sextracted
> Line: 3075  Objects:     6302 detected /     1697 sextracted
> Line: 3100  Objects:     6321 detected /     1752 sextracted
> Line: 3125  Objects:     6332 detected /     1800 sextracted
> Line: 3150  Objects:     6351 detected /     1854 sextracted
> Line: 3175  Objects:     6368 detected /     1906 sextracted
> Line: 3200  Objects:     6381 detected /     1942 sextracted
> Line: 3225  Objects:     6399 detected /     1976 sextracted
> Line: 3247  Objects:     6411 detected /     2000 sextracted
> Line: 3250  Objects:     6414 detected /     2002 sextracted
> Line: 3275  Objects:     6428 detected /     2033 sextracted
> Line: 3300  Objects:     6440 detected /     2068 sextracted
> Line: 3325  Objects:     6452 detected /     2108 sextracted
> Line: 3350  Objects:     6465 detected /     2160 sextracted
> Line: 3375  Objects:     6481 detected /     2229 sextracted
> Line: 3400  Objects:     6497 detected /     2314 sextracted
> Line: 3417  Objects:     6503 detected /     2400 sextracted
> Line: 3425  Objects:     6505 detected /     2449 sextracted
> Line: 3450  Objects:     6520 detected /     2563 sextracted
> Line: 3475  Objects:     6532 detected /     2662 sextracted
> Line: 3500  Objects:     6545 detected /     2753 sextracted
> Line: 3518  Objects:     6552 detected /     2800 sextracted
> Line: 3525  Objects:     6555 detected /     2820 sextracted
> Line: 3550  Objects:     6564 detected /     2890 sextracted
> Line: 3575  Objects:     6573 detected /     2955 sextracted
> Line: 3600  Objects:     6581 detected /     3011 sextracted
> Line: 3625  Objects:     6592 detected /     3058 sextracted
> Line: 3650  Objects:     6601 detected /     3102 sextracted
> Line: 3675  Objects:     6613 detected /     3138 sextracted
> Line: 3700  Objects:     6622 detected /     3176 sextracted
> Line: 3716  Objects:     6631 detected /     3200 sextracted
> Line: 3725  Objects:     6632 detected /     3211 sextracted
> Line: 3750  Objects:     6650 detected /     3256 sextracted
> Line: 3775  Objects:     6661 detected /     3295 sextracted
> Line: 3800  Objects:     6667 detected /     3320 sextracted
> Line: 3825  Objects:     6678 detected /     3344 sextracted
> Line: 3850  Objects:     6688 detected /     3370 sextracted
> Line: 3875  Objects:     6694 detected /     3390 sextracted
> Line: 3900  Objects:     6699 detected /     3412 sextracted
> Line: 3925  Objects:     6711 detected /     3431 sextracted
> Line: 3950  Objects:     6714 detected /     3452 sextracted
> Line: 3975  Objects:     6720 detected /     3468 sextracted
> Line: 4000  Objects:     6728 detected /     3484 sextracted
> Line: 4025  Objects:     6738 detected /     3503 sextracted
> Line: 4050  Objects:     6744 detected /     3518 sextracted
> Line: 4075  Objects:     6748 detected /     3530 sextracted
> Line: 4100  Objects:     6752 detected /     3544 sextracted
> Line: 4125  Objects:     6756 detected /     3556 sextracted
> Line: 4150  Objects:     6757 detected /     3572 sextracted
> Line: 4175  Objects:     6759 detected /     3591 sextracted
> Line: 4192  Objects:     6759 detected /     3600 sextracted
> Line: 4200  Objects:     6760 detected /     3604 sextracted
      Objects: detected 6760     / sextracted 3971            

> Closing files
> 
> All done (in 2.9 s: 1437.4 lines/s , 1354.8 detections/s)

In [37]:
from astropy import units as u
from astropy.coordinates import SkyCoord

# Se cargan listas con RA y DEC para cada imagen
RA475 = np.loadtxt('f475w.cat',usecols=(3,))
DE475 = np.loadtxt('f475w.cat',usecols=(4,))

RA850 = np.loadtxt('f850lp.cat',usecols=(3,))
DE850 = np.loadtxt('f850lp.cat',usecols=(4,))

# Match por parte de astropy. El catalogo del filtro f850lp contiene mas objetos
catalog = SkyCoord(ra=RA475*u.degree, dec=DE475*u.degree)  
c = SkyCoord(ra=RA850*u.degree, dec=DE850*u.degree)  
idx = c.match_to_catalog_sky(catalog)

# Del catalogo f475w.cat se extraen las filas que indica el match
matches = list(idx[0])
f475w = np.loadtxt('f475w.cat')
f850lp = np.loadtxt('f850lp.cat')
out = []

j = 0
for i in matches:
	out.append(np.concatenate([f475w[i]- 0.454- 0.023*3.773,f850lp[j]- 0.584- 0.023*1.498]))
	j = j+1

# Salida a archivo
np.savetxt('m87_match_f475w_f850lp.cat',out,
	fmt='%d\t%.4f\t%.4f\t%.7f\t%.7f\t%d\t%.4f\t%.4f\t%.7f\t%.7f',
	header='f475wN\tf475wMAG\tf475wMAGERR\tf475wALPHA\tf475wDELTA\tf850lpN\tf850lpMAG\tf850lpMAGERR\tf850lpALPHA\tf814wDELTA')

In [39]:
import matplotlib
import matplotlib.pyplot as plt
%matplotlib inline
from astropy.io import ascii

tbl = ascii.read('m87_match_f475w_f850lp.cat')

plt.figure(figsize=(10,10))
plt.hist(tbl["f475wMAG"] - tbl["f850lpMAG"])
plt.xlabel("$g-z$", fontsize=20)
plt.ylabel("Counts", fontsize=20)

plt.show()
plt.close()

plt.figure(figsize=(10,10))
plt.hist(tbl["f475wMAG"], histtype = 'step', color='b')
plt.hist(tbl["f850lpMAG"], histtype = 'step', color='r')
#plt.xlim(-0.5, 1.5)
#plt.ylim(24,16)
plt.xlabel("$g, z$", fontsize=20)
plt.ylabel("Counts", fontsize=20)

plt.show()
plt.close()



In [ ]: