[Top][All Lists]

[Date Prev][Date Next][Thread Prev][Thread Next][Date Index][Thread Index]

[task #15637] Match RA and Dec catalog to X and Y catalog to find WCS

From: Mohammad Akhlaghi
Subject: [task #15637] Match RA and Dec catalog to X and Y catalog to find WCS
Date: Mon, 11 May 2020 22:21:44 -0400 (EDT)
User-agent: Mozilla/5.0 (X11; Linux x86_64; rv:76.0) Gecko/20100101 Firefox/76.0

Follow-up Comment #1, task #15637 (project gnuastro):

To help in completing this task, I just uploaded some data to play with. A
link to each uploaded file is available at the bottom of this comment.

The main datasets are 8 short exposure (5 sec) images of Castor
<>, but each exposure is offset
compared to the others (this offset in individual exposures is called
dithering in astronomy). Please open these images and actually see how
Castor's position in the image changes in each of them.

A catalog of sources (actually "Clumps" in Segment
<>) is
generated for each image using NoiseChisel+Segment+MakeCatalog using the
script below (which is also available under that directory), the catalogs are
the actual inputs into this task and they have a `-x-y.fits' suffix. Each
catalog contains roughly 1500 clumps. Finally, there is also a reference
catalog with RA and Dec of 1179744 sources near Castor from ESO's Gaia

Once this task is complete, we should be able to have an accurate WCS
(including distortions) for each image.

Script to generate a catalog of sources in each image and their X and Y

# Base name of input files.
input="1881731715 3163025249 1887227377 3746244531 2923938133 4209692437
2966790193 595260509"

# Generate the catalog for each image.
for i in $input; do

    # First run NoiseChisel to separate signal from the background.
    astnoisechisel $i.fits --tilesize=50,50 --interpnumngb=21 \

    # Run Segment to get a label for each clump, don't bother with
    # detecting objects, they aren't relevant here.
    astsegment $i-nc.fits --onlyclumps --output=$i-seg.fits

    # Generate a catalog with the ID, X, Y, and magnitude (assuming a
    # zeropoint of 0).
    astmkcatalog $i-seg.fits --hdu=CLUMPS --ids --x --y --magnitude \

    # Clean up.
    rm $i-nc.fits $i-seg.fits

Reference Catalog with RA and Dec of many sources:

List of FITS images:

List of X-Y catalogs:

Script to generate catalog:


Reply to this item at:


  Message sent via Savannah

reply via email to

[Prev in Thread] Current Thread [Next in Thread]