Search

Deep neural network model extraction attack is the process of retraining a surrogate model based on the outputs of a target model with a given set of inputs. Such attacks are hard to defend for the sake of model owners' interest. Recently, some work propose model watermarking scheme for image processing networks, which is able to prove the intellectual property of deep models even after the model extraction attack. This scheme makes sure that, once the target model (an image processing network) is watermarked, we can extract the watermark from the output of the surrogate model. In this paper, we propose a new model extraction attack scheme to fight against the latest method. Instead of directly using the output images of a target model, we propose to use their reconstructed versions for model retraining, where an asymmetrical UNet is proposed for image reconstruction. To thoroughly remove the watermarking traces, we propose and incorporate a referenced subspace attention module in the asymmetrical UNet, which removes the watermark by projecting the outputs of the target model into the subspaces of the reference image. Various experiments demonstrate the effectiveness of our attack. [ABSTRACT FROM AUTHOR]

Motivation BISCAP is a state-of-the-art tool for automatically characterizing biofilm images obtained from Optical Coherence Tomography. Limited availability of other software tools is reported in the field. BISCAP's first version processes 2D images only. Processing 3D images is a problem of greater scientific relevance since it deals with the entire structure of biofilms instead of their 2D slices. Results Building on the image-processing principles and algorithms proposed earlier for 2D images, these were adapted to the 3D case, and a more general implementation of BISCAP was developed. The primary goal concerns the extension of the initial methodology to incorporate the depth axis in 3D images; multiple improvements were also made to boost computational performance. The calculation of structural properties and visual outputs was extended to offer new insights into the 3D structure of biofilms. BISCAP was tested using 3D images of biofilms with different morphologies, consistently delivering accurate characterizations of 3D structures in a few minutes using standard laptop machines. Low user dependency is required for image analysis. Availability and implementation BISCAP is available from https://github.com/diogonarciso/BISCAP. All images used in the tutorials and the validation examples are available from https://web.fe.up.pt/∼fgm/biscap3d. [ABSTRACT FROM AUTHOR]

We present a search for transient radio sources on time-scales of seconds to hours at 144 MHz using the LOFAR Two-metre Sky Survey (LoTSS). This search is conducted by examining short time-scale images derived from the LoTSS data. To allow imaging of LoTSS on short time-scales, a modern imaging procedure and fast filtering strategy are introduced. This includes sky model source subtraction, no cleaning or primary beam correction, a simple source finder, fast filtering schemes, and source catalogue matching. This new strategy is first tested by injecting simulated transients, with a range of flux densities and durations, into the data. We find the limiting sensitivity to be 113 and 6 mJy for 8 s and 1 h transients, respectively. The new imaging and filtering strategies are applied to 58 fields of the LoTSS survey, corresponding to LoTSS-DR1 (2 per cent of the survey). One transient source is identified in the 8 s and 2 min snapshot images. The source shows 1 min duration flare in the 8 h observation. Our method puts the most sensitive constraints on/estimates of the transient surface density at low frequencies at time-scales of seconds to hours; <4.0 × 10−4 deg−2 at 1 h at a sensitivity of 6.3 mJy; 5.7 × 10−7 deg−2 at 2 min at a sensitivity of 30 mJy; and 3.6 × 10−8 deg−2 at 8 s at a sensitivity of 113 mJy. In the future, we plan to apply the strategies presented in this paper to all LoTSS data. [ABSTRACT FROM AUTHOR]

We present extended point spread function (PSF) models for the Hyper Suprime-Cam Subaru Strategic Program Public Data Release 3 (HSC-SSP PDR3) in all g, r, i, Z , and Y -bands. Due to its 8.2 m primary mirror and long exposure periods, HSC combines deep images with wide-field coverage. Both properties make HSC one of the most suitable observing facilities for low-surface brightness (LSB) studies, which are particularly sensitive to the PSF. By applying a median stacking technique of point-like sources with different brightness, we show how to construct the HSC-SSP PDR3 PSF models to an extent of R ∼ 5.6 arcmin. These models are appropriate for the HSC-PDR3 intermediate-state data which do not have applied the final aggressive background subtraction. The intermediate-state data is especially stored for users interested in large extended objects, where our new PSFs provide them with a crucial tool to characterise LSB properties at large angles. We demonstrate that our HSC PSFs behave reasonably in two scenarios. In the first one, we generate 2D models of a bright star, showing no evidence of residual structures across the five bands. In the second scenario, we recreate the PSF-scattered light on mock images with special consideration of the effect of this additional flux on LSB measurements. Despite the well-behaved nature of the HSC-PDR3 PSFs, there is a non-negligible impact on the faint light present in the mock images. This impact could lead to incorrect LSB measurements if a proper star subtraction is not applied. The PSF models are publicly available to the astronomical community at github.com/luciagarate/HSC_PSFs (DOI: 10.5281/zenodo.11273118). [ABSTRACT FROM AUTHOR]

