3D Cube

ds9 can also display higher-order datasets (>2). The cube interface allows users to visualize slices through the data-cube; even playing the slices and saving as a movie.

The basic regions can be used to explore the data along the 3rd axis.

Data

RHO OPH CORE OBS_ID 635

download_chandra_obsid 635 evt2,fov

dmcopy "acisf00635N004_evt2.fits[sky=region(acisf00635_000N004_fov1.fits.gz)][bin sky=8,time=::#100]" time_cube.fits

This is a star forming region with many variable stars.

In [1]:
cd RegionFunctionality
dmlist time_cube.fits cols | grep "Axis#" --after-context=4
Group# Axis# 
   1   1,2    sky(x) = (+2791.0795) +(+8.0)* ((#1)-(+0.50))
                 (y)   (+2884.7055)  (+8.0)  ((#2) (+0.50))
   2   3      time                 = +72039163.853747 [s] +1019.7760 * (#3  -0.50)
 
--
Group# Axis# 
   1   1,2    EQPOS(RA ) = (+246.8247)[deg] +TAN[(-0.000136667)* (sky(x)-(+4096.50))]
                   (DEC)   (-24.5734 )           (+0.000136667)  (   (y) (+4096.50)) 
In [2]:
ds9 time_cube.fits -scale log &
[1] 28799
In [3]:
xpaset -p ds9 regions command "circle(3811,3456,100)"
In [4]:
xpaset -p ds9 cube 80

import -window "SAOImage ds9" ds9_04_regions_07.png 
display < ds9_04_regions_07.png

In [5]:
# Double click region, Plot 3D
# Close region info box
sleep 5

import -window "Circle" ds9_04_regions_08.png 
display < ds9_04_regions_08.png

In [6]:
xpaset -p ds9 quit
[1]+  Done                    ds9 time_cube.fits -scale log
In [ ]: