Enhanced resolution and sensitivity acoustic-resolution photoacoustic microscopy with semi/unsupervised GANs

Acoustic-resolution photoacoustic microscopy (AR-PAM) enables visualization of biological tissues at depths of several millimeters with superior optical absorption contrast. However, the lateral resolution and sensitivity of AR-PAM are generally lower than those of optical-resolution PAM (OR-PAM) owing to the intrinsic physical acoustic focusing mechanism. Here, we demonstrate a computational strategy with two generative adversarial networks (GANs) to perform semi/unsupervised reconstruction with high resolution and sensitivity in AR-PAM by maintaining its imaging capability at enhanced depths. The b-scan PAM images were prepared as paired (for semi-supervised conditional GAN) and unpaired (for unsupervised CycleGAN) groups for label-free reconstructed AR-PAM b-scan image generation and training. The semi/unsupervised GANs successfully improved resolution and sensitivity in a phantom and in vivo mouse ear test with ground truth. We also confirmed that GANs could enhance resolution and sensitivity of deep tissues without the ground truth.

et al. reported a dual orthogonal ultrasonic transducer that allowed dual-view signals 23 in nonlinear operations, which could affect quantitative measurements of the concentrations of molecules, such as blood oxygenation.M. Mozaffarzadeh et al. applied the double-stage delay-multiply-and-sum (DS-DMAS) algorithm to conventional AR-PAM combined with SAFT 24 , which incurred a higher computational complexity.Another popular solution was to use signal filtering techniques, such as FIR filtering 25 , empirical mode decomposition (EMD) 26 , and wavelet-based 27 filtering to enhance acoustic signals or integration with OR-PAM systems to show the overlay of blood vessels (captured by AR-PAM) and smaller microvessels (captured by OR-PAM) with well-quality lateral resolution while maintaining imaging depth 28 .However, the signal filtering techniques used in AR-PAM scanning require repeated recording from the same position or the use of Fourier transform that slows the speed of real-time monitoring.Integrating the OR-PAM configuration has shown that there is no actual enhancement of the resolution.
Because of these limitations, deep-learning methods that generally enhance PAI performance are regarded as new solutions [29][30][31] , such as enhancing contrast with reconstruction under low-illumination laser [32][33][34] , removing artifacts with recovery of the out-of-focus regions [35][36][37][38] , and super-resolution transformation 39 .J. Kim et al. introduced a 3D U-Net to reconstruct high-density super-resolution images from fewer 3D label-free localizations of OR-PAM 40 , but this approach required downscaling to sparse localization because of limitations with the GPU memory.Another research, Z. Zhang et al. presented a deep CNN network that provided model-based iterative enhancement algorithms to recover the latent clean image from its degraded counterparts 41,42 .Most studies use convolutional neural networks (CNNs) that learn to minimize the Euclidean distances between the predicted and ground-truth properties (U-Net, AlexNet, ResNet, etc.), but decision-making with CNNs requires expert knowledge 43 .To distinguish classification of the output image while simultaneously training a generative CNN to minimize loss, the recently proposed generative adversarial networks (GANs) can solve the generative modelling problem [44][45][46] , which specifically designed for synthesize new data instances that resemble a given dataset.In particular, several researchers have tried to use GANs: S. Cheng et al. introduced an approach with Wasserstein GAN and gradient penalty (WGAN-GP) as the training network to transform low-resolution AR-PAM images to match high-resolution OR-PAM images obtained at the same depth 47 .Unfortunately, only the 2D maximum amplitude projection (MAP) images were processed here, which limits the accuracy of the depth area information.Y. Zhou et al. reported application of a conditional adversarial network (cGAN) to suppress the side lobes and enhance the depth of focus of a Bessel beam from a Gaussian beam 48 .However, these reported models required manual pairing, which causes difficulties.
In this study, we used the dual AR/OR-PAM system for data mining with a customized agar phantom embedded with three hairs and in vivo mouse ear with manual labeling for conditional GAN (cGAN) and automatic labeling without paired samples for the cycle-consistent GAN (CycleGAN).These two-representative semi/ unsupervised GANs were compared for their performance enhancing abilities of low-resolution AR-PAM b-scan images to match the high-resolution OR-PAM b-scan images by measurement of the sensitivity parameters and differences in lateral/axial profile similarities.We further applied these GANs to a deep AR-PAM system to demonstrate quality enhancement with better lateral resolution at a greater depth without the validation from the OR-PAM configuration.

