im := image('my_image')
cs := im.coordsys()  # gets a copy of the coordinate system
cs.setrestfrequency(1420.4058e6)  # set the rest frequency here, in Hz)
im.setcoordsys(cs)  # and put that fixed copy back in the image
im.done()