Motion rejection and spectral unmixing for accurate estimation of in vivo oxygen saturation using multispectral optoacoustic tomography

Mitradeep Sarkar,Mailyn Pérez-Liva,Gilles Renault,Bertrand Tavitian,Jérôme Gateau
DOI: https://doi.org/10.1109/TUFFC.2023.3306592
2023-09-15
Abstract:Multispectral Optoacoustic Tomography (MSOT) uniquely enables spatial mapping in high resolution of oxygen saturation (SO$_2$), with potential applications in studying pathological complications and therapy efficacy. MSOT offers seamless integration with ultrasonography, by using a common ultrasound detector array. However, MSOT relies on multiple successive acquisitions of optoacoustic (OA) images at different optical wavelengths and the low frame rate of OA imaging makes the MSOT acquisition sensitive to body/respiratory motion. Moreover, estimation of SO$_2$ is highly sensitive to noise, and artefacts related to the respiratory motion of the animal were identified as the primary source of noise in <a class="link-external link-http" href="http://MSOT.In" rel="external noopener nofollow">this http URL</a> this work, we propose a two-step image processing method for SO$_2$ estimation in deep tissues. First, to mitigate motion artefacts, we propose a method of selection of OA images acquired only during the respiratory pause of the animal, using ultrafast ultrasound images (USIs) acquired immediately after each OA acquisition (USI acquisition duration of 1.4 ms and a total delay of 7 ms). We show that gating is more effective using USIs than OA images at different optical wavelengths. Secondly, we propose a novel method which can estimate directly the SO$_2$ value of a pixel and at the same time evaluate the amount of noise present in that pixel. Hence, the method can efficiently eliminate the pixels dominated by noise from the final SO$_2$ map. Our post-processing method is shown to outperform conventional methods for SO$_2$ estimation, and the method was validated by in vivo oxygen challenge experiments.
Medical Physics,Image and Video Processing,Optics
What problem does this paper attempt to address?
The problems that this paper attempts to solve mainly focus on two aspects: 1. **Reduction of motion artifacts**: Multispectral Optoacoustic Tomography (MSOT) has a low imaging frame rate when performing high - resolution oxygen saturation (\(SO_2\)) spatial mapping because it needs to continuously acquire optoacoustic (OA) images at different optical wavelengths. This low frame rate makes MSOT very sensitive to body movement or respiratory motion. Especially in small - animal imaging, the high respiratory rate under anesthesia can lead to significant motion artifacts. Therefore, the paper proposes a method to reduce the impact of these motion artifacts by selecting optoacoustic images acquired during animal apnea. 2. **Improvement of the accuracy of oxygen saturation estimation**: Traditional \(SO_2\) estimation methods usually involve first determining the absolute concentrations of deoxyhemoglobin (Hb) and oxyhemoglobin (HbO_2) in each pixel and then calculating their ratio. This method is easily affected by noise. Especially when there is a lot of noise in the data, it will seriously affect the final \(SO_2\) image quality. The paper proposes a new method that can directly estimate the \(SO_2\) value of each pixel and simultaneously evaluate the noise level in this pixel. By setting a threshold, those pixels dominated by noise can be effectively excluded, thereby improving the accuracy and reliability of \(SO_2\) estimation. Specifically, the solutions proposed in the paper include: - **Step 1**: Use ultrafast ultrasound images (USIs) to select optoacoustic images acquired during animal apnea. By calculating the autocorrelation coefficient of USIs, screen out those images acquired during apnea, thereby reducing the impact of motion artifacts. - **Step 2**: Develop a variance - based method to directly estimate the \(SO_2\) value of each pixel and evaluate the noise level in this pixel. Determine the most appropriate \(SO_2\) value by minimizing the Residual (Res) and exclude pixels with large noise by setting a threshold. Through these methods, the paper demonstrates its effectiveness in in - vivo oxygen challenge experiments and proves the superior performance of this method in \(SO_2\) estimation.