Adaptive optics (AO) systems correct optical phase aberrations of the incoming light generated by the atmosphere. To do so simultaneous estimators of the atmospheric turbulence parameters are required. For the family of wide-field AO (WFAO) systems, this information must be stratified in altitude. Among these vectorized estimations, wind profiling in altitude is needed for the reduction of temporal errors in AO loops or for the estimation of turbulence coherence time. This paper proposes a turbulence wind profiler called image processing based peak tracking algorithm (ipta). ipta is an image-processing based approach that automatically and reliably estimates wind velocity for several turbulent layers along the line of sight. The estimation of each wind layer is achieved by tracking peaks produced in cross-correlation maps from pairs of wavefront sensors (WFSs) slopes using the technique known as SLODAR (slope detection and ranging). Results for simulated and on-sky WFS data sets demonstrate that ipta outperforms one of the state of the art wind profiler methods [the profiler covariance parametrization of wind velocity (CAW)] in terms of accuracy and speed. Results also show that, in terms of execution time, our method scales better when the number of WFS lenslets is increased. Being an open source and reliable tool, we believe ipta can be a useful wind profiler for the AO community. [ABSTRACT FROM AUTHOR]

Seismic attributes are derived measures from seismic data that help characterize subsurface geological features and enhance the interpretation of subsurface structures: we propose to exploit the hidden layers of Long–Short Time Memory neural network predictions as possible new reflection seismic attributes. The idea is based on the inference process of a neural network, which in its hidden layers stores information related to different features embedded in the input data and which usually are not considered. Neural network applications typically ignore such intermediate steps because the main interest lies in the final output, which is considered as the exclusive exploitable feature of the process. On the contrary, here we analyse the possibility to exploit the intermediate prediction steps, hereafter referred as 'deep attributes' because they are produced by a deep learning algorithm, to highlight features and emphasize characteristics embedded in the data but neither recognizable by traditional interpretation, nor evident with classical attributes or multi-attribute approaches. Nowadays, classical signal attributes are numerous and used for different purposes; we here propose an original strategy to calculate attributes previously never exploited, which are potentially complementary or a good alternative to the classical ones. We tested the proposed procedure on synthetic and field 2-D and 3-D reflection seismic data sets to test and demonstrate the stability, affordability and versatility of the entire approach. Furthermore, we evaluated the performance of deep attributes on a 4-D seismic data set to assess the applicability and effectiveness for time-monitoring purposes and comparing them with the sweetness attribute. [ABSTRACT FROM AUTHOR]

The upcoming Nancy Grace Roman Space Telescope will carry out a wide-area survey in the near-infrared. A key science objective is the measurement of cosmic structure via weak gravitational lensing. Roman data will be undersampled, which introduces new challenges in the measurement of source galaxy shapes; a potential solution is to use linear algebra-based coaddition techniques such as imcom that combine multiple undersampled images to produce a single oversampled output mosaic with a desired 'target' point spread function (PSF). We present here an initial application of imcom to 0.64 square degrees of simulated Roman data, based on the Roman branch of the Legacy Survey of Space and Time (LSST) Dark Energy Science Collaboration (DESC) Data Challenge 2 (DC2) simulation. We show that imcom runs successfully on simulated data that includes features such as plate scale distortions, chip gaps, detector defects, and cosmic ray masks. We simultaneously propagate grids of injected sources and simulated noise fields as well as the full simulation. We quantify the residual deviations of the PSF from the target (the 'leakage'), as well as noise properties of the output images; we discuss how the overall tiling pattern as well as Moiré patterns appear in the final leakage and noise maps. We include appendices on interpolation algorithms and the interaction of undersampling with image processing operations that may be of broader applicability. The companion paper ('Paper II') explores the implications for weak lensing analyses. [ABSTRACT FROM AUTHOR]

