Search

When interpreting marine very high-resolution (VHR) single-channel seismic reflection data, the signal in the water column is generally considered as noise and is often eliminated by a water-mute application to focus on geological information under the seafloor. Alternatively, the signal in the water column can be used to study ocean currents or gas/fluid emissions. To provide images of the sedimentary formations and tectonic structures beneath the seafloor in shallow water regions, such as continental shelves and lakes, marine seismic reflection profiles are often acquired using a single-channel streamer and sparker-type source, providing VHR data, with limited penetration depth. To exploit the full potential of these single-channel data, we propose a simple algorithm, called REWARE (Recovery of Water-column Acoustic Reflectors). This algorithm allows to extract further geological information from the water-column data using open-source codes (Seismic Un*x), adding the coherent signal from the previous shots, recorded in the water column, to the previous traces. The record length becomes longer while maintaining a very high trace-to-trace consistency. To demonstrate its efficiency, we present two examples of the REWARE processing in two different geological contexts: the East Sardinia shelf (Italy) and the North Evia Gulf (Greece). This method provides deeper images than with original data for seismic data acquired across steep slopes, such as canyons or continental shelf breaks. Thus, depending on the seafloor geometry and subseafloor structures, it is possible to image or map sediment layers and tectonic structures at depth, keeping a very high structural resolution. [ABSTRACT FROM AUTHOR]

We show that convolution neural networks (CNNs) trained to find strong gravitational lens systems are biased towards systems with larger Einstein radii and large concentrated sources. This selection function is key to fully realizing the potential of the large samples of strong gravitational lens systems that will be found in upcoming wide-field surveys. In this paper, we use a CNN and three training data sets to quantify the network selection function and its implication for the many scientific applications of strong gravitational lensing. We use CNNs with similar architecture as is commonly found in the literature. The networks preferentially select systems with larger Einstein radii and larger sources with more concentrated source-light distributions. Increasing the detection significance threshold to 12 |$\sigma$| from 8 |$\sigma$| results in 50 per cent of the selected strong lens systems having Einstein radii |$\theta _\mathrm{E}$| |$\ge$| 1.04 arcsec from |$\theta _\mathrm{E}$| |$\ge$| 0.879 arcsec, source radii |$R_S$| |$\ge$| 0.194 arcsec from |$R_S$| |$\ge$| 0.178 arcsec, and source Sérsic indices |$n_{\mathrm{Sc}}^{\mathrm{S}}$| |$\ge$| 2.62 from |$n_{\mathrm{Sc}}^{\mathrm{S}}$| |$\ge$| 2.55. The model trained to find lensed quasars shows a stronger preference for higher lens ellipticities than those trained to find lensed galaxies. The selection function is independent of the slope of the power law of the mass profiles, hence measurements of this quantity will be unaffected. The lens finder selection function reinforces that of the lensing cross-section, and thus we expect our findings to be a general result for all galaxy–galaxy and galaxy–quasar lens finding neural networks. [ABSTRACT FROM AUTHOR]

In large-scale astronomical surveys, traditional supernova detection pipelines rely on complex and relatively inefficient image differencing techniques. This paper proposes an end-to-end deep-learning supernova detection network, the Real-Time SuperNova DEtection TRansformer (RT-SNDETR). This network partially replaces traditional pipelines by integrating image differencing, source detection, and Real-bogus classification, achieving a speed 51.49 times that of the fastest image differencing method, SFFT. Additionally, it remains competitive with methods like YOLO v8, offering a well-balanced trade-off between speed and accuracy. Experimental results highlight RT-SNDETR's superior performance, with an average precision(AP) of 96.30 per cent on synthetic samples and 76.60 per cent on real supernova samples. It significantly outperforms other detection networks, including RT-DETR (+5.6 per cent AP on synthetic/+5.1 per cent AP on real samples) and Cascade R-CNN (+8.9 per cent AP on synthetic/ +28.6 per cent AP on real samples). The incorporation of CycleGAN-based data generation methods plays a significant role in enhancing RT-SNDETR's performance. These methods simulate realistic PSF variations, enabling the object detection network to learn more robust features and improving its generalization to real supernovae data. Additionally, by integrating unsupervised domain adaptation techniques, RT-SNDETR achieves an AP of 81.70 per cent on real SDSS supernova survey samples. This study demonstrates RT-SNDETR's potential to significantly enhance both the speed and accuracy of supernova detection, making it a highly effective solution for large-scale astronomical surveys. [ABSTRACT FROM AUTHOR]

