Tags:
create new tag
, view all tags
N'hésitez pas à compléter cette page, à corriger les erreurs et à le signaler...

Pour tous les programmes (sinon, le signaler), il existe l'option -help qui donne le detail technique du fonctionnement.
Par exemple : find2mass -help
Il existe aussi la page web avec une description succinte de tous les programmes : http://cdsarc.u-strasbg.fr/doc/cdsclient.html

Sommaire

Programme gsc4sim : Pour retrouver les commandes de mises à jour d'un objet à partir de ses coo ou de son id dans GSC, 2MASS, UCAC, etc.
Programme astropos : Pour convertir des coordonnées
Programme ccd : Autre programme de conversion de coordonnées à partir d'offsets
Programmes findsdss et find2mass : Pour retrouver des objets par coordonnées ou id. tronqués (ou non) dans SDSS et 2MASS.
Programme vizquery : Permet d'interroger VizieR en ligne de commande pour retrouver des SDSS, Gaia, etc.
Programme viz4sim : Permet d'éditer directement les commandes de mise à jour SIMBAD à partir d'une table VizieR lorsqu'il y a une colonne Simbad ou SimbadName et qu'on souhaite faire un =v=.
Programme simref : Pour retrouver la fiche de référence bibliographique correspondant à un bibcode dans Simbad.
Programme simbib : Pour retrouver des références dans Simbad à partir de mots du titre
Programme oldsim : Pour voir l'historique des modifications pour un objet ou une référence SIMBAD
Programme findcat : Pour retrouver des catalogues de VizieR à partir de mots clefs
Programme jdate : Pour convertir une date en JD, MJD ou UT
Logiciels TOPCAT/Aladin : Comment insérer une table qui n'est pas dans VizieR dans Aladin via TopCat ou directement

Programme gsc4sim

Le programme gsc4sim permet de générer directement les commandes de MAJ Simbad pour un objet.
Par défaut, la recherche autour d'une position ou d'un nom d'objet est effectuée dans les catalogues suivant :
tycho2,gsc1.2,ucac2,ucac3,2mass,sdss7,usnob,gsc2.2,gsc2.3,denis,cmc14

  • Pour réduire la recherche à un ou plusieurs catalogues:
    Uniquement pour une interrogation par position, on peut limiter la recherche à un ou plusieurs catalogue avec l'option -source :
    gsc4sim -source=tycho2 -c="00 41 45.46 +89 58 37.3" Effectue la recherche dans Tycho2.
    gsc4sim -source=2MASS,tycho2 -c.rs=1.5 -c="00 41 45.46 +89 58 37.3" Effectue la recherche dans Tycho2 et 2MASS dans un rayon de recherche limité à 1.5".
N.B. : Si on fait une recherche avec un nom d'objet Tycho, le résultat liste tous les catalogues sans tenir compte de la source.

  • Exemples de recherche par identificateur:
    gsc4sim 1-145-1 Affiche la MAJ pour la tycho correspondante dans les différents catalogues.
    gsc4sim 2MASS J00414546+8958373 Même chose pour l'objet 2MASS.

    Pour le catalogue SDSS, on peut spécifier dans quelle version on souhaite chercher un objet; par défaut, il est recherché dans le SDSS8:
    gsc4sim SDSS7 J124039.96+134856.2
    gsc4sim SDSS8 J124039.96+134856.2
    (FW oct12)

  • Vérifier la correspondance entre le nombre d'identificateurs retrouvés et celui donnés en entrée:
    Il est important de vérifier que le même nombre d'objets données en entrée est retrouvé en sortie (grep '^o ' et grep '^a i ' |wc -l)
    Si un objet n'est pas retrouvé, on a un "Nothing found" qu'il est important de séparer du paragraphe précédent par une ligne blanche
    sed -e '/Nothing found/ i\ ' gsc4sim2MASS.txt par exemple.

Programme astropos voir help_astropos pour comprendre comment utiliser les options

  • Pour transformer des positions en degrés ou coordonnées B1950...
Exemple pour convertir une position en degrés à 5 décimales :
Commande : astropos -v -5d J2000 J2000
Input : 10 01 25.152 +02 08 52.34
Resultat : 
Input: 150.35480 +02.14787 
                FK5 Ep=J2000.00000 Eq=J2000.00000

