---
title: 'AutoSourceID: Modular Astrophysical Source Detection'
url: https://www.emergentmind.com/topics/autosourceid-asid
type: topic
---

# AutoSourceID: Modular Astrophysical Source Detection

AutoSourceID (ASID) is a modular deep-learning framework for astrophysical source detection that combines learned segmentation or representation modules with lightweight classical post-processing and downstream characterization networks. In its original gamma-ray formulation, ASID used a U-shaped convolutional network for source segmentation, alternative clustering stages for localization, and a downstream classifier for source typing [2103.11068]. In optical imaging, the best-known instantiation is AutoSourceID-Light (ASID-L), which localizes point sources in MeerLICHT images by chaining a U-Net with a single-scale Laplacian-of-Gaussian (LoG) filter [2202.00489]. The framework was subsequently extended with AutoSourceID-Classifier (ASID-C) for star-galaxy classification and AutoSourceID-FeatureExtractor (ASID-FE) for flux and centroid regression with uncertainty estimation [2307.14456], [2305.14495]. A later end-to-end gamma-ray pipeline generalized the same design pattern to multi-energy count maps from Fermi-LAT and CTAO, combining segmentation, LoG candidate extraction, classification, and ensemble regressors [2509.25128].

## 1. Origins and conceptual scope

The ASID line began with gamma-ray point-source identification in Fermi-LAT data. That system used a five-level U-shaped fully convolutional network on \(64\times 64\times 5\) energy-binned patches, followed by either \(k\)-means post-processing or Centroid-Net for localization, and then a three-class classifier that separated AGNs, PSRs, and FAKE sources [2103.11068]. In that setting, ASID was positioned as an alternative to traditional maximum-likelihood pipelines, with comparable flux thresholds and improved robustness to changes in the interstellar emission model.

ASID-L transferred the same hybrid philosophy to optical survey images. Its key design choice was to replace iterative or scale-search localization with a two-stage procedure in which a U-Net produces a soft mask of likely point sources and a fixed-scale LoG filter estimates source centers [2202.00489]. The method was motivated by the increasing survey speed of optical wide-field telescopes and the need for rapid and reliable transient-oriented source localization.

Later work made the broader scope of AutoSourceID explicit. ASID-C attached a star-galaxy classifier to ASID-L cutouts [2307.14456]; ASID-FE added a regression module for flux and sub-pixel position with calibrated uncertainties [2305.14495]; and the 2025 gamma-ray pipeline described ASID as a versatile framework already tested on Fermi-LAT simulated data, MeerLICHT optical data, CTAO simulated data, and transfer experiments on HST and WISE imagery [2509.25128]. This suggests that “ASID” is best understood not as a single network, but as a family of pipelines built around learned source proposal generation plus task-specific post-processing.

## 2. Optical localization pipeline in ASID-L

ASID-L converts a raw optical image into a catalog of source coordinates in two stages. First, the full-frame CCD image of \(10\,500\times 10\,500\) pixels is split into overlapping or non-overlapping \(256\times 256\)-pixel patches. Each patch is passed through a U-Net that outputs, for every pixel, a value in \([0,1]\) indicating the likelihood of belonging to a point source. The predicted mask is not a discrete label map; it is a soft, circular “blob” for each source [2202.00489].

Second, the continuous mask is thresholded at \(\tau=0.20\) to suppress background noise, and a single-scale LoG filter with \(\sigma=1.43\) is convolved with the mask. Local maxima in the LoG response are retained as source centers, subject to a maximum overlap rule \(\omega=0.8\). Pixel coordinates from all patches are then concatenated and transformed back to sky coordinates through the standard astrometric solution [2202.00489].

The LoG kernel used for localization is

\[
\mathrm{LoG}(x,y;\sigma)
\;=\;
-\frac{1}{\pi\sigma^4}\,
\Bigl[\,1 - \tfrac{x^2+y^2}{2\sigma^2}\Bigr]\,
\exp\!\Bigl(-\frac{x^2+y^2}{2\sigma^2}\Bigr).
\]

The response map is computed as

\[
R(x,y)\;=\;\bigl(\,m \;*\;\mathrm{LoG}\bigr)(x,y),
\]

