Redução de Imagens Polarimétricas Objetos Extensos
1. Os dados sao tomados com Lamina Retardadora + Polaroid HNP'B (BVR) ou HR (RI). Oito posicoes da lamina e' um bom numero.
2. Para o mesmo campo, tambem tomamos medidas com a Lamina + calcita.
3. Aplicamos correcao de Bias e Flat.
4. Imagens sao registradas por ex. com a primeira imagem da serie. Ver manual recente do Iraf a respeito para varias opcoes.
5. Para as imagens com a calcita, reducao usual com o phot e pccd.
- A partir delas, determinamos q,u para estrelas de campo.
6. Para imagens com o Polaroid, reduzimos (phot+pccd) as mesmas estrelas de campo.
- De (5)+(6), determinamos [q,u]polaroid vs. [q,u]calcita
7. Para as imagens com o Polaroid:
- Fazer subtracao de background, por ex. com a task imsurfit
Por ex., para subtrair uma constante:
xorder=1 yorder=1
regions=section
section=Rmon.sections
type=resid
- Binnar imagens usando a task blkavg:
blkavg RmonPV01B[400:912,250:762] RmonPVF01Bbl 16 16
- Usar a task imarith para fazer:
imarith Imagem01 - Imagem03 qnum
imarith Imagem01 + Imagem03 qden
imdivide qnum qden q
e analogamente para u usando as imagens 02 e 04 da sequencia.
- Listar os valores de q, u e da imagem usando a task listpix:
listpix q > q.dat listpix u > u.dat (listpix imagem > i.dat)
Por ex., q.dat vai conter i,j,q(i,j) para cada bin da imagem.
- rodar /users/mario/iraf/pexccd => p.dat = arquivo de x,u,polarizacao.
p.dat vai conter
i,j,p (0<p<1), theta (deg):
...
12.0 1.0 .60980 14.4
# 13.0 1.0 1.85774 20.0
# 14.0 1.0 1.15453 43.7
# 15.0 1.0 6.68155 42.1
# 16.0 1.0 3.55853 138.8
17.0 1.0 .64221 2.9
...
- fazer display de uma das imagens originais ou da soma delas e usar a
task fieldplot, com p.dat como entrada, para superpor os vetores
de polarizacao.