Output:150.35480 +02.14787 
                FK5 Ep=J2000.00000 Eq=J2000.00000

  • Convertir une position J2000 en B1950 avec en plus une "baisse" de la précision (Marianne) :
$ astropos -5 J2000 B1950

Input FK5 Position %x.y/rho.pa , mu1 mu2 [RV [plx]]: 10 01 25.152 +02 08 52.34

Input: 10 01 25.15 +02 08 52.3 
                FK5 Ep=J2000.00000 Eq=J2000.00000

Output:09 58 50.08 +02 23 20.6 
                FK4 Ep=J1949.99979 Eq=J1949.99979

  • Exemple précédent, avec un affichage plus compact (sur une seule ligne) (Marianne) :
$ astropos -5 -c1-25 J2000 B1950

Input FK5 Position %x.y/rho.pa , mu1 mu2 [RV [plx]]: 10 01 25.152 +02 08 52.34
10 01 25.152 +02 08 52.34;09 58 50.08 +02 23 20.6

  • Calcul des coordonnées d'une étoile à partir de sa position X,Y (X=-9.8 et Y=8.4) et de la position de l'étoile centrale (17 42 29.34 -28 59 12.9, B9150)
astropos B1950 J2000
Input FK4 Position %x.y/rho.pa , mu1 mu2 [RV [plx]]: 17 42 29.34 -28 59 12.9 %-9.8 8.4
Input: 17 42 28.59     -28 59 04.5     
      FK4 Ep=J1949.99979 Eq=J1949.99979

Output:17 45 39.32     -29 00 14.1     
      FK5 Ep=J2000.00000 Eq=J2000.00000

  • Calcul de coordonnées à partir d'une position de référence et d'offsets (02 20 16.603 -06 01 43.20%2.07 0.63) et affichage du résultat sur une seule ligne (-d% indique que le % est le séparateur de champ, -c1%2 indique que la position est dans le champ 1 et les offsets dans le champ 2) (Marianne):
$ astropos -5 -d% -c1%2 J2000 J2000

Input FK5 Position %x.y/rho.pa , mu1 mu2 [RV [plx]]: 02 20 16.603 -06 01 43.20%2.07 0.63 
02 20 16.603 -06 01 43.20%2.07 0.63 %02 20 16.74 -06 01 42.6

  • Calcul de coordonnées à partir d'une liste de positions de référence et d'offsets (ça marche aussi avec une liste de positions "simples"). Dans l'exemple ci-dessous, le fichier de position s'appelle listpos et contient deux lignes ; le fichier de sortie s'appelle listpos.out (Marianne) :
$ cat listpos
02 20 16.603 -06 01 43.20%2.07 0.63
02 20 16.603 -06 01 43.20%-4.76 0.3
$ astropos -5 -d% -c1%2 J2000 J2000 < listpos > listpos.out
$ cat listpos.out
02 20 16.603 -06 01 43.20%2.07 0.63%02 20 16.74 -06 01 42.6
02 20 16.603 -06 01 43.20%-4.76 0.3 %02 20 16.28 -06 01 42.9

Programme ccd