Optically observing and monitoring moving objects, both natural and artificial, is important to human space security. Non-sidereal tracking can improve the system's limiting magnitude for moving objects, which benefits the surveillance. However, images with non-sidereal tracking include complex background, as well as objects with different brightness and moving mode, posing a significant challenge for accurate multi-object detection in such images, especially in wide field-of-view telescope images. To achieve a higher detection precision in a higher speed, we proposed a novel object detection method, which combines the source feature extraction and the neural network. First, our method extracts object features from optical images such as centroid, shape, and flux. Then, it conducts a naive labelling based on those features to distinguish moving objects from stars. After balancing the labelled data, we employ it to train a neural network aimed at creating a classification model for point-like and streak-like objects. Ultimately, based on the neural network model's classification outcomes, moving objects whose motion modes consistent with the tracked objects are detected via track association, while objects with different motion modes are detected using morphological statistics. The validation, based on the space objects images captured in target tracking mode with the 1-m telescope at Nanshan, Xinjiang Astronomical Observatory, demonstrates that our method achieves 94.72 per cent detection accuracy with merely 5.02 per cent false alarm rate, and a processing time of 0.66 s per frame. Consequently, our method can rapidly and accurately detect objects with different motion modes from wide-field images with non-sidereal tracking. [ABSTRACT FROM AUTHOR]

Fault rupture dynamics is expected to be significantly affected by the geometry of fault system, especially for orthogonal faults. However, the rupture behaviours of orthogonal faults, especially the coseismic interactions, are far from fully understood. Here, we present experimental results from a series of laboratory earthquakes to elucidate the effect of the stress state and initiation location on the rupture behaviours of orthogonal faults. Our results reveal a phase diagram of rupture behaviours of orthogonal faults, which is collectively controlled by stress state and rupture initiation location. For events initiating from the main fault, the rupture cannot jump to the branch, which may be due to the clamping effect or the inhibited shear stress accumulation on the branch. On the contrary, events initiating from the branch can persistently trigger ruptures of the main fault. This difference highlights the directional effect associated with the rupture of orthogonal faults. Further, the rupture length of triggered ruptures on the main fault is controlled by the stress state of the fault system. With the increase of the ratio between the shear stress and normal stress, the rupture length of the main fault increases. Our results reproduce the rupture behaviours of orthogonal faults, which may provide insights into the rupture characteristics of natural earthquakes. [ABSTRACT FROM AUTHOR]

We present the joint astrometric and direct imaging discovery, mass measurement, and orbital analysis of HD 63754 B (HIP 38216 B), a companion near the stellar-substellar boundary orbiting |$\sim$| 20 au from its Sun-like host. HD 63754 was observed in our ongoing high-contrast imaging survey targeting stars with significant proper-motion accelerations between Hipparcos and Gaia consistent with wide-separation substellar companions. We utilized archival High Resolution Echelle Spectrograph and High Accuracy Radial Velocity Planet Searcher radial velocity (RV) data, together with the host star's astrometric acceleration extracted from the Hipparcos – Gaia Catalog of Accelerations, to predict the location of the candidate companion around HD 63754 A. We subsequently imaged HD 63754 B at its predicted location using the Near Infrared Camera 2 (NIRC2) in the |$L^{\prime }$| band at the W. M. Keck Observatory. We then jointly modelled the orbit of HD 63754 B with RVs, Hipparcos – Gaia accelerations, and our new relative astrometry, measuring a dynamical mass of |${81.9}_{-5.8}^{+6.4} {M_{\rm Jup}}$|⁠ , an eccentricity of |${0.260}_{-0.059}^{+0.065}$|⁠ , and a nearly face-on inclination of |$174.\!\!^\circ 81_{-0.50}^{+0.48}$|⁠. For HD 63754 B, we obtain an |$L^{\prime }$| -band absolute magnitude of |$L^{\prime } = 11.39\pm 0.06$| mag, from which we infer a bolometric luminosity of |$\rm{log(L_{\rm bol}/{\rm{L}_{\odot}})= -4.55 \pm 0.08}$| dex using a comparison sample of L and T dwarfs with measured luminosities. Although uncertainties linger in age and dynamical mass estimates, our analysis points towards HD 63754 B's identity as a brown dwarf on the L/T transition rather than a low-mass star, indicated by its inferred bolometric luminosity and model-estimated effective temperature. Future RV, spectroscopic, and astrometric data such as those from JWST and Gaia Data Release 4 will clarify HD 63754 B's mass, and enable spectral typing and atmospheric characterization. [ABSTRACT FROM AUTHOR]

