REDUÇÃO DADOS IV - (padrões, objetos individuais)

 

    Baixar imagens de teste: 2005jun25 - hd164740+flatoff+flaton (93 Mb)

 

  1. É precio ter instalado o PCCDPACK versão setembro/2007 ou posterior. As últimas atualizações podem ser lidas aqui.

     

     

  2. Editar cabeçalhos dos flats. Se for o caso, acrecentar aos flats (off e on) novo keyword imagetyp=flat usando hedit.

     

     

    pc> hedit (executa a rotina)

     

     

  3. Combinar os flatoff com flatcombine:

     

     

    pc> flatcombine (executa a rotina)

     

     

  4. Corrigir arquivos flaton com flatoff_med.fits usando imarith:

     

     

    pc> imarith (executa a rotina)

     

     

  5. Combinar arquivos flaton corrigidos usando flatcombine:

     

     

    pc> flatcombine (executa a rotina)

     

  1. Calcular o valor médio (MEAN) de flaton_med.fits usando imstat:

     

     

    pc> imstat flaton_med.fits (executa a rotina)

     

     

  2. Normalizar flaton_med.fits usando imarith:

     

     

    pc> imarith (executa a rotina)

     

  1. Corrigir imagens pontilhadas em cada seqüência da lâmina usando prepiv.

CASO 1: Pontilhamento feito em cada posição da lâmina e com NEXP=1. Se temos, por exemplo, cinco imagens pontilhadas para a 1a. posição da lâmina tiradas com NEXP=1:

 

hd164ch10001.fits (pontilhamento 1)

hd164ch10002.fits (pontilhamento 2)

hd164ch10003.fits (pontilhamento 3)

hd164ch10004.fits (pontilhamento 4)

hd164ch10005.fits (pontilhamento 5)

 

então, o prefixo hd164ch1 define a seqüência, e portanto, sequ=hd1641. O céu será calculado pela moda da seqüência e subtraído para cada imagem. Se flatcor=yes cada imagem (já corrigida pelo céu) será corrigida pelo flat (flatn). Se deseja preservar a imagem do ceú então delsky=no. Se flatcor=yes seram criadas cinco imagems com o prefixo flskyhd164ch1. Se flatcor=no, então o prefixo será somente skyhd164ch1.

 

 

 

pc> prepiv (executa a rotina)

 

O mesmo deve ser feito para as outras posições da lâmina.

 

 

CASO 2: Pontilhamento feito em cada posição da lâmina e com NEXP>1. Se temos, por exemplo, cinco seqüências de imagens pontilhadas para a 1a. posição da lâmina tiradas com NEXP=3:

 

hd164ch110001.fits, hd164ch110002.fits, hd164ch110003.fits, (pontilhamento 1)

hd164ch120001.fits, hd164ch120002.fits, hd164ch120003.fits, (pontilhamento 2)

hd164ch130001.fits, hd164ch130002.fits, hd164ch130003.fits, (pontilhamento 3)

hd164ch140001.fits, hd164ch140002.fits, hd164ch140003.fits, (pontilhamento 4)

hd164ch150001.fits, hd164ch150002.fits, hd164ch150003.fits, (pontilhamento 5)

 

então, o prefixo hd164ch1 define a seqüência completa e o procedimento a seguir é semelhante ao CASO 1. Neste exemplo o céu é calculado para a seqüência completa (5x3=15 imagens).

 

O mesmo deve ser feito para as outras posições da lâmina.

 

CASO 3: Pontilhamento feito entre as posições da lâmina e com NEXP>=1. Por exemplo, se o pontilhamento for feito entre quatro posições consecutivas da lâmina (digamos, as posições [1-4]) e com NEXP=3 por posição, então devem-se criar um diretório de trabalho paras essas primeiras quatro posições (digamos: pos1-4). Esse diretório deve conter, portanto, a seguinte seqüência de imagens:

 

pos1-4/

hd164ch10001.fits, hd164ch10002.fits, hd164ch10003.fits, (posição 1 da lâmina = pontilhamento 1)

hd164ch20001.fits, hd164ch20002.fits, hd164ch20003.fits, (posição 2 da lâmina = pontilhamento 2)

hd164ch30001.fits, hd164ch30002.fits, hd164ch30003.fits, (posição 3 da lâmina = pontilhamento 3)

hd164ch40001.fits, hd164ch40002.fits, hd164ch40003.fits, (posição 4 da lâmina = pontilhamento 4)

 