ccd [-v] [-p] [-f reference_file ] [-t#] [refcat=...  refmag=... ]
    [scale=min,max] [score=min] [sigma=min,max] [-c=...] [-c.r[b]=...]
    [asu-constraints...]

  • Ce programme calcule les coordonnées des étoiles qui ont des positions X,Y (Δα, Δδ peuvent être considérés comme des X/Y). Le calcul se fait par rapport à une liste des positions déjà connues.

  • exemple d'utilisation du programme "ccd" pour les étoiles sans coordonnées dans Simbad (BHA 1, BHA 2, BHA 4, BHA 5, BHA 7) de la réf 1984MNRAS.208..661B
 
      ccd <input_file.dat

------- - le fichier d'entrée input_file.dat contient sur chaque ligne:
--------- (a) RA Dec = X Y pour les étoiles pour lesquels on connait déjà les coordonnées RA Dec et la position X,Y
--------- (b) X Y ----------------pour les etoiles pour lesquels on connait que la position X,Y

------- - dans ce cas précis, la table 1 contient des étoiles BHA avec des positions X,Y et cross-id dans la colonne "BN number". Ce sont des étoiles :
--------- (a) qui ont déjà des coordonnées dans Simbad grâce aux cross-id. Dans le ficiher d'entrée on aura pour chaque de ces étoile une ligne avec les coordonnées RA,Dec de Simbad et la position X,Y de la table1
--------- (b) qui n'ont pas de cross-id et pas de coordonnées dans Simbad. Pour ces étoiles le programme ccd va calculer les coordonnées. Dans le fichier d'entrée on aura donc une ligne que avec les positions X et Y de la table 1

 
RA Dec = X Y
X Y


input_file.dat
17 45 39.5 -29 00 27 = -6.3 -5.0
17 45 39.389 -29 00 14.89 = -9.8 8.4
17 45 39.67 -29 00 32.4 = -3.9 -8.6
17 45 39.76 -29 00 29.8 = -3.4 -7.5
17 45 39.80 -29 00 23.9 = -3.0 -1.8
17 45 39.987 -29 00 22.24 = 0.0 0.0
17 45 39.9 -29 00 36 = 0.3 -14.5
-14.4 -17.3
-10.0 -17.3
-7.4 0.1
-6.9 -22.8
-4.1 -21.6

------- - voici la sortie du programme ccd, avec les coordonnées qui correspondent aux positions X,Y (ex. pour X,Y=-14.4 -17.3 le programme a calculé les coordonnées J2000 17 45 38.89 -29 00 40.9)

 ccd -v <input_file.dat 
.... Trial#  1: c=266.415546-29.007390 (sigma=0.492")
----Center=266.415546-29.007390 ( -3.49493  -4.16922)  1 iterations.
....O-C  1=     +0.42     -0.74 ( -6.30000  -5.00000)
....O-C  2=     +0.02     +0.05 ( -9.80000   8.40000)
....O-C  3=     +0.01     +1.35 ( -3.90000  -8.60000)
....O-C  4=     -0.62     -0.10 ( -3.40000  -7.50000)
....O-C  5=     -0.23     -0.30 ( -3.00000  -1.80000)
....O-C  6=     +0.19     +0.16 (  0.00000   0.00000)
....O-C  7=     +0.21     -0.43 (  0.30000 -14.50000)
....sigma =      0.31      0.62                0.940"/pix
-c=266.4155-29.0074 x0=-3.494932 y0=-4.169222 \
CD=4.3697e-06,4.67731e-07,5.43003e-07,4.81038e-06
266.4120-29.0113 = 17 45 38.89 -29 00 40.9 = -14.4 -17.3
266.4133-29.0112 = 17 45 39.19 -29 00 40.4 = -10.0 -17.3
266.4146-29.0063 = 17 45 39.49 -29 00 22.8 = -7.4 0.1
266.4140-29.0126 = 17 45 39.36 -29 00 45.5 = -6.9 -22.8
266.4148-29.0122 = 17 45 39.56 -29 00 44.0 = -4.1 -21.6

------- - si le programme n'a pas suffisament d'étoiles avec des coordonnées connues en entrée (exemple que 3), une erreur peut apparaître:

Can't transform pixels to position with only 3 point

Programmes findsdss et find2mass

  • Pour rechercher une liste d'objets du 2MASS :

  • find2mass JHHMM+DDMM affiche la liste des 2MASS correspondants.
  • détermination de la qualité pour les magnitudes 2MASS:
    flag_j2mass= {A,B,C} > qualité=E pour err>0.3, qualité=D pour err>0.1, sinon qualité=C
    flag_j2mass
    D > qualité=E pour err>0.3 , sinon qualité=C
    flag_j2mass
    {E,F} => qualité=E

  • Pour rechercher une liste d'objets du SDSS...

    findsdss -c "coord" -ebe -rs 1.5 -sr par exemple permet de faire une recherche dans la dernière version du SDSS sur une liste de positions à 1.5" de rayon de recherche (avec -r 1.5 le rayon est 1.5 arcmin !)
    -eb = edit basic-only (données par défaut (JHK+position pour 2MASS; un peu plus pour le SDSS) avec -ebe, on ajoute l'erreur sur la position pour le 2MASS)
    -sr = sort sur "r" => tri sur le paramètre r (séparation entre objet trouvé et position donnée en arcsec)
findsdss -i "SDSS JHHMMSS.ss+DDMMSS.s" permet une recherche sur un nom

findsdss -f liste_objid -i permet de rechercher tous les SDSS correspondants à une liste d'ObjId

En cas de non match, on a une ligne :
#SDSS 387312894892599296 [0x056003069d006800] ***NOT FOUND*** [in DR8?]

findsdss9 -f coo_sdss -eb -rs 1 -sr >| sdss_out utilise le fichier "coo_sdss" qui contient une position par ligne pour rechercher les SDSS de la DR9 correspondants dans 1".

Dans tous les cas, un espace après l'option -c ou -i est obligatoire.

Le programme fonctionne avec des SDSS au format tronqué.

Par exemple :
findsdss7 -z JHHMM+DDMM affiche la liste des SDSS, avec un redshift (option -z), qui peuvent correspondre au nom.

Le programme findsdss-sp permet une recherche dans les catalogues spectroscopiques du SDSS :

findsdss-sp "plate-mjd-fiber" permet une recherche via Plate, MJD et Fiber du SDSS. Dernière version modifiée par François le13/11/2014 avec les possibilités suivantes :
findsdss-sp 1759-53081-455 retourne le(s) Jname + versions du SDSS.
findsdss-sp 7 1759-53081-455 restreint à la recherche dans la DR7 (par exemple).
findsdss-sp 7 -f liste_Plate-MJD-Fiber permet de faire la manip' pour tout un fichier (comprenant la liste des plate-mjd-fiber) [ici limite sur DR7].
acut -c1-18 table1.dat | findsdss-sp 7 -f - si les colonnes 1 à 18 contiennent le Plate-MJD-Filber donne le même résultat que ci-dessus.

ATTENTION : La position JHHMM... donne la position photométrique (plus précise que celle en degrés qui est celle du catalogue spectro). MAIS, s'il y a une "*" après la (DRN), cela veut dire que l'objet est dans le catalogue spectro mais pas dans le photométrique, la position spectro est alors donnée avec HHMM...
Si l'objet n'est pas retrouvé, une ligne est conservée avec l'entrée et donne des "@@@@"... !

Exemple de sortie :
J001811.38+025640.0 (DR12); 004.547448+02.944454; 4300-55528-552 ; z=2.466049:0.000478; QSO BROADLINE
001811.90-004126.2 (DR12*); 004.54960-00.69062; 4219-55480-383 ; z=1.613566:0.000405; QSO BROADLINE
J001914.46+155555.9 (DR12); 004.810275+15.932201; 6185-56217-430 ; z=2.271479:0.000136; QSO BROADLINE
@@@@ (DR@@@@); @@@@ ; 4416-558281000 ; z=@@@@; @@@@
J001914.79+252113.1 (DR12); 004.811631+25.353654; 6279-56243-762 ; z=1.507042:0.000454; QSO BROADLINE

Programme vizquery

Le programme vizquery permet d'interroger les catalogues VizieR. Il y a une page d'aide complète ici : http://cdsarc.u-strasbg.fr/doc/vizquery.htx

Programme viz4sim

Le programme viz4sim permet d'obtenir un fichier de mise à jour SIMBAD à partir d'une colonne Simbad ou SimbadName de <nop>VizieR.
viz4sim "2019ApJ...881..154P" > paliya_maj
donne le fichier paliya_maj directement utilisable pour mettre à jour SIMBAD...

Programme simref

  • Pour afficher les commentaires et le bibcode d'une référence...
    simref GCN 2020, par exemple, affiche la fiche bibliographique de 2003GCN..2020....1G
    On obtient le même résultat avec simref 2003GCN..2020....1G ou simref 2003, GCN, 2020, 1
    De même, simref 2011, ApJ, 738 ou simref ApJ, 738 liste toutes les références du volume 738 de l'ApJ.

Programme simbib

  • Pour retrouver une référence de Simbad...
    simbib 2003 GREINER | less affiche les fiches bibliographiques de 2003 et des années suivantes qui ont GREINER en auteur.

Programme oldsim

  • Pour retrouver l'historique d'une référence ou d'un objet Simbad...
    oldsim -a "HD 42" affiche l'état de l'objet depuis sa création (minus une période de temps qui a été perdue).
    Voir oldsim -h pour l'aide.

Programme findcat

  • Pour retrouver un catalogue de VizieR...
    findcat Greiner | less affiche tous les catalogues de Vizier qui ont Greiner en auteur (ou dans le titre, mots-clefs et/ou abstract)

Commande jdate

Utiliser jdate -h pour voir l'ensemble des options.

Par exemple, jdate -a MJD57470.5 convertit la date MJD en date UT (DD-MMM-YYYY)

Charger une table qui n'est pas dans Vizier dans Aladin directement ou avec le logiciel TOPCAT (cf. =s=)

a) directement (voir un exemple tableegalsaladin1.odt)
Pour pouvoir charger dans Aladin une table qui n'est pas dans Vizier

  • on ajoute en tête du fichier une ligne qui contient la description des colonnes. Celle-ci doit contenir les coordonnées et l'ID (peu importe l'ordre : coordonndées ID ou ID coordonnées). Exemple RAJ2000 DEJ2000 ID ou ID Glon Glat si les coordonnées sont en galactique.
    exemple printf("RAJ2000 \t DEJ2000 \t ID \n") dans le bloc BEGIN {} du script
  • on ajoute un séparatuer (ex: tabulation \t, pipe | ) entre les colonnes (et entre les descriptions des colonnes)
    exemple printf("%s \t %s \t %s \n",RA,DE,ID)
  • dès que la table est prête on l'enregistre pour la charger ensuite dans Aladin (Aladin – Fichier - Charger un fichier local - mon_fichier) ou CTRL+o
