"""  
#  
print "\t----\t velocitytofrequency Ex 1 \t----"  
ia.fromshape(’hcn.cube’,[64,64,32,4], overwrite=true)  
csys = ia.coordsys()  
rtn = csys.findcoordinate(’spectral’)   # Find spectral axis  
pixel = csys.referencepixel();          # Use reference pixel for non-spectral  
pa = rtn[’pixel’]  
wa = rtn[’world’]  
nFreq = ia.shape()[pa]                  # Length of spectral axis  
freq = []  
for i in range(nFreq):  
  pixel[pa] = i;                        # Assign value for spectral axis of pixel coordinate  
  w = csys.toworld(value=pixel, format=’n’)# Convert pixel to world  
  freq.append(w[’numeric’][wa])         # Fish out frequency  
print "freq=", freq  
vel = csys.frequencytovelocity(value=freq, doppler=’optical’, velunit=’km/s’)  
freq2 = csys.velocitytofrequency(value=vel, doppler=’optical’, velunit=’km/s’)  
print "vel=",vel  
print "freq2=",freq2  
csys.done()  
#  
exit() # This is last example so exit casapy if you wish.  
#  
"""  
 
 
In this example, we find the optical velocity in km/s of every pixel  
along the spectral axis of our image.  First we  obtain the Coordinate  
System from the image.  Then we find which axis of the Coordinate System  
(image) pertain to the spectral coordinate.  Then we loop over each  
pixel of the spectral axis, and convert a pixel coordinate (one for each  
axis of the image) to world.  We obtain the value for the spectral axis  
from that world vector, and add it to the vector of frequencies.  Then  
we convert that vector of frequencies to velocity.  Then we convert it  
back to frequency.  They better agree.