You are viewing an old version of this page. View the current version.

Compare with Current View Page History

Version 1 Next »

 

 

Introduction

This starts with a use case by Sébastien Derrière (June 2015), enlarged with comments and implementation notes

 

Use case

Stéphane Erard comments are introduced by >

Summary : it is possible to convert planetary or satellite
to HiPS, in order to visualize them in Aladin, explore the
surface feature, and even overlay catalogues of features.

Introduction :
This is an extension of Use Case #28 of Aladin Beta Test page :
http://aladin.u-strasbg.fr/java/Demo/AladinDemo.gml

Step-by-step :
Retrieve image of planetary map, for example from here :
http://laps.noaa.gov/albers/sos/
or here
ftp://pdsimage2.wr.usgs.gov/pub/pigpen/

As example, work with Io
ftp://pdsimage2.wr.usgs.gov/pub/pigpen/io/io_global_images.zip
extract io_bjj_0dd.tif, convert to jpg
> convert io_bjj_0dd.tif io_bjj_0dd.jpg
    (in the terminal)

Load image in Aladin as local file.

Apply astrometric solution :
Image > Astrometrical calibration
Central pixel 1440.0 720.0, with pix ang size 7.5' and projection CARTESIAN

Convert to HiPS :
Tool > Convert current image to a HiPS
Apply grid, pan around, zoom in/out ...
    > Type Comm-G to overplot coordinate grid
    > Beware that  longitudes are handled as Righ Ascensions (inverted)


Adding catalogue :
Planetary nomenclature downloaded from
ftp://pdsimage2.wr.usgs.gov/pub/pigpen/nomenclature/nomenclature_all_feb2004.zip
Extract nomenclature_all_feb2004.dbf
convert to csv with http://dbfconv.com/

Little cleanup of csv file, and load it in TOPCAT.
    Virer le caractère de fin de ligne: <NULL> dans TextMate

    
Filter Io features only by creating subset (io, expression SA=="io")
-> 221 lines out of 8395 in io subset
    (Cliquer sur l'icone rouge/violette, ajouter un filtre
    + "Afficher le subset" io dans la fen principale)


Choose io Row subset in main TOPCAT window and broadcast to Aladin
via SAMP.

Right-click catalogue plane in Aladin, Column Information
Click on Coo for LATITUDE and LONG360, pretending they are DEC and RA,
respectively.

Create dedicated filter :
In advanced mode, use the following expression :

    (ds menu Catalogue / Créer nv filtre / advanced mode)

{
# scale is 113.3 arcsec/km
draw ellipse(113.3*${DIAM}, 113.3*${DIAM}, 0) rainbow(${AD},1979,2003)
draw ${NAME} rainbow(${AD},1979,2003)
}
Displays names and sizes of surface features, color-coded by year of discovery.

Remaining issues :
Defaut orientation of RA or longitudes in Aladin is different
from the one used in planetary science, therefore the surface
is viewed as from within the planetary body, not from outside.




• Autre pb chez moi :
long180GIS marche pas non plus ;
les cercles se déforment en arrivant au limbe ;
ds certains cas, les cercles s'affichent partiellement dans le ciel : disque trop petit ?

C'est décalé de 8h et N/S : Loki P. qui est à 20:30h +13N se retrouve à 4h -14S en healpix
Donc: l'image est mirorée en X et Y au moment du calcul de healpix (par rapport à la version calibrée en astrométrie). Pas la faute à java, ça part bien en v8 si je fais
    open /Applications/ApplisAstro/Aladin\ 8.0/Aladin.app/Contents/Resources/Java/Aladin.jar &


Loki Patera is located at 13°N 308.8°W[8]  ~ 20.6 h  - donc c'est la bonne soluton astrom, le pb st dans la projection Healpix.
(me les donne ds l'affichage catalogue après un select sur l'image)
Coordinates: 13°N 308.8°W[8].

Conclu:
- Pb avec calcul Hpx : l'image est retournée en X et Y. Ça se produit en java 1.6 (avec message d'info) ou 1.8
- Par contre, Catalogue OK, les objets sont aux bonnes coordonnées
- Se produit aussi en lançant avec 1Go de mem
- Les filtres sauvés ne sont pas relus.
- Il y a des fantomes à l'affichage en sphérique (limites de cratères déformées)
- L'affichage s mélange les pinceaux quand on passe d'une vue à deux (les limites en sphériques sont superposés à la carte cylindrique.


Soluce:
Le tuto marche avec la version 8 beta (d'après 8.162) dispo sur la page web mentionné plus haut.



Mars IRTM:
convert alb_.09.tif alb_.09.jpg
Taille pi * R / nb pix = pi * 3390. / 360. = 29.6'
Pole S mal placé (en bas de la carte, qui n'est pas complète)
+ fait pas le tour (mais devrait)
Faut faire une inversion pour avoir N/E correctement affichés
Non, en fait il faut prendre pi * 3390. *2 / 360. = 59.2'
Toujours inversé, et pôle mal placé
Un peu petit : prendre 60' ?
 en cyl, met l'E en biais
Pas mieux avec rotation.
Rotation = 180 ou inversion donnent toujours une image inversée sur la sphère.



En v9, Fev 2016 :
Sur nadja, en java 1.6 - attention, ça bloque de temps en temps au démarrage.
l'ex Mars est OK, image dans le bon sens.
Par contre, la sauvegarde image n'intègre pas l'image projetée !


  • No labels