1.1.1.0 image.coordmeasures

"""
#
print "\t----\t coordmeasures Ex 1 \t----"
ia.maketestimage(’myimage’,overwrite=true)
s = ia.shape()
for i in range(len(s)):
  s[i] = 0.5*s[i]
meas = ia.coordmeasures(s)
print meas.keys()                   # Get names of fields in record
#[’intensity’, ’mask’, ’measure’]
print meas[’intensity’]
#{’value’: 1.39924156665802, ’unit’: ’Jy/beam’}
print meas[’measure’][’direction’]
#{’type’: ’direction’,
# ’m1’: {’value’: 5.817764003289323e-05, ’unit’: ’rad’},
# ’m0’: {’value’: -5.8177644130875234e-05, ’unit’: ’rad’}, ’refer’: ’J2000’}
dir = meas[’measure’][’direction’]  # Get direction coordinate
me.doframe(me.observatory(’ATCA’))  # Set location on earth
me.doframe(me.epoch(’utc’,’16jun1999/12:30:20’))  # Set epoch
azel = me.measure(dir,’azel’)       # Convert to azimuth/elevation
print ’az,el=’, qa.angle(azel[’m0’]), qa.angle(azel[’m1’])  # Format nicely
#az,el= +105.15.47 -024.22.57
meas2=ia.coordmeasures()            # defaults to reference pixel
print meas2[’intensity’]
#{’value’: 2.5064315795898438, ’unit’: ’Jy/beam’}
print meas2[’measure’][’direction’]
#{’type’: ’direction’,
# ’m1’: {’value’: 0.0, ’unit’: ’rad’},
# ’m0’: {’value’: 0.0, ’unit’: ’rad’}, ’refer’: ’J2000’}
dir = meas2[’measure’][’direction’] # Get direction coordinate
me.doframe(me.observatory(’ATCA’))  # Set location on earth
me.doframe(me.epoch(’utc’,’16jun1999/12:30:20’))   # Set epoch
azel = me.measure(dir,’azel’)       # Convert to azimuth/elevation
print ’az,el=’, qa.angle(azel[’m0’]), qa.angle(azel[’m1’])
#az,el= +105.16.05 -024.23.00
#
"""


Inputs
pixel	Absolute pixel coordinate. Default is reference pixel.
	allowed:	doubleArray
	Default:	-1

image.coordmeasures - Function