...
Other Use cases with Aladin are available here (and pages below): VESPA-CDS surface meeting
3D map of surface features
(Updated Aug 2017 with Aladin v9.6 , SEand java 1.8, Stéphane Erard)
Use case by Sébastien Derrière, CDS (June 2015), enlarged with comments and implementation notes
...
Summary : Build a 3D spherical HiPS off of Io, and overplot surface features from USG USGS catalogue.
This is an extension of Use Case #28 of Aladin Beta Test page (requires Aladin v9 and up):
http://aladin.u-strasbg.fr/java/Demo/AladinDemo.gm
...
• extract io_bjj_0dd.tif, convert to jpg
In the terminal, type:
> convert io_bjj_0dd.tif io_bjj_0dd.jpg
• Load jpg image in Aladin as local file
Image size is 2880 x 1440
Set In the Frame drop menu (next to Location field) to Planet
• Apply astrometric solution:
, select Planet (will plot Lon-Lat coordinates, E-handed)
• Tie a coordinate frame to the image:
Click Image > Astrometrical calibration
Use And use parameters:
Central Central pixel = 1440.0 720.0
Pixel Pixel angular size = 7.5' (= 360/2880*60)
Projection Projection = CARTESIAN
Rotation Rotation = 0°
RA RA symmetry : False
• Convert to HiPS = False (default value, this is essential)
• Build a HiPS (Hierarchical Progressive Survey) from the image:
Tool > Convert Generate a HiPS based on… > current image to a HiPS
Apply grid
It will plot after computation. Then overplot grid (Comm-G or bottom left icon), 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 out…
Projection is normally set to Spherical, but you can play with this also
Longitudes should be correctly set, ie increasing towards the east, but may display as xxE and xxW in the two hemispheres.
• Add catalogue features (but see alternative below):
Download USGS planetary nomenclature from
ftp://pdsimage2.wr.usgs.gov/pub/pigpen/nomenclature/nomenclature_all_feb2004.zip
Extract nomenclature_all_feb2004.dbf
and convert to csv with (e.g. using http://dbfconv.com/
Little cleanup of csv file, and load it in TOPCAT.
> )
Remove EOL character <NULL> from csv file in TextMate or another handy editor
• Filter Load it in TOPCAT as csv file
• Select Io features only by creating a subset:
Click on red/violet icon (subsets), add a filter (Plus icon) with parameters (name: io, expression: SA=="io")
-
Then select "Row subset" io in TOPCAT main window
=> 221 lines out of 8395 in io subset
> Click on red/violet icon, add a filter
& "plot subset" io in main window
Choose io Row subset in main TOPCAT window and broadcast Broadcast to Aladin via SAMP (Interop > Send table to…)
• Right-click the catalogue plane in Aladin, select Column Information
Click on Coord for LATITUDE and LONG360, pretending they are Information…
Define source for coordinates carefully by clicking the Coo column and selecting a value in the drop list:
LATITUDE (col 16) => DE
LONG360 (col 18) => RA
(we're telling Aladin to handle these columns as its usual DEC and RA, respectively.)
• Create a dedicated filter in Catalog > Create new filter:
In advanced Select Advanced mode, use the following expression :
> this is Catalogue menu item / create new filter / advanced modeand type / copy:
{
# 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
Clicking Apply will display names and sizes of surface features, color-coded by according to 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.
> Other pb on my side :
- long180GIS does not work either?
- Filters are not saved ?
- Graphic mix up when ploting 2 views simultaneously
- Comment: healpix only supports spherical projections, not ellipsoids
...
resolution = 1°/pixel
size= 2 * pi * R / nb pix = pi * 3390. *2 / 360. = 59.2'
Can be enlarged to 60' ?
...
Aladin 9 / java 1.6 / Mac OS: often blocking at launch, but not always
Mars image is now OK
"Save image" does not include the projected image (only the sphere frame)!
Play with the 3D model in Pan mode
• Double click on nomenclature in Aladin's stack to see the catalogue
Clicking an item will center it on the map
The feature under your pointer in the list or map in Aladin will be highlighted in the original TOPCAT table.
Alternative (from Trent Hare):
• Use more recent USGS catalogue from here:
https://planetarynames.wr.usgs.gov/SearchResults?target=IO
and click the CSV link at the bottom of the page to retrieve file.
You may have to remove initial blank lines before entering TOPCAT
• Adapt the filter to this catalogue names as:
{
# scale is 113.3 arcsec/km
draw ellipse(113.3*${Diameter}, 113.3*${Diameter}, 0) rainbow(${Approval_Date},1979,2007)
draw ${Feature_Name} rainbow(${Approval_Date},1979,2007)
}
Feb 2016, Aladin 9
• Application by Mikhail Minin with Aladin Lite:
http://epn1.epn-vespa.jacobs-university.de:8080/MARS/
...
2) for each file, you will associate a xxx.hhh calibration file (same prefix name but .hhh extension).
These calibration files will contain the WCS header (*)
...
4) your HiPS will be generated in yourFolderHiPS directory that you can load directly in Aladin,
and/or publish on your apache HTTP server.
With this method, your will be able to create and publish deep HiPS (at any resolution) rather that
just a low resolution planet map.
Notice that you can already display your HiPS in various projections (SINUS, AITOFF, MOLLWEIDE, AITOFF, CAR, ...)
via the Propertie button in Aladin.
I will keep you inform of our progress in the Aladin code (inversion of longitude, ellipsoide projection...).
Do not hesitate to signal the various problem that you see
(*) The WCS header is presently sky oriented and you will have to adapt WCS to your image collection.
For instance, if you take this unique image http://i.stack.imgur.com/ojwD8.jpg (the earth in cartesian projection),
you may use this following calibration file ojwD8.hhh:
NAXIS1 = 2048
NAXIS2 = 1024
...