where \(m\) is the U-Net soft mask. Peaks in \(R\) correspond to blob centers of radius \(\simeq \sqrt{2}\,\sigma\), and fixing \(\sigma\) avoids expensive multi-scale loops because the LoG is run only once per patch [2202.00489].

Within this division of labor, the U-Net is responsible for background estimation, artifact rejection, and simple deblending, while the LoG acts as a lightweight centroiding stage. The artifacts explicitly named in the description are cosmic rays, diffraction spikes, and satellite trails [2202.00489].

## 3. Network architecture and objective functions

The ASID-L segmentation network follows the canonical U-Net design, adapted to small-object segmentation. Its encoder has five convolutional levels down to \(8\times 8\) pixels. At each level, two \(3\times 3\) convolutions with padding \(=1\) and ReLU activations are followed by \(2\times 2\) max-pooling with stride 2. The number of channels doubles after each pooling step, for example \(32\rightarrow 64\rightarrow 128\rightarrow 256\rightarrow 512\). The bottleneck consists of two \(3\times 3\) convolution-plus-ReLU layers at \(8\times 8\) resolution with 1024 channels. The decoder mirrors the encoder: each of its five levels uses a \(2\times 2\) up-convolution to halve the channel count and double spatial size, concatenates the corresponding encoder feature map through a skip connection, and applies two \(3\times 3\) convolution-plus-ReLU layers. A final \(1\times 1\) convolution to one channel followed by a Sigmoid activation produces the \([0,1]\) mask. The model has approximately \(2\) million trainable parameters [2202.00489].

Training uses a “Combo” loss equal to the sum of binary cross-entropy and Dice loss:

\[
L_{\rm BCE}
= -\frac{1}{n} \sum_{i=1}^n
\left[y_i \log \hat y_i + (1-y_i)\log(1-\hat y_i)\right]
\]

\[
L_{\rm Dice}
= 1 -
\frac{2 \sum y_i \hat y_i}
{\sum y_i^2 + \sum \hat y_i^2}
\]

\[
\text{Total loss} = L_{\rm BCE} + L_{\rm Dice}.
\]

The stated purpose of this combination is to balance pixel-wise accuracy with overlap-based segmentation quality in a highly imbalanced setting with few source pixels and many background pixels [2202.00489].

The broader ASID family preserves this modular logic but changes the output heads and losses according to the task. In the 2025 gamma-ray formulation, the segmentation module is a Multi-Input U-Net operating on energy-binned count maps, the classifier is a VGG-like CNN that outputs \(p_{\rm src}=P(\text{True source})\), and two deep-ensemble regressors estimate flux \(F\) and a residual \((\Delta x,\Delta y)\) for location refinement [2509.25128]. That formulation uses binary cross-entropy for segmentation and classification, while flux and localization can use either mean-squared-error or a negative-log-Gaussian-likelihood [2509.25128].

## 4. Training data, supervision, and empirical performance on MeerLICHT

ASID-L was trained on three fields from the MeerLICHT \(q\)-band: the Omega Cen globular cluster, the Fornax galaxy cluster, and the Chandra Deep Field-South. These yielded \(3\times 1681 = 5043\) patches of size \(256\times 256\) pixels. Gaia EDR3 provided the reference catalog used to place “true” point-source masks; the mask geometry was a \(3\times 3\) central square plus four cardinal pixels approximating a small circle. Only Gaia sources with predicted \(S/N \ge 3\), estimated from local zero point, PSF, and sky plus read noise, were included in the masks [2202.00489].

The split was \(80\%\) train, \(10\%\) validation, and \(10\%\) test, corresponding to approximately 4034, 504, and 504 patches. In addition, the 165 most crowded Omega Cen patches formed an independent test set for detailed evaluation [2202.00489]. The paper does not explicitly describe standard augmentations such as rotations, flips, or photometric jitter. It also states that optimizer and learning-rate schedule are not spelled out, although convergence was assessed through stabilization of the combined loss on the validation set and training stopped when no improvement was seen after a predefined patience [2202.00489].