Interferometric synthetic aperture radar (InSAR) decorrelation that creates great challenges to phase unwrapping has been a critical issue for mapping large earthquake deformation. Some studies have proposed a 'remove-and-return model' solution to tackle this problem, but it has not been fully validated yet, and therefore has rarely been applied to real earthquake scenarios. In this study, we use the 2023 M w 7.8 and 7.6 earthquake doublet in Turkey and Syria as a case example to develop an iterative modelling method for InSAR-based coseismic mapping. We first derive surface deformation fields using Sentinel-1 offset tracking and Sentinel-2 optical image correlation, and invert them for an initial coseismic slip model, based on which we simulate InSAR coseismic phase measurements. We then remove the simulated phase from the actual Sentinel-1 phase and conduct unwrapping. The simulated phase is added back to the unwrapped phase to produce the final phase measurements. Comparing to the commonly used unwrapping method, our proposed approach can significantly improve coherence and reduce phase gradients, enabling accurate InSAR measurements. Combining InSAR, offset tracking and optical image correlation, we implement a joint inversion to obtain an optimal coseismic slip model. Our model shows that slip on the Çardak Fault is concentrated on a ∼100 km segment; to both ends, slip suddenly diminished. On the contrary, rupture on the East Anatolian Fault Zone propagated much longer as its geometry is fairly smooth. The iterative coseismic modelling method is proven efficient and can be easily applied to other continental earthquakes. [ABSTRACT FROM AUTHOR]

Muography is a promising and rapidly developing physical prospecting technique based on the attenuation of muon flux. The feasibility and potential of applying muography to mining were presented in a small number of previous case studies. In this work, we aimed to address three challenges that limit the applicability and efficiency of muography in mineral exploration: (1) application to low-density-contrast ore body exploration, (2) analysis of altitudinal impacts on measurements and (3) precise and efficient reconstruction. We conducted the first case of applying muography to the exploration for low-density-contrast ore bodies. Six muon imaging systems were placed underground to collect surviving muons for roughly half a year. We analysed the altitudinal impact on the data measurements and proposed a simplified method to eliminate it. We also developed the seed algorithm, a novel 3-D reconstruction algorithm, that can reconstruct anomalies located far away from the detectors and avoid their elongation along the observed directions. Benefitting from the seed algorithm, a low-density-contrast ore body and a limonitic siliceous slate structure were reconstructed, demonstrating the sensitivity of this technique in density distinction; discoveries of several mined-out areas are important for accident avoidance; and reconstruction of the stope and scarps served as strong circumstantial evidence of the reliability of the result. The success of this experiment shows the great value of muography in the economic, research and safety aspects of mineral exploration and inspection. Moreover, the overcoming of challenges is helpful for the development of muography, making it an effective and competitive technique in this field. [ABSTRACT FROM AUTHOR]

This paper presents a new methodology to assist geophysicists in determining the first-break event in a 3-D seismic data set using the well-known technique called dynamic time warping algorithm (DTW), which is usually used to find the optimal alignment between two time-series. We used the optimal path from the cost matrix to identify the first break in the seismogram using a few picks (seeds) made by an interpreter as a reference to perform this task. Furthermore, the data were pre-conditioned by the topographic and linear moveout to improve the method's accuracy. To demonstrate the technique's robustness, first, we applied the methodology in a synthetic seismic data. After demonstrating the efficiency of the algorithm, we applied the aforementioned methodology in the Polo-Miranga 3-D seismic cube located in the Recôncavo sedimentary basin, Bahia-Brazil, and in the seismic data acquired from the Blackfoot field in Alberta, Canada. The high-quality results showed consistency in determining the first break in all ranges of offsets, demonstrating an alternative way to accelerate this seismic processing step. Furthermore, we compared the results obtained by the proposed methodology with an algorithm based on comparing the short-time averages with long-time averages. Finally, we performed the static correction calculation to ensure that the time distortion resulting from the terrain and the low-velocity layer was mitigated in shoot gathers and in the stacked section. [ABSTRACT FROM AUTHOR]

Dense array observation and seismic interferometry have revolutionized the imaging schemes of the earth structure. It is becoming possible to directly obtain the lateral variation of the earth's structure by applying array-based methods such as the cross-correlation beamforming (CBF) of the ambient noise to the subsets of the dense array, without tomography. CBF has been proven to extract the azimuth-averaged multimode surface wave dispersion curves. However, the resolution of the dispersion image generated by conventional CBF is low at high frequencies in the frequency–velocity (f - v) domain. Moreover, the irregular array geometry and uneven source distribution would bias the result of CBF, especially for the estimation of azimuth-dependent velocity. In this paper, two beamforming (BF) es are suggested to improve the resolution of multimode dispersion images in the f - v domain. First, the geometrical spreading of the wavefield is corrected to enhance the amplitude at high frequency (or large distance) and thereby improve the resolution of the dispersion image at high frequency. We call this scheme weighted correlation beamforming (WCBF). The azimuth-averaged velocity can be estimated with sufficient resolution using WCBF by stacking the BF output at each azimuth. We show that WCBF is the 2-D Fourier transform of the spatial wavefield from the viewpoint of the wavefield transform. Secondly, a modified beamforming scheme (MCBF) is suggested to reduce the effect of uneven source and/or irregular array geometry. The delay and summation in MCBF are performed only for plane waves incident from the stationary phase region. The azimuth-dependent velocity can therefore be estimated by MCBF with less dependence on the array geometry, as well as on the uneven source distribution. In terms of the estimation of azimuth-averaged phase velocity, we show the F-J method, another array-based method for extracting multimode surface waves from ambient noise using the Fourier–Bessel transform, is the azimuth-averaged version of WCBF. The reliability of WCBF and MCBF is verified based on the synthetic and field data using the array with different geometry. The dispersion image of multimode Rayleigh wave phase velocity at local and regional scales can be generated by WCBF or MCBF with high resolution. In particular, multimode dispersion curves at the local scale can be measured by MCBF with sufficient accuracy using quite short recordings from hours to days. This offers the possibility of a rapid assessment of the media properties. [ABSTRACT FROM AUTHOR]