Results
Phantom study and quantitative evaluation.To demonstrate enhanced depth imaging with the proposed GANs and image quality improvement in a well-prepared phantom (Fig. 1a), we used AR-and OR-PAM modes to find the maximum imaging depth achieved by both systems.The illustration (right) in Fig. 1a shows the phantom containing three human hairs at different depths from the top surface (left) to 2 mm depth (right).In Fig. 1b, we present the MAPs of the original AR / generation by cGAN / generation by CycleGAN / original OR.The series shows the depth mapping (Fig. 1c) with the central hair cross-sectional image along the y-z plane (Fig. 1d).Three cross-sectional images along the x-z plane (Fig. 1e) named "b1", "b2", and "b3" are selected to export the axial profiles (Fig. 1f) of the center hair that clearly shows the differences among AR-PAM mode at the deepest part of the hair (up to 1.9 mm) in Fig. 1e(i).The OR-PAM mode only shows the hair at ~ 1 mm (Fig. 1d(iv),e(iv)), where the deeper part of the hair cannot be easily recognized because the signal at 1.9 mm depth is low (same as background).Applying both cGAN (Fig. 1d(ii),e(ii)) and CycleGAN (Fig. 1d(iii),e(iii)), we identified the hair at a greater depth than the original AR.To prove this, we also acquired the lateral profile of the hair signal at 500 μm depth, which showed the smaller diameters of both the original OR and generated OR (~ 150 μm) from the original AR (~ 700 μm), as shown in Fig. 1g.Quantitatively, in Fig. 1h, we measured the sensitivity by calculating the signal-to-noise ratio (SNR) and the contrast-to-noise ratio (CNR) from the MAP images.For SNR, the signal between the original OR (5.8 ± 1.7 dB) and original AR (5.5 ± 0.8 dB) had no significant difference (p-value > 0.1).Therefore, neither trained network generated enhanced signals with cGAN (6.4 ± 0.3 dB) and CycleGAN (4.9 ± 0.1 dB).The CNR of original OR showed the highest value (23.4 ± 3.1 dB), converse to that of the original AR (18.8 ± 2.6 dB).The measured contrast was observed in the generators of both cGAN (20.1 ± 1.7 dB) and CycleGAN (14.9 ± 0.5 dB).To estimating sensitivity at difference depth from 100 to 1600 µm, we measured the SNR at different 3 depth positions as shown in Supplementary Fig. S1a.At the position R3, although the original OR showed nothing (0.023 ± 0.015 dB) and the original AR showed a low signal (0.2 ± 0.1 dB), PA signals at the same position were recovered well by cGAN (3.1 ± 0.3 dB) and CycleGAN (4.4 ± 0.1 dB).By measuring the full width at half maximum (FWHM) of the hair diameter, the original AR mode showed a larger diameter of the hair (338 ± 7 μm) than the original OR mode (134 ± 5 μm).Using trained GANs, the FWHM of the original AR was enhanced by both cGAN (144 ± 12 μm) and CycleGAN (136 ± 20 μm).In short, this hair phantom experiment showed the deeper parts from the AR-mode with better FWHM (~ 80%) than the OR-mode, and the recovery using GANs was successful for enhancing imaging resolution.To check its lateral resolution variance by imaging depth changing, we also measured the FWHM value at 6 different depth positions from 100 to 1100 μm as shown in Supplementary Fig. S2.The original AR and original OR showed at d3 (500 ± 10 μm) with 338 ± 7 μm and 134 ± 5 μm, respectively.Both cGAN and CycleGAN could maintain their FWHM well with 144 ± 12 μm and 136 ± 20 μm at d4 (700 ± 10 μm), respectively.Especially, CycleGAN achieved its latera resolution ability at d5 (900 μm) with 170 ± 38 μm.The structural similarity index measure (SSIM) value of the b-scans between the generated cGAN and original OR was 0.97 and that between the generated CycleGAN and original OR was 0.94.For more general visualization, we added the movie that showed 3D movies of AR-PAM and OR-PAM, and the recovering process of AR-PAM by cGAN and CycleGAN (Supplementary Fig. S3a and Movie S1). the in vivo experiment, we treated a mouse following the process described in the method section.We used ARand OR-PAM modes in turn to quickly scan the mouse ear under anesthesia.In Fig. 2a, we present the MAPs of the original AR, generation by cGAN, generation by CycleGAN, and original OR with the same configuration as the previous experiment.The series shows the depth mapping (Fig. 2b).We specially marked four regions of interest (ROIs) (yellow rectangle marked in Fig. 2a and magnified in Fig. 2c).The lateral profiles of the original AR (blue) / original OR (green) / generated cGAN (red)/generated CycleGAN (purple) were cropped and plotted after Gaussian fitting.The diameter in the lateral profile based on FWHM is shown in the comparison.In Fig. 2d(i), we measured the sensitivity by calculating SNR and CNR.For SNR, the signal between the original OR (77.4 ± 0.4 dB) and original AR (14.9 ± 1.1 dB) showed large gaps in performance (original OR higher than original AR by ~ 519%).Therefore, both trained GANs generated enhanced signals, with cGAN (60.3 ± 0.1 dB) and CycleGAN (74.7 ± 0.1 dB).The CNR of the original OR was higher (53.9 ± 0.7 dB) than that of the original AR (19.6 ± 1.9 dB).The measured contrast in the cGAN generator was 20.2 ± 0.5 dB, which is worse than that of CycleGAN (27.6 ± 0.8 dB).We measured the lateral resolution by the vascular diameter (Fig. 2d(ii)), and the original AR mode showed a larger diameter (124 ± 5 μm) than the original OR mode (75 ± 12 μm).Using trained GANs, the lateral resolution of the original AR was enhanced by cGAN (76 ± 20 μm) better than CycleGAN (95 ± 12 μm).The SSIM value in MAP between the generated cGAN and original OR was 0.75 ± 0.02 and that between the generated CycleGAN and original OR was 0.74 ± 0.02.For more general visualization, we added the movie that showed 3D movies of AR-PAM and OR-PAM, and the recovering process of AR-PAM by cGAN and CycleGAN (Supplementary Fig. S3b and Movie S2), with B-scan profile at yellow line in Fig. 2a (Supplementary Fig. S4).