b) via TOPCAT (voir un exemple tableegalsaladin2.odt)
  • dans un éditeur vi on transfrome le fichier en une table en format ASCII qu'on va pouvoir ensuite ouvir dans Topcat pour la transferer dans Aladin. Plus sur le format ASCII dans la doc de TOPCAT: http://www.star.bristol.ac.uk/~mbt/topcat/sun253/inAscii.html
  • on ajoute en tête du fichier une ligne qui contient la description des colonnes; cette ligne doit commencer par #
    exemple # RECNO Source RA DE SNR Flux Noise Deboosted_flux
    1 MMJ "04 54 38.96" "-03 07 39.8" 6.5 8.3 1.3 6.4+1.3−1.4
  • Les colonnes doivent être séparées par des espaces. Si dans la même colonne (exemple: la colonne RA) il y a des données séparées par des espaces (RAh RAm RAs ...), il faut mettre le tout entre double quotes.
  • dès que la table est prête on l'enregistre pour la charger ensuite dans Topcat (Topcat – File-Load Table – Format (choisir ASCII) ). La table doit être visible dans la fenêtre "Table List". Ensuite pour la charger dans Aladin : Interop – Send table to – Aladin (si on ne voit pas les objets il faut bien vérifier que la position des objets est bien interprétée par Aladin: click droit sur la table – Information sur les colonnes – RA et DE de la colonne "Coo" doivent correspondre aux colonnes qui décrivent la position dans la table)
-- EmmanuellePerret - 30 Jun 2011</verbatim>

</verbatim>

Topic attachments
I Attachment Action Size Date Who Comment
Unknown file formatodt tableegalsaladin1.odt manage 553.8 K 2011-07-04 - 15:04 MihaelaBuga Charger une table (pas dans Vizier) dans Aladin 1
Unknown file formatodt tableegalsaladin2.odt manage 528.0 K 2011-07-04 - 15:14 MihaelaBuga Charger une table (pas dans Vizier) dans Aladin 2
Topic revision: r31 - 2021-12-02 - MihaelaBuga
 
This site is powered by the TWiki collaboration platform Powered by PerlCopyright © 2008-2024 by the contributing authors. All material on this collaboration platform is the property of the contributing authors.
Ideas, requests, problems regarding TWiki? Send feedback