Earthquake source parameters can be estimated using seismological observations, but the identification of the fault responsible is often complicated by location uncertainties and the inherent ambiguity between nodal planes. Satellite Interferometric Synthetic Aperture Radar (InSAR) can be used to observe ground deformation and model fault geometry but is limited by climate conditions (water vapour) and ground coverage (dense vegetation). In the tropics, the atmosphere is dynamic and most regions are densely vegetated, making detecting coseismic deformation challenging. Here, we perform a systematic inspection of coseismic interferograms from Sentinel-1 SAR images, to assess their suitability for detecting coseismic deformation in Costa Rica. Using data from the seismological network, we target seven earthquakes between 2016 and 2020 with depths |$\le \, 20$| km and magnitudes M w 5.3–6.2. For each event, we use the seismic parameters to compute line-of-sight displacements for ascending and descending geometries and for both nodal planes and generate 12- and 24-d coseismic interferograms where available. We obtain interferograms with coseismic displacement signals for three of the seven earthquakes. We invert the geodetic data to retrieve the earthquake source parameters but the lack of offshore geodetic coverage causes trade-offs between parameters and large uncertainties. The Jacó and Golfito earthquakes likely occurred on the subduction interface and the geodetic locations were 6–9 km closer to the coast than previous seismic estimates. The Burica earthquake occurred on a shallow steeply dipping thrust fault in the outer forearc. For the other earthquakes, no coseismic deformation was detected due to atmospheric noise or poor coherence. These results demonstrate the suitability of 12-d Sentinel-1 interferograms for monitoring shallow earthquakes of magnitude > M w 5.7 in Central America. This approach can be used to begin a surface deformation catalogue for the region, which will ultimately help improve the understanding of active deformation processes and improve hazard maps. [ABSTRACT FROM AUTHOR]

Gravity surveys constitute an important method for investigating the Earth's interior based on density contrasts related to Earth material differentials. Because lithology depends on the environment and the period of formation, there are generally clear boundaries between rocks with different lithologies. Inversions with convex functions for approximating the L 0 norm are used to detect boundaries in reconstructed models. Optimizations can easily be found because of the convex transformations; however, the volume of the reconstructed model depends on the weighting parameter and the density constraint rather than the model sparsity. To determine and adapt the modelling size, a novel non-convex framework for gravity inversion is proposed. The proposed optimization aims to directly reduce the L 0 norm of the density matrix. An improved iterative hard thresholding algorithm is developed to linearly reduce the L 0 penalty during the inner iteration. Accordingly, it is possible to determine the modelling scale during the iteration and achieve an expected scale for the reconstructed model. Both simple and complex model experiments demonstrate that the proposed method efficiently reconstructs models. In addition, granites formed during the Yanshanian and Indosinian periods in the Nanling region, China, are reconstructed according to the modelling size evaluated in agreement with the magnetotelluric profile and density statistics of rock samples. The known ores occur at the contact zones between the sedimentary rocks and the reconstructed Yanshanian granites. The ore-forming bodies, periods, and processes are identified, providing guidance for further deep resource exploration in the study area. [ABSTRACT FROM AUTHOR]

