[packages/gnuastro] - updated to 0.21
qboosh
qboosh at pld-linux.org
Sun Jan 14 13:40:03 CET 2024
commit c2c679f9d8b5530360bc136f40c2e07124bb87cc
Author: Jakub Bogusz <qboosh at pld-linux.org>
Date: Sun Jan 14 12:47:15 2024 +0100
- updated to 0.21
gnuastro-info.patch | 796 ++++++++++++++++++++++++++++++----------------------
gnuastro.spec | 16 +-
2 files changed, 465 insertions(+), 347 deletions(-)
---
diff --git a/gnuastro.spec b/gnuastro.spec
index e330281..3ee8934 100644
--- a/gnuastro.spec
+++ b/gnuastro.spec
@@ -8,12 +8,12 @@
Summary: GNU Astronomy Utilities
Summary(pl.UTF-8): Narzędzia astronomiczne GNU
Name: gnuastro
-Version: 0.19
-Release: 3
+Version: 0.21
+Release: 1
License: GPL v3+
Group: Applications/Science
Source0: https://ftp.gnu.org/gnu/gnuastro/%{name}-%{version}.tar.lz
-# Source0-md5: 8b65b9bfd1541a7f47e039e5ba62af21
+# Source0-md5: 905c9f52dcadbb3832c90da514a8f7f3
Patch0: %{name}-info.patch
Patch1: ac.patch
Patch2: %{name}-link.patch
@@ -22,7 +22,7 @@ URL: http://www.gnu.org/software/gnuastro/
BuildRequires: autoconf >= 2.69
BuildRequires: automake
BuildRequires: bzip2-devel
-BuildRequires: cfitsio-devel
+BuildRequires: cfitsio-devel >= 3.30
BuildRequires: curl-devel
BuildRequires: ghostscript >= 9.10
BuildRequires: gsl-devel >= 2.0
@@ -159,7 +159,6 @@ rm -rf $RPM_BUILD_ROOT
%attr(755,root,root) %{_bindir}/astfits
%attr(755,root,root) %{_bindir}/astmatch
%attr(755,root,root) %{_bindir}/astmkcatalog
-%attr(755,root,root) %{_bindir}/astmknoise
%attr(755,root,root) %{_bindir}/astmkprof
%attr(755,root,root) %{_bindir}/astnoisechisel
%attr(755,root,root) %{_bindir}/astquery
@@ -170,14 +169,16 @@ rm -rf $RPM_BUILD_ROOT
%attr(755,root,root) %{_bindir}/astscript-psf-stamp
%attr(755,root,root) %{_bindir}/astscript-psf-subtract
%attr(755,root,root) %{_bindir}/astscript-psf-unite
+%attr(755,root,root) %{_bindir}/astscript-pointing-simulate
%attr(755,root,root) %{_bindir}/astscript-radial-profile
%attr(755,root,root) %{_bindir}/astscript-sort-by-night
+%attr(755,root,root) %{_bindir}/astscript-zeropoint
%attr(755,root,root) %{_bindir}/astsegment
%attr(755,root,root) %{_bindir}/aststatistics
%attr(755,root,root) %{_bindir}/asttable
%attr(755,root,root) %{_bindir}/astwarp
%attr(755,root,root) %{_libdir}/libgnuastro.so.*.*.*
-%attr(755,root,root) %ghost %{_libdir}/libgnuastro.so.17
+%attr(755,root,root) %ghost %{_libdir}/libgnuastro.so.19
%dir %{_sysconfdir}/gnuastro
%config(noreplace) %verify(not md5 mtime size) %{_sysconfdir}/gnuastro/ast*.conf
%config(noreplace) %verify(not md5 mtime size) %{_sysconfdir}/gnuastro/gnuastro.conf
@@ -193,12 +194,12 @@ rm -rf $RPM_BUILD_ROOT
%{_mandir}/man1/astfits.1*
%{_mandir}/man1/astmatch.1*
%{_mandir}/man1/astmkcatalog.1*
-%{_mandir}/man1/astmknoise.1*
%{_mandir}/man1/astmkprof.1*
%{_mandir}/man1/astnoisechisel.1*
%{_mandir}/man1/astquery.1*
%{_mandir}/man1/astscript-ds9-region.1*
%{_mandir}/man1/astscript-fits-view.1*
+%{_mandir}/man1/astscript-pointing-simulate.1*
%{_mandir}/man1/astscript-psf-scale-factor.1*
%{_mandir}/man1/astscript-psf-select-stars.1*
%{_mandir}/man1/astscript-psf-stamp.1*
@@ -206,6 +207,7 @@ rm -rf $RPM_BUILD_ROOT
%{_mandir}/man1/astscript-psf-unite.1*
%{_mandir}/man1/astscript-radial-profile.1*
%{_mandir}/man1/astscript-sort-by-night.1*
+%{_mandir}/man1/astscript-zeropoint.1*
%{_mandir}/man1/astsegment.1*
%{_mandir}/man1/aststatistics.1*
%{_mandir}/man1/asttable.1*
diff --git a/gnuastro-info.patch b/gnuastro-info.patch
index 70f207a..8b6c5c2 100644
--- a/gnuastro-info.patch
+++ b/gnuastro-info.patch
@@ -1,8 +1,11 @@
---- gnuastro-0.19/doc/gnuastro.texi.orig 2022-10-24 02:00:42.000000000 +0200
-+++ gnuastro-0.19/doc/gnuastro.texi 2023-01-19 21:52:35.888478772 +0100
-@@ -50,74 +50,74 @@ A copy of the license is included in the
+--- gnuastro-0.21/doc/gnuastro.texi.orig 2023-10-20 17:11:46.000000000 +0200
++++ gnuastro-0.21/doc/gnuastro.texi 2024-01-14 08:00:13.150076828 +0100
+@@ -47,75 +47,75 @@ A copy of the license is included in the
+ @end copying
+
@c To include in the info directory.
- @dircategory Astronomy
+- at dircategory Astronomy
++ at dircategory Astronomy:
@direntry
-* Gnuastro: (gnuastro). GNU Astronomy Utilities.
-* libgnuastro: (gnuastro)Gnuastro library. Full Gnuastro library doc.
@@ -10,7 +13,7 @@
+* libgnuastro: (gnuastro)Gnuastro library. Full Gnuastro library doc
-* help-gnuastro: (gnuastro)help-gnuastro mailing list. Getting help.
-+* help-gnuastro: (gnuastro)help-gnuastro mailing list. Getting help
++* help-gnuastro: (gnuastro)help-gnuastro mailing list. Getting help
-* bug-gnuastro: (gnuastro)Report a bug. How to report bugs
+* bug-gnuastro: (gnuastro)Report a bug. How to report bugs
@@ -27,8 +30,8 @@
-* ConvertType: (gnuastro)ConvertType. Convert different file types.
-* astconvertt: (gnuastro)Invoking astconvertt. Options to ConvertType.
-+* ConvertType: (gnuastro)ConvertType. Convert different file types
-+* astconvertt: (gnuastro)astconvertt. Options to ConvertType
++* ConvertType: (gnuastro)ConvertType. Convert different file types.
++* astconvertt: (gnuastro)astconvertt. Options to ConvertType.
-* Convolve: (gnuastro)Convolve. Convolve an input file with kernel.
-* astconvolve: (gnuastro)Invoking astconvolve. Options to Convolve.
@@ -42,7 +45,7 @@
-* Crop: (gnuastro)Crop. Crop region(s) from image(s).
-* astcrop: (gnuastro)Invoking astcrop. Options to Crop.
-+* Crop: (gnuastro)Crop. Crop region(s) from image(s)
++* Crop: (gnuastro)Crop. Crop region(s) from image(s).
+* astcrop: (gnuastro)astcrop. Options to Crop
-* Fits: (gnuastro)Fits. View and manipulate FITS extensions and keywords.
@@ -55,11 +58,6 @@
+* MakeCatalog: (gnuastro)MakeCatalog. Make a catalog from labeled image
+* astmkcatalog: (gnuastro)astmkcatalog. Options to MakeCatalog
--* MakeNoise: (gnuastro)MakeNoise. Make (add) noise to an image.
--* astmknoise: (gnuastro)Invoking astmknoise. Options to MakeNoise.
-+* MakeNoise: (gnuastro)MakeNoise. Make (add) noise to an image
-+* astmknoise: (gnuastro)astmknoise. Options to MakeNoise
-
-* MakeProfiles: (gnuastro)MakeProfiles. Make mock profiles.
-* astmkprof: (gnuastro)Invoking astmkprof. Options to MakeProfiles.
+* MakeProfiles: (gnuastro)MakeProfiles. Make mock profiles
@@ -101,29 +99,33 @@
+* astwarp: (gnuastro)astwarp. Options to Warp
-* astscript: (gnuastro)Installed scripts. Gnuastro's installed scripts.
--* astscript-sort-by-night: (gnuastro)Invoking astscript-sort-by-night. Options to this script
--* astscript-radial-profile: (gnuastro)Invoking astscript-radial-profile. Options to this script
-* astscript-ds9-region: (gnuastro)Invoking astscript-ds9-region. Options to this script
-* astscript-fits-view: (gnuastro)Invoking astscript-fits-view. Options to this script
+-* astscript-pointing-simulate: (gnuastro)Invoking astscript-pointing-simulate. Options to this script
+-* astscript-psf-scale-factor: (gnuastro)Invoking astscript-psf-scale-factor. Options to this script
-* astscript-psf-select-stars: (gnuastro)Invoking astscript-psf-select-stars. Options to this script
-* astscript-psf-stamp: (gnuastro)Invoking astscript-psf-stamp. Options to this script
--* astscript-psf-unite: (gnuastro)Invoking astscript-psf-unite. Options to this script
--* astscript-psf-scale-factor: (gnuastro)Invoking astscript-psf-scale-factor. Options to this script
-* astscript-psf-subtract: (gnuastro)Invoking astscript-psf-subtract. Options to this script
+-* astscript-psf-unite: (gnuastro)Invoking astscript-psf-unite. Options to this script
+-* astscript-radial-profile: (gnuastro)Invoking astscript-radial-profile. Options to this script
+-* astscript-sort-by-night: (gnuastro)Invoking astscript-sort-by-night. Options to this script
+-* astscript-zeropoint: (gnuastro)Invoking astscript-zeropoint. Options to this script
+* astscript: (gnuastro)Installed scripts. Gnuastro's installed scripts
-+* astscript-sort-by-night: (gnuastro)astscript-sort-by-night. Options to this script
-+* astscript-radial-profile: (gnuastro)astscript-radial-profile. Options to this script
+* astscript-ds9-region: (gnuastro)astscript-ds9-region. Options to this script
+* astscript-fits-view: (gnuastro)astscript-fits-view. Options to this script
++* astscript-pointing-simulate: (gnuastro)astscript-pointing-simulate. Options to this script
++* astscript-psf-scale-factor: (gnuastro)astscript-psf-scale-factor. Options to this script
+* astscript-psf-select-stars: (gnuastro)astscript-psf-select-stars. Options to this script
+* astscript-psf-stamp: (gnuastro)astscript-psf-stamp. Options to this script
-+* astscript-psf-unite: (gnuastro)astscript-psf-unite. Options to this script
-+* astscript-psf-scale-factor: (gnuastro)astscript-psf-scale-factor. Options to this script
+* astscript-psf-subtract: (gnuastro)astscript-psf-subtract. Options to this script
++* astscript-psf-unite: (gnuastro)astscript-psf-unite. Options to this script
++* astscript-radial-profile: (gnuastro)astscript-radial-profile. Options to this script
++* astscript-sort-by-night: (gnuastro)astscript-sort-by-night. Options to this script
++* astscript-zeropoint: (gnuastro)astscript-zeropoint. Options to this script
@end direntry
-@@ -433,7 +433,7 @@ Data containers
+@@ -465,7 +465,7 @@ Data containers
Fits
@@ -132,17 +134,17 @@
Invoking Fits
-@@ -448,7 +448,7 @@ ConvertType
+@@ -479,7 +479,7 @@ ConvertType
+ * Recognized file formats:: Recognized file formats
* Color:: Some explanations on color.
- * Color channels in same pixel grid:: When the WCS slightly differs.
* Annotations for figure in paper:: Adding coordinates or physical scale.
-* Invoking astconvertt:: Options and arguments to ConvertType.
+* astconvertt:: Options and arguments to ConvertType.
Color
-@@ -470,12 +470,12 @@ Table
-
+@@ -503,12 +503,12 @@ Table
+ * Vector columns:: How to keep more than one value in each column.
* Column arithmetic:: How to do operations on table columns.
* Operation precedence in Table:: Order of running options in Table.
-* Invoking asttable:: Options and arguments to Table.
@@ -156,7 +158,7 @@
Data manipulation
-@@ -489,7 +489,7 @@ Crop
+@@ -522,7 +522,7 @@ Crop
* Crop modes:: Basic modes to define crop region.
* Crop section syntax:: How to define a section to crop.
* Blank pixels:: Pixels with no value.
@@ -165,16 +167,16 @@
Invoking Crop
-@@ -502,7 +502,7 @@ Arithmetic
- * Reverse polish notation:: The current notation style for Arithmetic
- * Integer benefits and pitfalls:: Integers have major benefits, but require care
- * Arithmetic operators:: List of operators known to Arithmetic
--* Invoking astarithmetic:: How to run Arithmetic: options and output
-+* astarithmetic:: How to run Arithmetic: options and output
+@@ -536,7 +536,7 @@ Arithmetic
+ * Integer benefits and pitfalls:: Integers have benefits, but require care.
+ * Noise basics:: Introduction various noise models.
+ * Arithmetic operators:: List of operators known to Arithmetic.
+-* Invoking astarithmetic:: How to run Arithmetic: options and output.
++* astarithmetic:: How to run Arithmetic: options and output.
- Arithmetic operators
+ Noise basics
-@@ -531,7 +531,7 @@ Convolve
+@@ -575,7 +575,7 @@ Convolve
* Frequency domain and Fourier operations:: Using frequencies in input.
* Spatial vs. Frequency domain:: When to use which?
* Convolution kernel:: How to specify the convolution kernel.
@@ -183,16 +185,16 @@
Spatial domain convolution
-@@ -557,7 +557,7 @@ Warp
+@@ -600,7 +600,7 @@ Warp
+ * Linear warping basics:: Basics of coordinate transformation.
* Merging multiple warpings:: How to merge multiple matrices.
* Resampling:: Warping an image is re-sampling it.
- * Moire pattern and its correction:: Spatial resonance of the grid pattern on output.
-* Invoking astwarp:: Arguments and options for Warp.
+* astwarp:: Arguments and options for Warp.
Invoking Warp
-@@ -579,7 +579,7 @@ Statistics
+@@ -622,7 +622,7 @@ Statistics
* Sigma clipping:: Definition of @mymath{\sigma}-clipping.
* Least squares fitting:: Fitting with various parametric functions.
* Sky value:: Definition and derivation of the Sky value.
@@ -201,7 +203,7 @@
2D Histograms
-@@ -604,7 +604,7 @@ Invoking Statistics
+@@ -647,7 +647,7 @@ Invoking Statistics
NoiseChisel
* NoiseChisel changes after publication:: Updates since published papers.
@@ -210,7 +212,7 @@
Invoking NoiseChisel
-@@ -614,7 +614,7 @@ Invoking NoiseChisel
+@@ -657,7 +657,7 @@ Invoking NoiseChisel
Segment
@@ -219,16 +221,16 @@
Invoking Segment
-@@ -629,7 +629,7 @@ MakeCatalog
- * Quantifying measurement limits:: For comparing different catalogs.
+@@ -673,7 +673,7 @@ MakeCatalog
* Measuring elliptical parameters:: Estimating elliptical parameters.
* Adding new columns to MakeCatalog:: How to add new columns.
+ * MakeCatalog measurements:: List of all the measurements/columns by MakeCatalog.
-* Invoking astmkcatalog:: Options and arguments to MakeCatalog.
+* astmkcatalog:: Options and arguments to MakeCatalog.
Quantifying measurement limits
-@@ -651,7 +651,7 @@ Invoking MakeCatalog
+@@ -706,7 +706,7 @@ Invoking MakeCatalog
Match
* Matching algorithms:: Different ways to find the match
@@ -237,7 +239,7 @@
Data modeling
-@@ -663,7 +663,7 @@ MakeProfiles
+@@ -717,7 +717,7 @@ MakeProfiles
* Modeling basics:: Astronomical modeling basics.
* If convolving afterwards:: Considerations for convolving later.
* Profile magnitude:: Definition of total profile magnitude.
@@ -246,25 +248,16 @@
Modeling basics
-@@ -684,7 +684,7 @@ Invoking MakeProfiles
- MakeNoise
-
- * Noise basics:: Noise concepts and definitions.
--* Invoking astmknoise:: Options and arguments to MakeNoise.
-+* astmknoise:: Options and arguments to MakeNoise.
+@@ -743,7 +743,7 @@ CosmicCalculator
- Noise basics
-
-@@ -701,7 +701,7 @@ CosmicCalculator
-
- * Distance on a 2D curved space:: Distances in 2D for simplicity
+ * Distance on a 2D curved space:: Distances in 2D for simplicity.
* Extending distance concepts to 3D:: Going to 3D (our real universe).
--* Invoking astcosmiccal:: How to run CosmicCalculator
-+* astcosmiccal:: How to run CosmicCalculator
+-* Invoking astcosmiccal:: How to run CosmicCalculator.
++* astcosmiccal:: How to run CosmicCalculator.
Invoking CosmicCalculator
-@@ -719,28 +719,28 @@ Installed scripts
+@@ -763,23 +763,23 @@ Installed scripts
Sort FITS files by night
@@ -286,6 +279,20 @@
-* Invoking astscript-fits-view:: How to call this script
+* astscript-fits-view:: How to call this script
+ Zero point estimation
+
+-* Invoking astscript-zeropoint:: How to call the script
++* astscript-zeropoint:: How to call the script
+
+ Invoking astscript-zeropoint
+
+@@ -788,16 +788,16 @@ Invoking astscript-zeropoint
+
+ Pointing pattern simulation
+
+-* Invoking astscript-pointing-simulate:: Options and running mode.
++* astscript-pointing-simulate:: Options and running mode.
+
PSF construction and subtraction
* Overview of the PSF scripts:: Summary of concepts and methods
@@ -294,15 +301,15 @@
-* Invoking astscript-psf-unite:: Merge stacks of different regions of PSF.
-* Invoking astscript-psf-scale-factor:: Calculate factor to scale PSF to star.
-* Invoking astscript-psf-subtract:: Put the PSF in the image to subtract.
-+* astscript-psf-select-stars:: Select good starts within an image.
-+* astscript-psf-stamp:: Make a stamp of each star to stack.
-+* astscript-psf-unite:: Merge stacks of different regions of PSF.
-+* astscript-psf-scale-factor:: Calculate factor to scale PSF to star.
-+* astscript-psf-subtract:: Put the PSF in the image to subtract.
++* astscript-psf-select-stars:: Select good starts within an image.
++* astscript-psf-stamp:: Make a stamp of each star to stack.
++* astscript-psf-unite:: Merge stacks of different regions of PSF.
++* astscript-psf-scale-factor:: Calculate factor to scale PSF to star.
++* astscript-psf-subtract:: Put the PSF in the image to subtract.
Makefile extensions (for GNU Make)
-@@ -762,7 +762,7 @@ Review of library fundamentals
+@@ -819,7 +819,7 @@ Review of library fundamentals
BuildProgram
@@ -311,7 +318,7 @@
Gnuastro library
-@@ -988,7 +988,7 @@ $ sudo make install
+@@ -1047,7 +1047,7 @@ $ sudo make install
@c The last command is to enable Gnuastro's custom TAB completion in Bash.
@c For more on this useful feature, see @ref{Shell TAB completion}).
@@ -320,7 +327,35 @@
In @ref{Tutorials}, we have prepared some complete tutorials with common Gnuastro usage scenarios in astronomical research.
They even contain links to download the necessary data, and thoroughly describe every step of the process (the science, statistics and optimal usage of the command-line).
-@@ -1812,7 +1812,7 @@ But they need to be as realistic as poss
+@@ -1164,22 +1164,22 @@ The following scripts are used to estima
+
+ @table @code
+ @item astscript-psf-select-stars
+-(see @ref{Invoking astscript-psf-select-stars}) Find all the stars within an image that are suitable for constructing an extended PSF.
++(see @ref{astscript-psf-select-stars}) Find all the stars within an image that are suitable for constructing an extended PSF.
+ If the image has WCS, this script can automatically query Gaia to find the good stars.
+
+ @item astscript-psf-stamp
+-(see @ref{Invoking astscript-psf-stamp}) build a crop (stamp) of a certain width around a star at a certain coordinate in a larger image.
++(see @ref{astscript-psf-stamp}) build a crop (stamp) of a certain width around a star at a certain coordinate in a larger image.
+ This script will do sub-pixel re-positioning to make sure the star is centered and can optionally mask all other background sources).
+
+ @item astscript-psf-scale-factor
+-(see @ref{Invoking astscript-psf-scale-factor}) Given a PSF model, and the central coordinates of a star in an image, find the scale factor that has to be multiplied by the PSF to scale it to that star.
++(see @ref{astscript-psf-scale-factor}) Given a PSF model, and the central coordinates of a star in an image, find the scale factor that has to be multiplied by the PSF to scale it to that star.
+
+ @item astscript-psf-unite
+-(see @ref{Invoking astscript-psf-unite}) Unite the various components of a PSF into one.
++(see @ref{astscript-psf-unite}) Unite the various components of a PSF into one.
+ Because of saturation and non-linearity, to get a good estimate of the extended PSF, it is necessary to construct various parts from different magnitude ranges.
+
+ @item astscript-psf-subtract
+-(see @ref{Invoking astscript-psf-subtract}) Given the model of a PSF and the central coordinates of a star in the image, do sub-pixel re-positioning of the PSF, scale it to the star and subtract it from the image.
++(see @ref{astscript-psf-subtract}) Given the model of a PSF and the central coordinates of a star in the image, do sub-pixel re-positioning of the PSF, scale it to the star and subtract it from the image.
+ @end table
+
+ @end table
+@@ -2031,7 +2031,7 @@ Finally, in @ref{Pointing pattern design
In these tutorials, we have intentionally avoided too many cross references to make it more easy to read.
For more information about a particular program, you can visit the section with the same name as the program in this book.
Each program section in the subsequent chapters starts by explaining the general concepts behind what it does, for example, see @ref{Convolve}.
@@ -329,8 +364,8 @@
For an explanation of the conventions we use in the example codes through the book, please see @ref{Conventions}.
@menu
-@@ -1927,7 +1927,7 @@ To return to the previous page, press @k
- If you are searching for a specific phrase in the whole book (for example, an option name), press @key{s} and type your search phrase and end it with an @key{<ENTER>}.
+@@ -2152,7 +2152,7 @@ If you are searching for a specific phra
+ Finally, you can return to the command line and quit Info by pressing the @key{q} key.
You do not need to start from the top of the manual every time.
-For example, to get to @ref{Invoking astnoisechisel}, run the following command.
@@ -338,7 +373,7 @@
In general, all programs have such an ``Invoking ProgramName'' section in this book.
These sections are specifically for the description of inputs, outputs and configuration options of each program.
You can access them directly for each program by giving its executable name to Info.
-@@ -2036,7 +2036,7 @@ $ astscript-fits-view \
+@@ -2261,7 +2261,7 @@ $ astscript-fits-view \
After running this command, you will see that the DS9 window fully covers the height of your monitor, it is showing the whole image, using a more clear color-map, and many more useful things.
In fact, you see the DS9 command that is used in your terminal at footnote{When comparing DS9's command-line options to Gnuastro's, you will notice how SAO DS9 does not follow the GNU style of options where ``long'' and ``short'' options are preceded by @option{--} and @option{-} respectively (for example, @option{--width} and @option{-w}, see @ref{Options}).}.
On GNU/Linux operating systems (like Ubuntu, and Fedora), you can also set your graphics user interface to use this script for opening FITS files when you click on them.
@@ -347,16 +382,16 @@
As you hover your mouse over the image, notice how the ``Value'' and positional fields on the top of the ds9 window get updated.
The first thing you might notice is that when you hover the mouse over the regions with no data, they have a value of zero.
-@@ -2667,7 +2667,7 @@ As you see in the matrix that is printed
+@@ -2904,7 +2904,7 @@ As you see in the matrix that is printed
However, if you run Warp multiple times, the pixels will be mixed multiple times, creating a strong artificial blur/smoothing, or stronger correlated noise.
Recall that the merging of multiple warps is done through matrix multiplication, therefore order matters in the separate operations.
-At a lower level, through Warp's @option{--matrix} option, you can directly request your desired final warp and do not have to break it up into different warps like above (see @ref{Invoking astwarp}).
+At a lower level, through Warp's @option{--matrix} option, you can directly request your desired final warp and do not have to break it up into different warps like above (see @ref{astwarp}).
- Fortunately these datasets are already aligned to the same pixel grid, so you do not actually need the files that were just generated.You can safely delete them all with the following command.
- Here, you see why we put the processed outputs that we need later into a separate directory.
-@@ -3864,7 +3864,7 @@ $ cat apertures.txt
+ Fortunately these datasets are already aligned to the same pixel grid, so you do not actually need the files that were just generated.
+ You can safely delete them all with the following command.
+@@ -4125,7 +4125,7 @@ $ cat apertures.txt
We can now feed this catalog into MakeProfiles using the command below to build the apertures over the image.
The most important option for this particular job is @option{--mforflatpix}, it tells MakeProfiles that the values in the magnitude column should be used for each pixel of a flat profile.
Without it, MakeProfiles would build the profiles such that the @emph{sum} of the pixels of each profile would have a @emph{magnitude} (in log-scale) of the value given in that column (what you would expect when simulating a galaxy for example).
@@ -365,7 +400,25 @@
@example
$ astmkprof apertures.txt --background=flat-ir/xdf-f160w.fits \
-@@ -6673,7 +6673,7 @@ $ echo $scale
+@@ -6607,7 +6607,7 @@ Let's use @file{astscript-psf-select-sta
+ The advantage of using this script (instead of a simple @option{--range} in Table), is that it will also check distances to nearby stars and reject those that are too close (and not good for constructing the PSF).
+ Since we have very bright stars in this very wide-field image, we will also increase the distance to nearby neighbors with brighter or similar magnitudes (the default value is 1 arcmin).
+ To do this, we will set @option{--mindistdeg=0.02}, which corresponds to 1.2 arcmin.
+-The details of the options for this script are discussed in @ref{Invoking astscript-psf-select-stars}.
++The details of the options for this script are discussed in @ref{astscript-psf-select-stars}.
+
+ @example
+ $ mkdir outer
+@@ -6628,7 +6628,7 @@ $ astscript-ds9-region outer/67510-6-10.
+ Now that the catalog of good stars is ready, it is time to construct the individual stamps from the catalog above.
+ To create stamps, first, we need to crop a fixed-size box around each isolated star in the catalog.
+ The contaminant objects in the crop should be masked and finally, the fluxes in these cropped images should be normalized.
+-To do these, we will use @file{astscript-psf-stamp} (for more on this script see @ref{Invoking astscript-psf-stamp}).
++To do these, we will use @file{astscript-psf-stamp} (for more on this script see @ref{astscript-psf-stamp}).
+
+ One of the most important parameters for this script is the normalization radii @code{--normradii}.
+ This parameter defines a ring for the flux normalization of each star stamp.
+@@ -6941,7 +6941,7 @@ $ echo $scale
@end example
Now that we know the scaling factor, we are ready to unite the outer and the inner part of the PSF.
@@ -374,7 +427,16 @@
The basic parameters are the inner part of the PSF (given to @option{--inner}), the inner part's scale factor (@option{--scale}), and the junction radius (@option{--radius}).
The inner part is first scaled, and all the pixels of the outer image within the given radius are replaced with the pixels of the inner image.
Since the flux factor was computed for a ring of pixels between 10 and 15 pixels, let's set the junction radius to be 12 pixels (roughly in between 10 and 15):
-@@ -9567,7 +9567,7 @@ Afterwards, in @ref{Common options}, we
+@@ -8949,7 +8949,7 @@ $ aststatistics jplus-nc.fits -hSKY_STD
+ @end example
+
+ We are now ready to find the zero point!
+-First, let's run the @command{astscript-zeropoint} with @option{--help} to see the option names (recall that you can see more details of each option in @ref{Invoking astscript-zeropoint}).
++First, let's run the @command{astscript-zeropoint} with @option{--help} to see the option names (recall that you can see more details of each option in @ref{astscript-zeropoint}).
+ For the first time, let's use the script in the most simple state possible.
+ We will keep only the essential options: the names of the input and reference images (and their HDUs), the name of the output, and also two apertures with radii of 3 arcsec to start with:
+
+@@ -12233,7 +12233,7 @@ Afterwards, in @ref{Common options}, we
@cindex Arguments to programs
@cindex Command-line arguments
When you type a command on the command-line, it is passed onto the shell (a generic name for the program that manages the command-line) as a string of characters.
@@ -383,7 +445,7 @@
The shell then brakes up your string into separate @emph{tokens} or @emph{words} using any @emph{metacharacters} (like white-space, tab, @command{|}, @command{>} or @command{;}) that are in the string.
On the command-line, the first thing you usually enter is the name of the program you want to run.
-@@ -9679,7 +9679,7 @@ Command-line options allow configuring t
+@@ -12345,7 +12345,7 @@ Command-line options allow configuring t
A single option can be called in two ways: @emph{long} or @emph{short}.
All options in Gnuastro accept the long format which has two hyphens an can have many characters (for example, @option{--hdu}).
Short options only have one hyphen (@key{-}) followed by one character (for example, @option{-h}).
@@ -392,7 +454,7 @@
Both formats are shown for those which support both.
First the short is shown then the long.
-@@ -9753,7 +9753,7 @@ If you are satisfied with the change, yo
+@@ -12419,7 +12419,7 @@ If you are satisfied with the change, yo
If the change was not satisfactory, you can remove the one you just added and not worry about forgetting the original value.
Without this capability, you would have to memorize or save the original value somewhere else, run the command and then change the value again which is not at all convenient and is potentially cause lots of bugs.
@@ -401,7 +463,7 @@
In these cases, the order of stored values is the same order that you specified on the command-line.
@cindex Configuration files
-@@ -10224,7 +10224,7 @@ As another example, if an option needs i
+@@ -12934,7 +12934,7 @@ As another example, if an option needs i
Just note that the file name has to be already given on the command-line before reaching such options (that look into the contents of a file).
But TAB completion is not limited to file types or contents.
@@ -410,7 +472,7 @@
As another example, the option @option{--numthreads} option (to specify the number of threads to use by the program), will find the number of available threads on the system, and suggest the possible numbers with a TAB!
To activate Gnuastro's custom TAB completion in Bash, you need to put the following line in one of your Bash startup files (for example, @file{~/.bashrc}).
-@@ -10344,7 +10344,7 @@ All the values must either be stored in
+@@ -13170,7 +13170,7 @@ All the values must either be stored in
In case the necessary parameters are not given through any of these methods, the program will print a missing option error and abort.
The only exception to this is @option{--numthreads}, whose default value is determined at run-time using the number of threads available to your system, see @ref{Multi-threaded operations}.
Of course, you can still provide a default value for the number of threads at any of the levels below, but if you do not, the program will not abort.
@@ -419,7 +481,7 @@
-@@ -11153,7 +11153,7 @@ Each column in the table contains the va
+@@ -13981,7 +13981,7 @@ Each column in the table contains the va
For example, let's assume you have just ran MakeCatalog (see @ref{MakeCatalog}) on an image to measure some properties for the labeled regions (which might be detected galaxies for example) in the image.
For each labeled region (detected galaxy), there will be a @emph{row} which groups its measured properties as @emph{columns}, one column for each property.
One such property can be the object's magnitude, which is the sum of pixels with that label, or its center can be defined as the light-weighted average value of those pixels.
@@ -428,7 +490,7 @@
As a summary, for each labeled region (or, galaxy) we have one @emph{row} and for each measured property we have one @emph{column}.
This high-level structure is usually the first step for higher-level analysis, for example, finding the stellar mass or photometric redshift from magnitudes in multiple colors.
-@@ -11775,10 +11775,10 @@ for example, you can copy or cut (copy a
+@@ -14607,10 +14607,10 @@ For example, you can copy or cut (copy a
It also has features to delete, add, or edit meta-data keywords within one HDU.
@menu
@@ -441,7 +503,7 @@
@subsection Invoking Fits
Fits can print or manipulate the FITS file HDUs (extensions), meta-data keywords in a given HDU.
-@@ -11874,7 +11874,7 @@ These three classes of options cannot be
+@@ -14706,7 +14706,7 @@ These three classes of options cannot be
@@ -450,7 +512,7 @@
@subsubsection HDU information and manipulation
Each FITS file header data unit, or HDU (also known as an extension) is an independent dataset (data + meta-data).
Multiple HDUs can be stored in one FITS file, see @ref{Fits}.
-@@ -11892,7 +11892,7 @@ Note that this option must be called alo
+@@ -14724,7 +14724,7 @@ Note that this option must be called alo
It is thus useful in scripts, for example, when you need to do check the number of extensions in a FITS file.
For a complete list of basic meta-data on the extensions in a FITS file, do not use any of the options in this section or in @ref{Keyword inspection and manipulation}.
@@ -459,7 +521,7 @@
@item --hastablehdu
Print @code{1} (on standard output) if at least one table HDU (ASCII or binary) exists in the FITS file.
-@@ -12038,7 +12038,7 @@ If we had not used @option{--primaryimgh
+@@ -14870,7 +14870,7 @@ If we had not used @option{--primaryimgh
@end table
@@ -468,7 +530,7 @@
@subsubsection Keyword inspection and manipulation
The meta-data in each header data unit, or HDU (also known as extension, see @ref{Fits}) is stored as ``keyword''s.
Each keyword consists of a name, value, unit, and comments.
-@@ -12071,7 +12071,7 @@ $ astfits image-a.fits --keyvalue=NAXIS,
+@@ -14903,7 +14903,7 @@ $ astfits image-a.fits --keyvalue=NAXIS,
@end example
The output is internally stored (and finally printed) as a table (with one column per keyword).
@@ -477,7 +539,7 @@
The keyword metadata (comments and units) are extracted from the comments and units of the keyword in the input files (first file that has a comment or unit).
Hence if the keyword does not have units or comments in any of the input files, they will be empty.
For more on Gnuastro's plain-text metadata format, see @ref{Gnuastro text table format}.
-@@ -12513,7 +12513,7 @@ $ astfits in.fits --wcsdistortion=SIP --
+@@ -15355,7 +15355,7 @@ $ astfits in.fits --wcsdistortion=SIP --
@end table
@@ -486,7 +548,7 @@
@subsubsection Pixel information images
In @ref{Keyword inspection and manipulation} options like @option{--pixelscale} were introduced for information on the pixels from the keywords.
But that only provides a single value for all the pixels!
-@@ -12614,7 +12614,7 @@ The conversion is not only one way (from
+@@ -15456,7 +15456,7 @@ The conversion is not only one way (from
So you can also convert a JPEG image or text file into a FITS image.
Basically, other than EPS/PDF, you can use any of the recognized formats as different color channel inputs to get any of the recognized outputs.
@@ -495,16 +557,16 @@
In ConvertType, vector graphics are used to add markers over your originally rasterized data, producing high quality images, ready to be used in your exciting papers.
We will continue with a description of the recognized files types in @ref{Recognized file formats}, followed a short introduction to digital color in @ref{Color}.
A tutorial on how to add markers over an image is then given in @ref{Marking objects for publication} and we conclude with a @LaTeX{} based solution to add coordinates over an image.
-@@ -12625,7 +12625,7 @@ A tutorial on how to add markers over an
+@@ -15466,7 +15466,7 @@ A tutorial on how to add markers over an
+ * Recognized file formats:: Recognized file formats
* Color:: Some explanations on color.
- * Color channels in same pixel grid:: When the WCS slightly differs.
* Annotations for figure in paper:: Adding coordinates or physical scale.
-* Invoking astconvertt:: Options and arguments to ConvertType.
+* astconvertt:: Options and arguments to ConvertType.
@end menu
@node Raster and Vector graphics, Recognized file formats, ConvertType, ConvertType
-@@ -12775,7 +12775,7 @@ Most programs also support input as plai
+@@ -15612,7 +15612,7 @@ Most programs also support input as plai
As input, each plain text file is considered to contain one color channel.
In ConvertType, the recognized extensions for plain text files are @file{.txt} and @file{.dat}.
@@ -513,7 +575,7 @@
Besides these, when the format of a file cannot be recognized from its name, ConvertType will fall back to plain text mode.
So you can use any name (even without an extension) for a plain text input or output.
Just note that when the suffix is not recognized, automatic output will not be preformed.
-@@ -12863,8 +12863,8 @@ As a result, there is a lot of freedom i
+@@ -15700,8 +15700,8 @@ As a result, there is a lot of freedom i
The mapping of single-channel values to multi-channel colors is called called a ``color map''.
Since more information can be put in multiple channels, this usually results in better visualizing the dynamic range of your single-channel data.
@@ -524,16 +586,16 @@
@itemize
@item
-@@ -13000,7 +13000,7 @@ This shows how green and red channels ha
- If you don't want to have those, or if you want the outer parts of the final image (where there was no data) to be white, some more complex commands are necessary.
- We'll leave those as an exercise for you to try your self using @ref{Warp} and/or @ref{Crop} to pre-process the inputs before converting it to a color image.
+@@ -15769,7 +15769,7 @@ $ astconvertt --listcolors
+ @caption{Recognized color names in Gnuastro, shown with their numerical identifiers.}
+ @end float
-- at node Annotations for figure in paper, Invoking astconvertt, Color channels in same pixel grid, ConvertType
-+ at node Annotations for figure in paper, astconvertt, Color channels in same pixel grid, ConvertType
+- at node Annotations for figure in paper, Invoking astconvertt, Color, ConvertType
++ at node Annotations for figure in paper, astconvertt, Color, ConvertType
@subsection Annotations for figure in paper
@cindex Image annotation
-@@ -13407,7 +13407,7 @@ pdflatex -shell-escape -halt-on-error re
+@@ -16176,7 +16176,7 @@ pdflatex -shell-escape -halt-on-error re
@@ -542,7 +604,7 @@
@subsection Invoking ConvertType
ConvertType will convert any recognized input file type to any specified output type.
-@@ -13459,7 +13459,7 @@ The common options are not repeated here
+@@ -16228,7 +16228,7 @@ The common options are not repeated here
* Drawing with vector graphics:: Adding marks in many shapes and colors over the pixels.
@end menu
@@ -551,7 +613,7 @@
@subsubsection ConvertType input and output
@cindex Standard input
-@@ -13565,7 +13565,7 @@ If it is in color mode, some degradation
+@@ -16334,7 +16334,7 @@ If it is in color mode, some degradation
While the JPEG standard does support loss-less graphics, it is not commonly supported.
@end table
@@ -560,7 +622,7 @@
@subsubsection Pixel visualization
The main goal of ConvertType is to visualize pixels to/from print or web friendly formats.
-@@ -13697,7 +13697,7 @@ The result will be that a large are of t
+@@ -16466,7 +16466,7 @@ The result will be that a large are of t
Note that this behavior is ideal for gray-scale images, if you want a color image, the colors are going to be mixed up.
@end table
@@ -569,7 +631,7 @@
@subsubsection Drawing with vector graphics
With the options described in this section, you can draw marks over your to-be-published images (for example, in PDF).
-@@ -14000,18 +14000,18 @@ Table has a large set of operations that
+@@ -16769,20 +16769,20 @@ Table has a large set of operations that
For operations that Table does not do internally, FITS tables (ASCII or binary) are directly accessible to the users of Unix-like operating systems (in particular those working the command-line or shell, see @ref{Command-line interface}).
With Table, a FITS table (in binary or ASCII formats) is only one command away from AWK (or any other tool you want to use).
Just like a plain text file that you read with the @command{cat} command.
@@ -584,24 +646,44 @@
+Finally, in @ref{asttable}, we give some examples and describe each option in Table.
@menu
+ * Printing floating point numbers:: Optimal storage of floating point types.
+ * Vector columns:: How to keep more than one value in each column.
* Column arithmetic:: How to do operations on table columns.
* Operation precedence in Table:: Order of running options in Table.
-* Invoking asttable:: Options and arguments to Table.
+* asttable:: Options and arguments to Table.
@end menu
- @node Column arithmetic, Operation precedence in Table, Table, Table
-@@ -14074,7 +14074,7 @@ Column arithmetic changes the values of
- So the old column meta data cannot be used any more.
+ @node Printing floating point numbers, Vector columns, Table, Table
+@@ -16842,7 +16842,7 @@ However, if you choose to store the outp
+ @end cartouche
+
+ To customize how columns containing floating point values are printed (in a plain-text output file, or in the standard output in your terminal), Table has four options for the two different types: @option{--txtf32format}, @option{--txtf32precision}, @option{--txtf64format} and @option{--txtf64precision}.
+-They are fully described in @ref{Invoking asttable}.
++They are fully described in @ref{asttable}.
+
+ @cartouche
+ @noindent
+@@ -16906,7 +16906,7 @@ With the command below, let's have a loo
+ When you run it, you will have a different random number generator seed, so the numbers will be slightly different.
+ For making reproducible random numbers, see @ref{Generating random numbers}.
+ The @option{-Y} option is used for more easily readable numbers (without it, floating point numbers are written in scientific notation, for more see @ref{Printing floating point numbers}) and with the @option{-O} we are asking Table to also print the metadata.
+-For more on Table's options, see @ref{Invoking asttable} and for seeing how the short options can be merged (such that @option{-Y -O} is identical to @option{-YO}), see @ref{Options}.
++For more on Table's options, see @ref{asttable} and for seeing how the short options can be merged (such that @option{-Y -O} is identical to @option{-YO}), see @ref{Options}.
+
+ @example
+ $ asttable table.fits -YO
+@@ -17119,7 +17119,7 @@ Column arithmetic changes the values of
+ So the old column metadata cannot be used any more.
By default the output column of the arithmetic operation will be given a generic metadata (for example, its name will be @code{ARITH_1}, which is hardly useful!).
- But meta data are critically important and it is good practice to always have short, but descriptive, names for each columns, units and also some comments for more explanation.
+ But metadata are critically important and it is good practice to always have short, but descriptive, names for each columns, units and also some comments for more explanation.
-To add metadata to a column, you can use the @option{--colmetadata} option that is described in @ref{Invoking asttable} and @ref{Operation precedence in Table}.
+To add metadata to a column, you can use the @option{--colmetadata} option that is described in @ref{asttable} and @ref{Operation precedence in Table}.
Since the arithmetic expressions are a value to @option{--column}, it does not necessarily have to be a separate option, so the commands above are also identical to the command below (note that this only has one @option{-c} option).
Just be very careful with the quoting!
-@@ -14230,13 +14230,13 @@ Other than the units of the output, this
- See the description of that operator for an example.
+@@ -17371,13 +17371,13 @@ Such intervals can be useful in scenario
+
@end table
- at node Operation precedence in Table, Invoking asttable, Column arithmetic, Table
@@ -613,19 +695,37 @@
Knowing this precedence table is important to avoid confusion when you ask for more than one operation.
-For a description of each option, please see @ref{Invoking asttable}.
+For a description of each option, please see @ref{asttable}.
+ By default, column-based operations will be done first.
+ You can ask for switching to row-based operations to be done first, using the @option{--rowfirst} option.
- @table @asis
- @item Column information (@option{--information} or @option{-i})
-@@ -14381,7 +14381,7 @@ asttable table.fits -cRA,DEC --noblanken
+@@ -17474,7 +17474,7 @@ For more on type conversion, see @ref{Nu
+
+ @item Row selection by value in a column
+ The following operations select rows based on the values in them.
+-A more complete description of each of these options is given in @ref{Invoking asttable}.
++A more complete description of each of these options is given in @ref{asttable}.
+
+ @itemize
+ @item
+@@ -17526,7 +17526,7 @@ These options limit/select rows based on
+ @item Transpose vector columns (@option{--transpose})
+ Transposing vector columns will not affect the number or metadata of columns, it will just re-arrange them in their 2D structure.
+ As a result, after transposing, the number of rows changes, as well as the number of elements in each vector column.
+-See the description of this option in @ref{Invoking asttable} for more (with an example).
++See the description of this option in @ref{asttable} for more (with an example).
+ @end table
+
+ @item Column metadata (@option{--colmetadata})
+@@ -17586,7 +17586,7 @@ asttable table.fits -cRA,DEC --noblanken
+ | asttable --sort=MULTIP --output=selected.fits
@end example
- @end cartouche
- at node Invoking asttable, , Operation precedence in Table, Table
+ at node asttable, , Operation precedence in Table, Table
@subsection Invoking Table
Table will read/write, select, modify, or show the information of the rows and columns in recognized Table formats (including FITS binary, FITS ASCII, and plain text table files, see @ref{Tables}).
-@@ -14591,7 +14591,7 @@ $ asttable a.fits --catrowfile=b.fits --
+@@ -17870,7 +17870,7 @@ For a working example of this scenario,
@strong{How to avoid repetition when adding rows:} this option will simply add the rows of multiple tables into one, it does not check their contents!
Therefore if you use this option on multiple catalogs that may have some shared physical objects in some of their rows, those rows/objects will be repeated in the final table.
In such scenarios, to avoid potential repetition, it is better to use @ref{Match} (with @option{--notmatched} and @option{--outcols=AAA,BBB}) instead of Table.
@@ -634,7 +734,7 @@
@end cartouche
@item -X STR
-@@ -14623,7 +14623,7 @@ The chosen column does not have to be in
+@@ -17902,7 +17902,7 @@ The chosen column does not have to be in
This is good when you just want to select using one column's values, but do not need that column anymore afterwards.
For one example of using this option, see the example under
@@ -643,7 +743,7 @@
@item --inpolygon=STR1,STR2
Only return rows where the given coordinates are inside the polygon specified by the @option{--polygon} option.
-@@ -14938,10 +14938,10 @@ We will try to add high-level interfaces
+@@ -18299,10 +18299,10 @@ We will try to add high-level interfaces
@menu
* Available databases:: List of available databases to Query.
@@ -656,7 +756,7 @@
@subsection Available databases
The current list of databases supported by Query are listed at the end of this section.
-@@ -15166,7 +15166,7 @@ For details on each dataset with necessa
+@@ -18527,7 +18527,7 @@ For details on each dataset with necessa
@@ -665,7 +765,7 @@
@subsection Invoking Query
Query provides a high-level interface to downloading subsets of data from databases.
-@@ -15269,7 +15269,7 @@ When @option{--dataset} is specified, th
+@@ -18630,7 +18630,7 @@ When @option{--dataset} is specified, th
Some databases (like VizieR) contain tens of thousands of datasets, so you can limit the downloaded and printed information for available databases with the @option{--limitinfo} option (described below).
Dataset descriptions are often large and contain a lot of text (unlike column descriptions).
Therefore when printing the information of all datasets within a database, the information (e.g., database name) will be printed on separate lines before the description.
@@ -674,7 +774,7 @@
Important note to consider: the printed order of the datasets or columns is just for displaying in the printed output.
You cannot ask for datasets or columns based on the printed order, you need to use dataset or column names.
-@@ -15439,14 +15439,14 @@ So when you simply crop the image of suc
+@@ -18800,14 +18800,14 @@ So when you simply crop the image of suc
Therefore in its WCS mode, Crop will stitch parts of the tiles that are relevant for a target (with the given width) from all the input images that cover that region into the output.
Of course, the tiles have to be present in the list of input files.
@@ -691,7 +791,7 @@
@end menu
@node Crop modes, Crop section syntax, Crop, Crop
-@@ -15468,13 +15468,13 @@ Irrespective of how the crop region is d
+@@ -18829,13 +18829,13 @@ Irrespective of how the crop region is d
All coordinates are read as floating point numbers (not integers, except for the @option{--section} option, see below).
By setting the @emph{mode} in Crop, you define the standard that the given coordinates must be interpreted.
Here, the different ways to specify the crop region are discussed within each standard.
@@ -707,7 +807,7 @@
@table @asis
@item Image coordinates
-@@ -15501,7 +15501,7 @@ The former is lower-level (does not acce
+@@ -18862,7 +18862,7 @@ The former is lower-level (does not acce
Please see @ref{Crop section syntax} for a full description of this method.
The latter option (@option{--polygon}) is a higher-level method to define any polygon (with any number of vertices) with floating point values.
@@ -716,7 +816,7 @@
@end table
@item WCS coordinates
-@@ -15533,7 +15533,7 @@ If it exists in the input images, it wil
+@@ -18894,7 +18894,7 @@ If it exists in the input images, it wil
@item Vertices of a single crop
The @option{--polygon} option is a high-level method to define any convex polygon (with any number of vertices).
@@ -725,7 +825,7 @@
@end table
@cartouche
-@@ -15587,7 +15587,7 @@ If you forget the quotes, anything after
+@@ -18948,7 +18948,7 @@ If you forget the quotes, anything after
See @ref{Command-line} for a description of how the command-line works.
@@ -734,7 +834,7 @@
@subsection Blank pixels
@cindex Blank pixel
-@@ -15612,7 +15612,7 @@ So by default, when dealing with float o
+@@ -18973,7 +18973,7 @@ So by default, when dealing with float o
This can be turned off with the @option{--zeroisnotblank} option.
@@ -743,7 +843,7 @@
@subsection Invoking Crop
Crop will crop a region from an image.
-@@ -15681,7 +15681,7 @@ This has no effect on each output, see @
+@@ -19052,7 +19052,7 @@ This has no effect on each output, see @
* Crop known issues:: Known issues in running Crop.
@end menu
@@ -752,7 +852,7 @@
@subsubsection Crop options
The options can be classified into the following contexts: Input, Output and operating mode options.
-@@ -15947,7 +15947,7 @@ The value must either be @option{img} or
+@@ -19263,7 +19263,7 @@ When this column is not given, the row n
@@ -761,7 +861,7 @@
@subsubsection Crop output
The string given to @option{--output} option will be interpreted depending on how many crops were requested, see @ref{Crop modes}:
-@@ -15988,7 +15988,7 @@ The cropped image file name for that row
+@@ -19307,7 +19307,7 @@ The cropped image file name for that row
The number of input images that were used to create that image.
@item
A @code{0} if the central few pixels (value to the @option{--checkcenter} option) are blank and @code{1} if they are not.
@@ -770,7 +870,7 @@
@end enumerate
If the output crop(s) have a single element (pixel in an image) and @option{--oneelemstdout} has been called, no output file will be produced!
-@@ -16018,7 +16018,7 @@ In other words, the file in the first co
+@@ -19391,7 +19391,7 @@ You can disable this behavior with this
@end table
@@ -779,7 +879,7 @@
@subsubsection Crop known issues
When running Crop, you may encounter strange errors and bugs.
-@@ -16068,14 +16068,14 @@ For example, in the reduction of raw dat
+@@ -19441,7 +19441,7 @@ For example, in the reduction of raw dat
Later (once the images as warped into a single grid using Warp for example, see @ref{Warp}), the images are co-added (the output pixel grid is the average of the pixels of the individual input images).
Arithmetic is Gnuastro's program for such operations on your datasets directly from the command-line.
It currently uses the reverse polish or post-fix notation, see @ref{Reverse polish notation} and will work on the native data types of the input images/data to reduce CPU and RAM resources, see @ref{Numeric data types}.
@@ -788,15 +888,16 @@
@menu
- * Reverse polish notation:: The current notation style for Arithmetic
- * Integer benefits and pitfalls:: Integers have major benefits, but require care
- * Arithmetic operators:: List of operators known to Arithmetic
--* Invoking astarithmetic:: How to run Arithmetic: options and output
-+* astarithmetic:: How to run Arithmetic: options and output
+@@ -19449,7 +19449,7 @@ For more information on how to run Arith
+ * Integer benefits and pitfalls:: Integers have benefits, but require care.
+ * Noise basics:: Introduction various noise models.
+ * Arithmetic operators:: List of operators known to Arithmetic.
+-* Invoking astarithmetic:: How to run Arithmetic: options and output.
++* astarithmetic:: How to run Arithmetic: options and output.
@end menu
@node Reverse polish notation, Integer benefits and pitfalls, Arithmetic, Arithmetic
-@@ -16153,7 +16153,7 @@ There are no more operands or operators,
+@@ -19527,7 +19527,7 @@ There are no more operands or operators,
In the kitchen metaphor, you see that your recipe has no more steps, so you just pick up the remaining dish and take it to the dining room to enjoy a good dinner.
@end enumerate
@@ -805,16 +906,34 @@
In Table's column arithmetic, they can be any column in the table (a series of numbers in an array) or a single number (see @ref{Column arithmetic}).
With this notation, very complicated procedures can be created without the need for parenthesis or worrying about precedence.
-@@ -16251,7 +16251,7 @@ It is possible to do internal checks in
- However, we have not opted for this solution because all those checks will consume significant resources and slow down the program (especially with large datasets where RAM, storage and running time become important).
- To be optimal, we therefore trust that you (the wise Gnuastro user!) make the appropriate type conversion in your commands where necessary (recall that the operators are available in @ref{Numerical type conversion operators}).
+@@ -19849,7 +19849,7 @@ Gnuastro's programs will report the actu
+ In the example above, did you notice how the random number generator seed above is different between GSL and MakeProfiles?
+ However, if @option{--envseed} was given, both printed seeds would be the same.
-- at node Arithmetic operators, Invoking astarithmetic, Integer benefits and pitfalls, Arithmetic
-+ at node Arithmetic operators, astarithmetic, Integer benefits and pitfalls, Arithmetic
+- at node Arithmetic operators, Invoking astarithmetic, Noise basics, Arithmetic
++ at node Arithmetic operators, astarithmetic, Noise basics, Arithmetic
@subsection Arithmetic operators
In this section, list of recognized operators in Arithmetic (and the Table program's @ref{Column arithmetic}) and discussed in detail with examples.
-@@ -18112,7 +18112,7 @@ Similar to the @code{tofile} operator, w
+@@ -20988,7 +20988,7 @@ This operator only takes a single operan
+ Therefore it does not change the layout of the array contents in memory and is very fast.
+
+ If no further operation is requested on the 1D array, recall that Arithmetic will write a 1D array as a table column by default.
+-In case you want the output to be saved as a 1D image, or to see it on the standard output, please use the @code{--onedasimage} or @code{--onedonstdout} options respectively (see @ref{Invoking astarithmetic}).
++In case you want the output to be saved as a 1D image, or to see it on the standard output, please use the @code{--onedasimage} or @code{--onedonstdout} options respectively (see @ref{astarithmetic}).
+
+ This operator is useful in scenarios where after some operations on a 2D image or 3D cube, the dimensionality is no longer relevant for you and you just care about the values.
+ In the example below, we will first make a simple 2D image from a plain-text file, then convert it to a 1D array:
+@@ -22421,7 +22421,7 @@ The number of elements is read from the
+ The second popped operand is only used for its number of elements, its numeric data type, or its values are fully ignored and it is later freed.
+
+ @cindex Provenance
+-Here is one useful scenario for this operator in tables: you want to merge the objects/rows of some catalogs together, but you first want to give each source catalog a label/counter that distinguishes between the source of each rows in the merged/final catalog (using @ref{Invoking asttable}).
++Here is one useful scenario for this operator in tables: you want to merge the objects/rows of some catalogs together, but you first want to give each source catalog a label/counter that distinguishes between the source of each rows in the merged/final catalog (using @ref{asttable}).
+ The steps below show the the usage of this.
+
+ @example
+@@ -22555,7 +22555,7 @@ Similar to the @code{tofile} operator, w
@end table
@@ -823,7 +942,7 @@
@subsection Invoking Arithmetic
Arithmetic will do pixel to pixel arithmetic operations on the individual pixels of input data and/or numbers.
-@@ -18337,7 +18337,7 @@ $ echo "" | awk '@{print (10.32-3.84)^2.
+@@ -22787,7 +22787,7 @@ $ echo "" | awk '@{print (10.32-3.84)^2.
@cindex Average, weighted
@cindex Kernel, convolution
On an image, convolution can be thought of as a process to blur or remove the contrast in an image.
@@ -832,7 +951,7 @@
There are generally two methods to convolve an image.
The first and more intuitive one is in the ``spatial domain'' or using the actual image pixel values, see @ref{Spatial domain convolution}.
-@@ -18368,7 +18368,7 @@ However this text is written for an unde
+@@ -22818,7 +22818,7 @@ However this text is written for an unde
* Frequency domain and Fourier operations:: Using frequencies in input.
* Spatial vs. Frequency domain:: When to use which?
* Convolution kernel:: How to specify the convolution kernel.
@@ -841,7 +960,7 @@
@end menu
@node Spatial domain convolution, Frequency domain and Fourier operations, Convolve, Convolve
-@@ -19112,7 +19112,7 @@ The reason is that if you apply a freque
+@@ -23562,7 +23562,7 @@ The reason is that if you apply a freque
But when you have made the profiles in the image yourself, you can just make a larger input image and crop the central parts to completely remove the edge effect, see @ref{If convolving afterwards}.
Also due to oversampling, both the kernels and the images can become very large and the speed boost of frequency domain convolution will significantly improve the processing time, see @ref{Oversampling}.
@@ -850,7 +969,7 @@
@subsection Convolution kernel
All the programs that need convolution will need to be given a convolution kernel file and extension.
-@@ -19130,7 +19130,7 @@ By default MakeProfiles will make the Ga
+@@ -23580,7 +23580,7 @@ By default MakeProfiles will make the Ga
ConvertType: You can write your own desired kernel into a text file table and convert it to a FITS file with ConvertType, see @ref{ConvertType}.
Just be careful that the kernel has to have an odd number of pixels along its two axes, see @ref{Convolution process}.
All the programs that do convolution will normalize the kernel internally, so if you choose this option, you do not have to worry about normalizing the kernel.
@@ -859,7 +978,7 @@
@end itemize
-@@ -19170,7 +19170,7 @@ By default, the system configuration fil
+@@ -23620,7 +23620,7 @@ By default, the system configuration fil
@end table
@@ -868,34 +987,34 @@
@subsection Invoking Convolve
Convolve an input dataset (2D image or 1D spectrum for example) with a known kernel, or make the kernel necessary to match two PSFs.
-@@ -19395,7 +19395,7 @@ It is therefore necessary to warp the im
+@@ -23854,7 +23854,7 @@ It is therefore necessary to warp the im
+ * Linear warping basics:: Basics of coordinate transformation.
* Merging multiple warpings:: How to merge multiple matrices.
* Resampling:: Warping an image is re-sampling it.
- * Moire pattern and its correction:: Spatial resonance of the grid pattern on output.
-* Invoking astwarp:: Arguments and options for Warp.
+* astwarp:: Arguments and options for Warp.
@end menu
@node Linear warping basics, Merging multiple warpings, Warp, Warp
-@@ -19616,7 +19616,7 @@ However, when a non-linear distortion (f
- To account for such cases (which can only happen when correcting for non-linear distortions), Warp has the @option{--edgesampling} option to sample the output pixel over more vertices.
- For more, see the description of this option in @ref{Align pixels with WCS considering distortions}.
+@@ -23997,7 +23997,7 @@ These three operations can be merged in
+
+
+
+- at node Resampling, Invoking astwarp, Merging multiple warpings, Warp
++ at node Resampling, astwarp, Merging multiple warpings, Warp
+ @subsection Resampling
-- at node Moire pattern and its correction, Invoking astwarp, Resampling, Warp
-+ at node Moire pattern and its correction, astwarp, Resampling, Warp
- @subsection Moir@'e pattern and its correction
+ @cindex Pixel
+@@ -24078,7 +24078,7 @@ For more, see the description of this op
- @cindex Moir@'e pattern or fringes
-@@ -19838,7 +19838,7 @@ For example if you only have two dither
- Ideally, many more dither points should be chosen when you are planning your observation (not just for the Moir@'e pattern, but also for all the other reasons mentioned above).
- Based on the dithering pattern, you want to select the increased resolution such that the maximum @code{MAX-FRAC} values fall on every different pixel of the output grid for each exposure.
-- at node Invoking astwarp, , Moire pattern and its correction, Warp
-+ at node astwarp, , Moire pattern and its correction, Warp
+
+- at node Invoking astwarp, , Resampling, Warp
++ at node astwarp, , Resampling, Warp
@subsection Invoking Warp
- Warp an input image into a new pixel grid by pixel mixing (see @ref{Resampling}).
-@@ -19911,10 +19911,10 @@ The following options are shared between
+ Warp will warp an input image into a new pixel grid by pixel mixing (see @ref{Resampling}).
+@@ -24151,10 +24151,10 @@ The following options are shared between
@table @option
@item --hstartwcs=INT
@@ -908,7 +1027,7 @@
@item -C FLT
@itemx --coveredfrac=FLT
-@@ -19930,7 +19930,7 @@ As a result, with @option{--coveredfrac=
+@@ -24170,7 +24170,7 @@ As a result, with @option{--coveredfrac=
* Linear warps to be called explicitly:: Other warps.
@end menu
@@ -917,7 +1036,7 @@
@subsubsection Align pixels with WCS considering distortions
@cindex Resampling
-@@ -20198,7 +20198,7 @@ On the other hand, 0 means that the pixe
+@@ -24454,7 +24454,7 @@ On the other hand, 0 means that the pixe
@@ -926,7 +1045,7 @@
@subsubsection Linear warps to be called explicitly
Linear warps include operations like rotation, scaling, sheer, etc.
-@@ -20362,7 +20362,7 @@ The Statistics program is designed for s
+@@ -24618,7 +24618,7 @@ The Statistics program is designed for s
* Sigma clipping:: Definition of @mymath{\sigma}-clipping.
* Least squares fitting:: Fitting with various parametric functions.
* Sky value:: Definition and derivation of the Sky value.
@@ -935,7 +1054,7 @@
@end menu
-@@ -20376,7 +20376,7 @@ So on the horizontal axis we have the bi
+@@ -24632,7 +24632,7 @@ So on the horizontal axis we have the bi
You can use it to get a general view of the distribution: which values have been repeated the most? how close/far are the most significant bins? Are there more values in the larger part of the range of the dataset, or in the lower part? Similarly, many very important properties about the dataset can be deduced from a visual inspection of the histogram.
In the Statistics program, the histogram can be either output to a table to plot with your favorite plotting program at footnote{
We recommend @url{http://pgfplots.sourceforge.net/,PGFPlots} which generates your plots directly within @TeX{} (the same tool that generates your document).},
@@ -944,7 +1063,7 @@
@cindex Intervals, histogram
@cindex Bin width, histogram
-@@ -20397,7 +20397,7 @@ Normalizing a cumulative frequency plot
+@@ -24653,7 +24653,7 @@ Normalizing a cumulative frequency plot
Unlike the histogram which has a limited number of bins, ideally the cumulative frequency plot should have one point for every data element.
Even in small datasets (for example, a @mymath{200\times200} image) this will result in an unreasonably large number of points to plot (40000)! As a result, for practical reasons, it is common to only store its value on a certain number of points (intervals) in the input range rather than the whole dataset, so you should determine the number of bins you want when asking for a cumulative frequency plot.
In Gnuastro (and thus the Statistics program), the number reported for each bin is the total number of data points until the larger interval value for that bin.
@@ -953,16 +1072,16 @@
So as a summary, both the histogram and cumulative frequency plot in Statistics will work with bins.
Within each bin/interval, the lower value is considered to be within then bin (it is inclusive), but its larger value is not (it is exclusive).
-@@ -20967,7 +20967,7 @@ GSL provides many methods to reject outl
- For their full list, see the description of @option{--fitrobust} in @ref{Fitting options}.
- For a description of the outlier rejection methods, see the @url{https://www.gnu.org/software/gsl/doc/html/lls.html#c.gsl_multifit_robust_workspace, GSL manual}.
+@@ -25255,7 +25255,7 @@ You may have noticed that unlike the cas
+ This is becasue @option{--quiet} and @option{--fitestimate} were called together.
+ In this case, because all the fitting parameters are written as FITS keywords, because of the @option{--quiet} option, they are no longer printed on standard output.
- at node Sky value, Invoking aststatistics, Least squares fitting, Statistics
+ at node Sky value, aststatistics, Least squares fitting, Statistics
@subsection Sky value
@cindex Sky
-@@ -21225,7 +21225,7 @@ The same goes for the @option{--checksky
+@@ -25504,7 +25504,7 @@ The same goes for the @option{--checksky
@@ -971,7 +1090,7 @@
@subsection Invoking Statistics
Statistics will print statistical measures of an input dataset (table column or image).
-@@ -21345,7 +21345,7 @@ If you want another statistical paramete
+@@ -25624,7 +25624,7 @@ If you want another statistical paramete
* Statistics on tiles:: Possible to do single-valued measurements on tiles.
@end menu
@@ -980,7 +1099,7 @@
@subsubsection Input to Statistics
The following set of options are for specifying the input/outputs of Statistics.
-@@ -21386,7 +21386,7 @@ The quantile of each horizontal axis val
+@@ -25665,7 +25665,7 @@ The quantile of each horizontal axis val
@end table
@@ -989,7 +1108,7 @@
@subsubsection Single value measurements
@table @option
-@@ -21543,7 +21543,7 @@ Standard deviation after applying @mymat
+@@ -25822,7 +25822,7 @@ Standard deviation after applying @mymat
@end table
@@ -998,7 +1117,7 @@
@subsubsection Generating histograms and cumulative freq.
The list of options below are for those statistical operations that output more than one value.
-@@ -21711,7 +21711,7 @@ Similar to @option{--onebinstart}, but f
+@@ -25990,7 +25990,7 @@ Similar to @option{--onebinstart}, but f
@end table
@@ -1007,7 +1126,7 @@
@subsubsection Fitting options
With the options below, you can customize the least squares fitting features of Statistics.
-@@ -21879,7 +21879,7 @@ See @ref{Selecting table columns}.
+@@ -26163,7 +26163,7 @@ See @ref{Selecting table columns}.
@@ -1016,7 +1135,7 @@
@subsubsection Contour options
Contours are useful to highlight the 2D shape of a certain flux level over an image.
-@@ -21901,7 +21901,7 @@ If the image has World Coordinate System
+@@ -26185,7 +26185,7 @@ If the image has World Coordinate System
Note that currently, this is a very crude/simple implementation, please let us know if you find problematic situations so we can fix it.
@end table
@@ -1025,7 +1144,7 @@
@subsubsection Statistics on tiles
All the options described until now were from the first class of operations discussed above: those that treat the whole dataset as one.
-@@ -22055,7 +22055,7 @@ You can then directly feed NoiseChisel's
+@@ -26344,7 +26344,7 @@ You can then directly feed NoiseChisel's
Thanks to the published papers mentioned above, there is no need to provide a more complete introduction to NoiseChisel in this book.
However, published papers cannot be updated any more, but the software has evolved/changed.
The changes since publication are documented in @ref{NoiseChisel changes after publication}.
@@ -1034,7 +1153,7 @@
As discussed above, detection is one of the most important steps for your scientific result.
It is therefore very important to obtain a good understanding of NoiseChisel (and afterwards @ref{Segment} and @ref{MakeCatalog}).
-@@ -22067,17 +22067,17 @@ In the meantime, they also show the modu
+@@ -26356,17 +26356,17 @@ In the meantime, they also show the modu
Defining colors is a very common process in most science-cases.
Therefore it is also recommended to (patiently) complete that tutorial for optimal usage of NoiseChisel in conjunction with all the other Gnuastro programs.
@ref{Detecting large extended targets} shows you can optimize NoiseChisel's settings for very extended objects to successfully carve out to signal-to-noise ratio levels of below 1/10.
@@ -1056,7 +1175,7 @@
@subsection NoiseChisel changes after publication
NoiseChisel was initially introduced in @url{https://arxiv.org/abs/1505.01664, Akhlaghi and Ichikawa [2015]} and updates after the first four years were published in @url{https://arxiv.org/abs/1909.11230, Akhlaghi [2019]}.
-@@ -22105,7 +22105,7 @@ In our tests, this gave a much improved
+@@ -26395,7 +26395,7 @@ In our tests, this gave a much improved
@@ -1065,7 +1184,7 @@
@subsection Invoking NoiseChisel
NoiseChisel will detect signal in noise producing a multi-extension dataset containing a binary detection map which is the same size as the input.
-@@ -22232,7 +22232,7 @@ Finally, in @ref{NoiseChisel output} the
+@@ -26528,7 +26528,7 @@ Finally, in @ref{NoiseChisel output} the
* NoiseChisel output:: NoiseChisel's output options and format.
@end menu
@@ -1074,7 +1193,7 @@
@subsubsection NoiseChisel input
The options here can be used to configure the inputs and output of NoiseChisel, along with some general processing options.
-@@ -22336,7 +22336,7 @@ Except for the tile size, all the other
+@@ -26632,7 +26632,7 @@ Except for the tile size, all the other
The format is identical to that of the @option{--tilesize} option that is discussed in that section.
@end table
@@ -1083,7 +1202,7 @@
@subsubsection Detection options
Detection is the process of separating the pixels in the image into two groups: 1) Signal, and 2) Noise.
-@@ -22679,7 +22679,7 @@ By default the output will have the same
+@@ -26986,7 +26986,7 @@ By default the output will have the same
@@ -1092,8 +1211,8 @@
@subsubsection NoiseChisel output
NoiseChisel's output is a multi-extension FITS file.
-@@ -22850,14 +22850,14 @@ Those papers cannot be updated any more
- for example, Segment became a separate program (from NoiseChisel) in 2018 (after those papers were published).
+@@ -27157,14 +27157,14 @@ Those papers cannot be updated any more
+ For example, Segment became a separate program (from NoiseChisel) in 2018 (after those papers were published).
Therefore this book is the definitive reference.
@c To help in the transition from those papers to the software you are using, see @ref{Segment changes after publication}.
-Finally, in @ref{Invoking astsegment}, we will discuss Segment's inputs, outputs and configuration options.
@@ -1110,7 +1229,7 @@
@c @subsection Segment changes after publication
@c Segment's main algorithm and working strategy were initially defined and introduced in Section 3.2 of @url{https://arxiv.org/abs/1505.01664, Akhlaghi and Ichikawa [2015]} and @url{https://arxiv.org/abs/1909.11230, Akhlaghi [2019]}.
-@@ -22869,7 +22869,7 @@ Finally, in @ref{Invoking astsegment}, w
+@@ -27176,7 +27176,7 @@ Finally, in @ref{Invoking astsegment}, w
@c The aim of this section is to make the transition from the paper to your installed version, as smooth as possible through the list below.
@c For a more detailed list of changes in previous Gnuastro releases/versions, please follow the @file{NEWS} file at footnote{The @file{NEWS} file is present in the released Gnuastro tarball, see @ref{Release tarball}.}.
@@ -1119,7 +1238,7 @@
@subsection Invoking Segment
Segment will identify substructure within the detected regions of an input image.
-@@ -22922,7 +22922,7 @@ Finally, in @ref{Segment output}, we wil
+@@ -27229,7 +27229,7 @@ Finally, in @ref{Segment output}, we wil
* Segment output:: Outputs of Segment
@end menu
@@ -1128,7 +1247,7 @@
@subsubsection Segment input
Besides the input dataset (for example, astronomical image), Segment also needs to know the Sky standard deviation and the regions of the dataset that it should segment.
-@@ -23064,7 +23064,7 @@ It is important for them to be much larg
+@@ -27371,7 +27371,7 @@ It is important for them to be much larg
@end table
@@ -1137,7 +1256,7 @@
@subsubsection Segmentation options
The options below can be used to configure every step of the segmentation process in the Segment program.
-@@ -23178,7 +23178,7 @@ This behavior can be disabled with @opti
+@@ -27485,7 +27485,7 @@ This behavior can be disabled with @opti
@end table
@@ -1146,7 +1265,7 @@
@subsubsection Segment output
The main output of Segment are two label datasets (with integer types, separating the dataset's elements into different classes).
-@@ -23322,7 +23322,7 @@ Similarly, the sum of all these pixels w
+@@ -27629,7 +27629,7 @@ Similarly, the sum of all these pixels w
Pixels with labels equal to, or smaller than, zero will be ignored by MakeCatalog.
In other words, the number of rows in MakeCatalog's output is already known before running it (the maximum value of the labeled dataset).
@@ -1155,25 +1274,16 @@
A very important factor in any measurement is understanding its validity range, or limits.
Therefore in @ref{Quantifying measurement limits}, we will discuss how to estimate the reliability of the detection and basic measurements.
This section will continue with a derivation of elliptical parameters from the labeled datasets in @ref{Measuring elliptical parameters}.
-@@ -23334,7 +23334,7 @@ For those who feel MakeCatalog's existin
- * Quantifying measurement limits:: For comparing different catalogs.
+@@ -27642,7 +27642,7 @@ For those who feel MakeCatalog's existin
* Measuring elliptical parameters:: Estimating elliptical parameters.
* Adding new columns to MakeCatalog:: How to add new columns.
+ * MakeCatalog measurements:: List of all the measurements/columns by MakeCatalog.
-* Invoking astmkcatalog:: Options and arguments to MakeCatalog.
+* astmkcatalog:: Options and arguments to MakeCatalog.
@end menu
@node Detection and catalog production, Brightness flux magnitude, MakeCatalog, MakeCatalog
-@@ -23926,7 +23926,7 @@ When the object cannot be represented as
-
-
-
-- at node Adding new columns to MakeCatalog, Invoking astmkcatalog, Measuring elliptical parameters, MakeCatalog
-+ at node Adding new columns to MakeCatalog, astmkcatalog, Measuring elliptical parameters, MakeCatalog
- @subsection Adding new columns to MakeCatalog
-
- MakeCatalog is designed to allow easy addition of different measurements over a labeled image (see @url{https://arxiv.org/abs/1611.06387v1, Akhlaghi [2016]}).
-@@ -23942,7 +23942,7 @@ These two passes are designed to be raw
+@@ -28422,7 +28422,7 @@ These two passes are designed to be raw
This will allow parallel processing and simplicity/clarity.
So if your new calculation, needs new raw information from the pixels, then you will need to also modify the respective @code{mkcatalog_first_pass} and @code{mkcatalog_second_pass} functions (both in @file{bin/mkcatalog/mkcatalog.c}) and define new raw table columns in @file{main.h} (hopefully the comments in the code are clear enough).
@@ -1182,16 +1292,34 @@
This allows a particular column/option to be easily found in all steps.
Therefore in adding your new option, be sure to keep it in the same relative place in the list in all the separate places (it does not necessarily have to be in the end), and near conceptually similar options.
-@@ -23994,7 +23994,7 @@ Update this manual and add a description
+@@ -28473,7 +28473,7 @@ Update this manual and add a description
-- at node Invoking astmkcatalog, , Adding new columns to MakeCatalog, MakeCatalog
-+ at node astmkcatalog, , Adding new columns to MakeCatalog, MakeCatalog
+- at node MakeCatalog measurements, Invoking astmkcatalog, Adding new columns to MakeCatalog, MakeCatalog
++ at node MakeCatalog measurements, astmkcatalog, Adding new columns to MakeCatalog, MakeCatalog
+ @subsection MakeCatalog measurements
+
+ MakeCatalog's output measurements/columns can be specified using command-line options (@ref{Options}).
+@@ -28897,7 +28897,7 @@ To know when this measurement is reasona
+
+ @item --sigclip-mean-sb
+ Surface brightness (over 1 pixel's area in arcsec at mymath{^2}) of the sigma-clipped mean value of the pixel values distribution associated to each label (object or clump).
+-This is useful in scenarios where your labels have approximately @emph{constant} surface brightness values @emph{after} after removing outliers: for example in a radial profile, see @ref{Invoking astscript-radial-profile}).
++This is useful in scenarios where your labels have approximately @emph{constant} surface brightness values @emph{after} after removing outliers: for example in a radial profile, see @ref{astscript-radial-profile}).
+
+ In other scenarios it should be used with extreme care.
+ For example over the full area of a galaxy/star the pixel distribution is not constant (or symmetric after adding noise), their pixel distributions are inherently skewed (with fewer pixels in the center, having a very large value and many pixels in the outer parts having lower values).
+@@ -29178,7 +29178,7 @@ These measurements contain @code{-other-
+
+
+
+- at node Invoking astmkcatalog, , MakeCatalog measurements, MakeCatalog
++ at node astmkcatalog, , MakeCatalog measurements, MakeCatalog
@subsection Invoking MakeCatalog
MakeCatalog will do measurements and produce a catalog from a labeled dataset and optional values dataset(s).
-@@ -24048,7 +24048,7 @@ Finally, in @ref{MakeCatalog output} the
+@@ -29231,7 +29231,7 @@ Finally, in @ref{MakeCatalog output} the
* MakeCatalog output:: File names of MakeCatalog's output table.
@end menu
@@ -1200,7 +1328,7 @@
@subsubsection MakeCatalog inputs and basic settings
MakeCatalog works by using a localized/labeled dataset (see @ref{MakeCatalog}).
-@@ -24102,7 +24102,7 @@ In such cases, you can use Gnuastro's Ar
+@@ -29275,7 +29275,7 @@ In such cases, you can use Gnuastro's Ar
@end example
To summarize: if the input file to MakeCatalog is the default/full output of Segment (see @ref{Segment output}) you do not have to worry about any of the @option{--*file} options below.
@@ -1209,34 +1337,25 @@
To feed NoiseChisel's output into MakeCatalog, just change the labeled dataset's header (with @option{--hdu=DETECTIONS}).
The full list of input dataset options and general setting options are described below.
-@@ -24201,7 +24201,7 @@ But in practice real images my be over-s
+@@ -29381,7 +29381,7 @@ If you use this option, it will have 13
-- at node Upper-limit settings, MakeCatalog measurements, MakeCatalog inputs and basic settings, Invoking astmkcatalog
-+ at node Upper-limit settings, MakeCatalog measurements, MakeCatalog inputs and basic settings, astmkcatalog
+- at node Upper-limit settings, MakeCatalog output, MakeCatalog inputs and basic settings, Invoking astmkcatalog
++ at node Upper-limit settings, MakeCatalog output, MakeCatalog inputs and basic settings, astmkcatalog
@subsubsection Upper-limit settings
The upper-limit magnitude was discussed in @ref{Quantifying measurement limits}.
-@@ -24288,7 +24288,7 @@ The total number of rows is thus unknown
+@@ -29476,7 +29476,7 @@ For more on these positional columns, se
@end table
-- at node MakeCatalog measurements, MakeCatalog output, Upper-limit settings, Invoking astmkcatalog
-+ at node MakeCatalog measurements, MakeCatalog output, Upper-limit settings, astmkcatalog
- @subsubsection MakeCatalog measurements
-
- The final group of options particular to MakeCatalog are those that specify which measurements/columns should be written into the final output table.
-@@ -24849,7 +24849,7 @@ If you use this option, it will have 13
-
-
-
-- at node MakeCatalog output, , MakeCatalog measurements, Invoking astmkcatalog
-+ at node MakeCatalog output, , MakeCatalog measurements, astmkcatalog
+- at node MakeCatalog output, , Upper-limit settings, Invoking astmkcatalog
++ at node MakeCatalog output, , Upper-limit settings, astmkcatalog
@subsubsection MakeCatalog output
After it has completed all the requested measurements (see @ref{MakeCatalog measurements}), MakeCatalog will store its measurements in table(s).
If an output filename is given (see @option{--output} in @ref{Input output options}), the format of the table will be deduced from the name.
-@@ -24974,10 +24974,10 @@ The aperture can be a circle or an ellip
+@@ -29601,10 +29601,10 @@ The aperture can be a circle or an ellip
@menu
* Matching algorithms:: Different ways to find the match
@@ -1249,7 +1368,7 @@
@subsection Matching algorithms
Matching involves two catalogs, let's call them catalog A (with N rows) and catalog B (with M rows).
-@@ -25090,7 +25090,7 @@ If you always need to match against one
+@@ -29717,7 +29717,7 @@ If you always need to match against one
Therefore you can save its k-d tree into a file and simply give it to later calls, like the example given in the description of the k-d algorithm mentioned above.
@end itemize
@@ -1258,7 +1377,7 @@
@subsection Invoking Match
When given two catalogs, Match finds the rows that are nearest to each other within an input aperture.
-@@ -25476,7 +25476,7 @@ After all the transformations are applie
+@@ -30102,7 +30102,7 @@ After all the transformations are applie
* Modeling basics:: Astronomical modeling basics.
* If convolving afterwards:: Considerations for convolving later.
* Profile magnitude:: Definition of total profile magnitude.
@@ -1267,7 +1386,7 @@
@end menu
-@@ -25806,7 +25806,7 @@ You also have to shift all the pixel pos
+@@ -30432,7 +30432,7 @@ You also have to shift all the pixel pos
After convolution, you can crop the outer @mymath{n} pixels with the section crop box specification of Crop: @option{--section=n+1:*-n,n+1:*-n} (according to the FITS standard, counting is from 1 so we use @code{n+1}) assuming your PSF is a square, see @ref{Crop section syntax}.
This will also remove all discrete Fourier transform artifacts (blurred sides) from the final image.
@@ -1276,7 +1395,7 @@
-@@ -25817,7 +25817,7 @@ To facilitate this shift, MakeProfiles h
+@@ -30443,7 +30443,7 @@ To facilitate this shift, MakeProfiles h
@@ -1284,8 +1403,8 @@
+ at node Profile magnitude, astmkprof, If convolving afterwards, MakeProfiles
@subsection Profile magnitude
- @cindex Brightness
-@@ -25843,7 +25843,7 @@ If not, only the overlapping pixels of t
+ @cindex Truncation radius
+@@ -30468,7 +30468,7 @@ If not, only the overlapping pixels of t
@@ -1294,7 +1413,7 @@
@subsection Invoking MakeProfiles
MakeProfiles will make any number of profiles specified in a catalog either individually or in one image.
-@@ -25928,7 +25928,7 @@ Please see @ref{Sufi simulates a detecti
+@@ -30553,7 +30553,7 @@ Please see @ref{Sufi simulates a detecti
* MakeProfiles log file:: A description of the optional log file.
@end menu
@@ -1303,7 +1422,7 @@
@subsubsection MakeProfiles catalog
The catalog containing information about each profile can be in the FITS ASCII, FITS binary, or plain text formats (see @ref{Tables}).
The latter can also be provided using standard input (see @ref{Standard input}).
-@@ -26068,7 +26068,7 @@ If @option{--tunitinp} is given, this va
+@@ -30693,7 +30693,7 @@ If @option{--tunitinp} is given, this va
@end table
@@ -1312,7 +1431,7 @@
@subsubsection MakeProfiles profile settings
The profile parameters that differ between each created profile are specified through the columns in the input catalog and described in @ref{MakeProfiles catalog}.
-@@ -26278,7 +26278,7 @@ Because the maximum operator is independ
+@@ -30929,7 +30929,7 @@ Because the maximum operator is independ
@end table
@@ -1321,7 +1440,7 @@
@subsubsection MakeProfiles output dataset
MakeProfiles takes an input catalog uses basic properties that are defined there to build a dataset, for example, a 2D image containing the profiles in the catalog.
In @ref{MakeProfiles catalog} and @ref{MakeProfiles profile settings}, the catalog and profile settings were discussed.
-@@ -26456,7 +26456,7 @@ MakeProfiles will not complain if you us
+@@ -31107,7 +31107,7 @@ MakeProfiles will not complain if you us
@end table
@@ -1330,36 +1449,7 @@
@subsubsection MakeProfiles log file
Besides the final merged dataset of all the profiles, or the individual datasets (see @ref{MakeProfiles output dataset}), if the @option{--log} option is called MakeProfiles will also create a log file in the current directory (where you run MockProfiles).
-@@ -26499,16 +26499,16 @@ If an individual image was created, this
- @cindex Noise
- Real data are always buried in noise, therefore to finalize a simulation of real data (for example, to test our observational algorithms) it is essential to add noise to the mock profiles created with MakeProfiles, see @ref{MakeProfiles}.
- Below, the general principles and concepts to help understand how noise is quantified is discussed.
--MakeNoise options and argument are then discussed in @ref{Invoking astmknoise}.
-+MakeNoise options and argument are then discussed in @ref{astmknoise}.
-
- @menu
- * Noise basics:: Noise concepts and definitions.
--* Invoking astmknoise:: Options and arguments to MakeNoise.
-+* astmknoise:: Options and arguments to MakeNoise.
- @end menu
-
-
-
-- at node Noise basics, Invoking astmknoise, MakeNoise, MakeNoise
-+ at node Noise basics, astmknoise, MakeNoise, MakeNoise
- @subsection Noise basics
-
- @cindex Noise
-@@ -26730,7 +26730,7 @@ In the example above, did you notice how
- However, if @option{--envseed} was given, both printed seeds would be the same.
-
-
-- at node Invoking astmknoise, , Noise basics, MakeNoise
-+ at node astmknoise, , Noise basics, MakeNoise
- @subsection Invoking MakeNoise
-
- MakeNoise will add noise to an existing image.
-@@ -26848,7 +26848,7 @@ Even higher-level analysis is still need
+@@ -31172,7 +31172,7 @@ Even higher-level analysis is still need
To derive higher-level information regarding our sources in extra-galactic astronomy, cosmological calculations are necessary.
In Gnuastro, CosmicCalculator is in charge of such calculations.
@@ -1368,16 +1458,16 @@
In @ref{Distance on a 2D curved space} the basic idea of understanding distances in a curved and expanding 2D universe (which we can visualize) are reviewed.
Having solidified the concepts there, in @ref{Extending distance concepts to 3D}, the formalism is extended to the 3D universe we are trying to study in our research.
-@@ -26858,7 +26858,7 @@ There are many books thoroughly deriving
+@@ -31182,7 +31182,7 @@ There are many books thoroughly deriving
@menu
- * Distance on a 2D curved space:: Distances in 2D for simplicity
+ * Distance on a 2D curved space:: Distances in 2D for simplicity.
* Extending distance concepts to 3D:: Going to 3D (our real universe).
--* Invoking astcosmiccal:: How to run CosmicCalculator
-+* astcosmiccal:: How to run CosmicCalculator
+-* Invoking astcosmiccal:: How to run CosmicCalculator.
++* astcosmiccal:: How to run CosmicCalculator.
@end menu
@node Distance on a 2D curved space, Extending distance concepts to 3D, CosmicCalculator, CosmicCalculator
-@@ -27005,7 +27005,7 @@ We can thus parameterize the change in d
+@@ -31329,7 +31329,7 @@ We can thus parameterize the change in d
@dispmath{ds^2=c^2dt^2-a^2(t)ds_s^2 = c^2dt^2-a^2(t)(d\chi^2+r^2d\phi^2).}
@@ -1386,7 +1476,7 @@
@subsection Extending distance concepts to 3D
The concepts of @ref{Distance on a 2D curved space} are here extended to a 3D space that @emph{might} be curved.
-@@ -27054,7 +27054,7 @@ In a non-static universe (with a scale f
+@@ -31378,7 +31378,7 @@ In a non-static universe (with a scale f
@@ -1395,7 +1485,7 @@
@subsection Invoking CosmicCalculator
CosmicCalculator will calculate cosmological variables based on the input parameters.
-@@ -27112,7 +27112,7 @@ CosmicCalculator also has features to he
+@@ -31436,7 +31436,7 @@ CosmicCalculator also has features to he
* CosmicCalculator spectral line calculations:: How they get affected by redshift.
@end menu
@@ -1404,7 +1494,7 @@
@subsubsection CosmicCalculator input options
The inputs to CosmicCalculator can be specified with the following options:
-@@ -27330,9 +27330,9 @@ You can get this list on the command-lin
+@@ -31495,9 +31495,9 @@ $ astcosmiccal --listlines
@@ -1416,7 +1506,7 @@
The full list of calculations can be useful when you do not want any specific value, but just a general view.
In other contexts (for example, in a batch script or during a discussion), you know exactly what you want and do not want to be distracted by all the extra information.
-@@ -27447,7 +27447,7 @@ The comoving volume in Megaparsecs cube
+@@ -31612,7 +31612,7 @@ The comoving volume in Megaparsecs cube
@@ -1425,7 +1515,7 @@
@subsubsection CosmicCalculator spectral line calculations
@cindex Rest frame wavelength
-@@ -27638,10 +27638,10 @@ But we will leave the implementation of
+@@ -31827,10 +31827,10 @@ But we will leave the implementation of
@end itemize
@menu
@@ -1438,7 +1528,7 @@
@subsection Invoking astscript-sort-by-night
This installed script will read a FITS date formatted value from the given keyword, and classify the input FITS files into individual nights.
-@@ -27793,10 +27793,10 @@ Gnuastro's @file{astscript-radial-profil
+@@ -31982,10 +31982,10 @@ Gnuastro's @file{astscript-radial-profil
This script uses @ref{MakeProfiles} to generate elliptical apertures with the values equal to the distance from the center of the object and @ref{MakeCatalog} for measuring the values over the apertures.
@menu
@@ -1451,7 +1541,16 @@
@subsection Invoking astscript-radial-profile
This installed script will measure the radial profile of an object within an image.
-@@ -28003,7 +28003,7 @@ This is necessary when you request measu
+@@ -32240,7 +32240,7 @@ This option can also be called multiple
+ @cartouche
+ @noindent
+ @strong{Masking background/foreground objects:} For crude rejection of outliers, you can use sigma-clipping using MakeCatalog measurements like @option{--sigclip-mean} or @option{--sigclip-mean-sb} (see @ref{MakeCatalog measurements}).
+-To properly mask the effect of background/foreground objects from your target object's radial profile, you can use @command{astscript-psf-stamp} script, see @ref{Invoking astscript-psf-stamp}, and feed it the output of @ref{Segment}.
++To properly mask the effect of background/foreground objects from your target object's radial profile, you can use @command{astscript-psf-stamp} script, see @ref{astscript-psf-stamp}, and feed it the output of @ref{Segment}.
+ This script will mask unwanted objects from the image that is later used to measure the radial profile.
+ @end cartouche
+
+@@ -32273,7 +32273,7 @@ This is necessary when you request measu
@itemx --zeroisnotblank
Account for zero-valued pixels in the profile.
By default, such pixels are not considered (when this script crops the necessary region of the image before generating the profile).
@@ -1459,8 +1558,8 @@
+The long format of this option is identical to a similarly named option in Crop (see @ref{astcrop}).
When this option is called, it is passed directly to Crop, therefore the zero-valued pixels are not considered as blank and used in the profile creation.
- @item -v INT
-@@ -28075,10 +28075,10 @@ In this section we describe a simple ins
+ @item -i FLT/STR
+@@ -32332,10 +32332,10 @@ In this section we describe a simple ins
SAO DS9 at footnote{@url{http://ds9.si.edu}} is one of the most common FITS image visualization tools in astronomy and is free software.
@menu
@@ -1473,7 +1572,7 @@
@subsection Invoking astscript-ds9-region
This installed script will read two positional columns within an input table and generate an SAO DS9 region file to visualize the position of the given objects over an image.
-@@ -28239,13 +28239,13 @@ If @file{foo.fits} only consists of one
+@@ -32496,13 +32496,13 @@ If @file{foo.fits} only consists of one
On the other hand, for visualizing the contents of tables (that are also commonly stored in the FITS format), you need to call a different software (most commonly, people use TOPCAT, see @ref{TOPCAT}).
And to make things more inconvenient, by default both of these are only installed as command-line software, so while you are navigating in your GUI, you need to open a terminal there, and run these commands.
@@ -1490,7 +1589,61 @@
@subsection Invoking astscript-fits-view
Given any number of FITS files, this script will either open SAO DS9 (for images or cubes) or TOPCAT (for tables) to visualize their contents in a graphic user interface (GUI).
-@@ -28378,14 +28378,14 @@ The tutorial uses a real dataset and inc
+@@ -32659,13 +32659,13 @@ Since a zero point of 0 was used for the
+
+ In the ``Tutorials'' chapter of this Gnuastro book, there are two tutorials dedicated to the usage of this script.
+ The first uses an image as a reference (@ref{Zero point tutorial with reference image}) and the second uses a catalog (@ref{Zero point tutorial with reference catalog}).
+-For the full set of options an a detailed description of each, see @ref{Invoking astscript-zeropoint}.
++For the full set of options an a detailed description of each, see @ref{astscript-zeropoint}.
+
+ @menu
+-* Invoking astscript-zeropoint:: How to call the script
++* astscript-zeropoint:: How to call the script
+ @end menu
+
+- at node Invoking astscript-zeropoint, , Zero point estimation, Zero point estimation
++ at node astscript-zeropoint, , Zero point estimation, Zero point estimation
+ @subsection Invoking astscript-zeropoint
+ This installed script will calculate the zero point of an input image to calibrate it.
+ The reference can be an image or catalog (which have been previously calibrated)
+@@ -32716,7 +32716,7 @@ The @option{--keeptmp} option can be use
+ * zero point options:: List and details of options.
+ @end menu
+
+- at node zero point output, zero point options, Invoking astscript-zeropoint, Invoking astscript-zeropoint
++ at node zero point output, zero point options, astscript-zeropoint, astscript-zeropoint
+ @subsubsection astscript-zeropoint output
+
+ The output will be a multi-extension FITS table.
+@@ -32753,7 +32753,7 @@ $ astfits out.fits -h1 | grep ^ZP
+ $ astfits jplus-zeropoint.fits -h1 --keyvalue=ZPVALUE
+ @end example
+
+- at node zero point options, , zero point output, Invoking astscript-zeropoint
++ at node zero point options, , zero point output, astscript-zeropoint
+ @subsubsection astscript-zeropoint options
+
+ All the operating phases of the this script can be customized through the options below.
+@@ -32936,15 +32936,15 @@ For example in the XDF survey is focused
+ These are very small objects and within that small footprint (of just 1 arcmin) we have thousands of them.
+ However, the LIGHTS survey is focused on the halos of large nearby galaxies (that can be more than 10 arcminutes wide!).
+
+-In @ref{Invoking astscript-pointing-simulate} of Gnuastro's @ref{Installed scripts} is described in detail.
++In @ref{astscript-pointing-simulate} of Gnuastro's @ref{Installed scripts} is described in detail.
+ This script is designed to simplify the process of selecting the best pointing pattern for your observation strategy.
+ For a practical tutorial on using this script, see @ref{Pointing pattern design}.
+
+ @menu
+-* Invoking astscript-pointing-simulate:: Options and running mode.
++* astscript-pointing-simulate:: Options and running mode.
+ @end menu
+
+- at node Invoking astscript-pointing-simulate, , Pointing pattern simulation, Pointing pattern simulation
++ at node astscript-pointing-simulate, , Pointing pattern simulation, Pointing pattern simulation
+ @subsection Invoking astscript-pointing-simulate
+
+ This installed script will simulate a final stacked image from a certain pointing pattern (given as a table).
+@@ -33147,14 +33147,14 @@ The tutorial uses a real dataset and inc
@menu
* Overview of the PSF scripts:: Summary of concepts and methods
@@ -1511,8 +1664,8 @@
@subsection Overview of the PSF scripts
To obtain an extended and non-parametric PSF, several steps are necessary and we will go through them here.
-@@ -28413,11 +28413,11 @@ But you cannot blindly use all the stars
- for example, we do not want contamination from other bright, and nearby objects.
+@@ -33182,11 +33182,11 @@ But you cannot blindly use all the stars
+ For example, we do not want contamination from other bright, and nearby objects.
The first script below is therefore designed for selecting only good star candidates in your image.
It will use different criteria, for example, good parallax (where available, to avoid confusion with galaxies), not being near to bright stars, axis ratio, etc.
-For more on this script, see @ref{Invoking astscript-psf-select-stars}.
@@ -1525,14 +1678,14 @@
After obtaining a set of star stamps, they can be stacked for obtaining the combined PSF from many stars (for example, with @ref{Stacking operators}).
In the combined PSF, the masked background objects of each star's image will be covered and the signal-to-noise ratio will increase, giving a very nice view of the ``clean'' PSF.
-@@ -28433,16 +28433,16 @@ But in the example dataset from the tuto
+@@ -33202,16 +33202,16 @@ But in the example dataset from the tuto
Once clean stacks of different parts of the PSF have been constructed through the steps above, it is therefore necessary to blend them all into one.
This is done by finding a common radial region in both, and scaling the inner region by a factor to add with the outer region.
-This is not trivial, therefore, a third script is in charge of it, see @ref{Invoking astscript-psf-unite}.
+This is not trivial, therefore, a third script is in charge of it, see @ref{astscript-psf-unite}.
- Having constructed the PSF as described above (or by any other procedure), it can be scaled to the brightness of the various stars in the image to get subtracted (and thus remove the extended/bright wings; better showing the background objects of interest).
+ Having constructed the PSF as described above (or by any other procedure), it can be scaled to the magnitude of the various stars in the image to get subtracted (and thus remove the extended/bright wings; better showing the background objects of interest).
Note that the absolute flux of a PSF is meaningless (and in fact, it is usually normalized to have a total sum of unity!), so it should be scaled.
We therefore have another script that will calculate the scale (multiplication) factor of the PSF for each star.
-For more on the scaling script, see @ref{Invoking astscript-psf-scale-factor}.
@@ -1545,7 +1698,7 @@
As mentioned above, in the following sections, each script has its own documentation and list of options for very detailed customization (if necessary).
But if you are new to these scripts, before continuing, we recommend that you do the tutorial @ref{Building the extended PSF}.
-@@ -28452,7 +28452,7 @@ Just do not forget to run every command,
+@@ -33221,7 +33221,7 @@ Just do not forget to run every command,
@@ -1553,8 +1706,8 @@
+ at node astscript-psf-select-stars, astscript-psf-stamp, Overview of the PSF scripts, PSF construction and subtraction
@subsection Invoking astscript-psf-select-stars
This installed script will select good star candidates for constructing a PSF.
- It will consider stars within a given range of brightness without nearby contaminant objects.
-@@ -28608,7 +28608,7 @@ The output name of the final catalog con
+ It will consider stars within a given range of magnitudes without nearby contaminant objects.
+@@ -33378,7 +33378,7 @@ The output name of the final catalog con
@@ -1563,7 +1716,7 @@
@subsection Invoking astscript-psf-stamp
This installed script will generate a stamp of fixed size, centered at the provided coordinates (performing sub-pixel re-gridding if necessary) and normalized at a certain normalization radius.
Optionally, it will also mask all the other background sources.
-@@ -28793,7 +28793,7 @@ The main reason of setting this name is
+@@ -33566,7 +33566,7 @@ The main reason of setting this name is
@@ -1572,7 +1725,7 @@
@subsection Invoking astscript-psf-unite
This installed script will join two PSF images at a given radius.
This operation is commonly used when merging (uniting) the inner and outer parts of the PSF.
-@@ -28853,7 +28853,7 @@ The HDU/extension of the inner PSF (opti
+@@ -33626,7 +33626,7 @@ The HDU/extension of the inner PSF (opti
@itemx --scale=FLT
Factor by which the inner part (@option{--inner}) is multiplied.
This factor is necessary to put the two different parts of the PSF at the same flux level.
@@ -1581,7 +1734,7 @@
There is also a full tutorial on using all the @command{astscript-psf-*} installed scripts together, see @ref{Building the extended PSF}.
We recommend doing that tutorial before starting to work on your own datasets.
-@@ -28900,7 +28900,7 @@ This option is useful for debugging and
+@@ -33673,7 +33673,7 @@ This option is useful for debugging and
@@ -1590,7 +1743,7 @@
@subsection Invoking astscript-psf-scale-factor
This installed script will compute the multiplicative factor (scale) that is necessary to match the PSF to a given star.
The match in flux is done within a ring of pixels.
-@@ -28947,9 +28947,9 @@ To fully understand the process on first
+@@ -33720,9 +33720,9 @@ To fully understand the process on first
The most common use-cases of this scale factor are:
@enumerate
@item
@@ -1602,7 +1755,7 @@
@end enumerate
For a full tutorial on how to use this script along with the other @command{astscript-psf-*} scripts in Gnuastro, please see @ref{Building the extended PSF}.
-@@ -29030,7 +29030,7 @@ This option is useful for debugging and
+@@ -33803,7 +33803,7 @@ This option is useful for debugging and
@@ -1611,7 +1764,7 @@
@subsection Invoking astscript-psf-subtract
This installed script will put the provided PSF into a given position within the input image (implementing sub-pixel adjustments where necessary), and then it will subtract it.
It is aimed at modeling and subtracting the scattered light field of an input image.
-@@ -29114,7 +29114,7 @@ If the central position does not fall in
+@@ -33887,7 +33887,7 @@ If the central position does not fall in
@itemx --scale=FLT
Factor by which the PSF (@option{--psf}) is multiplied.
This factor is necessary to put the PSF with the desired flux level.
@@ -1620,7 +1773,7 @@
For a full tutorial on using the @command{astscript-psf-*} scripts together, see @ref{Building the extended PSF}.
@item -t
-@@ -29759,10 +29759,10 @@ Please see @ref{Optional dependencies} f
+@@ -34532,10 +34532,10 @@ Please see @ref{Optional dependencies} f
@end cartouche
@menu
@@ -1633,7 +1786,16 @@
@subsection Invoking BuildProgram
BuildProgram will compile and link a C source program with Gnuastro's library and all its dependencies, greatly facilitating the compilation and running of small programs that use Gnuastro's library.
-@@ -37430,7 +37430,7 @@ The @code{warp} library functions are @e
+@@ -37243,7 +37243,7 @@ This function is just for column informa
+
+ @deftypefun void gal_table_print_info (gal_data_t @code{*allcols}, size_t @code{numcols}, size_t @code{numrows}, char @code{*hdu_option_name})
+ Print the column information for all the columns (output of @code{gal_table_info}) to standard output.
+-The output is in the same format as this command with Gnuastro Table program (see @ref{Invoking asttable}):
++The output is in the same format as this command with Gnuastro Table program (see @ref{asttable}):
+ @example
+ $ asttable --info table.fits
+ @end example
+@@ -41782,7 +41782,7 @@ The @code{warp} library functions are @e
Acceptable fraction of output pixel that is covered by input pixels.
The value should be between 0 and 1 (inclusive).
If the area of an output pixel is covered by less than this fraction, its value will be @code{NaN}.
@@ -1642,7 +1804,7 @@
@item size_t edgesampling
Set the number of extra vertices along each edge of the output pixel's polygon to account for potential curvature due to projection or distortion.
-@@ -38442,7 +38442,7 @@ main(void)
+@@ -43040,7 +43040,7 @@ main(void)
@subsection Library demo - Warp to another image
Gnuastro's warp library (that you can access by including @file{gnuastro/warp.h}) allows you to resample an image from a grid to another entirely using the WCSLIB (while accounting for distortions if necessary; see @ref{Warp library}).
The Warp library uses a pixel-mixing or area-based resampling approach which is fully described in @ref{Resampling}.
@@ -1651,7 +1813,7 @@
For a related demo (where the output grid and WCS are constructed from scratch), see @ref{Library demo - Warp to new grid}.
In the example below, we are warping the @code{input.fits} file to the same pixel grid and WCS as @code{reference.fits} image (assuming it is in hdu @code{0}).
-@@ -38537,7 +38537,7 @@ main(void)
+@@ -43163,7 +43163,7 @@ main(void)
@subsection Library demo - Warp to new grid
Gnuastro's warp library (that you can access by including @file{gnuastro/warp.h}) allows you to resample an image from a grid to another entirely using the WCSLIB (while accounting for distortions if necessary; see @ref{Warp library}).
The Warp library uses a pixel-mixing or area-based resampling approach which is fully described in @ref{Resampling}.
@@ -1660,49 +1822,3 @@
For a related demo (where the output grid and WCS are imported from another file), see @ref{Library demo - Warp to another image}.
In the example below, we'll assume you have the SDSS image downloaded in @ref{Downloading and validating input data}.
-@@ -39509,7 +39509,7 @@ philosophy}.
-
- As you have already noticed for every program/library, it is very important
- that the basics of the science and technique be explained in separate
--subsections prior to the `Invoking Programname' subsection. If you are
-+subsections prior to the `Programname' subsection. If you are
- writing a new program or your addition to an existing program involves a
- new concept, also include such subsections and explain the concepts so a
- person completely unfamiliar with the concepts can get a general initial
-@@ -39527,7 +39527,7 @@ after reading a long text, nothing else.
- relevant/interesting for the reader, there is no page number limit/cost.
-
- It might also help if you start discussing the usage of your idea in the
--`Invoking ProgramName' subsection (explaining the options and arguments you
-+`ProgramName' subsection (explaining the options and arguments you
- have in mind) at this stage too. Actually starting to write it here will
- really help you later when you are coding.
-
-@@ -40629,22 +40629,22 @@ The following scripts are used to estima
-
- @table @code
- @item astscript-psf-select-stars
--(see @ref{Invoking astscript-psf-select-stars}) Find all the stars within an image that are suitable for constructing an extended PSF.
-+(see @ref{astscript-psf-select-stars}) Find all the stars within an image that are suitable for constructing an extended PSF.
- If the image has WCS, this script can automatically query Gaia to find the good stars.
-
- @item astscript-psf-stamp
--(see @ref{Invoking astscript-psf-stamp}) build a crop (stamp) of a certain width around a star at a certain coordinate in a larger image.
-+(see @ref{astscript-psf-stamp}) build a crop (stamp) of a certain width around a star at a certain coordinate in a larger image.
- This script will do sub-pixel re-positioning to make sure the star is centered and can optionally mask all other background sources).
-
- @item astscript-psf-scale-factor
--(see @ref{Invoking astscript-psf-scale-factor}) Given a PSF model, and the central coordinates of a star in an image, find the scale factor that has to be multiplied by the PSF to scale it to that star.
-+(see @ref{astscript-psf-scale-factor}) Given a PSF model, and the central coordinates of a star in an image, find the scale factor that has to be multiplied by the PSF to scale it to that star.
-
- @item astscript-psf-unite
--(see @ref{Invoking astscript-psf-unite}) Unite the various components of a PSF into one.
-+(see @ref{astscript-psf-unite}) Unite the various components of a PSF into one.
- Because of saturation and non-linearity, to get a good estimate of the extended PSF, it is necessary to construct various parts from different magnitude ranges.
-
- @item astscript-psf-subtract
--(see @ref{Invoking astscript-psf-subtract}) Given the model of a PSF and the central coordinates of a star in the image, do sub-pixel re-positioning of the PSF, scale it to the star and subtract it from the image.
-+(see @ref{astscript-psf-subtract}) Given the model of a PSF and the central coordinates of a star in the image, do sub-pixel re-positioning of the PSF, scale it to the star and subtract it from the image.
- @end table
-
- @end table
================================================================
---- gitweb:
http://git.pld-linux.org/gitweb.cgi/packages/gnuastro.git/commitdiff/c2c679f9d8b5530360bc136f40c2e07124bb87cc
More information about the pld-cvs-commit
mailing list