Problème sur l'astrométrie d'un NEO 1998 OR2 (52768)

[Marc DELDEM]

Bonsoir Cyril,

En ce moment, le NEO 1998 OR2 (52768) est très proche de la Terre, 0.087UA et magnitude 13.


J'ai fais une série hier soir pour compléter la CdR.
Le problème c'est que malgré une mise à jour de la base des astéroïdes, Prism n'indique pas la bonne position de l'astéroïde.

On récupère les coordonnnées Alpha/Delta de l'astéroïde via:

Code : Tout sélectionner

GETASTERCOORD NomAster$ jd[nImage] Latitude Longitude Altitude Alpha Delta SpeedRA SpeedDEC SpeedComb ApproxMagn GeoDistance Err 

Puis on transforme les coordonnées en X,Y avec:

Code : Tout sélectionner

GETIMG_RADDEC_TO_XY Img Alpha Delta PosRechX PosRechY ErrXY	

Mais si ensuite on essaye de faire le centroïde sur les coordonnées X,Y retournées, cela ne marche pas.
Il y a un fort décalage.
J'ai donc du introduire un offset sur X et Y pour pouvoir faire le centroïde et trouver la position X,Y du centroïde.

Voici les résultats sur les 10 premières images de la séquence:

Sélection des images CPA
Besoin d'un offset 1998OR2 / Image:T_1998OR2-R'-1.CPA Position théorique: (X=1066,Y=714) Position Centroïde: (X=1062,Y=722) AbsDelta: (X=5,Y=8)
Besoin d'un offset 1998OR2 / Image:T_1998OR2-R'-2.CPA Position théorique: (X=1064,Y=711) Position Centroïde: (X=1060,Y=719) AbsDelta: (X=5,Y=8)
Besoin d'un offset 1998OR2 / Image:T_1998OR2-R'-3.CPA Position théorique: (X=1062,Y=709) Position Centroïde: (X=1057,Y=717) AbsDelta: (X=5,Y=8)
Besoin d'un offset 1998OR2 / Image:T_1998OR2-R'-4.CPA Position théorique: (X=1062,Y=707) Position Centroïde: (X=1057,Y=715) AbsDelta: (X=5,Y=8)
Besoin d'un offset 1998OR2 / Image:T_1998OR2-R'-5.CPA Position théorique: (X=1061,Y=705) Position Centroïde: (X=1056,Y=713) AbsDelta: (X=5,Y=8)
Besoin d'un offset 1998OR2 / Image:T_1998OR2-R'-6.CPA Position théorique: (X=1059,Y=703) Position Centroïde: (X=1054,Y=711) AbsDelta: (X=5,Y=8)
Besoin d'un offset 1998OR2 / Image:T_1998OR2-R'-7.CPA Position théorique: (X=1058,Y=701) Position Centroïde: (X=1053,Y=708) AbsDelta: (X=5,Y=8)
Besoin d'un offset 1998OR2 / Image:T_1998OR2-R'-8.CPA Position théorique: (X=1057,Y=698) Position Centroïde: (X=1052,Y=706) AbsDelta: (X=5,Y=8)
Besoin d'un offset 1998OR2 / Image:T_1998OR2-R'-9.CPA Position théorique: (X=1056,Y=696) Position Centroïde: (X=1051,Y=704) AbsDelta: (X=5,Y=8)
Besoin d'un offset 1998OR2 / Image:T_1998OR2-R'-10.CPA Position théorique: (X=1045,Y=677) Position Centroïde: (X=1040,Y=685) AbsDelta: X=5,Y=8)
Script terminé ...

Le script et les images sont disponibles pendant 7 jours sur WeTransfer:
https://wetransfer.com/downloads/f3b23c ... 440/bd42f5

