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.