Motivation Cryo-electron microscopy (cryo-EM) is a powerful technique for determining the structures of large protein complexes. Picking single protein particles from cryo-EM micrographs (images) is a crucial step in reconstructing protein structures from them. However, the widely used template-based particle picking process requires some manual particle picking and is labor-intensive and time-consuming. Though machine learning and artificial intelligence (AI) can potentially automate particle picking, the current AI methods pick particles with low precision or low recall. The erroneously picked particles can severely reduce the quality of reconstructed protein structures, especially for the micrographs with low signal-to-noise ratio. Results To address these shortcomings, we devised CryoTransformer based on transformers, residual networks, and image processing techniques to accurately pick protein particles from cryo-EM micrographs. CryoTransformer was trained and tested on the largest labeled cryo-EM protein particle dataset—CryoPPP. It outperforms the current state-of-the-art machine learning methods of particle picking in terms of the resolution of 3D density maps reconstructed from the picked particles as well as F1-score, and is poised to facilitate the automation of the cryo-EM protein particle picking. Availability and implementation The source code and data for CryoTransformer are openly available at: https://github.com/jianlin-cheng/CryoTransformer. [ABSTRACT FROM AUTHOR]

The low surface brightness (LSB) regime (μg  ≳  26 mag arcsec−2) comprises a vast, mostly unexplored discovery space, from dwarf galaxies to the diffuse interstellar medium. Accessing this regime requires precisely removing instrumental signatures and light contamination, including, most critically, night sky emission. This is not trivial, as faint astrophysical and instrumental contamination can bias sky models at the precision needed to characterize LSB structures. Using idealized synthetic images, we assess how this bias impacts two common LSB-oriented sky-estimation algorithms: (1) masking and parametric modelling, and (2) stacking and smoothing dithered exposures. Undetected flux limits both methods by imposing a pedestal offset to all derived sky models. Careful, deep masking of fixed sources can mitigate this, but source density always imposes a fundamental limit. Stellar scattered light can contribute ∼28–29 mag arcsec−2 of background flux even in low-density fields; its removal is critical prior to sky estimation. For complex skies, image combining is an effective non-parametric approach, although it strongly depends on observing strategy and adds noise to images on the smoothing kernel scale. Preemptive subtraction of fixed sources may be the only practical approach for robust sky estimation. We thus tested a third algorithm, subtracting a preliminary sky-subtracted coadd from exposures to isolate sky emission. Unfortunately, initial errors in sky estimation propagate through all subsequent sky models, making the method impractical. For large-scale surveys like Legacy Survey of Space and Time, where key science goals constrain observing strategy, masking and modelling remain the optimal sky estimation approach, assuming stellar scattered light is removed first. [ABSTRACT FROM AUTHOR]

Disruptions to dopamine and noradrenergic neurotransmission are noted in several neurodegenerative and psychiatric disorders. Neuromelanin-sensitive (NM)-MRI offers a non-invasive approach to visualize and quantify the structural and functional integrity of the substantia nigra and locus coeruleus. This method may aid in the diagnosis and quantification of longitudinal changes of disease and could provide a stratification tool for predicting treatment success of pharmacological interventions targeting the dopaminergic and noradrenergic systems. Given the growing clinical interest in NM-MRI, understanding the contrast mechanisms that generate this signal is crucial for appropriate interpretation of NM-MRI outcomes and for the continued development of quantitative MRI biomarkers that assess disease severity and progression. To date, most studies associate NM-MRI measurements to the content of the neuromelanin pigment and/or density of neuromelanin-containing neurons, while recent studies suggest that the main source of the NM-MRI contrast is not the presence of neuromelanin but the high-water content in the dopaminergic and noradrenergic neurons. In this review, we consider the biological and physical basis for the NM-MRI contrast and discuss a wide range of interpretations of NM-MRI. We describe different acquisition and image processing approaches and discuss how these methods could be improved and standardized to facilitate large-scale multisite studies and translation into clinical use. We review the potential clinical applications in neurological and psychiatric disorders and the promise of NM-MRI as a biomarker of disease, and finally, we discuss the current limitations of NM-MRI that need to be addressed before this technique can be utilized as a biomarker and translated into clinical practice and offer suggestions for future research. [ABSTRACT FROM AUTHOR]

We report on three optically bright, ∼16th mag, point sources within 10 arcsec of each other that vanished within 1 h, based on two consecutive exposures at Palomar Observatory on 1952 July 19 (POSS I Red and Blue). The three point sources have continued to be absent in telescope exposures during 71 yr with detection thresholds of ∼21st mag. We obtained two deep exposures with the 10.4-m Gran Telescopio Canarias on 2023 April 25 and 27 in r and g band, both reaching magnitude 25.5 (3σ). The three point sources are still absent, implying they have dimmed by more than 10 mag within an hour back in 1952. When bright in 1952, the most isolated transient source has a profile nearly the same as comparison stars, implying the sources are subarcsec in angular size and they exhibit no elongation due to movement. This triple transient has observed properties similar to other cases where groups of transients ('multiple transients') have appeared and vanished in a small region within a plate exposure. The explanation for these three transients and the previously reported cases remains unclear. Models involving background objects that are optically luminous for less than 1 h coupled with foreground gravitational lensing seem plausible. If so, a significant population of massive objects with structure serving as the lenses, to produce three images, are required to explain the subhour transients. [ABSTRACT FROM AUTHOR]