Surgical robotics application in the field of minimally invasive surgery has developed rapidly and has been attracting increasingly more research attention in recent years. A common consensus has been reached that surgical procedures are to become less traumatic and with the implementation of more intelligence and higher autonomy, which is a serious challenge faced by the environmental sensing capabilities of robotic systems. One of the main sources of environmental information for robots are images, which are the basis of robot vision. In this review article, we divide clinical image into direct and indirect based on the object of information acquisition, and into continuous, intermittent continuous, and discontinuous according to the target-tracking frequency. The characteristics and applications of the existing surgical robots in each category are introduced based on these two dimensions. Our purpose in conducting this review was to analyze, summarize, and discuss the current evidence on the general rules on the application of image technologies for medical purposes. Our analysis gives insight and provides guidance conducive to the development of more advanced surgical robotics systems in the future. [ABSTRACT FROM AUTHOR]

Advancements in imaging technologies have revolutionized our ability to deeply profile pathological tissue architectures, generating large volumes of imaging data with unparalleled spatial resolution. This type of data collection, namely, spatial proteomics, offers invaluable insights into various human diseases. Simultaneously, computational algorithms have evolved to manage the increasing dimensionality of spatial proteomics inherent in this progress. Numerous imaging-based computational frameworks, such as computational pathology, have been proposed for research and clinical applications. However, the development of these fields demands diverse domain expertise, creating barriers to their integration and further application. This review seeks to bridge this divide by presenting a comprehensive guideline. We consolidate prevailing computational methods and outline a roadmap from image processing to data-driven, statistics-informed biomarker discovery. Additionally, we explore future perspectives as the field moves toward interfacing with other quantitative domains, holding significant promise for precision care in immuno-oncology. [ABSTRACT FROM AUTHOR]

SUMMARY: The resolution of velocity models obtained by tomography varies due to multiple factors and variables, such as the inversion approach, ray coverage, data quality, etc. Combining velocity models with different resolutions can enable more accurate ground motion simulations. Toward this goal, we present a novel methodology to fuse multiresolution seismic velocity maps with probabilistic graphical models (PGMs). The PGMs provide segmentation results, corresponding to various velocity intervals, in seismic velocity models with different resolutions. Further, by considering physical information (such as ray path density), we introduce physics-informed probabilistic graphical models (PIPGMs). These models provide data-driven relations between subdomains with low (LR) and high (HR) resolutions. Transferring (segmented) distribution information from the HR regions enhances the details in the LR regions by solving a maximum likelihood problem with prior knowledge from HR models. When updating areas bordering HR and LR regions, a patch-scanning policy is adopted to consider local patterns and avoid sharp boundaries. To evaluate the efficacy of the proposed PGM fusion method, we tested the fusion approach on both a synthetic checkerboard model and a fault zone structure imaged from the 2019 Ridgecrest, CA, earthquake sequence. The Ridgecrest fault zone image consists of a shallow (top 1 km) high-resolution shear-wave velocity model obtained from ambient noise tomography, which is embedded into the coarser Statewide California Earthquake Center Community Velocity Model version S4.26-M01. The model efficacy is underscored by the deviation between observed and calculated traveltimes along the boundaries between HR and LR regions, 38 per cent less than obtained by conventional Gaussian interpolation. The proposed PGM fusion method can merge any gridded multiresolution velocity model, a valuable tool for computational seismology and ground motion estimation. [ABSTRACT FROM AUTHOR]

In a strong gravitational lensing system, the distorted light from a source is analysed to infer the properties of the lens. However, light emitted by the lens itself can contaminate the image of the source, introducing systematic errors in the analysis. We present a simple and efficient lens light model based on the well-tested multi-Gaussian expansion (MGE) method for representing galaxy surface brightness profiles, which we combine with a semi-linear inversion scheme for pixelized source modelling. Testing it against realistic mock lensing images, we show that our scheme can fit the lensed images to the noise level, with relative differences between the true input and best-fitting lens light model remaining below 5 per cent. We apply the MGE lens light model to 38 lenses from the HST SLACS sample. We find that the new scheme provides a good fit for the majority of the sample with only 3 exceptions – these show clear asymmetric residuals in the lens light. We examine the radial dependence of the ellipticity and position angles and confirm that it is common for a typical lens galaxy to exhibit twisting non-elliptical isophotes and boxy / disky isophotes. Our MGE lens light model will be a valuable tool for understanding the hidden complexity of the lens mass distribution. [ABSTRACT FROM AUTHOR]

