Analyse de la Version du 12 mars 2010 du Catalogue CDS
Etude statistique du catalogue pour détecter toute éventuelle aberration
La version analysée du catalogue CDS date du 12/03/2010 et contient 38541 calibrateurs filtrés pour 21417 VOTables.
Il aura fallu 2 jours, 16 heures et 8 minutes pour l'obtenir, avec la version sclsvr 3.6.4 (avec limite de recherche CDS à 1000).
Discrétisation de V-R et V-I
Après découverte par
GillesDuvert puis étude par
XavierDelfosse sur sollicitation de
SylvainLafrasse et
DanielBonneau, il s'avère bien que cette discrétisation soit due à la méthode de calcul de R et I à partir du type spectral qui est lui-même discrétisé, et de B et V.
Déficit d'étoiles et de VOTable par rapport à la version de février
Il manquait 151 calibrateurs et 116 VOTables, sans explications apparentes. Elles ont toutes été re-générés sans problème, gagnant pour le coup :
- 61 nouveaux calibrateurs;
- 21530 VOTables;
- 38742 diamètres.
63 autres VOTables contiennent aussi un nombre supérieur d'étoiles (entre 1 et 174 de plus !), plus 1 avec moins de sources (box_RA+07:12:00.000_DEC+51:00:00.000, qui après relance retourne le même nombre de sources qu'avant).
%ACTION{ closed="2010-03-22" closer="Main.SylvainLafrasse" created="2010-03-22" creator="Main.SylvainLafrasse" due="2010-03-26" notify="" state="closed" uid="000705" who="Main.SylvainLafrasse" }% Etudier les différences sur cette VOTable, pour savoir si ces changements viennent de chez nous (éventuellement de la mise à jour de la
LibC) ou du CDS (quels catalogues ont changé depuis?)
Non-filtrage de CAPELLA
- DanielBonneau a découvert que CAPELLA était présente de le catalogue final, malgré son référencement dans SBC9 et WDS.
- SylvainLafrasse conclut que Vizier ne retourne pas les bonnes coordonnées de CAPELLA (par contre HIP26408 répond):
- 2MASS me donne pour Capella : 05:16:51360 / +45:59:52.78
- WDS et SBC9 eux donnent : 05:16:41.30 / +45 59 56.5
- Donc les coordonnées retournées par SBC9 et WDS sont fausses de plus d'une arcseconde, interdisant tout cross-matching dans sclsvrServer.
- GillesDuvert pense que c'est un problème CDS sur les coordonnées en époque J1991.25 au lieu de J2000. Mais il n'y a, selon lui, qu'une douzaine d'étoiles dérivant d'au moins 2 arcsecondes en 10 ans, ce qui limiterai grandement l'impact de ce problème sur notre catalogue.
%ACTION{ closed="2010-05-06" closer="Main.SylvainLafrasse" created="2010-04-27" creator="Main.SylvainLafrasse" due="2010-04-30" state="closed" uid="000798" who="Main.DanielBonneau" }% Remonter le bug des époques de coordonnées sur SBC9 et WDS au CDS
Déficit d'étoiles de magnitude inférieur à 2
DanielBonneau a constaté dans son étude que 11 étoiles brillantes manquent potentiellement à l'appel:
- HIP15863
- HIP67301
- HIP26311
- HIP50583
- HIP86228
- HIP90185
- HIP97649
- HIP30438
- HIP7588
- HIP45238
- HIP61084
SylvainLafrasse a retracé pour HIP15863. Il s'avère qu'elle est filtrée car sa magnitude en J est absente.
Elimination des etoiles WDS et SBC9 par interrogation Simbad
Simbad peut indiquer pour un identifier donné l'equivalent dans d'autres catalogues. Cette technique a été utilisée pour récupéré les identifiants WDS et SB9 afin de faire l'identification croisée sans passer par les coordonnées.
Interrogation par script:
votable vot1 {
MAIN_ID
id(SBC9)
id(WDS)
}
votable open vot1
---8<---
.....une ligne par etoile....
---8<---
votable close
Il suffit donc de récuperer la colonne Name du JSDC en sortie CSV pour fabriquer le fichier.
On constate au passage qu'il y a encore certain identifiants TYCO non reconnus:
TBD Guillaume
Seule 2TYC avait un id WDS/SBC9 TBD
Ensuite on recupere une votable (après suppression manuelle des premieres lignes...) que l'on va pouvoir croiser sur le critere 'meme Index'. mais il ne faut conserver que les wds a separation =< 2 arcseconds. Il faut aussi creer une nouvelle colonne WDS_id avec seulement l'id WDS : matchGroup(ID_2b,"WDS J([0-9+-]*)")
--+++ Recuperation des WDS 'proches'
Depuis vizier, recuperer en votable les index, sep1 et sep2. Faire un subset 'sep1<=2 || sep2 <=2'
Comptabilisation finale:
Il faut ensuite croiser en mode 'All from 1' WDS_id == WDS. On recupere ainsi une table avec
38736 etoiles au total dont
38472 a garder (subset : NULL_ID_2 && NULL_WDS).
Au passage, on note la disparition de Capella qui est WDS(sep1+sep2=0 et SBC9)
--
SylvainLafrasse - 15 Mar 2010