Source conditions are a key tool in regularisation theory that are needed to derive error estimates and convergence rates for ill-posed inverse problems. In this paper, we provide a recipe to practically compute source condition elements as the solution of convex minimisation problems that can be solved with first-order algorithms. We demonstrate the validity of our approach by testing it on two inverse problem case studies in machine learning and image processing: sparse coefficient estimation of a polynomial via LASSO regression and recovering an image from a subset of the coefficients of its discrete Fourier transform. We further demonstrate that the proposed approach can easily be modified to solve the machine learning task of identifying the optimal sampling pattern in the Fourier domain for a given image and variational regularisation method, which has applications in the context of sparsity promoting reconstruction from magnetic resonance imaging data. [ABSTRACT FROM AUTHOR]

We introduce a novel technique to mitigate the adverse effects of atmospheric turbulence on astronomical imaging. Utilizing a video-to-image neural network trained on simulated data, our method processes a sliding sequence of short-exposure (∼0.2 s) stellar field images to reconstruct an image devoid of both turbulence and noise. We demonstrate the method with simulated and observed stellar fields, and show that the brief exposure sequence allows the network to accurately associate speckles to their originating stars and effectively disentangle light from adjacent sources across a range of seeing conditions, all while preserving flux to a lower signal-to-noise ratio than an average stack. This approach results in a marked improvement in angular resolution without compromising the astrometric stability of the final image. [ABSTRACT FROM AUTHOR]

The physical mechanisms controlling deep-focus earthquakes, or those observed at depths greater than 300 km, remain enigmatic. The leading processes by which deep-focus earthquakes are thought to occur include transformational faulting, thermal runaway and dehydration embrittlement, but distinguishing observations in support of one or more prevailing mechanisms are needed. In this study, we use a modified back-projection method, data recorded by the Hi-net array in Japan and a 3-D velocity model to produce source images of 19 deep-focus earthquakes within the Izu-Bonin subduction zone. We find that the rupture properties and fault plane orientations of imaged events separate according to reported moment magnitude, indicating the distinct operation of two moment-dependent causal mechanisms of deep-focus earthquakes in this region. We discuss these results in the context of previous observational, laboratory and numerical studies and emphasize the importance of continued research to validate the dual-mechanism hypothesis both in and outside Izu-Bonin. Such work may not only improve our understanding of the nucleation and propagation of deep-focus earthquakes, but also help clarify slab structure and subduction zone dynamics. [ABSTRACT FROM AUTHOR]

We develop a new deconvolution method to recover the precise Crab Nebula image taken by the Hitomi HXT, suppressing the artifact due to the bright Crab pulsar. We extend the Richardson–Lucy method, introducing two components corresponding to the nebula and pulsar with regularization for smoothness and flux, respectively, and performing simultaneous deconvolution of multi-pulse-phase images. The structures, including the torus and jets, seen in the deconvolved nebula image at the lowest energy band of 3.6–15 keV appear consistent with those identified in the high-resolution Chandra X-ray image. Above 15 keV, we confirm NuSTAR's findings that the nebula size decreases in higher energy bands. We find that the north-east side of the nebula is fainter in higher energy bands. Our deconvolution method is applicable for any telescope images of faint diffuse objects containing a bright point source. [ABSTRACT FROM AUTHOR]

Muography is an imaging technique that relies on the attenuation of the muon flux traversing geological or anthropogenic structures. Several simulation frameworks help to perform muography studies by combining specialized codes: for muon generation through muon transport to muon detector performance. This methodology is precise but requires significant computational resources and time. We present an end-to-end python-based MUographY Simulation Code, which implements a muography simulation framework capable of rapidly estimating muograms of any geological structure worldwide. This framework considers the generated muon flux as the observation point; the energy loss of muons passing through the geological target; the integrated muon flux detected by the telescope and estimates the 3-D density distribution of the target using algebraic reconstruction techniques. The simulations ignore the relatively small muon flux variance caused by geomagnetic effects, solar modulation and atmospheric conditions. We validate the code performance by comparing our simulation results with data from other frameworks. [ABSTRACT FROM AUTHOR]