SUMMARY: To address complex subsurface structures, a high-resolution velocity model must be constructed. Conventionally, algorithms such as full waveform inversion (FWI) have been used to derive accurate high-resolution velocity structures, but obstacles such as high computational costs remain. Therefore, we propose a high-resolution U-NET (HR U-NET) machine learning model to derive a high-resolution velocity model from a low-resolution velocity model. The low-resolution velocity model and migration data obtained through the corresponding velocity information were used as input data for training. In addition, we tried to improve the accuracy of the high-resolution velocity model by using prior information containing accurate velocity values. A prior model generated through geophysical logging data and a weight model including the reliability information of the prior model were also utilized. Therefore, the HR U-NET model was trained using the low-resolution velocity model, the migration data, the prior model and the weight model. Numerical experiments conducted using synthetic and field data demonstrated that the proposed model could accurately construct a high-resolution velocity model and verified that the prior model and weight model play important roles in the training process. Additionally, we confirmed that the proposed method derived almost similar results using only 8.2 percent of the computational cost of the conventional inversion method. In other words, there is an advantage that it is possible to predict high-resolution velocity information more efficiently in terms of computational cost. [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]

Summary perox-per-cell automates cumbersome, image-based data collection tasks often encountered in peroxisome research. The software processes microscopy images to quantify peroxisome features in yeast cells. It uses off-the-shelf image processing tools to automatically segment cells and peroxisomes and then outputs quantitative metrics including peroxisome counts per cell and spatial areas. In validation tests, we found that perox-per-cell output agrees well with manually quantified peroxisomal counts and cell instances, thereby enabling high-throughput quantification of peroxisomal characteristics. Availability and implementation The software is coded in Python. Compiled executables and source code are available at https://github.com/AitchisonLab/perox-per-cell. [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]

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 growing demand for high-quality industrial products has led to a significant emphasis on image anomaly detection (AD). AD in industrial goods presents a formidable research challenge that demands the application of sophisticated techniques to identify and address deviations from the expected norm accurately. Manufacturers increasingly recognize the significance of employing intelligent systems to detect flaws and defects in product parts. However, industrial settings pose several challenges: diverse categories, limited abnormal samples and vagueness. Hence, there is a growing demand for advanced image AD techniques within industrial product manufacturing. In this paper, an intelligent industrial defective chips detection framework is proposed which mainly consists of three core components. First, the convolutional features of the efficient backbone model is effectively utilized to balance the computational complexity and performance of industrial resource-constrained devices. Secondly, a novel inverse feature matching followed by masking method is proposed to enhance the explanability that localizes the abnormal regions of the abnormal chips. Finally, to evaluate our proposed method a comprehensive ablation study is conducted, where different machine learning and deep learning algorithms are analysed to claim the superiority of our method. Furthermore, to help the research community, a benchmark dataset is collected from real-world industry manufacturing for defective chip detection. The empirical results from the dataset demonstrate the strength and effectiveness of the proposed model compared to the other models. [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]

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]

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]

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 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]

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]

The main advantage of full-waveform inversion (FWI) is the ability to obtain useful subsurface structure information, such as velocity and density, from complex seismic data. We have developed a novel inversion algorithm to improve the capability of FWI to achieve high-resolution imaging, even under complex conditions caused by random noise contamination, initial model dependence, or the selection of parameters to be estimated. Our algorithm considers an effective image processing and dimension reduction tool, randomized singular value decomposition-weighted truncated nuclear norm regularization, for embedding FWI to achieve high-resolution imaging results. This algorithm obtains a truncated matrix approximating the original matrix by reducing the rank of the velocity increment matrix, thus achieving the truncation of noisy data, with the truncation range controlled by weighted truncated nuclear norm regularization. Subsequently, we employ an inexact augmented Lagrangian method algorithm in the optimization to compress the solution space range, thus relaxing the dependence of FWI and randomized singular value decomposition-weighted truncated nuclear norm regularization on the initial model and accelerating the convergence rate of the objective function. We tested on one set of synthetic data, and the results show that, compared with traditional FWI, our method can more effectively suppress the impact of random noise, thus obtaining higher resolution and more accurate subsurface model information. This work indicates that the combination of randomized singular value decomposition-weighted truncated nuclear norm regularization and FWI is an effective imaging strategy that can help to solve the challenges faced by traditional FWI. [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 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]

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]

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]

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]

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]

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]

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]

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]