Assim, o prefixo hd164ch definirá a seqüência completa e o procedimento a seguir é semelhante ao CASO 1. Neste exemplo o céu é calculado para a seqüência completa (4x3=12 imagens).

 

O mesmo deve ser feito para as outras quatro posições da lâmina consecutivas. Assim, deve-se criar o respectivo diretório de trabalho (digamos, pos5-8 para as posições consecutivas [5-8]).

 

  1. Registrar imagens pontilhadas em cada seqüência da lâmina usando align (cada CASO é com respeito ao passo anterior).

     

    CASO 1: Se o resultado do paso anterior são imagens com o prefixo flskyhd164ch1, a lista de parâmetros para align deve ser:

     

     

    pc> align (executa a rotina)

     

    Neste paso sumfiles=yes permite que as imagens registradas sejam somadas criando o arquivo filesum (neste exemplo, filesum=sum1, pois refere-se à 1a. posição da lâmina). A rotina deve ser executada com o ds9 aberto. Inicialmente será mostrada a 1a. imagen da sequencia, onde deve-se marcar com o cursor (digitando 'r' para verificar o perfil e em seguida 'a' para anotar o centroide) as posições ordinária e extraordinária do objeto a analisar. Uma confirmação será pedida mostrando as coordenadas achadas. Na sequencia, seram mostradas cada una das imagens a registrar (incluindo a 1a. imagem novamente) onde deve-se marcar com o cursor (digitando 'r' para verificar o perfil e em seguida 'a' para anotar o centroide) únicamente a imagem ordinária do objeto a analizar. A confirmação de novo será pedida.

     

    O mesmo deve ser feito para as outras posições da lâmina. Assim, para uma seqüência de 8 posições da lâmina, no final teremos 8 imagens soma (sum1.fits, sum2.fits, ...., sum8.fits).

     

     

    CASO 2: Neste caso devem registrarse e somarse primeiro as imagens de cada posição de pontilhamento. Para o pontilhamehto 1, os parâmetros do align são:

     

    imgref   = flskyhd164ch110001.fits

    images   = flskyhd164ch11*.fits

    shifts   = flskyhd164ch110001

    confirm  = no

    sumfiles = sum11

     

    O parâmetro confirm=no pois neste caso nao é preciso confirmar as imagens individuais do pontilhamanto. A imagem soma parcial será sumfiles=sum11 neste caso.

     

    Isto deve repetirse para as outras posições de pontilhamento de modo a criar imagens soma parciais: sum11.fits, sum12.fits e sum13.fits.

     

    Finalmente, deve ser rodado mais uma vez o align de modo a registrar as imagens soma parciais para cada posição da lâmina. Os parâmetros para a 1a. posiçãao da lâmina devem ser:

     

    imgref   = sum11.fits

    images   = sum11*.fits

    shifts   = sum11

    confirm  = yes

    sumfiles = sumtot1

     

    Neste caso confirm=yes, pois é preciso comfirmar as coordenadas como no caso CASO 1. A imagem soma total para a 1a. posição da lâmina será sumfiles=sum1.

     

    O mesmo deve ser feito para as outras posições da lâmina. Assim, para uma seqüência de 8 posições da lâmina, no final teremos 8 imagens soma totais (sumtot1.fits, sumtot2.fits, ...., sumtot8.fits).

     

     

    CASO 3: Neste caso devem registrarse e somarse as imagens em cada posição da lâmina. Os parametros do align são os mesmos que no CASO 1 exceptuando confirm=no, pois não é preciso confirmar as coordenadas. No exemplo, ao final teremos quatro imagens soma no diretório pos1-4 e quatro no diretório pos5-8. Assim, é recomendável mover as oito imagens a um diretório acima de modo a ficar reunidas em conjunto para facilitar o passo seguinte.

     

  1. Registrar as imagens soma do paso anterior. É similar ao passo anterior mas com sumfiles=no pois somente registraremos as imagens de modo a facilitar o passo posterior de cálculo da polarização. Para o CASO 2 do passo anterior, deve trocar o prefixo sum por sumtot na lista de parâmetros a seguir:

     

     

    pc> align (executa a rotina)

     

    Os arquivos de saida terão o prefixo shsum (shsum1.fits, shsum2.fits, ...., shsum8.fits). Para o CASO 2 do passo anterior shsumtot (shsumtot1.fits, shsumtot2.fits,..., shsumtot8.fits).

     

     

  2. Cálculo da polarização - pccdgen. Exemplo para solução lambda/4 (retar=87, zero=47).

     

     

    Neste passo deve desconsiderar-se verificar filename e fileout pois o pccdgen será usado internamente pela rotina aperquick (passo 13). A rotina NÂO deve ser executada, somente devem ser editados os parâmetros apropiados.

     

     

  3. Cálculo da polarização - pospars. Exemplo para 8 posições consecutivas ativadas.

     

     

    A rotina somente permite editar seus parâmetros. NÂO é de execução.

     

     

  4. Cálculo da polarizacao - aperquick. Verifica o resultado para varios aneis de ceu e é semelhante ao quickpol. Para o CASO 2 do passo 10 deve-se trocar o prefixo sum por sumtot na lista de parâmetros a seguir:

     

     

    pc> aperquick (executa a rotina)

     

    No final devem-se conferir os resultados para os varios aneis do ceu no arquivo (namefile) criado no diretorio de trabalho (subdir). A seguir exemplo da saida da rotina (arquivo namefile):

     

    annulus    V    SIGMAV      Q       U     SIGMA     P    THETA SIGMAth.   rms    APER. STAR

    20.     0.00103 0.00083  0.00003 0.02155 0.00143 0.02155  45.0 0.00019 0.0013570   6.   1   first

    20.     0.00110 0.00076 -0.00037 0.02299 0.00132 0.02299  45.5 0.00019 0.0012476   10.  1   full

    30.     0.00099 0.00057 -0.00147 0.02255 0.00100 0.02260  46.9 0.00019 0.0009446   7.   1   first(melhor resultado)

    30.     0.00099 0.00057 -0.00147 0.02255 0.00100 0.02260  46.9 0.00019 0.0009446   7.   1   full

    40.     0.00110 0.00104 -0.00393 0.02185 0.00178 0.02220  50.1 0.00019 0.0017014   7.   1   first

    40.     0.00110 0.00104 -0.00393 0.02185 0.00178 0.02220  50.1 0.00019 0.0017014   7.   1   full

    50.     0.00106 0.00118 -0.00407 0.02228 0.00202 0.02265  50.2 0.00019 0.0019366   7.   1   first

    50.     0.00106 0.00118 -0.00407 0.02228 0.00202 0.02265  50.2 0.00019 0.0019366   7.   1   full

    60.     0.00124 0.00170 -0.00615 0.02181 0.00288 0.02266  52.9 0.00019 0.0027996   7.   1   first

    60.     0.00124 0.00170 -0.00615 0.02181 0.00288 0.02266  52.9 0.00019 0.0027996   7.   1   full

    70.     0.00103 0.00197 -0.00714 0.02127 0.00330 0.02243  54.3 0.00019 0.0032455   6.   1   first

    70.     0.00103 0.00197 -0.00714 0.02127 0.00330 0.02243  54.3 0.00019 0.0032455   6.   1   full

    80.     0.00105 0.00198 -0.00723 0.02136 0.00331 0.02255  54.3 0.00019 0.0032547   6.   1   first

    80.     0.00105 0.00198 -0.00723 0.02136 0.00331 0.02255  54.3 0.00019 0.0032547   6.   1   full

    90.     0.00105 0.00203 -0.00732 0.02124 0.00339 0.02247  54.5 0.00019 0.0033383   6.   1   first

    90.     0.00105 0.00203 -0.00732 0.02124 0.00339 0.02247  54.5 0.00019 0.0033383   6.   1   full

    ...

    ...

    280.    0.00104 0.00210 -0.00757 0.02131 0.00349 0.02261  54.8 0.00019 0.0034439   6.   1   first

    280.    0.00104 0.00210 -0.00757 0.02131 0.00349 0.02261  54.8 0.00019 0.0034439   6.   1   full

    290.    0.00101 0.00210 -0.00761 0.02131 0.00350 0.02262  54.8 0.00019 0.0034586   6.   1   first

    290.    0.00101 0.00210 -0.00761 0.02131 0.00350 0.02262  54.8 0.00019 0.0034586   6.   1   full

    300.    0.00102 0.00211 -0.00760 0.02126 0.00350 0.02257  54.8 0.00019 0.0034638   6.   1   first

    300.    0.00102 0.00211 -0.00760 0.02126 0.00350 0.02257  54.8 0.00019 0.0034638   6.   1   full

    310.    0.00101 0.00212 -0.00765 0.02130 0.00353 0.02264  54.9 0.00019 0.0034926   6.   1   first

    310.    0.00101 0.00212 -0.00765 0.02130 0.00353 0.02264  54.9 0.00019 0.0034926   6.   1   full

     

     

     

    FIM

     

    OUT/2007 AAPQ