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 ( 2 juin 2010 )
Simbad peut indiquer pour un identifiant donné l'équivalent 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.
Nous avons modifier la generation du nom pour le nom TYC en integrant les trois identifiants TYCO. Cela elimine 6 etoiles supplementaire (elles avaient un identifiant TYC2 mais pas de TYC1 ou TYC3)
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écupérer la colonne Name du JSDC en sortie CSV pour fabriquer le fichier.
On constate au passage qu'il y a encore certains identifiants TYCO non reconnus:
TYC 5274-1318-1
TYC 6425-1838-1
TYC 4328-1781-1
TYC 2880-1969-1
TYC 9161-1289-1
TYC 4751-968-1
TYC 114-1537-1
TYC 2933-1394-1
TYC 154-347-1
TYC 4803-540-1
TYC 8931-1806-1
TYC 8211-279-1
TYC 8976-2499-1
TYC 5619-1257-1
TYC 7384-264-1
TYC 8751-1567-1
TYC 6273-1114-1
TYC 3544-222-1
TYC 7940-1403-1
TYC 2166-91-1
TYC 3559-1552-1
TYC 1618-829-1
Aussi parmi les TYCXXX , seuls 2TYC avaientt un id WDS:
- HD 60683 ( TYC 6543-261-1) WDS J07347-2602 impossible a retrouver dans WDS???
- LTT12714 ( TYC 833-134-1 ) WDS J10084+1158BC avec min_sep1=4.0 min_sep2 = 2.5
En sortie du script on recupere une votable (après suppression manuelle des premieres lignes...) que l'on va pouvoir croiser sur le critère 'même Index' de laquelle il ne faut conserver que les wds à separation =< 2 arcseconds.
--+++ Recuperation des WDS 'proches'
Depuis vizier, recuperer en votable les index, sep1 et sep2. Faire un subset 'sep1<=2 || sep2 <=2'
Comptabilisation finale:
Afin de croiser les entrees WDS, il faut créer une nouvelle colonne WDS_id avec seulement l'id WDS : matchGroup(ID_2b,"WDS J([0-9+-]*)")
Enfin, on peut croiser en mode 'All from 1' WDS_id == WDS. On recupere ainsi une table avec
38736 etoiles au total dont
38472 a garder (soit 264 etoiles en moins avec le subset : NULL_ID_2 && NULL_WDS).
Au passage, on note la disparition de Capella qui est WDS(sep1+sep2=0 et SBC9)
Memo des commandes:
# aggregation des identifiants SBC9 et WDS:
stilts tmatch2 out='JSDC_v1.0b4TMP1.fits' matcher='exact' values1='$0' values2='$0' join='all1' in1='JSDC_v1.0b4.fits' in2='post_result.vot'
# calcul de l'identifiant WDS pour comparaison:
stilts tpipe in='JSDC_v1.0b4TMP1.fits' out='JSDC_v1.0b4TMP2.fits' cmd='addcol WDS_id matchGroup(ID_2a,\"WDS\ J([0-9+-]*)\")'
# selection des WDS a supprimer:
stilts tpipe in='WDS.vot' out='WDSToRemove.vot' cmd='select (sep1<=2\ ||\ sep2<=2)'
# croisement avec les WDS a supprimer:
stilts tmatch2 out='JSDC_v1.0b4TMP3.fits' matcher='exact' values1='WDS_id' values2='WDS' join='all1' in1='JSDC_v1.0b4TMP2.fits' in2='WDSToRemove.vot'
# selection finale:
stilts tpipe in='JSDC_v1.0b4TMP3.fits' out='JSDC_v1.0b4TMP4.fits' cmd='select (NULL_ID_2\ &&\ NULL_WDS)' cmd='keepcols Name\ RAJ2000\ DEJ2000\ pmRA\ pmDEC\ Bmag\ Vmag\ Rmag\ f_Rmag\ Imag\ f_Imag\ Jmag\ Hmag\ Kmag\ LDD\ e_LDD\ UDDB\ UDDV\ UDDR\ UDDI\ UDDJ\ UDDH\ UDDK\ plx\ e_plx\ SpType\ Teff_SpType\ logg_SpType'
# affichage du nombre d'etoiles:
stilts tcat in=JSDC_v1.0b4TMP4.fits omode=count
--
SylvainLafrasse - 15 Mar 2010