Performance was quantified through true positives, false negatives, false positives, and the Dice coefficient

\[
\mathrm{Dice} = \frac{2\,TP}{2\,TP + FP + FN}.
\]

On the Omega Cen test set, Gaia EDR3 ground truth corresponded to \(36.958\) sources arcmin\(^{-2}\). SExtractor, with default settings except \(\texttt{DETECT\_MINAREA}=3\) and \(\texttt{BACK\_SIZE}=60\), achieved \(TP=14.634\), \(FN=22.324\), \(FP=0.178\), and \(\mathrm{Dice}=0.5653\). ASID-L, with \(\sigma=1.43\), \(\tau=0.20\), and \(S/N\ge 3\), achieved \(TP=27.001\), \(FN=9.957\), \(FP=2.918\), and \(\mathrm{Dice}=0.8075\) [2202.00489].

The density dependence is central to the method’s positioning. ASID-L maintained \(\mathrm{Dice}\sim 0.75\)–\(0.9\) across \(10\)–\(300\) sources arcmin\(^{-2}\), whereas SExtractor’s Dice dropped sharply above approximately \(100\) sources arcmin\(^{-2}\) [2202.00489]. Computation time was measured on \(2560\times 2560\) pixel images, approximately \(580\) arcmin\(^2\), on an Intel i9 plus RTX 2080 system. For source densities of \(10\), \(75\), and \(250\) sources arcmin\(^{-2}\), SExtractor required \(1.97\) s, \(10.22\) s, and \(25.25\) s, while ASID-L required \(0.30\) s, \(0.28\) s, and \(0.30\) s, respectively [2202.00489]. The method therefore had essentially constant runtime with source density and was reported as giving a factor \(>7\) speed-up in sparse fields and \(>80\times\) in crowded fields.

## 5. Extension to classification and feature extraction

ASID-C extends the localization stage into a two-stage optical pipeline in which ASID-L first detects and centers sources and then a dedicated classifier labels the resulting cutouts as stars or galaxies [2307.14456]. ASID-L processes full MeerLICHT images from the \(65\) cm telescope and outputs \(34\times 34\) pixel cutouts with source centers. These cutouts are matched to DECaLS DR10 morphological labels and passed to a dual-branch CNN. One branch ingests a single-band \(32\times 32\) cutout, randomly shifted by \(\pm 1\) pixel from the original \(34\times 34\) stamp; the other ingests normalized pixel coordinates \([x,y]\) within the full frame. Their fused representation is mapped to a Sigmoid output that yields a raw score \(f\), and a modified Platt scaling stage then produces calibrated posterior probabilities \(p\in[0,1]\) for star versus galaxy [2307.14456].

The ASID-C training set comprises \(12\) million images from \(718\) “green” MeerLICHT exposures across SDSS \(ugriz+q\) filters, with a class imbalance of approximately \(10{:}1\) stars to galaxies. Splits are \(50\%\) training, \(20\%\) validation, \(20\%\) test, and \(10\%\) calibration. The network was trained with Adam, initial learning rate \(\eta_0=10^{-3}\), exponential decay \(\eta_t=\eta_0\cdot 0.99^{(t/8\ \mathrm{epochs})}\), and early stopping after \(10\) epochs without validation improvement [2307.14456]. Quantitatively, ASID-C maintained \(>95\%\) ROC-AUC down to \(S/N\approx 4\), produced sharply bimodal well-calibrated predictions, and processed approximately \(28\)k cutouts per second on a GeForce RTX2080, corresponding to about \(36\,\mu\)s per source [2307.14456].

ASID-FE addresses a different downstream problem: continuous feature regression with uncertainty characterization. It takes \(32\times 32\) cutouts around localized sources from ASID-L or other catalogs and estimates flux and sub-pixel center coordinates together with their uncertainties [2305.14495]. The method uses a two-step mean-variance estimation (TS-MVE) architecture. Part I predicts \(\mu_x\), \(\mu_y\), and \(\mu_{\rm flux}\) through three convolutional blocks and three independent MLP tails; Part II reuses the mean branch, adds a parallel variance branch, and outputs refined means and \(\sigma_x\), \(\sigma_y\), and \(\sigma_{\rm flux}\), assuming diagonal covariance [2305.14495].