Quantitative evaluation of in vivo deep mouse-back experiment (w/o ground-truth).
By applying the deep AR-PAM system alone, we conducted the same experimental process without the ground truth.We used the deep AR-PAM system to scan a side of the mouse's back skin under anesthesia.In Fig. 3a, we present the MAP of the original AR, generation by cGAN, and generation by CycleGAN without original OR.

Discussion
By demonstrating the proposed unsupervised and semi-supervised GANs for enhanced depth and image quality improvement in a customized phantom (Fig. 1a), we proved that the enhanced depth signal with the AR-PAM system was better than that with OR-PAM under 1 mm.The proposed GANs enhanced the source AR-PAM b-scan images and showed lateral resolutions similar to the OR-PAM b-scan images on the surface through FWHM measurements, as shown in Fig. 1g (close to ~ 80%).From the series showing the depth mappings (Fig. 1c) for the central hair cross-sectional image at the y-z plane (Fig. 1d), we chose three b-scan images at "b1", "b2", and "b3" to prove that both cGAN and CycleGAN work well.In conjunction with Supplementary Fig. S1, we also observed that both cGAN and CycleGAN were successful in effectively restoring the high PA signal in regions located at depths beyond 1400 µm when compared to the original AR and the original OR images.The measured lateral profiles of FWHM in both cGAN and CycleGAN could keep and recover well the original OR lateral resolution when original AR out-of-focus deeper than 700 µm (both GANs model are shallower 4 times than original AR), as shown in Supplementary Fig. S2.Despite the attenuation of the PA signal with imaging depth, the FWHM of the axial profiles in Fig. 1f remains relatively similar.The generated axial profile's FWHM exhibits approximately 10% lower values compared to the originals.This can be attributed to the utilization of the same measurement transducer, resulting in comparable axial resolution.The quality measurements via SNR and CNR showed insignificant differences between the original OR and original AR (not significant), and we consider that both cGAN and CycleGAN did not improve SNR in the phantom experiment.
From the in vivo experiment with mouse ear monitoring by the dual AR/OR-PAM system, as shown in Fig. 2a, we present the MAP of the original AR-PAM with more blurring and loss of most of the small vascular area than the original OR-PAM.By applying the proposed GAN pretrained datasets, the generations by cGAN and CycleGAN could recover the missing vessels from the original AR b-scan images but not perfectly.The salt and pepper noise from the original OR caused unequal tissue distribution (skin, fat, blood, muscle, cartilage, and small hairs on mouse ear) as it was difficult to prepare a noise-free training dataset.The small gaps between the vasculature were hard to observe in the proposed GANs.In the specially marked four ROIs in Fig. 2c, the lateral profile of the original AR (blue) was significantly larger in diameter than the original OR (green), so the proposed generated cGAN (red) and generated CycleGAN (purple) successfully enhanced and recovered the unlined small vascular structures.Therefore, both trained GANs generated enhanced signals, with cGAN and CycleGAN having three times higher results than the original AR.For visualization, we showed the B-scan www.nature.com/scientificreports/image (yellow line in Fig. 2a) and their lateral profile in Supplementary Fig. S4.The CNR of CycleGAN was higher by a few points (~ 2.9 dB) than both the cGAN generator, and no significant difference between original AR the original OR.In Fig. 2d(ii), we cropped the lateral profile to show single small vascular structures, where the proposed GANs converted the original AR to sharp images better by cGAN than CycleGAN.This in vivo experiment proved that better lateral resolution could be repeatedly recovered using GANs with similar structures as the original OR (~ 0.75 by SSIM).Even with unstable vessel structures, the sensitivity was successfully improved from the low-level original AR.The absorption properties in the custom human hair phantom (having the most amount of melanin) are stronger than those of soft biomaterials (blood, vessel borders, skin, etc.) 49 .These indicate that the GANs recovered complex patterns with blurring and noise under the skin and produced good effects by increasing the SNR.
To improve the process without the ground truth images, we used the deep AR-PAM system without OR-PAM to scan a mouse's side back that showed the twisted vascular structures and changes based on depth.Here, the proposed GANs still produced better reconstruction and lateral resolution than the original AR (~ 190% sharper).We measured the lateral profiles at three ROIs (in Fig. 3c-e), which showed similar trends of sharp vascular structures and diameter from both the pretrained cGAN and CycleGAN.The measured sensitivity by SNR and CNR (Fig. 3f(i)) showed significantly different when compared with the original AR (decrease ~ 34% noise and ~ 18% enhancing contrast for the proposed GANs) (p < 0.1).We also checked the sensitivity at specific ROIs of different depth images as shown in Supplementary Fig. S5.The great visualization that showed the well recovery from saturating PA signal from near surface AR-PAM image (enhancing contrast and reducing noise) but be poorer at lower depth it showed unclear structures.Fortunately, the generated PAM images showed the detail vasculature map at more deep depth.Due to the use of different transducer specification in preparing the deep AR-PAM data for the testing dataset, both cGAN and CycleGAN models were unable to significantly enhance the lateral resolution as effectively as the OR-PAM system.Subsequent to the enhancement of the deep AR-PAM b-scan images, no thresholding or bandpass filter was employed.Nevertheless, notable contrast enhancement was observed in the synthetic b-scan images.That could be noticeable that thresholding and bandpass filter configuration in dual OR/AR-PAM was learned and applied by GANs models, so the manual denoising was able to pass temporarily.
In most experiments with GANs, researchers have only selected the best results they could show because the GANs are difficult to control while maintaining good conditions.To compare the differences between the cGAN and CycleGAN based on training time (epochs), we chose in vivo experiments that showed significant bias efficiency over the rest for the original OR (SNR, CNR, lateral resolution, and SSIM for p < 0.1 with the best weight).In Supplementary Fig. S6, we showed the epochs 10-50-100-150-200 as the marks and separately use cGAN and CycleGAN to generate results from the original AR in vivo dataset.The SNR (Supplementary Fig. S6a) shows balance during training time in both GANs.The CNR (Supplementary Fig. S6b) shows a similar value for cGAN as compared to the original AR (similar in the phantom experiment).Hence, the CNR was not improved by cGAN in both the phantom and in vivo experiments.However, the CycleGAN showed the highest score for 10 epochs and then decreased close to that of the original OR.This was opposite to the result of the phantom experiment, where the CNR of CycleGAN was lowest (~ 1.5 times less than original OR).The lateral resolution of cGAN was smaller than that of the original AR and had an earlier balance and similar diameter as the original OR.The generated signal by CycleGAN slowly decreased at 150 epochs and was close to that of the original OR.These observations were also proved by the SSIM via comparison of the original AR, generated cGAN, generated CycleGAN separately with the original OR as initial comparison (Supplementary Fig. S6d).Because the SNR and CNR values of the original AR were the smallest, and the lateral resolution was highest, the SSIM was also the smallest (0.14 compared with the original OR).The value generated by cGAN slowly increased (from ~ 0.68 to 0.8), following the same trend as the value generated by CycleGAN (from ~ 0.63 to 0.77), whose learning was best after 200 epochs.We conducted the custom hair phantom experiment with a simple condition (fixed diameter, simple structure) and the in vivo mouse ear experiment with a more complex condition (dynamic vascular structure based on time, flexible structure).These proved the efficiencies of both cGAN and CycleGAN with well-prepared samples.The problem with the mouse ear prepared dataset was that it did not contain multilayer information.However, enhancement for depth changes and noise-free lateral resolution are still efficient features of the deep AR-PAM scanning, as in the previous experiments with MAP 47 , b-scan image training 48 , or 3D training 40 .These results were verified by in/ex vivo animal experiments with random samples.
This study still has some major limitations.First, the model has a lot of parameters that should be carefully optimized cause the stability of the GAN is highly dependent on network architecture (our generator was based on ResNet-9 blocks, which can be modified to U-Net or another generator networks).Second, there were many factors that affected the choice of GPU (price was too high for the experiments, memory of GPU was not enough for full-scale b-scan, training time was too long to optimize the parameters, and training failure caused random corrupted memory during training).Lastly, there is a limitation with the dataset in that it does not cover all scales of biostructures (only proven with vascular shapes and not organs, lymph nodes, bones, etc.) and the dataset was prepared on the limited system specification with our own experiment configuration.To deal with current issues, we are going to design a translation standard which can combine our current GANs following different transducers configuration in simple GANs models.The GANs can be used to replace manual enhancing and reduce the processing time to be able for research and commercial PAM system.

