Quote by: fitz
We don't (yet) have a a generalized cutout tool, the closest you might find in distribution is the ACECUTOUTS task in the NFEXTERN.ACE package (see the help page for details).
In general, the problem involves converting the RA/DEC coords to image coords and then simply copying out the pixels. This can be done in some cases with a script using tasks like WCSCTRAN and IMCOPY or with a compiled task, but it's starts to get messy depending on whether you need to deal with registration, rotation, image extensions and so on.
So I've got the X/Y coords of the objects (converted from RA/DEC using the header) and have aligned the frames together with WCS. I'm currently trying a pyraf script to imcopy a 150 by 150 pixel region. Not being very familiar with pyraf or python, so far I've got this simple script:
Are there any quick changes/additions you could suggest through which I could feed in a text file containing all the X and Y of the objects and get the regions out? I seem to run into a brick wall when I try to put in a general defined variable within imcopy rather than a number like this:
I'm not quite sure why this gives an error even though its essentially doing the same as above. The other issue is that I haven't yet found a way to tell pyraf within the script what to name the output file. It always prompts me to name the file on the terminal after I've run the script. Sorry for the long question and thank you for the help.