Le script est un extrait du script utilisé pour faire la réduction photométrie des astéroïdes et variables, je ne pense pas que le problème vienne de là.
Si tu veux jeter un oeil, cela peut être intéressant de comprendre pourquoi il y a une erreur aussi importante.
Sachant que j'ai un échantillonnage de 1.13"/pixel, on a donc une erreur de plus de 10" sur la position réelle de l'astéroïde.
La vitesse de la bestiole au cours des acquisitions est de 4.1"/minute (info de Prism).
Mes images sont calibrées en astrométrie avec un polynôme de degré 2. Je n'ai jamais rencontré ce problème jusqu'à ce jour (mais c'est mon premier NEO).
Et je pense que les paramètres orbitaux de l'astéroïde sont bien définis (en tout cas faut l'espérer pour un NEO...)

Marc

[Cyril Cavadore]

10" ... et tu affiche seul cet astéroïde ? Parce que si tu fais un affichage avec les 950 000 autres, la precision (pour etre assez rapide) ne sera pas aussi bonne.

[Marc DELDEM]

OK alors je compare les éphémérides de Prism et celles du MPC (https://www.minorplanetcenter.net/iau/MPEph/MPEph.html):

PRISM:

Site d'observation:
Longitude=5.0666667°, Latitude=43.6833333°, Altitude=125m

Ephémérides pour 1998 OR2:
14/04/2020 19:54:00 ; 08 29 36.233 +27 33 09.28
14/04/2020 19:55:00 ; 08 29 36.388 +27 33 05.80
14/04/2020 19:56:00 ; 08 29 36.544 +27 33 02.33
14/04/2020 19:57:00 ; 08 29 36.699 +27 32 58.85
14/04/2020 19:58:00 ; 08 29 36.855 +27 32 55.38
14/04/2020 19:59:00 ; 08 29 37.011 +27 32 51.90
14/04/2020 20:00:00 ; 08 29 37.167 +27 32 48.42
14/04/2020 20:01:00 ; 08 29 37.322 +27 32 44.94
14/04/2020 20:02:00 ; 08 29 37.478 +27 32 41.46
14/04/2020 20:03:00 ; 08 29 37.634 +27 32 37.98
14/04/2020 20:04:00 ; 08 29 37.790 +27 32 34.50

MPC:

2020 04 14 195400 08 29 35.8 +27 33 13
2020 04 14 195500 08 29 36.0 +27 33 09
2020 04 14 195600 08 29 36.2 +27 33 06
2020 04 14 195700 08 29 36.3 +27 33 02
2020 04 14 195800 08 29 36.5 +27 32 59
2020 04 14 195900 08 29 36.6 +27 32 55
2020 04 14 200000 08 29 36.8 +27 32 52
2020 04 14 200100 08 29 36.9 +27 32 48
2020 04 14 200200 08 29 37.1 +27 32 45
2020 04 14 200300 08 29 37.2 +27 32 41

Comparaison:

Code : Tout sélectionner

Mon image 14/04/2020 - 19:56:57.043  08 29 36.119 +27 33 06.18
MPC       14/04/2020 - 19:57:00.000  08 29 36.300 +27 33 02.00
PRISM     14/04/2020 - 19:57:00.000  08 29 36.699 +27 32 58.85

Code : Tout sélectionner

Mon image 14/04/2020 - 19:59:00.743  08 29 36.446 +27 32 58.78
MPC       14/04/2020 - 19:59:00.000  08 29 36.600 +27 32 55.00
PRISM     14/04/2020 - 19:59:00.000  08 29 37.011 +27 32 51.90

Marc

[Cyril Cavadore]

et l'ecart est de combien ?

[Marc DELDEM]

En gros 5 pixels en X/RA soit 5.6" et 8 pixels en Y/DEC soit 9.0".

Mais bon on voit qu'il y a aussi un petit écart entre mes mesures par centroïdes et les éphémérides du MPC.

Marc

[Cyril Cavadore]

Tu peux en parler a Raoul Behrend, c'est l'expert du sujet.
Mais je pense que le passage pres de la terre torture le modele et le calcul doit moins bien marcher.
A quelle distance de la terre c'etait ?

[Marc DELDEM]

C'était à une distance de 0.087UA.
OK je vais en parler à Raoul.

Merci.
Marc