Materials and methods
Dual AR/OR-PAM.We prepared a dual-mode AR/OR-PAM system as shown in Fig. 4a.A pulsed laser with a central wavelength of 532 nm (SPOT-1010-200, 10 kHz repetition rate, Elforlight, France) served as the main source for the dual AR/OR-PAM system.The beam was first delivered to a collimator (C, F280APC-A, Thorlabs, USA) and separately connected to two optical fibers (i.e., single-mode fiber for OR-PAM and multimode fiber for AR-PAM) depending on the PAM operating mode.To achieve OR-PAM mode, a 5 μm single-core optical fiber (SMF, P1-405BPM-FC-1, Thorlabs, USA) was used, and the delivered beam was focused using a doublet lens (OB, AC254-060-A, Thorlabs, USA).To demonstrate AR-PAM mode, a 400 μm multicore optical fiber (MMF, M74L01, Thorlabs, USA) was used, where the focused beam was first passed through a homemade beam combiner (BC) coated by thin aluminum-sheet.The focused laser beam was then reflected and directed with a one-axis MEMS scanner (MEMS, OpitchoMS-001, 50 Hz, Opticho Inc., Ltd., Korea).The laser beam was absorbed by the biological tissue, and through subsequent thermal expansion, the acoustic signal was exposed to 360 degrees.The generated acoustic signal was reflected by the MEMS mirror and directed through the aluminum-prism to a high-frequency 50 MHz ultrasonic transducer (TR, V214-BC-RM, Olympus NDT, Japan).The whole transducer was aligned and immersed in the water bath.To amplify the acoustic signal, an amplifier (AMP, ZX60-3018G-S+, MiniCircuit, USA) was used, followed by crystal low-pass filtering (CLPFL-050, CRYTEK, USA).The analog signal was converted to digital signal with a high-speed digital board (DAQ, ATS9371, AlazarTech, Canada) that digitized the 10 μm depth-resolved analog signal to 16-bits digital A-scans before storage.To scan a large-field surface, two linear stepper-motorized stages (L-509-10SD00, Physik Instrumente, Germany) along the x-y plane and a MEMS scanner, with a 0.127 μm x-step and 0.386 μm y-step, were used to obtain data.The measured lateral resolutions in OR-PAM and AR-PAM were 12 μm and 85 μm, respectively.The axial resolution was 27 μm, which was fixed for the same transducer.For the synchronous control stages, a custom LabVIEW program was developed to manage the hybrid scanning system with a data acquisition board (PCIe-6321, NI Instruments, TX, USA).The laser was controlled at approximately 5 mJ/cm2 below the American National Standards Institute (ANSI) limit (20 mJ/cm2 under visible light).All data analyses and reconstructions were performed using MATLAB (R2020b, The MathWorks, MA, USA).The 3D visualization in Supplementary Fig. S3 was created by movie maker tools in 3D PHOVIS 50 .
Deep AR-PAM.We configured the deep AR-PAM system 7 as shown in Fig. 4b.A pulsed laser beam source with a central wavelength of 700 nm (Surelite OPO PLUS, Continuum, USA) was controlled by a Q-switched Nd:YAG laser (5 ns pulse duration and 10 Hz pulse repetition rate, SLII-10, Continuum, USA).The diverged ring-shaped beam configuration formed a spherical conical lens (AX2540-B, Thorlabs, USA) and a handmade optical condenser.Using this imaging setup, the dark-field light illumination was coaxially aligned with ultrasound focusing on water.The incident laser pulse energy was limited to below the ANSI safety limit of 20 mJ/ cm 2 .The generated ultrasound waves were detected with a single-element 5 MHz ultrasonic transducer (V308, Olympus NDT, Japan) which immersed in the water bath.The axial and lateral directions were 150 μm and 590 μm, respectively.

