Package: sep-pjw
SEP-PJWPython and C library for Source Extraction and Photometry, forked fromkbarbary/septo provide additional features and bug fixes."... [it's] an SEP: Somebody Else's Problem."
"Oh, good. I can relax then."SEP-PJW vs. SEPThe original release ofsepby Kyle Barbary,kbarbary/sep, no longer appears to be
maintained. Whilst the package is incredibly useful, there are a few
outstanding bugs, and support for the latest versions of Python
(python>=3.11) will be limited. The aim ofsep-pjwis to offer a
version ofsepthat resolves these issues, whilst maintaining
compatibility as much as is feasibly possible. Any fixes or updates to
the original library are documented inCHANGES.md. For
existing workflows, the only necessary update will be to change the import
toimportsep_pjwassepAboutSource Extractor(Bertin & Arnouts 1996) is a widely used
command-line program for segmentation and analysis of astronomical
images. It reads in FITS format files, performs a configurable series
of tasks, including background estimation, source detection,
deblending and a wide array of source measurements, and finally
outputs a FITS format catalog file.While Source Extractor is highly useful, the fact that it can only be
used as an executable can limit its applicability or lead to awkward
workflows. There is often a desire to have programmatic access to
perform one or more of the above tasks on in-memory images as part of
a larger custom analysis.SEP makes the core algorithms of Source Extractor available as a
library of stand-alone functions and classes.These operate directly
on in-memory arrays (no FITS files or configuration files).  The code
is derived from the Source Extractor code base (written in C) and aims
to produce results compatible with Source Extractor whenever possible.
SEP consists of a C library with no dependencies outside the standard
library, and a Python module that wraps the C library in a Pythonic
API. The Python wrapper operates on NumPy arrays with NumPy as its
only dependency. See below for language-specfic build and usage
instructions.PythonDocumentation:http://sep-pjw.readthedocs.ioRequirements:Python 3.9+numpy 1.23+Install release version:SEP-PJW can be installed withpip:python -m pip install sep-pjwIf you get an error about permissions, you are probably using your
system Python. In this case, we recommend usingpip's "user
install"option to install sep-pjw into your user directory:python -m pip install --user sep-pjwDonotinstall sep or other third-party Python packages usingsudounless you are fully aware of the risks.Install development version:Building the development version (from github) requires Cython.
Build and install in the usual place:python -m pip install --editable .Run tests:Tests require thepytestPython
package. To run the tests, execute./test.pyin the top-level
directory. Some tests require a FITS reader (either fitsio or astropy)
and will be skipped if neither is present.C LibraryNote: The build process only works on Linux and OS X.CMake:To build using CMake, enter these commands:cd sep
mkdir -p build
cd build
cmake -DCMAKE_INSTALL_PREFIX=/usr ../and follow next steps from the build folderBuild:To build the C library from source:makeRun tests:make testInstallThe static library and header can be installed withmake install
make PREFIX=/path/to/prefix installThis will install the shared and static library in/path/to/prefix/liband header file in/path/to/prefix/include. The default prefix is/usr/local.API:The C library API is documented in the header filesep.h.Rust bindings:Low-level Rust wrapper for the C library can be found athttps://crates.io/crates/sep-sysContributingReport a bug or documentation issue:http://github.com/PJ-Watson/sep-pjw/issuesAsk (or answer) a question:https://github.com/PJ-Watson/sep-pjw/discussions/categories/q-aDevelopment of SEP-PJW takes place on GitHub atPJ-Watson/sep-pjw.  Contributions of bug fixes,
documentation improvements and minor feature additions are welcome via
GitHub pull requests. For major features, it is best to discuss the change first
viaGitHub Discussions.CitationIf you use SEP in a publication, please cite the following article in the Journal of Open Source Software:Please also cite the original Source Extractor paper (Bertin & Arnouts
1996).The DOI for the sep v1.0.0 code release is:LicenseThe license for all parts of the code derived from Source Extractor is
LGPLv3. The license for code derived from photutils (src/overlap.h)
is BSD 3-clause. Finally, the license for the Python wrapper
(sep.pyx) is MIT. The license for the library as a whole is
therefore LGPLv3. The license for each file is explicitly stated at
the top of the file and the full text of each license can be found inlicenses.FAQWhy isn't the C library part of Source Extractor?Source Extractor isnotdesigned as a library with an
executable built on top of the library. In Source Extractor, background
estimation, object detection and photometry are deeply integrated into the
Source Extractor executable. Many changes to the code were necessary in
order to put the functionality in stand-alone C functions. It's too much
to ask of the Source Extractor developer to rewrite large parts of the
core of the Source Extractor program with little gain for the executable.What sort of changes?Source Extractor reads in only a small portion of each image at a
time.  This allows it to keep its memory footprint extremely low and
to operate on images that are much larger than the system's physical
memory. It also means that a FITS reader is deeply integrated into
the code.  SEP operates on images in memory, so all the FITS I/O
machinery in Source Extractor is not used here.Error handling: When it encounters a problem, Source Extractor
immediately exits with an error message. This is fine for an
executable, but a library function doesn't have that luxury. Instead
it must ensure that allocated memory is freed and return an error
code.Options: Source Extractor has many options that affect its behavior. These
are stored in a global structure used throughout the executable. In SEP,
options for a particular function are passed as function parameters.Array types: Source Extractor can operate on FITS images containing various
types of data (float, double, int, etc). Internally, it does this by
converting all data tofloatimmediately when reading from disk.
SEP does something similar, but in memory: SEP functions typically convert
input arrays to float on the fly within each function, then perform
all operations as floating point.Is SEP as fast as Source Extractor?It's fast. It should be similar to Source Extractor as a lot of the code
is identical. Source Extractor has the advantage of doing all the
operations (detection and analysis) simultaneously on each image
section, which may confer CPU cache advantages, but this hasn't been
tested at all. On the other hand, depending on your usage SEP might
let you avoid writing files to disk, which is likely to be a bigger
win.What happens when Source Extractor is updated in the future?SEP can be considered a fork of the Source Extractor code base: it's
development will not track that of Source Extractor in any automated
way. However, the algorithms implemented so far in SEP are stable in
Source Extractor: the SEP code was forked from v2.18.11, yet it is tested
against the results of v2.8.6. This indicates that the algorithms have
not changed in Source Extractor over the last few years.In the Python interface, why do I have to byte swap data when using
astropy.io.fits?This occurs because FITS files have big-endianbyte
order, whereas most widely
used CPUs have little-endian byte order. In order for the CPU to
operate on the data, it must be byte swapped atsome point.Some
FITS readers such asfitsiodo
the byte swap immediately when reading the data from disk to memory,
returning numpy arrays in native (little-endian) byte order. However,
astropy.io.fits does not (for reasons having to do with memory
mapping). Most of the time you never notice this because when you do
any numpy operations on such arrays, numpy uses an intermediate buffer
to byte swap the array behind the scenes and returns the result as a
native byte order array. Internally, SEP is not using numpy
operations; it's just getting a pointer to the data in the array and
passing it to C code. As the C code does not include functionality to
do buffered byte swapping, the input array must already be in native
byte order.It would be possible to add buffered byte swapping capability to the
SEP code, but it would increase the code complexity. A simpler
alternative would be to make a byte swapped copy of the entire input
array, whenever necessary. However, this would significantly increase
memory use, and would have to be done repeatedly in multiple SEP
functions:Background,extract,sum_circle, etc. Each would make
a copy of the entire data array. Given these considerations, it seemed
best to just explicitly tell the user to do the byte swap operation
themselves so they could just do it once, immediately after reading in
the data.I have more questions!Open a discussion on theGitHub Discussions page!