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 [ ]:
Content source: ffpenaloza/AstroExp
Similar notebooks: