[packages/gnuastro] - updated to 0.17 - updated info patch
qboosh
qboosh at pld-linux.org
Sun Jun 26 22:22:44 CEST 2022
commit f89484e7f80d96436908549ade8eacf6c222a566
Author: Jakub Bogusz <qboosh at pld-linux.org>
Date: Sun Jun 26 22:24:13 2022 +0200
- updated to 0.17
- updated info patch
gnuastro-info.patch | 586 ++++++++++++++++++++++++++++++++++++----------------
gnuastro.spec | 24 ++-
2 files changed, 430 insertions(+), 180 deletions(-)
---
diff --git a/gnuastro.spec b/gnuastro.spec
index 871b9ad..6008534 100644
--- a/gnuastro.spec
+++ b/gnuastro.spec
@@ -1,3 +1,6 @@
+# TODO: ds9, topcat (for Suggests)
+# https://sites.google.com/cfa.harvard.edu/saoimageds9
+# http://www.star.bris.ac.uk/~mbt/topcat/
#
# Conditional build:
%bcond_without static_libs # static library
@@ -5,12 +8,12 @@
Summary: GNU Astronomy Utilities
Summary(pl.UTF-8): Narzędzia astronomiczne GNU
Name: gnuastro
-Version: 0.16
-Release: 2
+Version: 0.17
+Release: 1
License: GPL v3+
Group: Applications/Science
Source0: https://ftp.gnu.org/gnu/gnuastro/%{name}-%{version}.tar.lz
-# Source0-md5: 2685beb70cd4fddfe7ddc903d8ed661a
+# Source0-md5: dc1a9d0dba6550cc8eded53ecfd7535a
Patch0: %{name}-info.patch
Patch1: ac.patch
URL: http://www.gnu.org/software/gnuastro/
@@ -155,6 +158,12 @@ rm -rf $RPM_BUILD_ROOT
%attr(755,root,root) %{_bindir}/astnoisechisel
%attr(755,root,root) %{_bindir}/astquery
%attr(755,root,root) %{_bindir}/astscript-ds9-region
+%attr(755,root,root) %{_bindir}/astscript-fits-view
+%attr(755,root,root) %{_bindir}/astscript-psf-scale-factor
+%attr(755,root,root) %{_bindir}/astscript-psf-select-stars
+%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-radial-profile
%attr(755,root,root) %{_bindir}/astscript-sort-by-night
%attr(755,root,root) %{_bindir}/astsegment
@@ -162,10 +171,11 @@ rm -rf $RPM_BUILD_ROOT
%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.14
+%attr(755,root,root) %ghost %{_libdir}/libgnuastro.so.15
%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
+%{_datadir}/gnuastro
%{_infodir}/gnuastro.info*
%{_infodir}/gnuastro-figures
%{_mandir}/man1/astarithmetic.1*
@@ -182,6 +192,12 @@ rm -rf $RPM_BUILD_ROOT
%{_mandir}/man1/astnoisechisel.1*
%{_mandir}/man1/astquery.1*
%{_mandir}/man1/astscript-ds9-region.1*
+%{_mandir}/man1/astscript-fits-view.1*
+%{_mandir}/man1/astscript-psf-scale-factor.1*
+%{_mandir}/man1/astscript-psf-select-stars.1*
+%{_mandir}/man1/astscript-psf-stamp.1*
+%{_mandir}/man1/astscript-psf-subtract.1*
+%{_mandir}/man1/astscript-psf-unite.1*
%{_mandir}/man1/astscript-radial-profile.1*
%{_mandir}/man1/astscript-sort-by-night.1*
%{_mandir}/man1/astsegment.1*
diff --git a/gnuastro-info.patch b/gnuastro-info.patch
index 85f60c9..53ae86f 100644
--- a/gnuastro-info.patch
+++ b/gnuastro-info.patch
@@ -1,117 +1,129 @@
---- gnuastro-0.16/doc/gnuastro.texi.orig 2021-10-10 20:14:28.000000000 +0200
-+++ gnuastro-0.16/doc/gnuastro.texi 2021-11-27 16:39:48.326765277 +0100
-@@ -47,68 +47,68 @@ A copy of the license is included in the
+--- gnuastro-0.17/doc/gnuastro.texi.orig 2022-03-20 01:46:15.000000000 +0100
++++ gnuastro-0.17/doc/gnuastro.texi 2022-06-26 20:52:42.165083932 +0200
+@@ -50,74 +50,74 @@ A copy of the license is included in the
@c To include in the info directory.
@dircategory Astronomy
@direntry
-* Gnuastro: (gnuastro). GNU Astronomy Utilities.
-* libgnuastro: (gnuastro)Gnuastro library. Full Gnuastro library doc.
-+* Gnuastro: (gnuastro). GNU Astronomy Utilities
++* Gnuastro: (gnuastro). GNU Astronomy Utilities
+* 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
-* bug-gnuastro: (gnuastro)Report a bug. How to report bugs
-+* bug-gnuastro: (gnuastro)Report a bug. How to report bugs
++* bug-gnuastro: (gnuastro)Report a bug. How to report bugs
-* Arithmetic: (gnuastro)Arithmetic. Arithmetic operations on pixels.
-* astarithmetic: (gnuastro)Invoking astarithmetic. Options to Arithmetic.
-+* Arithmetic: (gnuastro)Arithmetic. Arithmetic operations on pixels
++* Arithmetic: (gnuastro)Arithmetic. Arithmetic operations on pixels
+* astarithmetic: (gnuastro)astarithmetic. Options to Arithmetic
-* BuildProgram: (gnuastro)BuildProgram. Compile and run programs using Gnuastro's library.
-* astbuildprog: (gnuastro)Invoking astbuildprog. Options to BuildProgram.
-+* BuildProgram: (gnuastro)BuildProgram. Compile and run programs using Gnuastro's library
-+* astbuildprog: (gnuastro)astbuildprog. Options to BuildProgram
++* BuildProgram: (gnuastro)BuildProgram. Compile and run programs using Gnuastro's library
++* astbuildprog: (gnuastro)astbuildprog. Options to BuildProgram
-* 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.
-+* Convolve: (gnuastro)Convolve. Convolve an input file with kernel
-+* astconvolve: (gnuastro)astconvolve. Options to Convolve
++* Convolve: (gnuastro)Convolve. Convolve an input file with kernel
++* astconvolve: (gnuastro)astconvolve. Options to Convolve
-* CosmicCalculator: (gnuastro)CosmicCalculator. For cosmological params.
-* astcosmiccal: (gnuastro)Invoking astcosmiccal. Options to CosmicCalculator.
+* CosmicCalculator: (gnuastro)CosmicCalculator. For cosmological params
-+* astcosmiccal: (gnuastro)astcosmiccal. Options to CosmicCalculator
++* astcosmiccal: (gnuastro)astcosmiccal. Options to CosmicCalculator
-* 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)
-+* astcrop: (gnuastro)astcrop. Options to Crop
++* 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.
-* astfits: (gnuastro)Invoking astfits. Options to Fits.
-+* Fits: (gnuastro)Fits. View and manipulate FITS extensions and keywords
-+* astfits: (gnuastro)astfits. Options to Fits
++* Fits: (gnuastro)Fits. View and manipulate FITS extensions and keywords
++* astfits: (gnuastro)astfits. Options to Fits
-* MakeCatalog: (gnuastro)MakeCatalog. Make a catalog from labeled image.
-* astmkcatalog: (gnuastro)Invoking astmkcatalog. Options to MakeCatalog.
-+* MakeCatalog: (gnuastro)MakeCatalog. Make a catalog from labeled image
-+* astmkcatalog: (gnuastro)astmkcatalog. Options to MakeCatalog
++* 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
++* 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
-+* astmkprof: (gnuastro)astmkprof. Options to MakeProfiles
++* MakeProfiles: (gnuastro)MakeProfiles. Make mock profiles
++* astmkprof: (gnuastro)astmkprof. Options to MakeProfiles
-* Match: (gnuastro)Match. Match two separate catalogs.
-* astmatch: (gnuastro)Invoking astmatch. Options to Match.
-+* Match: (gnuastro)Match. Match two separate catalogs
-+* astmatch: (gnuastro)astmatch. Options to Match
++* Match: (gnuastro)Match. Match two separate catalogs
++* astmatch: (gnuastro)astmatch. Options to Match
-* NoiseChisel: (gnuastro)NoiseChisel. Detect signal in noise.
-* astnoisechisel: (gnuastro)Invoking astnoisechisel. Options to NoiseChisel.
-+* NoiseChisel: (gnuastro)NoiseChisel. Detect signal in noise
++* NoiseChisel: (gnuastro)NoiseChisel. Detect signal in noise
+* astnoisechisel: (gnuastro)astnoisechisel. Options to NoiseChisel
-* Segment: (gnuastro)Segment. Segment detections based on signal structure.
-* astsegment: (gnuastro)Invoking astsegment. Options to Segment.
-+* Segment: (gnuastro)Segment. Segment detections based on signal structure
-+* astsegment: (gnuastro)astsegment. Options to Segment
++* Segment: (gnuastro)Segment. Segment detections based on signal structure
++* astsegment: (gnuastro)astsegment. Options to Segment
-* Query: (gnuastro)Query. Access remote databases for downloading data.
-* astquery: (gnuastro)Invoking astquery. Options to Query.
-+* Query: (gnuastro)Query. Access remote databases for downloading data
-+* astquery: (gnuastro)astquery. Options to Query
++* Query: (gnuastro)Query. Access remote databases for downloading data
++* astquery: (gnuastro)astquery. Options to Query
-* Statistics: (gnuastro)Statistics. Get image Statistics.
-* aststatistics: (gnuastro)Invoking aststatistics. Options to Statistics.
-+* Statistics: (gnuastro)Statistics. Get image Statistics
++* Statistics: (gnuastro)Statistics. Get image Statistics
+* aststatistics: (gnuastro)aststatistics. Options to Statistics
-* Table: (gnuastro)Table. Read and write FITS binary or ASCII tables.
-* asttable: (gnuastro)Invoking asttable. Options to Table.
-+* Table: (gnuastro)Table. Read and write FITS binary or ASCII tables
-+* asttable: (gnuastro)asttable. Options to Table
++* Table: (gnuastro)Table. Read and write FITS binary or ASCII tables
++* asttable: (gnuastro)asttable. Options to Table
-* Warp: (gnuastro)Warp. Warp a dataset to a new grid.
-* astwarp: (gnuastro)Invoking astwarp. Options to Warp.
-+* Warp: (gnuastro)Warp. Warp a dataset to a new grid
-+* astwarp: (gnuastro)astwarp. Options to Warp
++* Warp: (gnuastro)Warp. Warp a dataset to a new grid
++* 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: (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)Invoking astscript-fits-view. 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: (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-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
@end direntry
-@@ -392,7 +392,7 @@ Data containers
+
+@@ -412,7 +412,7 @@ Data containers
Fits
@@ -120,7 +132,7 @@
Invoking Fits
-@@ -405,7 +405,7 @@ ConvertType
+@@ -425,7 +425,7 @@ ConvertType
* Color:: Some explanations on color.
* Aligning images with small WCS offsets:: When the WCS slightly differs.
* Annotations for figure in paper:: Adding coordinates or physical scale.
@@ -129,7 +141,7 @@
Annotations for figure in paper
-@@ -415,12 +415,12 @@ Table
+@@ -435,12 +435,12 @@ Table
* Column arithmetic:: How to do operations on table columns.
* Operation precedence in Table:: Order of running options in Table.
@@ -144,7 +156,7 @@
Data manipulation
-@@ -434,7 +434,7 @@ Crop
+@@ -454,7 +454,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.
@@ -153,16 +165,16 @@
Invoking Crop
-@@ -446,7 +446,7 @@ Arithmetic
-
+@@ -467,7 +467,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
Arithmetic operators
-@@ -473,7 +473,7 @@ Convolve
+@@ -494,7 +494,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.
@@ -171,7 +183,7 @@
Spatial domain convolution
-@@ -498,7 +498,7 @@ Warp
+@@ -519,7 +519,7 @@ Warp
* Warping basics:: Basics of coordinate transformation.
* Merging multiple warpings:: How to merge multiple matrices.
* Resampling:: Warping an image is re-sampling it.
@@ -180,7 +192,7 @@
Data analysis
-@@ -514,7 +514,7 @@ Statistics
+@@ -535,7 +535,7 @@ Statistics
* 2D Histograms:: Plotting the distribution of two variables.
* Sigma clipping:: Definition of @mymath{\sigma}-clipping.
* Sky value:: Definition and derivation of the Sky value.
@@ -189,7 +201,7 @@
2D Histograms
-@@ -530,7 +530,7 @@ Sky value
+@@ -551,7 +551,7 @@ Sky value
NoiseChisel
* NoiseChisel changes after publication:: Updates since published papers.
@@ -198,7 +210,7 @@
Invoking NoiseChisel
-@@ -540,7 +540,7 @@ Invoking NoiseChisel
+@@ -561,7 +561,7 @@ Invoking NoiseChisel
Segment
@@ -207,7 +219,7 @@
Invoking Segment
-@@ -555,7 +555,7 @@ MakeCatalog
+@@ -576,7 +576,7 @@ MakeCatalog
* Quantifying measurement limits:: For comparing different catalogs.
* Measuring elliptical parameters:: Estimating elliptical parameters.
* Adding new columns to MakeCatalog:: How to add new columns.
@@ -216,16 +228,16 @@
Quantifying measurement limits
-@@ -575,7 +575,7 @@ Invoking MakeCatalog
-
+@@ -598,7 +598,7 @@ Invoking MakeCatalog
Match
+ * Matching algorithms:: Different ways to find the match
-* Invoking astmatch:: Inputs, outputs and options of Match
+* astmatch:: Inputs, outputs and options of Match
Modeling and fitting
-@@ -587,7 +587,7 @@ MakeProfiles
+@@ -610,7 +610,7 @@ MakeProfiles
* Modeling basics:: Astronomical modeling basics.
* If convolving afterwards:: Considerations for convolving later.
* Profile magnitude:: Definition of total profile magnitude.
@@ -234,7 +246,7 @@
Modeling basics
-@@ -608,7 +608,7 @@ Invoking MakeProfiles
+@@ -631,7 +631,7 @@ Invoking MakeProfiles
MakeNoise
* Noise basics:: Noise concepts and definitions.
@@ -243,7 +255,7 @@
Noise basics
-@@ -625,7 +625,7 @@ CosmicCalculator
+@@ -648,7 +648,7 @@ CosmicCalculator
* Distance on a 2D curved space:: Distances in 2D for simplicity
* Extending distance concepts to 3D:: Going to 3D (our real universe).
@@ -252,26 +264,45 @@
Invoking CosmicCalculator
-@@ -641,15 +641,15 @@ Installed scripts
+@@ -666,28 +666,28 @@ Installed scripts
Sort FITS files by night
-* Invoking astscript-sort-by-night:: Inputs and outputs to this script.
-+* astscript-sort-by-night:: Inputs and outputs to this script.
++* astscript-sort-by-night:: Inputs and outputs to this script.
Generate radial profile
-* Invoking astscript-radial-profile:: How to call astscript-radial-profile
-+* astscript-radial-profile:: How to call astscript-radial-profile
++* astscript-radial-profile:: How to call astscript-radial-profile
SAO DS9 region files from table
-* Invoking astscript-ds9-region:: How to call astscript-ds9-region
-+* astscript-ds9-region:: How to call astscript-ds9-region
++* astscript-ds9-region:: How to call astscript-ds9-region
+
+ Viewing FITS file contents with DS9 or TOPCAT
+
+-* Invoking astscript-fits-view:: How to call this script
++* astscript-fits-view:: How to call this script
+
+ PSF construction and subtraction
+
+ * Overview of the PSF scripts:: Summary of concepts and methods
+-* Invoking astscript-psf-select-stars:: Select good starts within an image.
+-* Invoking astscript-psf-stamp:: Make a stamp of each star to stack.
+-* 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.
Library
-@@ -666,7 +666,7 @@ Review of library fundamentals
+@@ -704,7 +704,7 @@ Review of library fundamentals
BuildProgram
@@ -280,7 +311,7 @@
Gnuastro library
-@@ -884,7 +884,7 @@ $ echo "source /usr/local/share/gnuastro
+@@ -919,7 +919,7 @@ $ echo "source /usr/local/share/gnuastro
The last command is to enable Gnuastro's custom TAB completion in Bash.
For more on this useful feature, see @ref{Shell TAB completion}).
@@ -289,7 +320,7 @@
Some complete Tutorials are also available in this book 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).
-@@ -1639,7 +1639,7 @@ Because bright/large galaxies and stars@
+@@ -1695,7 +1695,7 @@ The tutorial concludes with subtracting
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}.
@@ -298,7 +329,7 @@
For an explanation of the conventions we use in the example codes through the book, please see @ref{Conventions}.
@menu
-@@ -2120,8 +2120,8 @@ To return to the previous page, press @k
+@@ -2179,8 +2179,8 @@ 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>}.
You don't need to start from the top of the manual every time.
@@ -309,7 +340,7 @@
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.
-@@ -2819,7 +2819,7 @@ As you see in the matrix that is printed
+@@ -2882,7 +2882,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.
@@ -318,7 +349,7 @@
Fortunately these datasets are already aligned to the same pixel grid, so you don't 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.
-@@ -3762,7 +3762,7 @@ $ cat apertures.txt
+@@ -3832,7 +3832,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).
@@ -327,7 +358,16 @@
@example
$ astmkprof apertures.txt --background=flat-ir/xdf-f160w.fits \
-@@ -6929,7 +6929,7 @@ Afterwards, in @ref{Common options}, we'
+@@ -6202,7 +6202,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.
+-To do that, we will use the script @file{astscript-psf-unite} with the command below (for more on this script, see @ref{Invoking astscript-psf-unite}).
++To do that, we will use the script @file{astscript-psf-unite} with the command below (for more on this script, see @ref{astscript-psf-unite}).
+ 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):
+@@ -8406,7 +8406,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.
@@ -336,7 +376,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.
-@@ -6955,7 +6955,7 @@ Hence, when there multiple arguments, th
+@@ -8432,7 +8432,7 @@ Hence, when there multiple arguments, th
The outputs of @option{--usage} and @option{--help} shows which arguments are optional and which are mandatory, see @ref{--usage}.
As their name suggests, @emph{options} can be considered to be optional and most of the time, you don't have to worry about what order you specify them in.
@@ -345,7 +385,7 @@
@cindex Metacharacters on the command-line In case your arguments or option values contain any of the shell's meta-characters, you have to quote them.
If there is only one such character, you can use a backslash (@command{\}) before it.
-@@ -6985,8 +6985,8 @@ $ astcrop --hdu="3; images(exposure > 10
+@@ -8462,8 +8462,8 @@ $ astcrop --hdu="3; images(exposure > 10
@node Arguments, Options, Arguments and options, Arguments and options
@subsubsection Arguments
In Gnuastro, arguments are almost exclusively used as the input data file names.
@@ -356,7 +396,7 @@
@cindex Filename suffix
@cindex Suffix (filename)
-@@ -7043,7 +7043,7 @@ Command-line options allow configuring t
+@@ -8520,7 +8520,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}).
@@ -365,7 +405,7 @@
Both formats are shown for those which support both.
First the short is shown then the long.
-@@ -7117,7 +7117,7 @@ If you are satisfied with the change, yo
+@@ -8594,7 +8594,7 @@ If you are satisfied with the change, yo
If the change wasn't 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.
@@ -374,7 +414,7 @@
In these cases, the order of stored values is the same order that you specified on the command-line.
@cindex Configuration files
-@@ -7588,7 +7588,7 @@ As another example, if an option needs i
+@@ -9065,7 +9065,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.
@@ -383,7 +423,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}).
-@@ -7708,7 +7708,7 @@ All the values must either be stored in
+@@ -9185,7 +9185,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 don't, the program will not abort.
@@ -392,7 +432,7 @@
-@@ -8083,7 +8083,7 @@ $ info astprogramname
+@@ -9560,7 +9560,7 @@ $ info astprogramname
@end example
@noindent
@@ -401,7 +441,7 @@
Finally, if you run Info with the official program name, for example Crop or NoiseChisel:
@example
-@@ -8514,7 +8514,7 @@ Each column in the table contains the va
+@@ -9991,7 +9991,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.
@@ -410,7 +450,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.
-@@ -9084,10 +9084,10 @@ For example you can copy or cut (copy an
+@@ -10561,10 +10561,10 @@ For example you can copy or cut (copy an
It also has features to delete, add, or edit meta-data keywords within one HDU.
@menu
@@ -423,7 +463,7 @@
@subsection Invoking Fits
Fits can print or manipulate the FITS file HDUs (extensions), meta-data keywords in a given HDU.
-@@ -9176,7 +9176,7 @@ These two classes of options cannot be c
+@@ -10653,7 +10653,7 @@ These two classes of options cannot be c
@@ -432,7 +472,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}.
-@@ -9194,7 +9194,7 @@ Note that this option must be called alo
+@@ -10671,7 +10671,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, don't use any of the options in this section or in @ref{Keyword inspection and manipulation}.
@@ -441,7 +481,7 @@
@item --hastablehdu
Print @code{1} (on standard output) if at least one table HDU (ASCII or binary) exists in the FITS file.
-@@ -9330,7 +9330,7 @@ If we hadn't used @option{--primaryimghd
+@@ -10807,7 +10807,7 @@ If we hadn't used @option{--primaryimghd
@end table
@@ -450,7 +490,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.
-@@ -9363,7 +9363,7 @@ $ astfits image-a.fits --keyvalue=NAXIS,
+@@ -10840,7 +10840,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).
@@ -459,7 +499,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 doesn't 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}.
-@@ -9824,14 +9824,14 @@ The various types will increase with fut
+@@ -11309,14 +11309,14 @@ The various types will increase with fut
The conversion is not only one way (from FITS to other formats), but two ways (except the EPS and PDF formats at footnote{Because EPS and PDF are vector, not raster/pixelated formats}).
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.
@@ -476,7 +516,7 @@
@end menu
@node Recognized file formats, Color, ConvertType, ConvertType
-@@ -9945,7 +9945,7 @@ Most programs also support input as plai
+@@ -11430,7 +11430,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}.
@@ -485,7 +525,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.
-@@ -9990,7 +9990,7 @@ To make measurements in different filter
+@@ -11475,7 +11475,7 @@ To make measurements in different filter
Therefore, the FITS format that is used to store astronomical datasets is inherently a mono-channel format (see @ref{Recognized file formats} or @ref{Fits}).
When a subject has been imaged in multiple filters, you can feed each different filter into the red, green and blue channels and obtain a colored visualization.
@@ -494,7 +534,7 @@
@cindex Grayscale
@cindex Visualization
-@@ -10004,7 +10004,7 @@ The most basic is to use shades of black
+@@ -11489,7 +11489,7 @@ The most basic is to use shades of black
This scheme is called grayscale.
To help in visualization, more complex mappings can be defined.
For example, the values can be scaled to a range of 0 to 360 and used as the ``Hue'' term of the @url{https://en.wikipedia.org/wiki/HSL_and_HSV, Hue-Saturation-Value} (HSV) color space (while fixing the ``Saturation'' and ``Value'' terms).
@@ -503,7 +543,7 @@
Since grayscale is a commonly used mapping of single-valued datasets, we'll continue with a closer look at how it is stored.
One way to represent a gray-scale image in different color spaces is to use the same proportions of the primary colors in each pixel.
-@@ -10094,7 +10094,7 @@ $ ./my-align.sh
+@@ -11579,7 +11579,7 @@ $ ./my-align.sh
@noindent
Of course, feel free to hack it and modify it to fit your datasets, like the rest of Gnuastro, this script is released under GNU GPLv.3 and above, see @ref{Your rights}.
@@ -512,7 +552,7 @@
@subsection Annotations for figure in paper
@cindex Image annotation
-@@ -10503,7 +10503,7 @@ pdflatex -shell-escape -halt-on-error re
+@@ -11988,7 +11988,7 @@ pdflatex -shell-escape -halt-on-error re
@@ -521,7 +561,7 @@
@subsection Invoking ConvertType
ConvertType will convert any recognized input file type to any specified output type.
-@@ -10791,18 +10791,18 @@ Table has a large set of operations that
+@@ -12276,18 +12276,18 @@ Table has a large set of operations that
For operations that Table doesn't 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.
@@ -543,7 +583,7 @@
@end menu
@node Column arithmetic, Operation precedence in Table, Table, Table
-@@ -10865,7 +10865,7 @@ Column arithmetic changes the values of
+@@ -12350,7 +12350,7 @@ Column arithmetic changes the values of
So the old column meta data can't 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.
@@ -552,8 +592,8 @@
Finally, since the arithmetic expressions are a value to @option{--column}, it doesn't 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!
-@@ -11001,13 +11001,13 @@ If you don't need to see the Unix-second
- For more on @option{--keyvalue}, see @ref{Keyword inspection and manipulation}.
+@@ -12497,13 +12497,13 @@ Other than the units of the output, this
+ See the description of that operator for an example.
@end table
- at node Operation precedence in Table, Invoking asttable, Column arithmetic, Table
@@ -568,7 +608,7 @@
@table @asis
@item Column information (@option{--information} or @option{-i})
-@@ -11131,7 +11131,7 @@ asttable table.fits -cRA,DEC --noblank=M
+@@ -12648,7 +12648,7 @@ asttable table.fits -cRA,DEC --noblanken
@end example
@end cartouche
@@ -577,7 +617,16 @@
@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}).
-@@ -11333,7 +11333,7 @@ The chosen column doesn't have to be in
+@@ -12846,7 +12846,7 @@ $ asttable a.fits --catrowfile=b.fits --
+ @strong{How to avoid repetition when adding rows:} this option will simply add the rows of multiple tables into one, it doesn't 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.
+-For more on using Match for this scenario, see the description of @option{--outcols} in @ref{Invoking astmatch}.
++For more on using Match for this scenario, see the description of @option{--outcols} in @ref{astmatch}.
+ @end cartouche
+
+ @item -X STR
+@@ -12878,7 +12878,7 @@ The chosen column doesn't have to be in
This is good when you just want to select using one column's values, but don't need that column anymore afterwards.
For one example of using this option, see the example under
@@ -586,7 +635,7 @@
@item --inpolygon=STR1,STR2
Only return rows where the given coordinates are inside the polygon specified by the @option{--polygon} option.
-@@ -11570,10 +11570,10 @@ We will try to add high-level interfaces
+@@ -13160,10 +13160,10 @@ We will try to add high-level interfaces
@menu
* Available databases:: List of available databases to Query.
@@ -599,7 +648,7 @@
@subsection Available databases
The current list of databases supported by Query are listed at the end of this section.
-@@ -11794,7 +11794,7 @@ For details on each dataset with necessa
+@@ -13384,7 +13384,7 @@ For details on each dataset with necessa
@@ -608,7 +657,7 @@
@subsection Invoking Query
Query provides a high-level interface to downloading subsets of data from databases.
-@@ -11896,7 +11896,7 @@ When @option{--dataset} is specified, th
+@@ -13487,7 +13487,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.
@@ -617,7 +666,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.
-@@ -12064,14 +12064,14 @@ So when you simply crop the image of suc
+@@ -13657,14 +13657,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.
@@ -634,7 +683,7 @@
@end menu
@node Crop modes, Crop section syntax, Crop, Crop
-@@ -12093,13 +12093,13 @@ Irrespective of how the crop region is d
+@@ -13686,13 +13686,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.
@@ -650,7 +699,7 @@
@table @asis
@item Image coordinates
-@@ -12126,7 +12126,7 @@ The former is lower-level (doesn't accep
+@@ -13719,7 +13719,7 @@ The former is lower-level (doesn't accep
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.
@@ -659,7 +708,7 @@
@end table
@item WCS coordinates
-@@ -12156,7 +12156,7 @@ If it exists in the input images, it wil
+@@ -13751,7 +13751,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).
@@ -668,7 +717,7 @@
@end table
@cartouche
-@@ -12210,7 +12210,7 @@ If you forget the quotes, anything after
+@@ -13805,7 +13805,7 @@ If you forget the quotes, anything after
See @ref{Command-line} for a description of how the command-line works.
@@ -677,7 +726,7 @@
@subsection Blank pixels
@cindex Blank pixel
-@@ -12235,7 +12235,7 @@ So by default, when dealing with float o
+@@ -13830,7 +13830,7 @@ So by default, when dealing with float o
This can be turned off with the @option{--zeroisnotblank} option.
@@ -686,7 +735,7 @@
@subsection Invoking Crop
Crop will crop a region from an image.
-@@ -12300,7 +12300,7 @@ This has no effect on each output, see @
+@@ -13899,7 +13899,7 @@ This has no effect on each output, see @
* Crop known issues:: Known issues in running Crop.
@end menu
@@ -695,7 +744,7 @@
@subsubsection Crop options
The options can be classified into the following contexts: Input, Output and operating mode options.
-@@ -12556,7 +12556,7 @@ The value must either be @option{img} or
+@@ -14165,7 +14165,7 @@ The value must either be @option{img} or
@@ -704,7 +753,7 @@
@subsubsection Crop output
The string given to @option{--output} option will be interpreted depending
-@@ -12598,10 +12598,10 @@ The cropped image file name for that row
+@@ -14207,10 +14207,10 @@ 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 aren't.
@@ -717,7 +766,7 @@
@subsubsection Crop known issues
When running Crop, you may encounter strange errors and bugs.
-@@ -12651,13 +12651,13 @@ For example, in the reduction of raw dat
+@@ -14260,14 +14260,14 @@ 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}.
@@ -727,13 +776,14 @@
@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
@end menu
- @node Reverse polish notation, Arithmetic operators, Arithmetic, Arithmetic
-@@ -12735,7 +12735,7 @@ There are no more operands or operators,
+ @node Reverse polish notation, Integer benefits and pitfalls, Arithmetic, Arithmetic
+@@ -14345,7 +14345,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
@@ -742,16 +792,16 @@
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.
-@@ -12746,7 +12746,7 @@ This is a very powerful notation and is
-
+@@ -14443,7 +14443,7 @@ It is possible to do internal checks in
+ However, we haven't 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}).
-
-- at node Arithmetic operators, Invoking astarithmetic, Reverse polish notation, Arithmetic
-+ at node Arithmetic operators, astarithmetic, Reverse polish notation, Arithmetic
+- at node Arithmetic operators, Invoking astarithmetic, Integer benefits and pitfalls, Arithmetic
++ at node Arithmetic operators, astarithmetic, Integer benefits and pitfalls, 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.
-@@ -14110,7 +14110,7 @@ Similar to the @code{tofile} operator, w
+@@ -15833,7 +15833,7 @@ Similar to the @code{tofile} operator, w
@end table
@@ -760,7 +810,7 @@
@subsection Invoking Arithmetic
Arithmetic will do pixel to pixel arithmetic operations on the individual pixels of input data and/or numbers.
-@@ -14311,7 +14311,7 @@ $ echo "" | awk '@{print (10.32-3.84)^2.
+@@ -16034,7 +16034,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.
@@ -769,7 +819,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}.
-@@ -14342,7 +14342,7 @@ However this text is written for an unde
+@@ -16065,7 +16065,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.
@@ -778,7 +828,7 @@
@end menu
@node Spatial domain convolution, Frequency domain and Fourier operations, Convolve, Convolve
-@@ -15086,7 +15086,7 @@ The reason is that if you apply a freque
+@@ -16809,7 +16809,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}.
@@ -787,7 +837,7 @@
@subsection Convolution kernel
All the programs that need convolution will need to be given a convolution kernel file and extension.
-@@ -15104,7 +15104,7 @@ By default MakeProfiles will make the Ga
+@@ -16827,7 +16827,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 don't have to worry about normalizing the kernel.
@@ -796,7 +846,7 @@
@end itemize
-@@ -15144,7 +15144,7 @@ By default, the system configuration fil
+@@ -16867,7 +16867,7 @@ By default, the system configuration fil
@end table
@@ -805,7 +855,7 @@
@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.
-@@ -15354,7 +15354,7 @@ It is therefore necessary to warp the im
+@@ -17077,7 +17077,7 @@ It is therefore necessary to warp the im
* Warping basics:: Basics of coordinate transformation.
* Merging multiple warpings:: How to merge multiple matrices.
* Resampling:: Warping an image is re-sampling it.
@@ -814,7 +864,7 @@
@end menu
@node Warping basics, Merging multiple warpings, Warp, Warp
-@@ -15497,7 +15497,7 @@ These three operations can be merged in
+@@ -17220,7 +17220,7 @@ These three operations can be merged in
@@ -823,7 +873,7 @@
@subsection Resampling
@cindex Pixel
-@@ -15559,7 +15559,7 @@ So if the input image has fringes, they
+@@ -17282,7 +17282,7 @@ So if the input image has fringes, they
Because of the PSF no astronomical target has a sharpchange in the signal so this issue is less important for astronomical applications, see @ref{PSF}.
@@ -832,7 +882,7 @@
@subsection Invoking Warp
Warp an input dataset into a new grid.
-@@ -15699,10 +15699,10 @@ This option is applied after the final w
+@@ -17422,10 +17422,10 @@ This option is applied after the final w
See the explanation above for coordinates in the FITS standard to better understand this option and when it should be used.
@item --hstartwcs=INT
@@ -845,7 +895,7 @@
@item -k
@itemx --keepwcs
-@@ -15772,7 +15772,7 @@ The Statistics program is designed for s
+@@ -17495,7 +17495,7 @@ The Statistics program is designed for s
* 2D Histograms:: Plotting the distribution of two variables.
* Sigma clipping:: Definition of @mymath{\sigma}-clipping.
* Sky value:: Definition and derivation of the Sky value.
@@ -854,7 +904,7 @@
@end menu
-@@ -15786,7 +15786,7 @@ So on the horizontal axis we have the bi
+@@ -17509,7 +17509,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).},
@@ -863,7 +913,7 @@
@cindex Intervals, histogram
@cindex Bin width, histogram
-@@ -15807,7 +15807,7 @@ Normalizing a cumulative frequency plot
+@@ -17530,7 +17530,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.
@@ -872,7 +922,7 @@
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).
-@@ -16117,7 +16117,7 @@ The ASCII histogram that is printed on t
+@@ -17840,7 +17840,7 @@ The ASCII histogram that is printed on t
@@ -881,7 +931,7 @@
@subsection Sky value
@cindex Sky
-@@ -16374,7 +16374,7 @@ The same goes for the @option{--checksky
+@@ -18098,7 +18098,7 @@ The same goes for the @option{--checksky
@@ -890,7 +940,7 @@
@subsection Invoking Statistics
Statistics will print statistical measures of an input dataset (table column or image).
-@@ -17007,7 +17007,7 @@ You can then directly feed NoiseChisel's
+@@ -18731,7 +18731,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}.
@@ -899,7 +949,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}).
-@@ -17019,17 +17019,17 @@ In the meantime, they also show the modu
+@@ -18743,17 +18743,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.
@@ -921,7 +971,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]}.
-@@ -17057,7 +17057,7 @@ In our tests, this gave a much improved
+@@ -18781,7 +18781,7 @@ In our tests, this gave a much improved
@@ -930,7 +980,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.
-@@ -17184,7 +17184,7 @@ Finally, in @ref{NoiseChisel output} the
+@@ -18908,7 +18908,7 @@ Finally, in @ref{NoiseChisel output} the
* NoiseChisel output:: NoiseChisel's output options and format.
@end menu
@@ -939,7 +989,7 @@
@subsubsection NoiseChisel input
The options here can be used to configure the inputs and output of NoiseChisel, along with some general processing options.
-@@ -17288,7 +17288,7 @@ Except for the tile size, all the other
+@@ -19012,7 +19012,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
@@ -948,7 +998,7 @@
@subsubsection Detection options
Detection is the process of separating the pixels in the image into two groups: 1) Signal, and 2) Noise.
-@@ -17566,7 +17566,7 @@ By default the output will have the same
+@@ -19355,7 +19355,7 @@ By default the output will have the same
@@ -957,7 +1007,7 @@
@subsubsection NoiseChisel output
NoiseChisel's output is a multi-extension FITS file.
-@@ -17737,14 +17737,14 @@ Those papers cannot be updated any more
+@@ -19526,14 +19526,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}.
@@ -975,7 +1025,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]}.
-@@ -17756,7 +17756,7 @@ Finally, in @ref{Invoking astsegment}, w
+@@ -19545,7 +19545,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}.}.
@@ -984,7 +1034,7 @@
@subsection Invoking Segment
Segment will identify substructure within the detected regions of an input image.
-@@ -17809,7 +17809,7 @@ Finally, in @ref{Segment output}, we'll
+@@ -19598,7 +19598,7 @@ Finally, in @ref{Segment output}, we'll
* Segment output:: Outputs of Segment
@end menu
@@ -993,7 +1043,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.
-@@ -17951,7 +17951,7 @@ It is important for them to be much larg
+@@ -19740,7 +19740,7 @@ It is important for them to be much larg
@end table
@@ -1002,7 +1052,7 @@
@subsubsection Segmentation options
The options below can be used to configure every step of the segmentation process in the Segment program.
-@@ -18065,7 +18065,7 @@ This behavior can be disabled with @opti
+@@ -19854,7 +19854,7 @@ This behavior can be disabled with @opti
@end table
@@ -1011,7 +1061,7 @@
@subsubsection Segment output
The main output of Segment are two label datasets (with integer types, separating the dataset's elements into different classes).
-@@ -18213,7 +18213,7 @@ Similarly, the sum of all these pixels w
+@@ -19998,7 +19998,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).
@@ -1020,7 +1070,7 @@
A very important factor in any measurement is understanding its validity range, or limits.
Therefore in @ref{Quantifying measurement limits}, we'll 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}.
-@@ -18225,7 +18225,7 @@ For those who feel MakeCatalog's existin
+@@ -20010,7 +20010,7 @@ For those who feel MakeCatalog's existin
* Quantifying measurement limits:: For comparing different catalogs.
* Measuring elliptical parameters:: Estimating elliptical parameters.
* Adding new columns to MakeCatalog:: How to add new columns.
@@ -1029,7 +1079,7 @@
@end menu
@node Detection and catalog production, Brightness flux magnitude, MakeCatalog, MakeCatalog
-@@ -18770,7 +18770,7 @@ When the object cannot be represented as
+@@ -20573,7 +20573,7 @@ When the object cannot be represented as
@@ -1038,7 +1088,7 @@
@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]}).
-@@ -18786,7 +18786,7 @@ These two passes are designed to be raw
+@@ -20589,7 +20589,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).
@@ -1047,7 +1097,7 @@
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 doesn't necessarily have to be in the end), and near conceptually similar options.
-@@ -18838,7 +18838,7 @@ Update this manual and add a description
+@@ -20641,7 +20641,7 @@ Update this manual and add a description
@@ -1056,7 +1106,7 @@
@subsection Invoking MakeCatalog
MakeCatalog will do measurements and produce a catalog from a labeled dataset and optional values dataset(s).
-@@ -18892,7 +18892,7 @@ Finally, in @ref{MakeCatalog output} the
+@@ -20695,7 +20695,7 @@ Finally, in @ref{MakeCatalog output} the
* MakeCatalog output:: File names of MakeCatalog's output table.
@end menu
@@ -1065,7 +1115,7 @@
@subsubsection MakeCatalog inputs and basic settings
MakeCatalog works by using a localized/labeled dataset (see @ref{MakeCatalog}).
-@@ -18946,7 +18946,7 @@ In such cases, you can use Gnuastro's Ar
+@@ -20749,7 +20749,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 don't have to worry about any of the @option{--*file} options below.
@@ -1074,7 +1124,7 @@
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.
-@@ -19045,7 +19045,7 @@ But in practice real images my be over-s
+@@ -20848,7 +20848,7 @@ But in practice real images my be over-s
@@ -1083,7 +1133,7 @@
@subsubsection Upper-limit settings
The upper-limit magnitude was discussed in @ref{Quantifying measurement limits}.
-@@ -19132,7 +19132,7 @@ The total number of rows is thus unknown
+@@ -20935,7 +20935,7 @@ The total number of rows is thus unknown
@end table
@@ -1092,7 +1142,7 @@
@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.
-@@ -19689,7 +19689,7 @@ If you use this option, it will have 13
+@@ -21492,7 +21492,7 @@ If you use this option, it will have 13
@@ -1101,20 +1151,29 @@
@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.
-@@ -19813,10 +19813,10 @@ The nearest objects in the two catalogs,
- The aperture can be a circle or an ellipse with any orientation.
+@@ -21617,10 +21617,10 @@ The aperture can be a circle or an ellip
@menu
+ * Matching algorithms:: Different ways to find the match
-* Invoking astmatch:: Inputs, outputs and options of Match
+* astmatch:: Inputs, outputs and options of Match
@end menu
-- at node Invoking astmatch, , Match, Match
-+ at node astmatch, , Match, Match
+- at node Matching algorithms, Invoking astmatch, Match, Match
++ at node Matching algorithms, astmatch, Match, Match
+ @subsection Matching algorithms
+
+ Matching involves two catalogs, let's call them catalog A (with N rows) and catalog B (with M rows).
+@@ -21733,7 +21733,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
+
+- at node Invoking astmatch, , Matching algorithms, Match
++ at node astmatch, , Matching algorithms, Match
@subsection Invoking Match
When given two catalogs, Match finds the rows that are nearest to each other within an input aperture.
-@@ -20154,7 +20154,7 @@ After all the transformations are applie
+@@ -22119,7 +22119,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.
@@ -1123,7 +1182,7 @@
@end menu
-@@ -20480,7 +20480,7 @@ You also have to shift all the pixel pos
+@@ -22449,7 +22449,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:*-n,n:*-n} 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.
@@ -1132,7 +1191,7 @@
-@@ -20491,7 +20491,7 @@ To facilitate this shift, MakeProfiles h
+@@ -22460,7 +22460,7 @@ To facilitate this shift, MakeProfiles h
@@ -1141,7 +1200,7 @@
@subsection Profile magnitude
@cindex Brightness
-@@ -20517,7 +20517,7 @@ If not, only the overlapping pixels of t
+@@ -22486,7 +22486,7 @@ If not, only the overlapping pixels of t
@@ -1150,7 +1209,7 @@
@subsection Invoking MakeProfiles
MakeProfiles will make any number of profiles specified in a catalog either individually or in one image.
-@@ -20602,7 +20602,7 @@ Please see @ref{Sufi simulates a detecti
+@@ -22571,7 +22571,7 @@ Please see @ref{Sufi simulates a detecti
* MakeProfiles log file:: A description of the optional log file.
@end menu
@@ -1159,7 +1218,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}).
-@@ -20730,7 +20730,7 @@ If @option{--tunitinp} is given, this va
+@@ -22699,7 +22699,7 @@ If @option{--tunitinp} is given, this va
@end table
@@ -1168,7 +1227,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}.
-@@ -20897,7 +20897,7 @@ Because the maximum operator is independ
+@@ -22866,7 +22866,7 @@ Because the maximum operator is independ
@end table
@@ -1177,7 +1236,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.
-@@ -21068,7 +21068,7 @@ MakeProfiles won't complain if you use n
+@@ -23039,7 +23039,7 @@ MakeProfiles won't complain if you use n
@end table
@@ -1186,7 +1245,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).
-@@ -21111,16 +21111,16 @@ If an individual image was created, this
+@@ -23082,16 +23082,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.
@@ -1206,7 +1265,7 @@
@subsection Noise basics
@cindex Noise
-@@ -21317,7 +21317,7 @@ In the example above, did you notice how
+@@ -23313,7 +23313,7 @@ In the example above, did you notice how
However, if @option{--envseed} was given, both printed seeds would be the same.
@@ -1215,7 +1274,7 @@
@subsection Invoking MakeNoise
MakeNoise will add noise to an existing image.
-@@ -21434,7 +21434,7 @@ Even higher-level analysis is still need
+@@ -23430,7 +23430,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.
@@ -1224,7 +1283,7 @@
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.
-@@ -21444,7 +21444,7 @@ There are many books thoroughly deriving
+@@ -23440,7 +23440,7 @@ There are many books thoroughly deriving
@menu
* Distance on a 2D curved space:: Distances in 2D for simplicity
* Extending distance concepts to 3D:: Going to 3D (our real universe).
@@ -1233,7 +1292,7 @@
@end menu
@node Distance on a 2D curved space, Extending distance concepts to 3D, CosmicCalculator, CosmicCalculator
-@@ -21592,7 +21592,7 @@ We can thus parameterize the change in d
+@@ -23588,7 +23588,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).}
@@ -1242,7 +1301,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.
-@@ -21641,7 +21641,7 @@ In a non-static universe (with a scale f
+@@ -23637,7 +23637,7 @@ In a non-static universe (with a scale f
@@ -1251,7 +1310,7 @@
@subsection Invoking CosmicCalculator
CosmicCalculator will calculate cosmological variables based on the input parameters.
-@@ -21699,7 +21699,7 @@ CosmicCalculator also has features to he
+@@ -23695,7 +23695,7 @@ CosmicCalculator also has features to he
* CosmicCalculator spectral line calculations:: How they get affected by redshift.
@end menu
@@ -1260,7 +1319,7 @@
@subsubsection CosmicCalculator input options
The inputs to CosmicCalculator can be specified with the following options:
-@@ -21917,9 +21917,9 @@ You can get this list on the command-lin
+@@ -23913,9 +23913,9 @@ You can get this list on the command-lin
@@ -1272,7 +1331,7 @@
The full list of calculations can be useful when you don't 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 don't want to be distracted by all the extra information.
-@@ -22034,7 +22034,7 @@ The comoving volume in Megaparsecs cube
+@@ -24030,7 +24030,7 @@ The comoving volume in Megaparsecs cube
@@ -1281,7 +1340,7 @@
@subsubsection CosmicCalculator spectral line calculations
@cindex Rest frame wavelength
-@@ -22223,10 +22223,10 @@ But we'll leave the implementation of th
+@@ -24221,10 +24221,10 @@ But we'll leave the implementation of th
@end itemize
@menu
@@ -1294,7 +1353,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.
-@@ -22378,10 +22378,10 @@ Gnuastro's @file{astscript-radial-profil
+@@ -24376,10 +24376,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
@@ -1307,7 +1366,7 @@
@subsection Invoking astscript-radial-profile
This installed script will measure the radial profile of an object within an image.
-@@ -22624,10 +22624,10 @@ In this section we describe a simple ins
+@@ -24628,10 +24628,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
@@ -1320,7 +1379,154 @@
@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.
-@@ -23196,10 +23196,10 @@ Please see @ref{Optional dependencies} f
+@@ -24792,13 +24792,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.
+-All of the issues above are the founding purpose of the installed script that is introduced in @ref{Invoking astscript-fits-view}.
++All of the issues above are the founding purpose of the installed script that is introduced in @ref{astscript-fits-view}.
+
+ @menu
+-* Invoking astscript-fits-view:: How to call this script
++* astscript-fits-view:: How to call this script
+ @end menu
+
+- at node Invoking astscript-fits-view, , Viewing FITS file contents with DS9 or TOPCAT, Viewing FITS file contents with DS9 or TOPCAT
++ at node astscript-fits-view, , Viewing FITS file contents with DS9 or TOPCAT, Viewing FITS file contents with DS9 or TOPCAT
+ @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).
+@@ -24924,14 +24924,14 @@ The tutorial uses a real dataset and inc
+
+ @menu
+ * Overview of the PSF scripts:: Summary of concepts and methods
+-* Invoking astscript-psf-select-stars:: Select good starts within an image.
+-* Invoking astscript-psf-stamp:: Make a stamp of each star to stack.
+-* 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.
+ @end menu
+
+- at node Overview of the PSF scripts, Invoking astscript-psf-select-stars, PSF construction and subtraction, PSF construction and subtraction
++ at node Overview of the PSF scripts, astscript-psf-select-stars, PSF construction and subtraction, PSF construction and subtraction
+ @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.
+@@ -24959,11 +24959,11 @@ But you can't blindly use all the stars
+ For example we don't 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}.
++For more on this script, see @ref{astscript-psf-select-stars}.
+
+ Once the catalog of stars is constructed, another script is in charge of making appropriate stamps of the stars.
+ Each stamp is a cropped image of the star with the desired size, normalization of the flux, and mask of the contaminant objects.
+-For more on this script, see @ref{Invoking astscript-psf-stamp}
++For more on this script, see @ref{astscript-psf-stamp}
+ 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.
+@@ -24979,16 +24979,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).
+ 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}.
++For more on the scaling script, see @ref{astscript-psf-scale-factor}.
+
+ Once the flux factor has been computed, a final script is in charge of placing the scaled PSF over the proper location in the image, and subtracting it.
+ It is also possible to only obtain the modeled star by the PSF.
+-For more on the scaling and positioning script, see @ref{Invoking astscript-psf-subtract}.
++For more on the scaling and positioning script, see @ref{astscript-psf-subtract}.
+
+ 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}.
+@@ -24998,7 +24998,7 @@ Just don't forget to run every command,
+
+
+
+- at node Invoking astscript-psf-select-stars, Invoking astscript-psf-stamp, Overview of the PSF scripts, PSF construction and subtraction
++ 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.
+@@ -25154,7 +25154,7 @@ The output name of the final catalog con
+
+
+
+- at node Invoking astscript-psf-stamp, Invoking astscript-psf-unite, Invoking astscript-psf-select-stars, PSF construction and subtraction
++ at node astscript-psf-stamp, astscript-psf-unite, astscript-psf-select-stars, PSF construction and subtraction
+ @subsection Invoking astscript-psf-stamp
+ This installed script will will generate a stamp of fixed size, centered at the provided coordinates (performing sub-pixel regridding if necessary) and normalized at a certain normalization radius.
+ Optionally, it will also mask all the other background sources.
+@@ -25296,7 +25296,7 @@ The main reason of setting this name is
+
+
+
+- at node Invoking astscript-psf-unite, Invoking astscript-psf-scale-factor, Invoking astscript-psf-stamp, PSF construction and subtraction
++ at node astscript-psf-unite, astscript-psf-scale-factor, astscript-psf-stamp, PSF construction and subtraction
+ @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.
+@@ -25356,7 +25356,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.
+-A convenient way of obtaining this value is by using the script @file{astscript-model-scale-factor}, see @ref{Invoking astscript-psf-scale-factor}.
++A convenient way of obtaining this value is by using the script @file{astscript-model-scale-factor}, see @ref{astscript-psf-scale-factor}.
+ 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.
+
+@@ -25403,7 +25403,7 @@ This option is useful for debugging and
+
+
+
+- at node Invoking astscript-psf-scale-factor, Invoking astscript-psf-subtract, Invoking astscript-psf-unite, PSF construction and subtraction
++ at node astscript-psf-scale-factor, astscript-psf-subtract, astscript-psf-unite, PSF construction and subtraction
+ @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.
+@@ -25450,9 +25450,9 @@ To fully understand the process on first
+ The most common use-cases of this scale factor are:
+ @enumerate
+ @item
+-To find the factor for joining two different parts of the same PSF, see @ref{Invoking astscript-psf-unite}.
++To find the factor for joining two different parts of the same PSF, see @ref{astscript-psf-unite}.
+ @item
+-When modeling a star in order to subtract it using the PSF, see @ref{Invoking astscript-psf-subtract}.
++When modeling a star in order to subtract it using the PSF, see @ref{astscript-psf-subtract}.
+ @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}.
+@@ -25544,7 +25544,7 @@ This option is useful for debugging and
+
+
+
+- at node Invoking astscript-psf-subtract, , Invoking astscript-psf-scale-factor, PSF construction and subtraction
++ at node astscript-psf-subtract, , astscript-psf-scale-factor, PSF construction and subtraction
+ @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.
+@@ -25628,7 +25628,7 @@ If the central position doesn't 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.
+-A convenient way of obtaining this value is by using the script @file{astscript-scale-factor}, see @ref{Invoking astscript-psf-scale-factor}.
++A convenient way of obtaining this value is by using the script @file{astscript-scale-factor}, see @ref{astscript-psf-scale-factor}.
+ For a full tutorial on using the @command{astscript-psf-*} scripts together, see @ref{Building the extended PSF}.
+
+ @item -t
+@@ -26102,10 +26102,10 @@ Please see @ref{Optional dependencies} f
@end cartouche
@menu
@@ -1333,7 +1539,7 @@
@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.
-@@ -32094,7 +32094,7 @@ philosophy}.
+@@ -35057,7 +35057,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
@@ -1342,7 +1548,7 @@
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
-@@ -32112,7 +32112,7 @@ after reading a long text, nothing else.
+@@ -35075,7 +35075,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
@@ -1351,3 +1557,31 @@
have in mind) at this stage too. Actually starting to write it here will
really help you later when you are coding.
+@@ -36230,22 +36230,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, its 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/f89484e7f80d96436908549ade8eacf6c222a566
More information about the pld-cvs-commit
mailing list