Skip to content

Problem with "compute_shift_list" #64

@FFalkenbergM

Description

@FFalkenbergM

Percebi que em algumas das minhas séries de imagens o "compute_shift_list" tem gerado uma lista de shifts incorreta.

Por exemplo:

Quando eu executo o comando:

"shift_list = compute_shift_list(hd94851_framesv,ref_image=0, algorithm='cross-correlation', upsample_factor=10)"

o shift list de um conjunto de 16 imagens de uma estrela padrão sai da seguinte forma:

image dx dy
1.0 0.0 0.0
2.0 60.1 58.9
3.0 0.9 -2.3
4.0 0.3 -2.1
5.0 0.9 -1.7
6.0 0.4 -1.8
7.0 1.2 -1.0
8.0 1.8 0.1
9.0 2.3 1.3
10.0 2.3 3.2
11.0 2.7 5.1
12.0 2.3 4.1
13.0 1.9 2.2
14.0 1.7 1.7
15.0 1.2 1.4
16.0 0.8 0.9

Percebi que na imagem n.2 o shifts estão deslocados (aparentemente o algoritmo está confundindo feixe ordinário com o extraordinário). Já verifiquei manualmente que esse shift da imagem 2 não é real.

Além disso, quando eu mudo a imagem de referência os shifts mudam de uma forma que não faz sentido pra mim:

"shift_list = compute_shift_list(hd94851_framesv,ref_image=1, algorithm='cross-correlation', upsample_factor=1000)"

image dx dy
1.0 -60.103 -58.881
2.0 0.0 0.0
3.0 -0.501 -0.6950000000000001
4.0 -0.33 -1.396
5.0 1.038 -1.058
6.0 -1.158 -1.147
7.0 0.501 1.296
8.0 1.685 2.531
9.0 2.833 3.84
10.0 -0.18 4.471
11.0 0.788 5.782
12.0 0.69 4.705
13.0 0.317 3.001
14.0 0.696 2.161
15.0 0.295 1.873
16.0 3.885 2.348

E também:

"shift_list = compute_shift_list(hd94851_framesv,ref_image=2, algorithm='cross-correlation', upsample_factor=1000)"

image dx dy
1.0 -0.864 2.29
2.0 0.501 0.6950000000000001
3.0 0.0 0.0
4.0 -0.139 0.088
5.0 0.009 0.237
6.0 0.012 0.215
7.0 0.353 1.541
8.0 1.117 2.605
9.0 1.555 3.928
10.0 1.295 5.494
11.0 1.824 6.404
12.0 1.674 5.849
13.0 0.812 3.94
14.0 1.197 3.851
15.0 0.509 2.802
16.0 0.486 3.162

Metadata

Metadata

Assignees

No one assigned

    Labels

    No labels
    No labels

    Type

    No type

    Projects

    No projects

    Milestone

    No milestone

    Relationships

    None yet

    Development

    No branches or pull requests

    Issue actions