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.
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.
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))
ds9 time_cube.fits -scale log &
[1] 28799
xpaset -p ds9 regions command "circle(3811,3456,100)"
xpaset -p ds9 cube 80
import -window "SAOImage ds9" ds9_04_regions_07.png
display < ds9_04_regions_07.png
# 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
xpaset -p ds9 quit
[1]+ Done ds9 time_cube.fits -scale log