The upcoming Square Kilometre Array Observatory will produce images of neutral hydrogen distribution during the epoch of reionization by observing the corresponding 21-cm signal. However, the 21-cm signal will be subject to instrumental limitations such as noise and galactic foreground contamination that pose a challenge for accurate detection. In this study, we present the SegU-Net v2 framework, an enhanced version of our convolutional neural network, built to identify neutral and ionized regions in the 21-cm signal contaminated with foreground emission. We trained our neural network on 21-cm image data processed by a foreground removal method based on Principal Component Analysis achieving an average classification accuracy of 71 per cent between redshift z  = 7 and 11. We tested SegU-Net v2 against various foreground removal methods, including Gaussian Process Regression, Polynomial Fitting, and Foreground-Wedge Removal. Results show comparable performance, highlighting SegU-Net v2 's independence on these pre-processing methods. Statistical analysis shows that a perfect classification score with |${\rm AUC}=95~{{\ \rm per\ cent}}$| is possible for 8 < z < 10. While the network prediction lacks the ability to correctly identify ionized regions at higher redshift and differentiate well the few remaining neutral regions at lower redshift due to low contrast between 21-cm signal, noise, and foreground residual in images. Moreover, as the photon sources driving reionization are expected to be located inside ionized regions, we show that SegU-Net v2 can be used to correctly identify and measure the volume of isolated bubbles with |$V_{\rm ion}\gt (10\, {\rm cMpc})^3$| at z > 9, for follow-up studies with infrared/optical telescopes to detect these sources. [ABSTRACT FROM AUTHOR]

STUDY QUESTION Can generative artificial intelligence (AI) models produce high-fidelity images of human blastocysts? SUMMARY ANSWER Generative AI models exhibit the capability to generate high-fidelity human blastocyst images, thereby providing substantial training datasets crucial for the development of robust AI models. WHAT IS KNOWN ALREADY The integration of AI into IVF procedures holds the potential to enhance objectivity and automate embryo selection for transfer. However, the effectiveness of AI is limited by data scarcity and ethical concerns related to patient data privacy. Generative adversarial networks (GAN) have emerged as a promising approach to alleviate data limitations by generating synthetic data that closely approximate real images. STUDY DESIGN, SIZE, DURATION Blastocyst images were included as training data from a public dataset of time-lapse microscopy (TLM) videos (n = 136). A style-based GAN was fine-tuned as the generative model. PARTICIPANTS/MATERIALS, SETTING, METHODS We curated a total of 972 blastocyst images as training data, where frames were captured within the time window of 110–120 h post-insemination at 1-h intervals from TLM videos. We configured the style-based GAN model with data augmentation (AUG) and pretrained weights (Pretrained-T: with translation equivariance; Pretrained-R: with translation and rotation equivariance) to compare their optimization on image synthesis. We then applied quantitative metrics including Fréchet Inception Distance (FID) and Kernel Inception Distance (KID) to assess the quality and fidelity of the generated images. Subsequently, we evaluated qualitative performance by measuring the intelligence behavior of the model through the visual Turing test. To this end, 60 individuals with diverse backgrounds and expertise in clinical embryology and IVF evaluated the quality of synthetic embryo images. MAIN RESULTS AND THE ROLE OF CHANCE During the training process, we observed consistent improvement of image quality that was measured by FID and KID scores. Pretrained and AUG + Pretrained initiated with remarkably lower FID and KID values compared to both Baseline and AUG + Baseline models. Following 5000 training iterations, the AUG + Pretrained-R model showed the highest performance of the evaluated five configurations with FID and KID scores of 15.2 and 0.004, respectively. Subsequently, we carried out the visual Turing test, such that IVF embryologists, IVF laboratory technicians, and non-experts evaluated the synthetic blastocyst-stage embryo images and obtained similar performance in specificity with marginal differences in accuracy and sensitivity. LIMITATIONS, REASONS FOR CAUTION In this study, we primarily focused the training data on blastocyst images as IVF embryos are primarily assessed in blastocyst stage. However, generation of an array of images in different preimplantation stages offers further insights into the development of preimplantation embryos and IVF success. In addition, we resized training images to a resolution of 256 × 256 pixels to moderate the computational costs of training the style-based GAN models. Further research is needed to involve a more extensive and diverse dataset from the formation of the zygote to the blastocyst stage, e.g. video generation, and the use of improved image resolution to facilitate the development of comprehensive AI algorithms and to produce higher-quality images. WIDER IMPLICATIONS OF THE FINDINGS Generative AI models hold promising potential in generating high-fidelity human blastocyst images, which allows the development of robust AI models as it can provide sufficient training datasets while safeguarding patient data privacy. Additionally, this may help to produce sufficient embryo imaging training data with different (rare) abnormal features, such as embryonic arrest, tripolar cell division to avoid class imbalances and reach to even datasets. Thus, generative models may offer a compelling opportunity to transform embryo selection procedures and substantially enhance IVF outcomes. STUDY FUNDING/COMPETING INTEREST(S) This study was supported by a Horizon 2020 innovation grant (ERIN, grant no. EU952516) and a Horizon Europe grant (NESTOR, grant no. 101120075) of the European Commission to A.S. and M.Z.E. the Estonian Research Council (grant no. PRG1076) to A.S. and the EVA (Erfelijkheid Voortplanting & Aanleg) specialty program (grant no. KP111513) of Maastricht University Medical Centre (MUMC+) to M.Z.E. TRIAL REGISTRATION NUMBER Not applicable. [ABSTRACT FROM AUTHOR]

