[Date Prev][Date Next][Thread Prev][Thread Next][Date Index][Thread Index]
[gnuastro-commits] master 87f5383 23/32: astscript-radial-profile: minor
From: |
Mohammad Akhlaghi |
Subject: |
[gnuastro-commits] master 87f5383 23/32: astscript-radial-profile: minor corrections of previous commits |
Date: |
Wed, 24 Feb 2021 22:36:19 -0500 (EST) |
branch: master
commit 87f53836407fdabc71bf3e9ef40052ece521e01e
Author: Raul Infante-Sainz <infantesainz@gmail.com>
Commit: Raul Infante-Sainz <infantesainz@gmail.com>
astscript-radial-profile: minor corrections of previous commits
With this commit, several minor corrections from previous commits have
been included. They are mostly related with comments and not actually
code itself.
---
bin/script/radial-profile.in | 43 +++++++++++++++++++++++++------------------
1 file changed, 25 insertions(+), 18 deletions(-)
diff --git a/bin/script/radial-profile.in b/bin/script/radial-profile.in
index e243369..38688aa 100644
--- a/bin/script/radial-profile.in
+++ b/bin/script/radial-profile.in
@@ -368,8 +368,8 @@ astmkprof $aperturestxt --background=$inputs --backhdu=$hdu
\
# Detection signal
# ----------------
#
-# To obtain signal to noise ratio for each radial profile; at first we have
-# to detect signal of noise.
+# To obtain signal to noise ratio radial profile it is necessary to have the
+# detection map. For now, it uses the default option of NoiseChisel.
detection=$(echo $output | sed -e"s|.fits|_detected.fits|g")
astnoisechisel $inputs -o$detection
@@ -380,15 +380,16 @@ astnoisechisel $inputs -o$detection
# Obtain the radial profile
# -------------------------
#
-# The radial profile is obtained using Catalog. In practice, what is done
-# is to obtain a catalogue using the segmentation image previously
-# generated (the elliptical apertures) and the original input image for
-# computing the values. Also using of SKY_STD extinction of noisechisel
-# output for computing signal to noise ratio.
+# The radial profile is obtained using Catalog. In practice, what is done is to
+# obtain a catalogue using the profile image previously generated (the
+# elliptical apertures) and the original input image for computing the values.
+# Also using of SKY_STD extension of NoiseChisel output for computing signal
+# to noise ratio.
fprofile=$(echo $output | sed -e"s|.fits|_catalogue.fits|g")
astmkcatalog $aperturesfits -h1 --valuesfile=$inputs --valueshdu=$hdu \
- --instd=$detection --stdhdu=SKY_STD --ids --$m --sn \
- --config=$c -o$fprofile
+ --instd=$detection --stdhdu=SKY_STD \
+ --config=$c -o$fprofile \
+ --ids --$m --sn
@@ -397,15 +398,13 @@ astmkcatalog $aperturesfits -h1 --valuesfile=$inputs
--valueshdu=$hdu \
# Binning data
# ------------
#
-# In order to increase the signal-to-noise ratio of the radial profile, it
-# is possible to bin the data. It is done in any case, if the user set the
-# default value of binning variable($b) except 1, because the binning will
-# be equal to 1 so the output binned will be the same as the input. To do
-# the binning of the data, a small Awk script is used.Because of the
-# behavior of signal to noise ratio isn't linear. We used the different
-# way to do binning the S/N column. Since the Awk script will print the
-# columns as float values, it is necessary to change the headers. To do
-# that, Sed is used to replace all ocurrences of i32 to f32.
+# In order to increase the signal-to-noise ratio of the radial profile, it is
+# possible to bin the data. If the user specify the binning parameter (`-b' or
+# `--binning') different of 1. To do the binning of the data, a small Awk
+# script is used. Because of the behavior of signal to noise ratio is not
+# linear, a different way of binning the S/N column is used. Since the Awk
+# script will print the columns as float values, it is necessary to change the
+# headers. To do that, Sed is used to replace all ocurrences of i32 to f32.
bprofile=$(echo $output | sed -e"s|.fits|_binned.fits|g")
if [ $b != 1 ]; then
@@ -434,6 +433,7 @@ fi
+
# Modify the radial profile
# -------------------------
#
@@ -487,7 +487,11 @@ asttable $bprofile \
cat $aprofile | asttable -o$output
+
+
+
# Calculate the Half-Light Radii
+# ------------------------------
#
# The half-light or 'effective' radius as the radius within which half of
# the object's luminosity is contained.
@@ -497,6 +501,9 @@ radius=$(asttable $output \
echo "Half-Light-Radii: $radius pixels from the center of the object"
+
+
+
# Remove temporal files
# ---------------------
#
- [gnuastro-commits] master aec3057 13/32: Book: describing most important options of radial profile script, (continued)
- [gnuastro-commits] master aec3057 13/32: Book: describing most important options of radial profile script, Mohammad Akhlaghi, 2021/02/24
- [gnuastro-commits] master 1b505d1 21/32: Binning the S/N column, Mohammad Akhlaghi, 2021/02/24
- [gnuastro-commits] master 4294506 32/32: astscript-radial-profile: polished the script and book, Mohammad Akhlaghi, 2021/02/24
- [gnuastro-commits] master ad7fb58 26/32: Central coordinates of the cropped image, Mohammad Akhlaghi, 2021/02/24
- [gnuastro-commits] master d085379 22/32: Radial-profile script documentation, Mohammad Akhlaghi, 2021/02/24
- [gnuastro-commits] master 8089ea3 27/32: astscript-radial-profile: Samaeh's changes included, conflicts fixed, Mohammad Akhlaghi, 2021/02/24
- [gnuastro-commits] master 04a8bec 29/32: Book: simplifying documentation of astscript-radial-profile, Mohammad Akhlaghi, 2021/02/24
- [gnuastro-commits] master 362819c 30/32: astscript-radial-profile: fixed bug when creating the temporal directory, Mohammad Akhlaghi, 2021/02/24
- [gnuastro-commits] master c9647b5 31/32: Imported work on the radial profile script, minor conflict fixed, Mohammad Akhlaghi, 2021/02/24
- [gnuastro-commits] master 1072b2b 02/32: astscript-radial-profile: added the option -b for binning the data, Mohammad Akhlaghi, 2021/02/24
- [gnuastro-commits] master 87f5383 23/32: astscript-radial-profile: minor corrections of previous commits,
Mohammad Akhlaghi <=
- [gnuastro-commits] master e41c2e4 24/32: astscript-radial-profile: sigma clip parameters for measurements added, Mohammad Akhlaghi, 2021/02/24
- [gnuastro-commits] master 337e46d 18/32: The half-light radii, Mohammad Akhlaghi, 2021/02/24
- [gnuastro-commits] master 2df5a22 19/32: Signal to noise ratio for each radius, Mohammad Akhlaghi, 2021/02/24
- [gnuastro-commits] master 846fc39 25/32: Radial profile script modifications, Mohammad Akhlaghi, 2021/02/24
- [gnuastro-commits] master 8e5ab9e 28/32: astscript-radial-profile: simplifying the script, Mohammad Akhlaghi, 2021/02/24