[task #16377] Log-scale and custom radial profiles

[Mohammad Akhlaghi]

URL:
  <https://savannah.gnu.org/task/?16377>

                 Summary: Log-scale and custom radial profiles
                   Group: GNU Astronomy Utilities
               Submitter: makhlaghi
               Submitted: Thu 08 Jun 2023 04:22:30 PM BST
         Should Start On: Thu 08 Jun 2023 12:00:00 AM BST
   Should be Finished on: Thu 08 Jun 2023 12:00:00 AM BST
                Category: Installed scripts
                Priority: 5 - Normal
              Item Group: Enhancement
                  Status: None
                 Privacy: Public
        Percent Complete: 0%
             Assigned to: None
             Open/Closed: Open
         Discussion Lock: Any
                  Effort: 0.00


    _______________________________________________________

Follow-up Comments:


-------------------------------------------------------
Date: Thu 08 Jun 2023 04:22:30 PM BST By: Mohammad Akhlaghi <makhlaghi>
As of version 0.20, 'astscript-radial-profile' will only produce a linear
profile (where the distances between radial bins are constant). But we can
implement custom-scale radial profiles also (where the spacing is different).
For example a log-scale radial profile!

As a demonstration, with the commands below, we'll generate a Sersic profile,
and generate its radial profile (while keeping the temporary directory). With
the following 'astarithmetic' command, we will generate an image where the
radial annuli have a log-scale distance between them:


echo "1 500 500 1 100 2.5 0 1 5 10" \
     | astmkprof --mergedsize=999,999 --oversample=1 \
                 --output=sersic.fits

astscript-radial-profile sersic.fits --keeptmp \
         --tmpdir=tmp-radial --output=profile.fits

astarithmetic tmp-radial/apertures.fits log10 10 x uint8 \
              --output=apertures-log.fits

astscript-fits-view --ds9scale=minmax \
                    tmp-radial/apertures.fits apertures-log.fits 


The extra 'astarithmetic' command can be added when the user give the option
'--scale=log' to the radial profile script. The 'apertures-log.fits' file can
then be given to MakeCatalog to generate the log-scaled radial profile.
Afterwards (on the output table), we can do the inverse of the operation above
(divide by 10, and put to the power of 10) to recover the original radius of
each annulus.

Further more, because MakeProfiles accepts custom profiles (see the
description of '--customtable' in
[https://www.gnu.org/software/gnuastro/manual//html_node/MakeProfiles-profile-settings.html
MakeProfiles profile settings], users can optionally give
'--scale=custom.fits' to define any custom method of binning they desire (not
just log-scale!).






    _______________________________________________________
File Attachments:


-------------------------------------------------------
Date: Thu 08 Jun 2023 04:22:30 PM BST  Name: radial-profile-linear-log.png 
Size: 75KiB   By: makhlaghi

<http://savannah.gnu.org/task/download.php?file_id=54824>

    _______________________________________________________

Reply to this item at:

  <https://savannah.gnu.org/task/?16377>

_______________________________________________
Message sent via Savannah
https://savannah.gnu.org/