iraf.net

I have a user who has stacked series of 2D images into a cube (e.g. a single
image from the cube is tt0066[1:12,1:12,4000]). He can display such individual
images OK but commands like imexam or fitpsf always seem to take only the first
image in the cube,i.e. tt0066[*,*,1]I have verified that, for example, trying to use :
imexam tt0066[*,*,4000]
and typing "r" to display a radial profile results in the output being for
image 1, not 4000. His only workaround is :
imcopy tt0066[*,*,4000] junk
imexam junk
which is not a suitable workaround.This is with iraf 2.10.3 beta under Solaris2.3. Do you have any suggestions? Stephen Meatheringham
Mt. Stromlo Observatory
Private Bag, Weston Creek P.O.
ACT 2611, AUSTRALIA
phone : +61 6 2490293
email : sjm@mso.anu.edu.au

Hello Stephen,What you report, that IMEXAMINE ignores any image section, is a
property of this task. This has to do with wanting to always use the
parent image coordinates and not the image section coordinates; I won't
go into the technical reasons why. It does this by removing any image
section. There are better ways to do this now and IMEXAM should
probably changed.The solution is as you noted, to copy the desired band out of the 3D
image. Rather than have user do this a simple script can be used.
Below I include such a script. A more sophisticated one could be made
with embellishments. The user would copy this to their IRAF home
directory. Then in either login.cl or loginuser.cl they can define the
task: task cubeexam=home$cubeexam.cl

The usage would then be: cl> cubeexam tt00066 7

to work on band 7. A temporary image will automatically be created
and then deleted upon exiting.Let us know if there are questions.Cheers,
Frank
=============== file cubeexam.cl =====================
# CUBEEXAM -- script to apply IMEXAM to a band of a 3D cube.
# Notes: All IMEXAM parameters are used with their current settings.
# Only one image and band may be specified at a time.
#
# usage: cubeexam <cube image> bandprocedure cubeexam (cube, band)string cube {prompt="cube to be examined"}
int band {prompt="band in cube to be examined"}begin
string im3, sec, im2 im3 = cube
im2 = mktemp (im3)
imcopy (im3 // "[*,*,"//band//"]", im2, verbose=no)
imexamine (im2)
imdelete (im2, verify=no)
end