We studied the benefits of a nine-component (9C) seismic survey over a landslide-prone area in southwest of Sweden to retrieve ultrahigh-resolution shear wave reflection images of the subsurface as well as crucial information on physical properties of the sediments. A complete, 1 m shot and receiver spacing, multicomponent 2-D seismic profile was acquired using three-component microelectromechanical-system-based landstreamer receivers, and a 5-kg sledgehammer strike in vertical and horizontal orientations as three-component seismic source. Given the rich number of shear wave reflections observed on all the 9C data, the processing work focused on their retrievals. It revealed three distinct reflections, two of which associated with coarse-grained materials and one with an extremely undulating bedrock surface. Given the extremely slow shear wave velocities on the order of 60–100 m s−1, we obtained ultrahigh-resolution shear wave sections avoiding temporal and spatial aliasing. Imaging results suggest vertical-source and horizontal-radial receiver (V–HR), and horizontal-transverse source–receiver orientations (HT–HT) provided the most optimum images of the subsurface. A non-hyperbolic algorithm was applied to the normal-moveout corrections justified by the traveltime differences of the bedrock reflection in different shear wave sections. The improved images by incorporating the anisotropy term suggest that the data set reveals moderate shear wave anisotropy along some portions of the profile. The Vp / Vs ratio obtained by using bedrock reflection in P - and S -wave sections suggests values ranging 10–16, which implies high water content. Areas with lower Vp / Vs coincides with greater anisotropic parameters and this can indicate disturbed clays or presence of sensitive clays. [ABSTRACT FROM AUTHOR]

Seismic interferometry (SI) is a technique that allows one to estimate the wavefields accounting for the wave propagation between seismometers, any of which can act as a virtual source (VS). Interferometry, particularly noise interferometry, has been applied to several geophysical disciplines such as passive monitoring and distributed acoustic sensing. In practice, one requires long recordings of seismic noise for noise interferometry. Additionally, one can have missing seismic interferometric traces because some receivers in seismic arrays may be absent or inoperative due to issues of receiver installation and malfunction. Thus, filling the gap of seismic interferometric profile requires wavefield reconstruction and regularization techniques. Compressive sensing (CS) is one such method that can reconstruct seismic interferometric wavefields and help mitigate the limitations by exploiting the sparsity of seismic waves. In our work, we use CS to reconstruct missing seismic interferometric wavefields. One can interpolate interferometric wavefields using correlograms provided by one VS. We call this method of reconstructing an individual VS gather single-source wavefield reconstruction. We propose an alternative technique called multi-source wavefield reconstruction , which applies CS to reconstruct multiple interferometric wavefields using a volume of VS gathers provided from all available VSs. Using numerical examples, we show that one can apply CS to recover interferometric wavefields resulting from interferometry of a linear seismic array. To exploit the sparsity of interferometric wavefields, we apply the Fourier and Curvelet transforms to the two reconstruction schemes. Using the signal-to-noise ratio (SNR) to compare reconstruction of interferometric wavefields, the Fourier multi-source method improves the recovery of interferometric wavefields by approximately 50 dB compared to the Fourier and Curvelet single-source wavefield reconstructions. [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]

We present a new method, called 'forced-spectrum fitting', for physically based spectral modelling of radio sources during deconvolution. This improves upon current common deconvolution fitting methods, which often produce inaccurate spectra. Our method uses any pre-existing spectral index map to assign spectral indices to each model component cleaned during the multifrequency deconvolution of wsclean , where the pre-determined spectrum is fitted. The component magnitude is evaluated by performing a modified weighted linear least-squares fit. We test this method on a simulated LOFAR high-band antenna (HBA) observation of the 3C 196 QSO and a real LOFAR HBA observation of the 4C+55.16 FRI galaxy. We compare the results from the forced-spectrum fitting with traditional joined-channel deconvolution using polynomial fitting. Because no prior spectral information was available for 4C+55.16, we demonstrate a method for extracting spectral indices in the observed frequency band using 'clustering'. The models generated by the forced-spectrum fitting are used to improve the calibration of the data sets. The final residuals are comparable to existing multifrequency deconvolution methods, but the output model agrees with the provided spectral index map, embedding correct spectral information. While forced-spectrum fitting does not solve the determination of the spectral information itself, it enables the construction of accurate multifrequency models that can be used for wide-band calibration and subtraction. [ABSTRACT FROM AUTHOR]

One of the fundamental questions in plant developmental biology is how cell proliferation and cell expansion coordinately determine organ growth and morphology. An amenable system to address this question is the Arabidopsis root tip, where cell proliferation and elongation occur in spatially separated domains, and cell morphologies can easily be observed using a confocal microscope. While past studies revealed numerous elements of root growth regulation including gene regulatory networks, hormone transport and signaling, cell mechanics and environmental perception, how cells divide and elongate under possible constraints from cell lineages and neighboring cell files has not been analyzed quantitatively. This is mainly due to the technical difficulties in capturing cell division and elongation dynamics at the tip of growing roots, as well as an extremely labor-intensive task of tracing the lineages of frequently dividing cells. Here, we developed a motion-tracking confocal microscope and an Artificial Intelligence (AI)-assisted image-processing pipeline that enables semi-automated quantification of cell division and elongation dynamics at the tip of vertically growing Arabidopsis roots. We also implemented a data sonification tool that facilitates human recognition of cell division synchrony. Using these tools, we revealed previously unnoted lineage-constrained dynamics of cell division and elongation, and their contribution to the root zonation boundaries. [ABSTRACT FROM AUTHOR]

