# Biomedical Photoacoustic Imaging Optimization with Deconvolution and EMD Reconstruction
## 3 Results
### Experiments  
The schematic of our experimental setup is shown in Figure 4a. A tunable optical parametric oscillator (OPO) system (Phocus Mobile, Opotek, CA, USA) is pumped by a Nd; YAG laser (Brilliant B, Bigsky, MT, USA) is employed as the laser source. The laser works at the wavelength of 720 nm with pulse repetition rate of 10 Hz and pulse duration of 5.5 ns. A Verasonics system (Verasonics Inc., Kirkland, WA, USA) is controlled by a computer to receive raw PA data. An ultrasound linear array (L7-4, Philips, Amsterdam, The Netherlands) with 128 elements, 5.208 MHz center frequency and 0.298 mm element pitch was used. Ultrasonic coupling agent is used as the acoustic transmission medium between phantom and ultrasound probe.  
The phantoms used in this experiment were made of porcine gelatin with 8% concentration as shown in Figure 4b,c. Figure 4b is the phantom used as the imaging target in our experiment which  
Figure 4: (**a**) Schematic of the experimental setup; (**b**) the phantom used as imaging target (two rubber threads); (**c**) the phantom to acquire the deconvolution kernel (one tungsten wire); (**d**) a photo of the experiment setup.  
embedded with two rubber threads of diameter 1 mm and Figure 4d shows the placement of phantom, probe and laser in this experiment. The arrows in Figure 4b,\(c\) both point to the position where a probe acquired a signal of each phantom while carrying out the experiment. During the process of signal acquisition, each set of signals contained approximately 200 frames for a continuous period to average.  
The deconvolution kernel is acquired from the middle transducer element of the probe, which detected a single tungsten wire of diameter 0.3 mm (as Figure 4c); approximately the smallest size L7-4 probe can measure. The normalized deconvolution kernel (as shown in Figure 5), averaged from 200 frame signals, was used as the deconvolution kernel in this experiment.  
The comparison of the reconstructed result is shown in Figure 6. We can observe that the two oblique rubber wires are too close to look separated in the image reconstructed with original PA signal, as shown in Figure 6a. Without knowing the imaging targets in advance, it is difficult to distinguish the structure visually. After applying our proposed method, as shown in Figure 6b, the two lines are completely separated, and the realistic inner structures of this phantom can be more clearly identified. Meanwhile, image resolution is enhanced, and the artifact is slightly reduced.  
To quantify the width of gaps in Figure 6 and compared to the realistic width of gap in Figure 4b, we measured the width of gaps in Figure 6. In Figure 4b, the diameter of the two rubber wires in the phantom are both 1 mm, we can observe that the realistic gap, where the arrow mark is pointing to, is approximately 0.4 mm. In respect of the reconstructed image of the original signal as shown in Figure 6a, the width of the gap is about 0.08 mm. For the reconstructed image of the processed signal as shown in Figure 6b, the width of the gap is about 0.3 mm. From the quantification, we can find that the image of processed signal more aptly indicates the accurate width and position of realistic targets.