iraf.net

Hi,I am trying to use imsurfit to fit a plane to data within a circular region of an image. Right now I haveimsurfit(infile, outfile, 2, 2, cross_terms="no", type_output="fit", region="circle", circle=?)The x and y coordinates and the radius of the circle need to be passed at "?". On the command line this can be done manually (e.g., circle="430 430 350") but our application requires these numbers to be read from a fits header and fed to imsurfit on the fly. The routine reads these parameters as integers but can't yet feed them to imsurfit. Please could you help?Many thanks,Gordon

The 'circle' parameter is a simple string: If you are able to construct this string within a script then there's no problem: The HSELECT task can extract keyword values, the '//' operator in the CL will concatenate strings, and if you're doing all of this from a script you should be able to easily construct a string of values derived from header values you can pass in to the 'circle' param of the task.If you post more info about what you're doing I can suggest some things, if you want the imsurfit.circle param to be able to read a keyword value, then it simply won't happen.-Mike

OK, thank you. I have something working along the lines that you suggested. I am now wondering how to record the coefficients of a polynomial surface fit using imsurfit.Thanks again,Gordon