Data preprocessing.
Figure 4c shows the data preprocessing flow chart.To prepare the dataset acquired with the dual AR/OR-PAM system, we saved all normalized b-scan images as stacks in a 3D volume array.The 3D volume was loaded as a single task during processing and separated into multiple frames with 256 × 256 pixels for each image.The random jitter function was disabled to avoid scaling errors.The photoacoustic amplitudes were normalized for the entire volume to ensure that the laser power remained identical.To reduce the background noise, we applied bandpass filtering based on difference ultrasound transducer properties to each envelope A-line dataset for dual AR/OR-PAM and deep AR-PAM system.

Experimental sample preparation.
As experimental materials, we prepared three samples: a customdesigned three hairs embedded phantom, a mouse ear, and a mouse back skin.These samples were prepared to monitor the sensitivity and structure of the objects and their effects on the training models.
A phantom embedded with three hairs was prepared to show different imaging depths and resolutions between the AR-and OR-PAM images.Using a 3D printer (da Vinci 1.0 Pro, XYZ printing Inc., USA) 51 , a white polylactic acid frame (30 × 50 mm) with three top and three bottom holes was designed to fix three long black hairs at linearly increasing depths from the surface to 5 mm.The frame was internally filled with 5 g of 1.5% concentration mold jelly (BD BACTO™ Agar, Becton-Dickinson and Company, USA) 52 mixed with 3 g of milk and 100 mL pure water.We cooked this mixture at 80 °C for 5 min and stirred well until the Agar dissolved completely.A solid formed stably after 15 min in the refrigerator and maintained under a wet towel.
We prepared two in vivo samples for this study: a mouse ear for use with dual AR/OR-PAM scanning and a side-back skin of a mouse for deep AR-PAM scanning.Two healthy six-week-old male BALB/c mice (Orientbio, Korea) were used according to the instructions of laboratory animal protocols approved by the Institutional Animal care and Use Committee of Chonnam National University Hwasun Hospital (CNU IACUC-H-2020-35).The mice each weighed ~ 20 g and were anesthetized by intraperitoneal injection of a cocktail of Ketamine/Xylazine at 70/15 ratio.After removing the surface soft hairs, we placed the experimental mouse on a small animal table.For better acoustic coupling in vivo animal experiment, we used the ultrasound gel (hypo-allergenic water soluble, Sono Jelly, Meditop Co., Ltd, Korea).An isoflurane system (Luna Vaporiser, NorVap International Ltd., Barrowford, UK) was used for gaseous anesthetization, and a temperature table was used for short-term in vivo observations to maintain stable body condition of the mouse 53 .The energy of the laser pulse was ~ 5 mJ/cm 2 .Each sample was prepared and imaged separately.We scanned the mouse ear vasculature in both AR-and OR-PAM modes because the size of the mouse ear was small enough (10 × 10 mm 2 ) for fast scanning.For the side-back skin of mouse scanning (100 × 50 mm 2 ), the deep AR-PAM system required a longer scanning time (~ 20 min) owing to the mechanical scanning speed and laser repetition rate 7 .Hence, we divided the volume into four vertical stacks and displayed each b-scan image with 300 × 512 pixels to acquire the 3D volume.
Image-to-image GANs.Typical GAN models have been used in various imaging tasks 54 .To clarify the opposing relationships between AR-PAM and OR-PAM b-scan images, we applied different pixelwise contrastive GAN models, cGAN and CycleGAN, to learn the mapping functions.This was to prove the effectiveness of unsupervised learning with different datasets (labeled in cGAN and label-free in CycleGAN) in AR-PAM for enhancement of the intended sample (phantom, mouse ear vascular with ground truth from OR-PAM scanning) and uncertain sample (mouse's back skin without ground truth).
The schematic of the GANs consists of competing between generators and discriminators to create and discern in the network.In both cGAN and CycleGAN, we used ResNet-9 blocks (Fig. 5a) in the generator to The output of the generator was 256 × 256 patches of b-scan images.The discriminator used was a Markovian discriminator (Fig. 5b) that relied on a term to force low-frequency accuracy even with large b-scans.By multiplying the generated patch b-scan images with the target patch OR-PAM b-scan images, the discriminator encoded multiple pool mask layers to a 1 × 1 decision vector by the sigmoid function.Finally, the L1 distance was updated loss distance in the generator.
To replicate a probability photoacoustic amplitude distribution, the GANs were subjected to get the minimax task (Eq. 1) between the generator and discriminator to discriminate the learning process.In our study, that minimax task performed distribution of the generated data (G) between the original/source AR-PAM and original/target OR-PAM data (O).This was aimed at distinguishing the actual b-scan image map (X, Z) with the discriminator (D).The final generator G * and discriminator D * were created using the minimum argument ( arg ) of loss function L(G, D) to minimize the generator min G and maximize the discriminator max D to ensure lower difference.P. Isola et al. proposed the modified cGAN 55 to learn features from AR-PAM b-scan images paired with OR-PAM b-scan images.In Fig. 6a, the cGAN is shown with one generator and one discriminator.By labeling the AR-PAM b-scan images with the OR-PAM b-scan images in a one-to-one manner, the generator encodes the AR-PAM b-scan images to a 1D vector as "features".The residual blocks bypass the nonlinear layers with an identity mapping to remove the AR-PAM's features.The decoder part of the generator adds the "learned features" that are updated from the previous epoch.Figure 7 shows that the loss function used for training both cGAN and CycleGAN is L1 loss.Both networks' weights assigned for loss are 0.01.For optimizing the network, the Adam optimization algorithm was used with a learning rate scheduler to reduce the learning rate by a factor of 0.5 when the initial learning rate was 0.0002.To obtain 200 iterations, cGAN only required around 6 h, which was less than that of CycleGAN that was trained for around 20 h (Table 1).The CW-SSIM was used to analyze 4000 paired b-scans, and the best value of 0.72 was observed at 200 epochs.
Quality comparison.We evaluated the effectiveness of the model and compared it with the target using quality assessment methods.Apart from perceptual quality, three concept metrics, including SNR, CNR, and

Figure 3 .
Figure 3.In vivo experiment of the side-back skin of mouse with deep AR-PAM imaging and generated OR-PAM images with GANs.(a) MAP images, (b) corresponding depth-resolved MAP images, (c-e) lateral profile comparisons of the original AR (blue), generated cGAN (red), and generated CycleGAN (purple).(f) Sensitivity and lateral resolution of mouse vascular mapping.Scale bar 1000 µm (*p-value: (*) < 0.1).

Figure 4 .
Figure 4. (a) System configuration of the dual AR/OR-PAM system: (i) Part of OR-PAM with optical focusing (green), (ii) Part of AR-PAM with acoustic focusing (gray), (b) system configuration of the deep AR-PAM system, (c) preprocessing workflow.M mirror, SMF single-mode fiber, MMF multimode fiber, C collimator, OB objective lens, BC beam combiner, TR transducer, AMP amplifier, MEMS MEMS scanner, DAQ data acquisition, PS prism, SCL conical lens, OPO optical parametric oscillator, WT water bath.
13:13423 | https://doi.org/10.1038/s41598-023-40583-xwww.nature.com/scientificreports/avoid overfitting and achieve lower training errors when the details increase in the b-scan images.The input to the generator was the 256 × 256 patches of b-scans from AR-PAM.Between the encoder and decoder layers, the generator used nine residual blocks (green block) containing the identify function as one of the elements.
By comparing the generated b-scan images with the paired OR-PAM b-scan images, the discriminator then encodes the multiple generated/OR to a sigmoid activation function that creates the updated decision for the next epoch.Using the loss function (Eq.2) for cGAN ( L cGAN ), cGAN drives the generator ( G ) from the expected values ( E ) of the original AR b-scan images ( A ) and paired pixelwise discriminator ( D ) from the original OR b-scan images ( O).
Based on image-to-image translation, J. Y. Zhu et al. introduced the unpaired learning method using CycleGAN56 .Inspired by the CycleGAN structure, we denoted the data distribution as two groups of PAM images with expected values of original AR-PAM b-scan images ( E A ) and original OR-PAM b-scan images ( E O ) without pairing.Similar to cGAN, the process with the generator/discriminator learned from the updated decision that affected the next epoch.In Fig.6b, by identifying the backward-looking generator (G) , the reverse is confirmed so that can the cGAN can be used in both directions (from AR to OR-L cGAN (A, O) and inverse-L cGAN (O, A)).