Digital rock physics (DRP) is important for characterizing the pore characteristics and petrophysical parameters of rocks from a microscopic perspective. Among the digital rock reconstruction methods, the most widely used is the computer tomography (CT) scan method. However, the pore structure of rocks has multiscale features, and CT scan images with a single resolution cannot completely describe the pore structure characteristics of rocks. In this paper, we analysed reconstructed multisource and multiscale 3-D digital rocks based on the cycle-consistent generative adversarial network (CycleGAN) method. This method fully integrates the advantages of the large field of view of low-resolution images and the high-precision features of high-resolution images. To investigate the performance of the method, two sequences of CT scan images of a sandstone (sample A) and a carbonate rock (sample B) collected from oilfields were studied. Moreover, due to the high resolution of scanning electron microscope (SEM) images, we also constructed 3-D digital rocks with different resolutions of the Berea sandstone based on the combination of CT scan images and SEM images. Finally, the statistical properties and absolute permeabilities were calculated to evaluate the accuracies of reconstructed multisource and multiscale 3-D digital rocks. The results show that the reconstructed multiscale digital rocks based on the CycleGAN method have good accuracy in terms of statistical properties and petrophysical properties. Considering the computational cost and computational accuracy, for high- and low-resolution CT scan images and CT scan images with SEM images as training image data sources, we suggest that the resolution of the reconstructed multiscale digital cores is 4–8 times and 4–16 times higher than that of the low-resolution CT scan images, respectively. The findings of our research will be helpful in gaining insight into the petrophysical properties of heterogeneous rocks. [ABSTRACT FROM AUTHOR]

Full-waveform inversion (FWI) is the current standard method to determine final and detailed model parameters to be used in the seismic imaging process. However, FWI is an ill-posed problem that easily achieves a local minimum, leading the model solution in the wrong direction. Recently, some works proposed integrating FWI with Convolutional Neural Networks (CNN). In this case, the CNN weights are updated following the FWI gradient, defining the process as a Physics-Informed Neural Network (PINN). FWI integrated with CNN has an important advantage. The CNN stabilizes the inversion, acting like a regularizer, avoiding local minima-related problems and sparing an initial velocity model in some cases. However, such a process, especially when not requiring an initial model, is computationally expensive due to the high number of iterations required until the convergence. In this work, we propose an approach which relies on combining supervised learning and physics-informed by using a previously trained CNN to start the DL-FWI inversion. Loading the pre-trained weights configures transfer learning. The pre-trained CNN is obtained using a supervised approach based on training with a reduced and simple data set to capture the main velocity trend at the initial FWI iterations. The proposed training process is different from the initial works on the area which obtained the velocity model from the shots in supervised learning tasks and that required a large amount of labelled data to ensure reasonable model predictions. We investigated in our approach two CNN architectures, obtaining more robust results and a reduced number of parameters when using a modified U-Net. The method was probed over three benchmark models, showing consistently that the pre-training phase reduces the process's uncertainties and accelerates the model convergence using minimal prior information. Besides, the final scores of the iterative process are better than the examples without transfer learning. Thus, transfer learning solved one main limitation of the previous PINN approaches: the unfeasible number of iterations when not using an initial model. Moreover, we tested the method using data with low-frequency band limitations, since the lack of low frequencies is a common issue within real seismic data. The inversion converges to reasonable results probing the method's robustness with restricted frequency content. [ABSTRACT FROM AUTHOR]

Accurate near-surface velocity models are necessary for land seismic imaging. First-arrival traveltime tomography (FTT) routinely used for estimating near-surface velocity models may fail in geological complex areas. Supervised deep learning (SDL) is capable of building accurate velocity models, based on tens of thousands of velocity model-shot gathers training pairs. It takes lots of time and memory space, which may be unaffordable for practical applications. We propose integrating the FTT and SDL to build near-surface velocity models. During the neural network training, the FTT-inverted models rather than the original seismic data are used as the network inputs and corresponding true models are the outputs. The FTT-inverted and true models are the same physical quantities and with the same dimensions. Their relationship is less non-linear than that between shot gathers and true models. Thus, the neural network of the proposed method can be trained well using only a small number of training samples, dramatically reducing the time and memory costs. Numerical tests demonstrate the feasibility and effectiveness of the proposed method. We applied the proposed method to a land data set obtained in mountainous areas in the west of China and obtained satisfactory near-surface velocity models and stacking images. [ABSTRACT FROM AUTHOR]

