-
以厚度1 mm的猪肉切片为散射介质,光源采用光纤卤钨灯和红外滤光片,对USAF-1951分辨率板进行成像,InGaAs探测器为单像素探测器进行信号采集,具体参数和光路如表1、图4所示。
Key device Parameter Light source Wave length: 260-2500 nm Lens 1 AR COATED 650-1050 nm, f=26 mm Lens 2 f=100 mm Infrared filter 800-2500 nm DMD Operating 350-2700 nm, 1024×768 pixel Projection lens f=35 mm Single point detector Operating 800-1750 nm Focusing lens AR COATED 650-1050 nm, f=60 mm Table 1. Key device parameters in the experiment
Figure 4. Layout of active single-pixel imaging system based on DMD. (a) When there is no scattering medium; (b) When the scattering medium in detection path; (c) When the scattering medium inimaging path
首先在无散射介质的情况下实施近红外单像素成像实验,采用分辨率为64×64的阿达玛矩阵,分别使用CGI算法和TVAL3算法在不同采样率下恢复。为了减少随机误差的影响,循环采集四次数据并取平均值作为实验结果的最终测量值。
CGI算法使用4096个测量值重构了目标图像,PSNR与SSIM分别为16.3051、0.6337。TVAL3算法在采样率为0.4时就可以重构出图像的轮廓和细节特征,PSNR与SSIM分别为16.9636、0.6326,如图5(b)所示,两种算法的重构效果相近。可以看出,在无散射介质时,基于压缩感知理论的TVAL3算法的成像质量优于基于关联成像理论的CGI算法。
Figure 5. Comparison of the reconstruction effects of the two algorithms with non-scattering medium. (a) The first image recovered by CGI, and the last four images are recovered by TVAL3, with a sampling rate of 0.4-0.7; (b) PSNR and SSIM comparison of reconstructed image
在图4所示的光路上加入散射介质后,介质的多重散射作用促使光束随机偏离原来的传播方向。成像光路上结构光的空间信息将被打乱[20-21],抵消了结构光照明可增加有效信息提取[22]的优势,造成成像质量下降。如图6 (a)、(b)所示,在结构光的传播路径上放置毛玻璃做散射介质。在d1距离处,光束经过毛玻璃后空间结构细节变得模糊,并且距离越远,图像越模糊,最终如图6 (b)所示,在d2距离处图像失去正方形形状,变成一个大的光斑。同样的,如果将散射介质放置在非成像光路中,光被混匀,可探测到的光束强度衰减。如图6(c)、(d)所示,均匀光经过物体和紧贴物体的毛玻璃后,光束的传播方向同样发生了变化,如同经过了一个光扩散板,并且随着距离变长,光束愈发扩散。
Figure 6. Performance of the light beam under the scattering of ground glass. (a)-(b) Intensity distribution comparison of structured light with or without scattering medium; (e) Uniform light passing through the object and scattering medium; The upper of (c1)-(c2) is a photo taken directly, and the below is a picture taken from the observation surface; (c3) Uniform light passing through the object and scattering medium; (d) Uniform light passing through the object
选择厚度约为1 mm的生物组织,分别固定在分辨率板的前面(图4 (c))与后面(图4 (b))。介质的散射作用会限制成像深度,使被介质包裹着的目标信号不容易被探测器接收到,并且探测到的微弱信号易与探测器噪声混杂。因此,在重构运算前,需要在测量值中减去无光和有光环境下探测器的暗噪声平均值。此外,在重复测量的数据中选择最大的数据作为测量值。经过预处理的测量值可以在一定限度内减小探测器噪声、环境光照、光源强度波动对测量值的影响。对于TVAL3算法,选取多个采样率进行图像重构并对比,调整了重构图像的灰度,增加非线性映射操作,增强了对比度。
图7的成像结果表明,生物组织散射介质对成像有一定程度的影响,重构图像质量下降。CGI算法得到的图像底噪多,其SSIM值略低于TVAL3算法,但PSNR值(10.1886、9.9831)均高于TVAL3算法重构图像的PSNR(9.170456)。当TVAL3算法的采样率减小时,图像结构越发不明显。
Figure 7. Reconstructed image of penetrating scattering medium. (a) When the medium in imaging path, using CGI algorithm and TVAL3 algorithm, the latter sampling rate is 1-0.6, its peak signal-to-noise ratio and similarity. Comparison of (a) is in (b); (c) When scattering medium in detection path, and the comparison is in (d)
当涉及到动态成像时,单像素成像对成像时间也有要求。实验选用了高速切换状态的DMD对光束进行编码[23],其照明时间与探测器的响应时间相对应,因此不同场景下的数据采集时间相同,在静态成像下对成像结果影响不大。图8为在有散射介质的情况下两种算法分别重构图像所需的时间,可以看出CGI算法所需的时间远小于TVAL3算法。
Figure 8. Comparison of reconstruction time between CGI algorithm and TVAL3 algorithm in penetrating scattering media imaging
从以上实验结果可以看出,在有散射介质的情况下,CGI算法重构的图像质量略优于欠采样率的TVAL3算法,而可压缩性正是压缩感知理论的特点,采样率越低时,基于压缩感知理论的重构算法优势逐渐降低。
Scattering interference suppression for single-pixel imaging reconstruction of biological tissues
doi: 10.3788/IRLA20210722
- Received Date: 2021-09-30
- Rev Recd Date: 2021-12-06
- Publish Date: 2022-04-07
-
Key words:
- single-pixel imaging /
- scattering medium imaging /
- compressed sensing /
- correlation algorithm
Abstract: In view of the influence of the scattering medium in the reconstruction of single-pixel imaging, the reconstructed image cannot achieve the best effect. The applicability of the correlation algorithm and compressed sensing algorithm for image reconstruction with or without the scattering medium was investigated. The influence of the spatial structure change of modulated information in the imaging path and the signal loss in detection path caused by the medium was analyzed, a near-infrared single-pixel imaging system was established, and the single-pixel imaging of penetrating the biological tissues scattering medium with the CGI algorithm and the TVAL3 algorithm was realized. It was found that the reconstruction time, peak signal-to-noise ratio and SSIM of the TVAL3 were better than CGI when there was no medium; while two of the three values of the CGI were better when there was medium, its maximum reconstruction time (0.304091 s) was smaller than the minimum (1.766299 s) of the TVAL3, and its minimum PSNR (9.9831dB) was higher than the maximum (9.170456 dB) of the TVAL3, and its SSIM (0.0982,0.1178) lay within the range of the SSIM of the TVAL3 (0.099258-0.497622). The results show that the CGI based on correlation imaging theory is more suitable for imaging scattering media, and the TVAL3 based on compressed perception theory is more suitable for imaging non-scattering media.