 |
| FIGURE 1: Processing of part of a PFC image of NGC 3664. The run script is in T_runs/momcat/ex03_n3664. |
Using and image and some corresponding pixel mask (like that from ccl4) compute image moments for every unique source identified in the mask.
% momcat.sh a.fits mask.fits 2350 arg1 - name of input fits image arg2 - name of mask image arg3 - background level to be subtracted
|
| FIGURE 1: Processing of part of a PFC image of NGC 3664. The run script is in T_runs/momcat/ex03_n3664. |
Sometime we may want to inspecte the objects that have been processed with momcat (i.e. those targets selected by a run of something like ccl4). The momcat routine will prepare a (and make executable) a script that can be used to prepare such stamps by creating a local file called "stamp.cut":
% ls RUN* S/ stamp.cut % pwd /home/sco/sco/projects/Test_Data_for_Codes/T_runs/momcat/ex03_n3664 % ls RUN* S/ stamp.cut % cat stamp.cut 100 % RUN 50 0 Running momcal check with N3664 image Row 1052 extracted t2.fits Writing FITS image: out.fits Writing FITS image: pass1.fits Writing FITS image: pass2.fits Number of stamps you can cut = 9 Name of the cutting script = Stamp_Cut.sh I am making Stamp_Cut.sh an executable! % ls Sta* Stamp_Cut.sh*So, because I put the value "100" into the file named "stamp.cut", I created an executable file called Stamp_Cut.sh that will cut stamps for 9 targets, each target having been composed of at least 100 pixels in the mask image fed to momcat.
After running Stamp_Cut.sh you will be able to run XYP, with the file momcat_xyp.dat as input. You can plot user-specified sets of data, but more importantly in thsi context, you can use the "I" command to see information about any point. If you have opened a ds9 gui prior to running XYP, then an image of the stamp linked to that point will be diplayed. This is a powerful way to quickly understand your image data. Near the bottom of the XYP page is an illustration of how this is done.