Exemplo: Modo Campo Denso (Automático)
Lista típica de parâmetros neste modo:
Tomar cuidado que intera = "no" e reimag = "yes", neste modo.
É preciso calcular a distancia em X e Y (parâmetros shiftx e shifty) das imagens ordinárias e extraordinárias.
Usamos as coordenadas da padrão do exemplo anterior (file: daoedit6748a)
Os valores calculados são: shiftx = 4.0 ; shifty = 36.6
O parâmetro side = right, neste caso (imagem superior à direita da inferior).
Rodando a rotina:
pc> quick
sub-directory to create for temporary files (field): Confirmação de diretório onde será gravada a informação criada
Running DAOFIND ... Rodando daofind na imagem de referencia para achar a maior quantidade de objetos com perfil estelar
FWHM of features in scale units (2.5) (CR or value): 3.4 Usamos o valor calculado no exemplo anterior
New FWHM of features: 3.4 scale units 3.4 pixels
Standard deviation of background in counts (INDEF) (CR or value): 3.2 Usamos o valor calculado no exemplo anterior
New standard deviation of background: 3.2 counts
Detection threshold in sigma (4.) (CR or value): O valor default (4) é um bom compromisso
New detection threshold: 4. sigma 12.8 counts
Minimum good data value (INDEF) (CR or value): 50 Valor mínimo: 50 ~ ( sky + 6*sigma )
New minimum good data value: 50. counts
Maximum good data value (INDEF) (CR or value): 31000 Valor massimo: algo abaixo do nível de saturação
New maximum good data value: 31000. counts
Image: ../hd110cv0001.fits fwhmpsf: 3.4 ratio: 1. theta: 0. nsigma: 1.5 Procurando objetos
738.60 1.92 -1.305 INDEF -0.874 0.244 1
736.95 3.00 -1.979 0.411 0.396 0.016 2
739.25 3.74 -1.377 INDEF 0.730 0.105 3
178.73 13.28 -0.034 INDEF -0.950 0.001 4
...
...
981.11 1012.31 -1.305 INDEF -0.032 0.062 240
17.92 1014.73 -1.243 0.315 -0.424 -0.091 241
578.17 1011.67 -1.412 INDEF -0.089 0.048 242
572.84 1016.59 -2.933 0.424 0.625 -0.020 243
916.33 1017.86 -0.129 0.405 0.663 -0.082 244 Foram achados 244 objetos
threshold: 12.8 relerr: 1.444 0.2 <= sharp <= 1. -1. <= round <= 1.
Running ORDEM ... Rodando ordem para ordenar em pares
PAIR 1
887.873 13.852 5
892.01 50.055 10
PAIR 2
185.204 18.274 6
188.987 54.887 13
PAIR 3
794.188 51.2 11
797.924 87.841 18
PAIR 4
546.769 88.792 19
549.589 125.262 28
PAIR 5
252.469 95.958 21
256.323 132.555 30
...
...
PAIR 58
734.944 936.427 212
738.778 972.219 228
PAIR 59
144.766 953.725 219
148.718 991.234 235
PAIR 60
971.162 970.812 229
975.222 1006.889 239
PAIR 61
574.431 974.697 232
578.17 1011.669 242
PAIR 62
568.9 980.584 233
572.838 1016.593 243 Foram ordenados 62 pares (124 objetos)
Running DAOFIND for bright stars ... Rodando daofind (novamente) na imagem de referencia para achar objetos brilhantes para o alinhamento de sequencia de imagens
Input image (../hd110cv0001.fits):
FWHM of features in scale units (2.5) (CR or value): 3.4 Mesmo valor usado na vez anterior
New FWHM of features: 3.4 scale units 3.4 pixels
Standard deviation of background in counts (INDEF) (CR or value): 3.2 Mesmo valor usado na vez anterior
New standard deviation of background: 3.2 counts
Detection threshold in sigma (4.) (CR or value): Mesmo valor usado na vez anterior
New detection threshold: 4. sigma 12.8 counts
Minimum good data value (INDEF) (CR or value): 5000 Valor alto para selecionar somente objetos brilhantes
New minimum good data value: 5000. counts
Maximum good data value (INDEF) (CR or value): 31000 Mesmo valor usando na vez anterior
New maximum good data value: 31000. counts
Image: ../hd110cv0001.fits fwhmpsf: 3.4 ratio: 1. theta: 0. nsigma: 1.5
189.04 54.84 -7.205 0.213 0.034 0.068 1
516.68 530.21 -6.816 0.367 -0.213 -0.140 2
519.17 566.79 -6.870 0.425 -0.091 -0.133 3
734.21 840.27 -6.774 0.443 0.210 0.264 4
738.16 876.83 -7.171 0.415 0.111 -0.017 5 Foram achados 5 objetos brilhantes
threshold: 12.8 relerr: 1.444 0.2 <= sharp <= 1. -1. <= round <= 1.
Running IMALIGN ... Rodando o imalign para alinhar imagens
Warning: bigbox must be odd, using 41
#Coords Image X-center Err Y-center Err Num
../hd110cv0001.fits 188.920 (0.027) 54.889 (0.028) 1
../hd110cv0001.fits 514.716 (0.023) 530.137 (0.023) 2
../hd110cv0001.fits 518.677 (0.024) 566.609 (0.022) 3
../hd110cv0001.fits 733.692 (0.020) 840.454 (0.021) 4
../hd110cv0001.fits 737.779 (0.021) 876.768 (0.021) 5
../hd110cv0002.fits 188.732 (0.026) 56.323 (0.028) 1
../hd110cv0002.fits 514.374 (0.021) 531.527 (0.022) 2
../hd110cv0002.fits 518.330 (0.023) 567.999 (0.022) 3
../hd110cv0002.fits 733.267 (0.019) 842.047 (0.021) 4
../hd110cv0002.fits 737.346 (0.020) 878.359 (0.020) 5
../hd110cv0003.fits 188.747 (0.030) 55.623 (0.030) 1
../hd110cv0003.fits 514.134 (0.025) 530.993 (0.025) 2
../hd110cv0003.fits 518.080 (0.026) 567.470 (0.024) 3
../hd110cv0003.fits 733.137 (0.022) 841.244 (0.024) 4
../hd110cv0003.fits 737.211 (0.023) 877.564 (0.023) 5
../hd110cv0004.fits 188.420 (0.026) 56.156 (0.027) 1
../hd110cv0004.fits 514.000 (0.022) 531.565 (0.023) 2
../hd110cv0004.fits 517.949 (0.022) 568.042 (0.020) 3
../hd110cv0004.fits 732.705 (0.019) 842.014 (0.021) 4
../hd110cv0004.fits 736.786 (0.020) 878.324 (0.020) 5
../hd110cv0005.fits 187.468 (0.028) 54.816 (0.030) 1
../hd110cv0005.fits 513.317 (0.022) 530.268 (0.023) 2
../hd110cv0005.fits 517.277 (0.023) 566.737 (0.022) 3
../hd110cv0005.fits 732.021 (0.019) 840.463 (0.021) 4
../hd110cv0005.fits 736.110 (0.021) 876.775 (0.021) 5
../hd110cv0006.fits 187.296 (0.036) 55.876 (0.036) 1
../hd110cv0006.fits 512.936 (0.027) 531.477 (0.027) 2
../hd110cv0006.fits 516.901 (0.029) 567.951 (0.027) 3
../hd110cv0006.fits 731.756 (0.026) 841.715 (0.026) 4
../hd110cv0006.fits 735.837 (0.027) 878.028 (0.025) 5
../hd110cv0007.fits 188.288 (0.029) 53.806 (0.029) 1
../hd110cv0007.fits 514.054 (0.023) 529.467 (0.022) 2
../hd110cv0007.fits 518.008 (0.024) 565.935 (0.022) 3
../hd110cv0007.fits 733.031 (0.020) 839.843 (0.021) 4
../hd110cv0007.fits 737.121 (0.021) 876.163 (0.020) 5
../hd110cv0008.fits 188.711 (0.026) 55.759 (0.028) 1
../hd110cv0008.fits 514.307 (0.022) 531.162 (0.022) 2
../hd110cv0008.fits 518.265 (0.022) 567.622 (0.021) 3
../hd110cv0008.fits 733.069 (0.018) 841.644 (0.020) 4
../hd110cv0008.fits 737.162 (0.020) 877.956 (0.020) 5
#Refcoords Reference X-center Err Y-center Err Num
../hd110cv0001.fits 188.920 (0.027) 54.889 (0.028) 1
../hd110cv0001.fits 514.716 (0.023) 530.137 (0.023) 2
../hd110cv0001.fits 518.677 (0.024) 566.609 (0.022) 3
../hd110cv0001.fits 733.692 (0.020) 840.454 (0.021) 4
../hd110cv0001.fits 737.779 (0.021) 876.768 (0.021) 5
#Shifts Image X-shift Err Y-shift Err N Internal
../hd110cv0001.fits 0.00 (0.01) 0.00 (0.01) 5 (0.00,0.00)
../hd110cv0002.fits 0.35 (0.01) -1.48 (0.01) 5 (0.04,0.05)
../hd110cv0003.fits 0.50 (0.02) -0.81 (0.02) 5 (0.08,0.02)
../hd110cv0004.fits 0.78 (0.01) -1.45 (0.01) 5 (0.09,0.05)
../hd110cv0005.fits 1.52 (0.01) -0.04 (0.01) 5 (0.06,0.04)
../hd110cv0006.fits 1.81 (0.02) -1.24 (0.02) 5 (0.06,0.07)
../hd110cv0007.fits 0.66 (0.01) 0.73 (0.01) 5 (0.01,0.09)
../hd110cv0008.fits 0.45 (0.01) -1.06 (0.01) 5 (0.08,0.06)
#Trim_Section = [3:1050,2:1048]
continue (yes|no)? Confirmar alinhamento correto
y
Running COORSHIFT ... Rodando o coorshift para calcular deslocamentos
x-minimum distance to the border (pixels)?: 30 Colocar valor de segurança no eixo-x
30.
y-minimum distance to the border (pixels)?: 30 Colocar valor de segurança no eixo-y
30.
Analyzing file...
shift (X,Y) : 0.45 -1.06
delta (X,Y) : 30. 30.
limits (X,Y) : 30.45 1018.94
filtered objects 62
creating coordinate file for ../hd110cv0001.fits xshift 0. yshift 0.
creating coordinate file for ../hd110cv0002.fits xshift 0.35 yshift -1.48
creating coordinate file for ../hd110cv0003.fits xshift 0.5 yshift -0.81
creating coordinate file for ../hd110cv0004.fits xshift 0.78 yshift -1.45
creating coordinate file for ../hd110cv0005.fits xshift 1.52 yshift -0.04
creating coordinate file for ../hd110cv0006.fits xshift 1.81 yshift -1.24
creating coordinate file for ../hd110cv0007.fits xshift 0.66 yshift 0.73
creating coordinate file for ../hd110cv0008.fits xshift 0.45 yshift -1.06
Running PHOT ... Rodando o phot para fazer a fotometria
FWHM of the PSF in scale units: 2.5
is it correct (yes|no)? Pede confirmar valor do fwhm
n Colocar "n" = "no"
FWHM of the PSF in scale units:
3.4 Colocar valor usado no daofind acima
Standard deviation of background in counts: INDEF
is it correct (yes|no)? Pede confirmar valor default do sigma
n Colocar "n" = "no"
Standard deviation of background in counts:
3.2 Colocar valor usado no daofind acima
Centering algorithm (centroid) (CR or value): none Indicar algoritmo para calculo do centroide
New centering algorithm: none Neste caso none não voltará a fazer o centroide mas pode ser usado centroid
Sky fitting algorithm (mode) (CR or value): Indicar algoritmo para calculo de céu
Sky fitting algorithm: mode Deixar valor default: mode
Inner radius of sky annulus in scale units (10.) (CR or value): 20 Indicar raio interno do anel do céu
New inner radius of sky annulus: 20. scale units 20. pixels Neste caso 20 ~ 6*fwhm
Width of the sky annulus in scale units (10.) (CR or value): 10 Indicar largura do anel de céu
New width of the sky annulus: 10. scale units 10. pixels Neste caso 10 ~ 2.5*fwhm
File/list of aperture radii in scale units (3:12:1) (CR or value): Confirmar aperturas usadas
Aperture radius 1: 3. scale units 3. pixels
Aperture radius 2: 4. scale units 4. pixels
Aperture radius 3: 5. scale units 5. pixels
Aperture radius 4: 6. scale units 6. pixels
Aperture radius 5: 7. scale units 7. pixels
Aperture radius 6: 8. scale units 8. pixels
Aperture radius 7: 9. scale units 9. pixels
Aperture radius 8: 10. scale units 10. pixels
Aperture radius 9: 11. scale units 11. pixels
Aperture radius 10: 12. scale units 12. pixels
Minimum good data value (INDEF) (CR or value): Valor mínimo de contagens para fazer fotometria
New minimum good data value: INDEF counts Manter valor default: INDEF
Maximum good data value (INDEF) (CR or value): Valor máximo de contagens para fazer fotometria
New maximum good data value: INDEF counts Manter valor default: INDEF
Running TXDUMP ... Rodando txdump para obter a fotometria
textfiles = "*.mag.1" Input apphot/daophot text database(s)
fields = "image,msky,nsky,rapert[1-10],sum[1-10],area[1-10]" Fields to be extracted
expr = "yes" Boolean expression for record selection
(headers = no) Print the field headers ?
(parameters = yes) Print the parameters if headers is yes ?
(mode = "ql") Mode of task
Running PCCD ... Rodando pccd para calcular informação polarimétrica
filename = "txdump6748md" input file (.dat)
(nstars = 62) number of stars (maximum 2000)
(nhw = 8) number of postions of wave-plate (maximum 16)
(nap = 10) number of apertures (maximum 10)
(calc = "c") analyser: calcite (c) / polaroid (p)
(readnoise = 0.82) CCD readnoise (adu)
(ganho = 5.) CCD gain (e/adu)
(deltatheta = 0.) correction in polarization angle (degrees)
(fileout = "pccd6748nd") output file (.log)
(fileexe = "/iraf/extern/pccdpack/pccd/pccd2000.exe") pccd execute file (.exe)
(mode = "ql")
Running MACROL ... Rodando macrol para obter apertura que minimiza o erro polarimétrico
File : pccd6748nd
Evaluating...star 1 of 62 , aperture 1. of 10.
Evaluating...star 1 of 62 , aperture 2. of 10.
Evaluating...star 1 of 62 , aperture 3. of 10.
Evaluating...star 1 of 62 , aperture 4. of 10.
Evaluating...star 1 of 62 , aperture 5. of 10.
...
...
Evaluating...star 62 of 62 , aperture 5. of 10.
Evaluating...star 62 of 62 , aperture 6. of 10.
Evaluating...star 62 of 62 , aperture 7. of 10.
Evaluating...star 62 of 62 , aperture 8. of 10.
Evaluating...star 62 of 62 , aperture 9. of 10.
Evaluating...star 62 of 62 , aperture 10. of 10.
Q U SIGMA P THETA SIGMAtheor. APERTURE STAR
0.01291 0.01646 0.00435 0.02092 25.9 0.00409 3. 1
0.01375 0.01807 0.00024 0.02271 26.4 0.00051 8. 2
0.00640 0.04037 0.00452 0.04087 40.5 0.00942 3. 3
0.00787 0.01326 0.00681 0.01542 29.7 0.00721 3. 4
0.01331 0.03005 0.00248 0.03287 33.1 0.00258 4. 5
...
...
-0.00634 0.04545 0.00817 0.04589 49.0 0.00396 9. 58
0.01844 0.03644 0.01130 0.04084 31.6 0.00941 3. 59
0.01178 0.00962 0.00178 0.01521 19.6 0.00220 5. 60
0.16252 -0.18449 0.09849 0.24587 155.7 0.06503 3. 61
0.00979 0.05214 0.00283 0.05305 39.7 0.00296 3. 62
Running SELECT ... Rodando select (ver)
XCENTER YCENTER P THETA SIGMA ID APERTURE STAR
892.010 50.055 0.02092 25.9 0.00435 10 3 1
188.987 54.887 0.02271 26.4 0.00024 13 8 2
797.924 87.841 0.04087 40.5 0.00452 18 3 3
549.589 125.262 0.01542 29.7 0.00681 28 3 4
256.323 132.555 0.03287 33.1 0.00248 30 4 5
...
...
212.021 970.195 0.01948 18.6 0.00368 227 3 57
738.778 972.219 0.04589 49.0 0.00817 228 9 58
148.718 991.234 0.04084 31.6 0.01130 235 3 59
975.222 1006.889 0.01521 19.6 0.00178 239 5 60
578.170 1011.669 0.24587 155.7 0.09849 242 3 61
572.838 1016.593 0.05305 39.7 0.00283 243 3 62
pc>
FIM
(Ver exemplo de uso das rotinas REFER e VECPLOT)