Seismic velocity analysis is the basis for seismic imaging and understanding complex subsurface geological structures. Although the performance of automatic velocity analysis methods based on Common Middle Point (CMP) data or Velocity Spectra (VS) is encouraging, particularly deep learning methods. However, most methods ignore the complementarity between CMP data and VS data, and only one of them is selected for velocity modelling. We propose a multimodal neural network (MMN) that combines the advantages of CMP data details representation and simplification of VS. MMN includes multilayer convolution structures and auto-encoder structures, which are used to extract time–space amplitude information from CMP gathers and energy groups features from VS data, respectively. This paper compared MMN with the CMP single-modal network (CSN) and the velocity spectra single-modal network (VSSN). Based on synthetic data, we investigated their differences in terms of continuity, accuracy, noise resistance and generalization. The MMN prediction results makes a trade-off between the overall continuity and local details. Visualization analysis of the intermediate feature maps explains the MMN velocity prediction mechanism, that is, the multi-angle representation and complementary fusion of velocity information. Finally, the performance of the proposed method is demonstrated using the braided river deposited field data example. [ABSTRACT FROM AUTHOR]

PAUCam is an innovative optical narrow-band imager mounted at the William Herschel Telescope built for the Physics of the Accelerating Universe Survey (PAUS). Its set of 40 filters results in images that are complex to calibrate, with specific instrumental signatures that cannot be processed with traditional data reduction techniques. In this paper, we present two pipelines developed by the PAUS data management team with the objective of producing science-ready catalogues from the uncalibrated raw images. The Nightly pipeline takes care of entire image processing, with bespoke algorithms for photometric calibration and scatter-light correction. The Multi-Epoch and Multi-Band Analysis pipeline performs forced photometry over a reference catalogue to optimize the photometric redshift (photo-z) performance. We verify against spectroscopic observations that the current approach delivers an inter-band photometric calibration of 0.8 per cent across the 40 narrow-band set. The large volume of data produced every night and the rapid survey strategy feedback constraints require operating both pipelines in the Port d'Informació Cientifica data centre with intense parallelization. While alternative algorithms for further improvements in photo-z performance are under investigation, the image calibration and photometry presented in this work already enable state-of-the-art photo-z down to i AB = 23.0. [ABSTRACT FROM AUTHOR]

Motivation Cell membrane segmentation in electron microscopy (EM) images is a crucial step in EM image processing. However, while popular approaches have achieved performance comparable to that of humans on low-resolution EM datasets, they have shown limited success when applied to high-resolution EM datasets. The human visual system, on the other hand, displays consistently excellent performance on both low and high resolutions. To better understand this limitation, we conducted eye movement and perceptual consistency experiments. Our data showed that human observers are more sensitive to the structure of the membrane while tolerating misalignment, contrary to commonly used evaluation criteria. Additionally, our results indicated that the human visual system processes images in both global–local and coarse-to-fine manners. Results Based on these observations, we propose a computational framework for membrane segmentation that incorporates these characteristics of human perception. This framework includes a novel evaluation metric, the perceptual Hausdorff distance (PHD), and an end-to-end network called the PHD-guided segmentation network (PS-Net) that is trained using adaptively tuned PHD loss functions and a multiscale architecture. Our subjective experiments showed that the PHD metric is more consistent with human perception than other criteria, and our proposed PS-Net outperformed state-of-the-art methods on both low- and high-resolution EM image datasets as well as other natural image datasets. Availability and implementation The code and dataset can be found at https://github.com/EmmaSRH/PS-Net. [ABSTRACT FROM AUTHOR]

SUMMARY: With the increasing advances in the oil and gas industry, seismic imaging near or under salt structures has become an important point in deep-water exploration. Detailed velocity models of these areas are particularly interesting not only to characterize hydrocarbon reservoirs but also to identify potential sites for hydrogen and carbon dioxide storage in offshore salt caverns. Thus, we study the full-waveform inversion for the salt reconstruction in acoustic media with constant density considering the time-harmonic wave propagation in a finite element formulation using the topology optimization (TO) method. This problem is challenging due to the strong velocity contrast between salt bodies and the sedimentary background, in addition to the lack of low-frequency data and the inherent ill-posedness of the inverse problem. In this context, we incorporate techniques from the TO field, usually used in design applications, to overcome or reduce these known problems. We initially defined the squared slowness as a combination of two fields, one related to the salt shape and the other to the background. An interpolation rule based on the solid isotropic material with penalization method, combined with filtering and projection schemes, is used to find the shape of the salt bodies with increased sharpness interfaces. A Helmholtz-type filter is applied to modify the gradient aimed to regularize the problem and provide a more stable way for the salt shape to evolve during the inversion process. In particular, we demonstrate that the proposed approach may be relevant for reconstructing media with salt bodies when a suitable starting model is unavailable, and sharp interfaces are required. In addition, we present inversion results from synthetic data generated by a variable density model to demonstrate the approach capability when subjected to a reconstruction application. [ABSTRACT FROM AUTHOR]