Spectral induced polarization (SIP) and electrical impedance tomography (EIT) are advanced approaches for structural characterisation of soils and rocks based on measurements of the frequency-dependent complex electrical resistivity. Accurate SIP and EIT measurements require the reduction of electrode effects, which have often been attributed to the polarization of metal electrodes in the current pathway. In this study, we demonstrate an alternative electrode effect that influences the accuracy of complex impedance measurements. This alternative electrode effect is caused by a heterogeneous distribution of the contact impedance, which can occur due to ageing or corrosion of the electrode. First, theoretical considerations are given to explain how this electrode effect influences the accuracy of complex impedance measurements. Next, we demonstrate using SIP and EIT measurements that this alternative mechanism better explains the observed electrode effects than previous explanations based only on the presence of metal. The SIP measurements were made on a water-filled container with a metal sphere between two non-polarizing electrodes, which showed that the mere presence of metal in the current pathway only leads to small electrode effects that are typically well below 1 mrad even for small electrode spacings. When the same metal sphere was used as a potential electrode, the observed electrode effect was about 2.5 times larger and varied with rotation. These changes can be explained by a heterogeneous contact impedance resulting in a varying electrode voltage. The EIT measurements were performed with a cylindrical sample holder filled with water. They are used to illustrate how the presented theoretical considerations can be used to investigate the magnitude of the electrode effects for different electrode configurations. We conclude that electrode effects are important to consider for laboratory SIP and EIT measurements and field measurement with closely spaced electrodes. [ABSTRACT FROM AUTHOR]

Seismic interrogation of the upper mantle from the base of the crust to the top of the mantle transition zone has revealed discontinuities that are variable in space, depth, lateral extent, amplitude and lack a unified explanation for their origin. Improved constraints on the detectability and properties of mantle discontinuities can be obtained with P-to-S receiver function (Ps-RF) where energy scatters from P to S as seismic waves propagate across discontinuities of interest. However, due to the interference of crustal multiples, uppermost mantle discontinuities are more commonly imaged with lower resolution S-to-P receiver function (Sp-RF). In this study, a new method called CRISP-RF (Clean Receiver-function Imaging using SParse Radon Filters) is proposed, which incorporates ideas from compressive sensing and model-based image reconstruction. The central idea involves applying a sparse Radon transform to effectively decompose the Ps-RF into its underlying wavefield contributions, that is direct conversions, multiples, and noise, based on the phase moveout and coherence. A masking filter is then designed and applied to create a multiple-free and denoised Ps-RF. We demonstrate, using synthetic experiment, that our implementation of the Radon transform using a sparsity-promoting regularization outperforms the conventional least-squares methods and can effectively isolate direct Ps conversions. We further apply the CRISP-RF workflow on real data, including single station data on cratons, common-conversion-point stack at continental margins and seismic data from ocean islands. The application of CRISP-RF to global data sets will advance our understanding of the enigmatic origins of the upper mantle discontinuities like the ubiquitous mid-lithospheric discontinuity and the elusive X-discontinuity. [ABSTRACT FROM AUTHOR]

The treatment of recurrent ovarian cancer has been based on systemic therapy. The role of secondary cytoreductive surgery has been addressed recently in several trials. Imaging plays a key role in helping the surgical team to decide which patients will have resectable disease and benefit from surgery. The role of staging laparoscopy and several imaging and clinical scores has been extensively debated in the field. In other surgical fields there have been reports of using 3D imaging software and 3D printed models to help surgeons better plan the surgical approach. To the best of our knowledge, we report the first case of a patient with recurrent ovarian cancer undergoing 3D modeling before secondary cytoreductive surgery. The 3D modeling was of most value to evaluate the extension of the disease in our patient who underwent a successful secondary cytoreductive surgery and is currently free of the disease. [ABSTRACT FROM AUTHOR]