ASID-FE was trained on synthetic cutouts derived from real MeerLICHT images. PSFEx provided position-dependent PSF models; Gaia DR2 sources were projected onto blank full-frame arrays; sky background and its \(\sigma\) were taken from real pipeline products; and approximately \(3.5\) million \(34\times 34\) cutouts were generated and randomly shifted before center-cropping to \(32\times 32\) [2305.14495]. On synthetic tests, ASID-FE achieved centroid uncertainty calibration \(\mathrm{WRSD}\simeq 1.00\pm 0.05\) across \(S/N\), flux uncertainty calibration near \(1.00\), and \(5\)–\(10\%\) lower flux MARE than SourceExtractor’s FLUX\_AUTO. Relative to a small deep-ensemble MVE without the two-step procedure, TS-MVE Part I already yielded approximately \(20\)–\(30\%\) lower MSE and MAE on the mean predictions, and Part II further lowered the negative log-likelihood by a similar margin [2305.14495].

## 6. Generalization to gamma-ray surveys, limitations, and future directions

The 2025 end-to-end ASID pipeline generalizes the AutoSourceID pattern from optical images and earlier gamma-ray work to multi-energy count maps from Fermi-LAT and CTAO [2509.25128]. Its first stage is a Multi-Input U-Net operating on six energy slices for Fermi-LAT or three for CTAO, each covering a \(10^\circ\times 10^\circ\) sky patch. A LoG clustering step extracts discrete source candidates from the segmentation map, a VGG-like CNN classifies each candidate as a true source or a fake fluctuation, and two deep-ensemble networks then estimate flux and refine the location [2509.25128].

In that formulation, the inputs are \(N_E\)-channel count maps with \(N_E=6\) for Fermi-LAT and \(N_E=3\) for CTAO, while the outputs for each detected source are sky location \((l_i,b_i)\), estimated flux \(\hat F_i\pm \sigma_{F,i}\), and classification score \(p_{{\rm src},i}\) [2509.25128]. Reported results include a Fermi-LAT flux sensitivity of approximately \(2\times 10^{-10}\,\mathrm{cm^{-2}\,s^{-1}}\) at \(90\%\) completeness, a \(98\%\) association rate on real high-latitude 4FGL-DR2 sources with \(\sigma>20\) and \(|b|>20^\circ\), approximately \(93\%\) association for two more difficult source subsets, and CTAO recall reaching \(90\%\) at \(F(>1\;\mathrm{TeV})\approx 2\times 10^{-14}\,\mathrm{cm^{-2}\,s^{-1}}\), consistent with standard likelihood analyses in Gammapy [2509.25128].

The same paper reports latent-space clustering in which Fermi-LAT and CTAO source patches form one cluster while background patches form another, and interprets this as evidence that a single “foundation” encoder could serve multiple instruments [2509.25128]. This suggests a path toward jointly trained, heterogeneous source-detection models, although that remains prospective rather than established.

The documented limitations of ASID-L remain specific and important. Current U-Net training is limited to point-like masks derived from Gaia, so extended sources such as galaxies and nebulae are not represented. The method has not been exhaustively tested on images with dramatically different PSFs or backgrounds, although preliminary transfer tests on HST and WISE were described as encouraging. ASID-L currently outputs only centroid coordinates; photometry, shape parameters, and classification are delegated to follow-up modules [2202.00489]. Proposed extensions include adding a refinement step for flux and PSF metrics, incorporating extended-source segmentation masks, developing a full real-time multi-band ASID that ingests all filters simultaneously, and carrying out more thorough transfer-learning experiments across telescopes and instruments [2202.00489].

Taken together, these developments define AutoSourceID as a layered astronomical analysis framework: learned localization through segmentation, lightweight centroid extraction, and modular downstream inference for classification or regression. The available evidence shows that this design is effective in both optical and gamma-ray point-source detection, especially in crowded or background-dominated regimes where classical heuristics or likelihood scans become slow, brittle, or both [2202.00489].

Source: https://www.emergentmind.com/topics/autosourceid-asid