SUMMARY: Large continental earthquakes are infrequent, but they are important for understanding active tectonics. Many pre-1990, magnitude greater than 7 earthquakes have not been well documented due to a lack of data. The declassification of KeyHole (KH) imagery provides opportunities for investigating those historical events. In this study, we present new geodetic observations of surface deformation from the 1976 November 24 Ms 7.3 Chaldiran earthquake in eastern Turkey. We use various historical KH imagery and modern Pleiades stereo imagery to determine the coseismic deformation, providing new constraints on the 1976 rupture. Based on the image measurements, we resolve the coseismic slip at depth. The inverted slip from a simple elastic dislocation model, ∼4.0 m at shallow depth of the Hidirmentes segment, is larger than the surface displacement, around 3 m from previous field surveys, suggesting that the coseismic deformation is a combination of slip on the primary fault with 30 ± 16 per cent off-fault deformation. We also measure cumulative offsets at two sites with well-preserved features, and find a constant ratio of horizontal to vertical offsets throughout the past few earthquake cycles. Both the observed constant slip ratio in this study and offset clusters from previous field measurements support that the Chaldiran fault exhibits a characteristic slip behaviour in the late Quaternary. By combining all the geodetic and geological observations in the broader region, we also present a brief analysis of the slip rate of the Chaldiran fault. [ABSTRACT FROM AUTHOR]

Sewing defect detection is an essential step in garment production quality control. Although sewing defects significantly influence the quality of clothing, they are yet to be studied widely compared to fabric defects. In this study, to address sewing defect detection and develop an appropriate method for small and labor-intensive garment companies, an on-machine broken stitch detection system is proposed. In hardware, a versatile mounting kit, including clamping, display, and adjustable linkage for a camera, is presented for easy installation on a typical industrial sewing machine and for placing the camera close to the sewing position. Additionally, a prototype is implemented using a low-cost single-board computer, Raspberry Pi 4 B, its camera, and Python language. For automated broken stitch detection, a method is proposed that includes removing the texture of the background fabric, image processing in the HSV color space, and edge detection for robust broken detection under various fabric and thread colors and lighting conditions. The proposed system demonstrates reasonable real-time detection accuracy. The maximum accuracy obtained on a sewing stitch dataset with 880 images and on-site tests of various industrial sewing machines is 82.5%, which is 12.1–34.6% higher than that of the two existing methods. [ABSTRACT FROM AUTHOR]

Reflector-normal angles and reflector-curvature parameters are the principal geometric attributes used in seismic interpretation for characterizing the orientations and shapes, respectively, of geological reflecting surfaces. Commonly, the input data set for their computation consists of fine 3-D grids of scalar fields representing either the seismic-driven reflectivities (e.g. amplitudes of 3-D seismic migrated volumes) or model-driven reflectivities, computed, for example, from the derived elastic impedance parameters. Conventionally, the computation of curvature parameters at each gridpoint is based on analysing the local change in the inline/crossline dips, considering the potential existence of a local quadratic reflecting surface in the vicinity of that point. This assumption breaks down for subsurface points in the vicinity of either complex reflecting surfaces (e.g. brittle/rough/tilted synclines/anticlines, ridges/troughs and saddles) and/or sharp, discontinuous geological features (e.g. fault edges/tips, pinch-outs, fracture systems, channels and small geobodies), where the values of the computed curvature become extremely high. However, while these high values can indicate the existence of non-reflecting objects, they do not deliver their specific geometric characteristics. In this study, we present a novel method that better characterizes the shapes of these complex geological features by extending the assumption of local surfaces (2-D surfaces in 3-D space) into local hypersurfaces (3-D hypersurfaces in 4-D space), with their corresponding (three rather than two) principal (and effective) curvature parameters. We demonstrate the advantages of our method by comparing the conventional dip-based surface curvature parameters with the hypersurface curvature parameters, using a synthetic model/image with different types and shapes of geological features and a seismic image of real data containing a complex fault and hidden buried channels. [ABSTRACT FROM AUTHOR]