## FITS I/O
im = readfits('input.img', hdr)
im = readfits('input.img')
im = mrdfits('input.img')
rawevent = mrdfits('input.img', 0)
rawimage = mrdfits('input.img', 0, header)
## Smoothing
sm = smooth(im, 5)
## Figures
tv, im
tvscl, im
surface, im
surface, im, ax=70, az=25
shade_surf, im
contour, im, nlevels=20
show3, im
live_surface, im
contour, im, nlevels=25, color=0, background=255, /fill
contour, im, nlevels=25, color=0, background=255, /cell_fill
xsurface, im
tvimage, im
## Circle and Box
tvcircle, 120, 768.5, 768.5, color=0, /data
tvbox, 1024, 768.5, 768.5, angle=100, color=0, /data
## Print EPS
set_plot, 'ps'
device, filename='output.eps', /port, /enc, /inches, xsize=7, ysize=7
contour, im
device, /close
set_plot, 'x'