搜索

x

留言板

尊敬的读者、作者、审稿人, 关于本刊的投稿、审稿、编辑和出版的任何问题, 您可以本页添加留言。我们将尽快给您答复。谢谢您的支持!

姓名
邮箱
手机号码
标题
留言内容
验证码

基于数字全息的血红细胞显微成像技术

张益溢 吴佳琛 郝然 金尚忠 曹良才

基于数字全息的血红细胞显微成像技术

张益溢, 吴佳琛, 郝然, 金尚忠, 曹良才
PDF
HTML
导出引用
  • 数字全息显微成像有别于传统光学显微成像, 可根据重建全息图获取细胞的生物学参数与形貌信息, 是一种有效的非接触无损三维成像技术. 随着像感器的发展与硬件计算能力的提升, 数字全息显微成像技术在活体生物细胞检测尤其在血红细胞检测领域取得了显著进展和突破. 本文介绍了同轴、离轴以及光镊辅助离轴的数字全息显微技术, 这些技术利用瑞利索末菲反向传播算法、清晰度量化算法、分水岭分割算法、数字重聚焦方法与热涨落方法等来实现血红细胞的形变、空间分布、三维体积信息的高精度提取, 有助于糖尿病、心血管疾病、帕金森氏疾病等病理研究. 数字全息显微成像技术实现了传统三维显微成像技术难以达到的实时性和定量化检测, 由于独有的非接触、无损性特点, 在细胞成像领域应用前景广阔.
      通信作者: 金尚忠, jinsz@cjlu.edu.cn ; 曹良才, clc@mail.tsinghua.edu.cn
    [1]

    Cuche E, Marquet P, Depeursinge C 1999 Appl. Opt. 38 6994

    [2]

    Ferraro P, Coppola G, De Nicola S, Finizio A, Pierattini G 2003 Opt. Lett. 28 1257

    [3]

    Ferraro P, Miccio L, Grilli S, Paturzo M, De Nicola S, Finizio A, Osellame R, Laporta P 2007 Opt. Express 15 14591

    [4]

    Charrière F, Marian A, Montfort F, Kuehn J, Colomb T, Cuche E, Marquet P, Depeursinge C 2006 Opt. Lett. 31 178

    [5]

    Kemper B, Von Bally G 2008 Appl. Opt. 47 A52

    [6]

    Potcoava M, Kim M 2008 Meas Sci Technol 19 074010

    [7]

    袁操今, 翟宏琛 2010 光子学报 39 893

    Yuan C J, Zhai H C 2010 Acta Phot. Sin. 39 893

    [8]

    Ma J, Yuan C, Situ G, Pedrini G, Osten W 2013 Chin. Opt. Lett. 11 090901

    [9]

    王华英, 刘飞飞, 成惠, 廖薇, 赵宝群, 于梦杰, 刘佐强 2013 强激光与粒子束 25 345

    Wang H Y, Liu F F, Cheng H, Liao W, Zhao B Q, Yu M J, Liu Z Q 2013 High Pow Las Part Beam 25 345

    [10]

    王华英, 刘飞飞, 宋修法, 廖微, 赵宝群, 于梦杰, 刘佐强 2013 物理学报 62 024207

    Wang H Y, Liu F F, Liao W, Song X F, Yu M J, Liu Z Q 2013 Acta Phys. Sin. 62 024207

    [11]

    Zheng J, Pedrini G, Gao P, Yao B, Osten W 2015 J. Opt. 17 085301

    [12]

    Lin Q, Wang D, Wang Y, Rong L, Chang S 2015 Opt. Lasers Eng. 66 279

    [13]

    Wang Y, Guo S, Wang D, Lin Q, Rong L, Zhao J 2016 Opt. Commun. 366 81

    [14]

    董可平, 钱晓凡, 张磊, 张永安 2007 光子学报 36 2013

    Dong K P, Qian X F, Zhang L, Zhang Y A 2007 Acta Phot. Sin. 36 2013

    [15]

    邸江磊, 赵建林, 范琦, 姜宏振, 孙伟伟 2008 光学学报 28 56

    Di J L, Zhao J L, Fan Q, Jiang H Z, Sun W W 2008 Acta Optic. Sin. 28 56

    [16]

    马利红, 王辉, 金洪震, 李勇 2012 中国激光 39 215

    Ma L H, Wang H, Jin H Z, Li Y 2012 Chin. J. Las. 39 215

    [17]

    Cao R, Xiao W, Wu X, Sun L, Pan F 2018 Biomed. Opt. Express 9 72

    [18]

    Zhang J, Chen Q, Sun J, Zuo C 2019 Seventh International Conference on Optical and Photonic Engineering (icOPEN 2019), Phuket, Thailand, October 16, 2019 p112050 C

    [19]

    Zhang J, Sun J, Chen Q, Zuo C 2020 IEEE Trans. Comput. Imaging. 6 697

    [20]

    Amos W, White J 2003 Biol. Cell 95 335

    [21]

    Subramaniam V, Kirsch A, Jovin T 1998 Cell. Mol. Biol. (Noisy-le-grand) 44 689

    [22]

    Huang D, Swanson E A, Lin C P, Schuman J S, Stinson W G, Chang W, Hee M R, Flotte T, Gregory K, Puliafito C A 1991 Science 254 1178

    [23]

    Kemper B, Carl D, Knoche S, Thien R, von Bally G 2004 Optical Metrology in Production Engineering Strasbourg, France, September 10, 2004 p581

    [24]

    Jin H, Wan H, Zhang Y, Li Y, Qiu P 2008 J. Mod. Opt. 55 2989

    [25]

    周海宪 2006 应用光学 27 I0007

    Zhou H X 2006 J. Appl. Opt. 27 I0007

    [26]

    Nehmetallah G T, Aylo R, Williams L 2015 Analog and Digital Holography with MATLAB (Bellingham: Society of Photo-Optical Instrumentation Engineers (SPIE)) pp137–143

    [27]

    Kreis T M, Adams M, Jüptner W P 1997 Optical Inspection and Micromeasurements II Munich, Germany, September 17, 1997 p224

    [28]

    Liu C, Wang D, Zhang Y 2009 Opt. Eng. 48 105802

    [29]

    Song X, Zhao B, Qin A 4th International Symposium on Advanced Optical Manufacturing Chengdu, China, November 19−21, 2008 p728340

    [30]

    Huaying W, Baoqun Z, Wei L 2009 IEEE International Conference on Intelligent Computing and Intelligent Systems, Shanghai, China, November 20−22, 2009 p838

    [31]

    王广俊, 王大勇, 王华英 2010 激光与光电子学进展 47 030901

    Wang G J, Wang D Y, Wang H Y 2010 Laser Optoelectron. Prog. 47 030901

    [32]

    Goodman J W 2005 Introduction To Fourier Optics (Placerville: Roberts and Company Publishers) pp298–306

    [33]

    Lai S, King B, Neifeld M A 2000 Opt. Commun. 173 155

    [34]

    Yamaguchi I, Zhang T 1997 Opt. Lett. 22 1268

    [35]

    Kim M-K 2010 J. Opt. Soc. Korea 14 77

    [36]

    Colomb T, Cuche E, Charrière F, Kühn J, Aspert N, Montfort F, Marquet P, Depeursinge C 2006 Appl. Opt. 45 851

    [37]

    惠梅, 王东生, 李庆祥, 邓年茂, 徐毓娴 2003 光学学报 23 1245

    Hui M, Wang D S, Li Q X, Deng N M, Xu Y X 2003 Acta Optic. Sin. 23 1245

    [38]

    葛宝臻, 崔鹏, 吕且妮, 魏耀林 2008 中国石油大学学报: 自然科学版 32 169

    Ge B Z, Cui P, Lv Q N, Wei Y L 2008 J. China U. Pet.: Nat. Sci. Ed. 32 169

    [39]

    朱娜娜, 金孝岠, 郭建荣 2019 中国糖尿病杂志 27 12

    Zhu N N, Jin X J, Guo J R 2019 Chin. J. Diabetes Mellitus 27 12

    [40]

    Lee S, Park H, Kim K, Sohn Y, Jang S, Park Y 2017 Sci. Rep. 7 1

    [41]

    丁宁, 陈长强, 樊绮诗 2019 检验医学 34 889

    Ding N, Chen C Q, Fan Q S 2019 Lab. Med. 34 889

    [42]

    Li N, Zhou H, Tang Q 2017 Dis. Markers 2017 1

    [43]

    Jenner P 2003 Ann. Neurol. 53 S26

    [44]

    Hald A, Lotharius J 2005 Exp. Neurol. 193 279

    [45]

    Spirito P, Maron B J, Bonow R O 1986 J. Am. Coll. Cardiol. 7 518

    [46]

    Harders A G, Gilsbach J M 1987 J. Neurosurg. 66 718

    [47]

    Sugii Y, Okuda R, Okamoto K, Madarame H 2005 Meas. Sci. Technol. 16 1126

    [48]

    Lima R, Wada S, Takeda M, Tsubota K-i, Yamaguchi T 2007 J. Biomech. 40 2752

    [49]

    Jia Y, Morrison J C, Tokayer J, Tan O, Lombardi L, Baumann B, Lu C D, Choi W, Fujimoto J G, Huang D 2012 Biomed. Opt. Express 3 3127

    [50]

    van den Berg P, Daoudi K, Steenbergen W 2015 Photoacoustics 3 89

    [51]

    Choi Y S, Lee S J 2009 Appl. Opt. 48 2983

    [52]

    Zakrisson J, Schedin S, Andersson M 2015 Appl. Opt. 54 7442

    [53]

    Rinehart M T, Park H S, Wax A 2015 Biomed. Opt. Express 6 2067

    [54]

    Yi F, Lee C-G, Moon I-K 2012 J. Opt. Soc. Korea 16 115

    [55]

    Fällman E, Schedin S, Jass J, Andersson M, Uhlin B E, Axner O 2004 Biosens. Bioelectron. 19 1429

    [56]

    Andersson M, Fällman E, Uhlin B E, Axner O 2006 Biophys. J. 91 2717

    [57]

    Groen F C, Young I T, Ligthart G 1985 Cytometry A 6 81

    [58]

    Sheng J, Malkiel E, Katz J 2006 Appl. Opt. 45 3893

    [59]

    Langehanenberg P, Kemper B, Dirksen D, Von Bally G 2008 Appl. Opt. 47 D176

    [60]

    Moon I, Daneshpanah M, Javidi B, Stern A 2009 Proc. IEEE Inst. Electr. Electron Eng. 97 990

    [61]

    Moon I, Javidi B, Yi F, Boss D, Marquet P 2012 Opt. Express 20 10295

    [62]

    Gonzalez R, Woods R E 2002 Digital Image Processing (New York: Prentice-Hall) pp588−596

    [63]

    Rencher A C 1998 Multivariate Statistical Inference and Applications (New York: Wiley) pp41−63

    [64]

    Masters B R 2012 J. Biomed. Opt. 17 029901

    [65]

    Mir M, Wang Z, Shen Z, Bednarz M, Bashir R, Golding I, Prasanth S G, Popescu G 2011 P NATL ACAD SCI USA 108 13124

    [66]

    Shaked N T, Finan J D, Guilak F, Wax A P 2010 J Biomed Opt 15 010505

    [67]

    Park Y, Diez-Silva M, Popescu G, Lykotrafitis G, Choi W, Feld M S, Suresh S 2008 Proc. Natl. Acad. Sci. 105 13730

    [68]

    Dubois F, Schockaert C, Callens N, Yourassowsky C 2006 Opt. Express 14 5895

    [69]

    Xu L, Mater M, Ni J 2011 Opt. Express 19 14779

    [70]

    Liebling M, Unser M 2004 J. Opt. Soc. Am. A. Opt. Image. Sci. Vis. 21 2424

    [71]

    Hwang O 2013 Exp. Neurobiol. 22 11

    [72]

    Liu J, Zhu L, Zhang F, Dong M, Qu X 2019 Appl. Opt. 58 4042

    [73]

    Felgner H, Müller O, Schliwa M 1995 Appl. Opt. 34 977

    [74]

    Florin E-L, Pralle A, Stelzer E, Hörber J 1998 Appl. Phys. A 66 S75

    [75]

    Schlenker T, Feranchak A P, Schwake L, Stremmel W, Roman R M, Fitz J G 2000 Gastroenterology 118 395

  • 图 1  数字全息原理示意图

    Fig. 1.  Optical layout of digital holography.

    图 2  数字全息显微重建原理图

    Fig. 2.  Optical layout of digital holographic reconstruction.

    图 3  同轴数字全息显微成像系统检测血红细胞形变[52]

    Fig. 3.  Inline digital holographic microscopic system for detecting micro-deformation of RBC[52].

    图 4  RBC形变与重建的复振幅实部信息之间联系的结果分析 (a)不同记录距离z下最佳反向重建距离z', 绿线代表实际实验, 红线代表模拟实验, 插图代表模拟实验所用RBC; (b)绿线代表实际实验RBC的Re(U), 红线与紫线分别代表模拟实验中初始RBC与变形后RBC的Re(U); (c)所用RBC比(b)所用小约20%, 绿线与红线分别代表实际实验与模拟实验所得Re(U)结果; (b)与(c)所用黑线表示RBC的记录距离[52]

    Fig. 4.  The relationship between the RBC deformation and the real part of reconstructed amplitude. (a) The reconstruction distance z' to the focus at different recording distances z of the simulated RBC (red) and the experiment using a real RBC (green). Inset shows the Cassini model of the RBC used in the simulations; (b) the experiment using a real RBC (green), the reconstructed Re(U) for the simulated unaltered RBC and a deformed RBC represented by the red and blue curves, respectively; (c) reconstructed data from a ~20% smaller RBC compared with the one used in (b). Red and green curves represent the Re(U) of simulation and experiment, respectively; Gray vertical line in (b) and (c) indicates position of the RBC[52].

    图 5  同轴数字全息显微成像系统追踪红细胞空间分布[51]

    Fig. 5.  Inline digital holographic microscopic system for tracking spatial distribution of RBCS[51].

    图 6  RBC的全息图与重建图像 (a) CMOS像感器拍摄FEP微管内RBC所得全息图; (b), (c), (d)分别表示不同重建深度下RBC重建图像, 箭头表示聚焦的RBC[51]

    Fig. 6.  The hologram and the reconstruction images of RBCs. (a) The RBC hologram obtained from CMOS; (b), (c), and (d) represent RBC reconstruction images at different reconstruction depths, respectively. Each focused RBC is shown by an arrow[51].

    图 7  离轴数字全息显微成像系统测量不同形状RBC三维体积[54]

    Fig. 7.  Off-axis digital holographic microscopy system for measuring RBCs’ three-dimensional volume of different shapes[54].

    图 8  RBC的相位重建图像 (a)重建的口腔形状RBC相位图像; (b)重建的盘状RBC相位图像; (c)重建后经分水岭算法分割的口腔形状RBC相位图像; (d)重建后经分水岭算法分割的盘状RBC相位图像; (e)分割的单个RBC相位图像; (f)(g)(h)分别经标记分水岭算法进一步分割得到的单个RBC的A、B、C部分[54]

    Fig. 8.  The reconstructed phase image for RBCs (a) The reconstructed phase image for RBCs having a stomatocyte shape; (b) the reconstructed RBCs phase image for RBCs having a discocyte shape; (c)the segmented phase image for RBCs having a stomatocyte shape; (d) the segmented phase image for RBCs having a discocyte shape; (e) the segmented phase image for single RBC(f), (g) and (h) represent the A, B and C parts by the marker-controlled watershed algorithm in RBC, respectively[54].

    图 9  离轴数字全息显微系统研究散焦现象对RBC三维体积测量影响[53]

    Fig. 9.  Off-axis digital holographic microscopy system for investigating the effect of defocus on RBC three-dimensional volume measurement[53].

    图 10  对单一RBC图像进行数字重聚焦与相应的RBC光学体积测量 (a)通过人工聚焦方法与数字重聚焦方法对单一RBC重建所得振幅图与相位图; (b)A—G的振幅方差分布; (c)RBC在人工聚焦方法所得光学体积(黑线)与数字重聚焦方法所得光学体积(蓝线), 光学体积OV表示为平均值 ± 标准差[53]

    Fig. 10.  Digital refocusing of a single red blood cell image and corresponding optical volume measurements. (a) The amplitude and phase images by the manually-focused method and digitally-refocused method from a single RBC; (b) amplitude variance metric of holograms A-G; (c) computed OV of RBC from manually-focused phase images(black) and digitally-refocused phase images(blue). OV reported as mean ± standard deviation[53].

    图 11  光镊辅助离轴数字全息显微系统测量RBC体积[72]

    Fig. 11.  Off-axis digital holographic microscopy system with optical tweezer for measuring RBCs’ three-dimensional volume[72].

    图 12  0−3 pN陷阱拉伸力变化下不同浓度(0−200 μmol/L)氧化应激下重建RBC的高度变化, 颜色深浅代表RBC高度的大小[72]

    Fig. 12.  Height change of reconstructed RBCs under different concentrations of oxidative stress (0−200 μmol/L). Four images in each group are corresponding to trap force varying from 0−3 pN. Color bar represents different thickness[72].

    图 13  RBC在不同浓度氧化应激下的性能 (a)不同浓度氧化应激下RBC最大高度H与陷阱拉伸力关系; (b)不同浓度氧化应激下RBC体积[72]

    Fig. 13.  Performance of RBC under different oxidative stress. (a) The relationship between the maximum height H of RBC and the trap tensile force under different oxidative stress; (b) the volume of RBC under different oxidative stress[72].

    表 1  RBC横向信息与轴向信息的测量精度[51]

    Table 1.  The lateral and axial measurement accuracy of RBC[51].

    坐标测量精度/μm均方根误差/μm
    $x$ ± 0.30.63
    $y$ ± 0.30.52
    $z$ ± 1.02.05
    下载: 导出CSV

    表 2  两种不同形状RBC的A、B部分的三维体积

    Table 2.  The different shapes of RBC’s three-dimensional volume of A and B parts.

    A部分B部分
    平均体积/μm3均方根
    误差
    平均体积/μm3均方根
    误差
    口腔状7.67.241.514.7
    盘状14.68.032.77.3
    下载: 导出CSV

    表 3  人工聚焦方法与数字重聚焦方法测得微球光学体积对比[53]

    Table 3.  Comparison of OV measured by manually-focused and digitally-refocused methods[53].

    测量方法光学体积/fL误差比/%
    人工聚焦方法281.75.90
    数字重聚焦方法266.60.22
    真实值266.0 ± 7.980.00
    下载: 导出CSV
  • [1]

    Cuche E, Marquet P, Depeursinge C 1999 Appl. Opt. 38 6994

    [2]

    Ferraro P, Coppola G, De Nicola S, Finizio A, Pierattini G 2003 Opt. Lett. 28 1257

    [3]

    Ferraro P, Miccio L, Grilli S, Paturzo M, De Nicola S, Finizio A, Osellame R, Laporta P 2007 Opt. Express 15 14591

    [4]

    Charrière F, Marian A, Montfort F, Kuehn J, Colomb T, Cuche E, Marquet P, Depeursinge C 2006 Opt. Lett. 31 178

    [5]

    Kemper B, Von Bally G 2008 Appl. Opt. 47 A52

    [6]

    Potcoava M, Kim M 2008 Meas Sci Technol 19 074010

    [7]

    袁操今, 翟宏琛 2010 光子学报 39 893

    Yuan C J, Zhai H C 2010 Acta Phot. Sin. 39 893

    [8]

    Ma J, Yuan C, Situ G, Pedrini G, Osten W 2013 Chin. Opt. Lett. 11 090901

    [9]

    王华英, 刘飞飞, 成惠, 廖薇, 赵宝群, 于梦杰, 刘佐强 2013 强激光与粒子束 25 345

    Wang H Y, Liu F F, Cheng H, Liao W, Zhao B Q, Yu M J, Liu Z Q 2013 High Pow Las Part Beam 25 345

    [10]

    王华英, 刘飞飞, 宋修法, 廖微, 赵宝群, 于梦杰, 刘佐强 2013 物理学报 62 024207

    Wang H Y, Liu F F, Liao W, Song X F, Yu M J, Liu Z Q 2013 Acta Phys. Sin. 62 024207

    [11]

    Zheng J, Pedrini G, Gao P, Yao B, Osten W 2015 J. Opt. 17 085301

    [12]

    Lin Q, Wang D, Wang Y, Rong L, Chang S 2015 Opt. Lasers Eng. 66 279

    [13]

    Wang Y, Guo S, Wang D, Lin Q, Rong L, Zhao J 2016 Opt. Commun. 366 81

    [14]

    董可平, 钱晓凡, 张磊, 张永安 2007 光子学报 36 2013

    Dong K P, Qian X F, Zhang L, Zhang Y A 2007 Acta Phot. Sin. 36 2013

    [15]

    邸江磊, 赵建林, 范琦, 姜宏振, 孙伟伟 2008 光学学报 28 56

    Di J L, Zhao J L, Fan Q, Jiang H Z, Sun W W 2008 Acta Optic. Sin. 28 56

    [16]

    马利红, 王辉, 金洪震, 李勇 2012 中国激光 39 215

    Ma L H, Wang H, Jin H Z, Li Y 2012 Chin. J. Las. 39 215

    [17]

    Cao R, Xiao W, Wu X, Sun L, Pan F 2018 Biomed. Opt. Express 9 72

    [18]

    Zhang J, Chen Q, Sun J, Zuo C 2019 Seventh International Conference on Optical and Photonic Engineering (icOPEN 2019), Phuket, Thailand, October 16, 2019 p112050 C

    [19]

    Zhang J, Sun J, Chen Q, Zuo C 2020 IEEE Trans. Comput. Imaging. 6 697

    [20]

    Amos W, White J 2003 Biol. Cell 95 335

    [21]

    Subramaniam V, Kirsch A, Jovin T 1998 Cell. Mol. Biol. (Noisy-le-grand) 44 689

    [22]

    Huang D, Swanson E A, Lin C P, Schuman J S, Stinson W G, Chang W, Hee M R, Flotte T, Gregory K, Puliafito C A 1991 Science 254 1178

    [23]

    Kemper B, Carl D, Knoche S, Thien R, von Bally G 2004 Optical Metrology in Production Engineering Strasbourg, France, September 10, 2004 p581

    [24]

    Jin H, Wan H, Zhang Y, Li Y, Qiu P 2008 J. Mod. Opt. 55 2989

    [25]

    周海宪 2006 应用光学 27 I0007

    Zhou H X 2006 J. Appl. Opt. 27 I0007

    [26]

    Nehmetallah G T, Aylo R, Williams L 2015 Analog and Digital Holography with MATLAB (Bellingham: Society of Photo-Optical Instrumentation Engineers (SPIE)) pp137–143

    [27]

    Kreis T M, Adams M, Jüptner W P 1997 Optical Inspection and Micromeasurements II Munich, Germany, September 17, 1997 p224

    [28]

    Liu C, Wang D, Zhang Y 2009 Opt. Eng. 48 105802

    [29]

    Song X, Zhao B, Qin A 4th International Symposium on Advanced Optical Manufacturing Chengdu, China, November 19−21, 2008 p728340

    [30]

    Huaying W, Baoqun Z, Wei L 2009 IEEE International Conference on Intelligent Computing and Intelligent Systems, Shanghai, China, November 20−22, 2009 p838

    [31]

    王广俊, 王大勇, 王华英 2010 激光与光电子学进展 47 030901

    Wang G J, Wang D Y, Wang H Y 2010 Laser Optoelectron. Prog. 47 030901

    [32]

    Goodman J W 2005 Introduction To Fourier Optics (Placerville: Roberts and Company Publishers) pp298–306

    [33]

    Lai S, King B, Neifeld M A 2000 Opt. Commun. 173 155

    [34]

    Yamaguchi I, Zhang T 1997 Opt. Lett. 22 1268

    [35]

    Kim M-K 2010 J. Opt. Soc. Korea 14 77

    [36]

    Colomb T, Cuche E, Charrière F, Kühn J, Aspert N, Montfort F, Marquet P, Depeursinge C 2006 Appl. Opt. 45 851

    [37]

    惠梅, 王东生, 李庆祥, 邓年茂, 徐毓娴 2003 光学学报 23 1245

    Hui M, Wang D S, Li Q X, Deng N M, Xu Y X 2003 Acta Optic. Sin. 23 1245

    [38]

    葛宝臻, 崔鹏, 吕且妮, 魏耀林 2008 中国石油大学学报: 自然科学版 32 169

    Ge B Z, Cui P, Lv Q N, Wei Y L 2008 J. China U. Pet.: Nat. Sci. Ed. 32 169

    [39]

    朱娜娜, 金孝岠, 郭建荣 2019 中国糖尿病杂志 27 12

    Zhu N N, Jin X J, Guo J R 2019 Chin. J. Diabetes Mellitus 27 12

    [40]

    Lee S, Park H, Kim K, Sohn Y, Jang S, Park Y 2017 Sci. Rep. 7 1

    [41]

    丁宁, 陈长强, 樊绮诗 2019 检验医学 34 889

    Ding N, Chen C Q, Fan Q S 2019 Lab. Med. 34 889

    [42]

    Li N, Zhou H, Tang Q 2017 Dis. Markers 2017 1

    [43]

    Jenner P 2003 Ann. Neurol. 53 S26

    [44]

    Hald A, Lotharius J 2005 Exp. Neurol. 193 279

    [45]

    Spirito P, Maron B J, Bonow R O 1986 J. Am. Coll. Cardiol. 7 518

    [46]

    Harders A G, Gilsbach J M 1987 J. Neurosurg. 66 718

    [47]

    Sugii Y, Okuda R, Okamoto K, Madarame H 2005 Meas. Sci. Technol. 16 1126

    [48]

    Lima R, Wada S, Takeda M, Tsubota K-i, Yamaguchi T 2007 J. Biomech. 40 2752

    [49]

    Jia Y, Morrison J C, Tokayer J, Tan O, Lombardi L, Baumann B, Lu C D, Choi W, Fujimoto J G, Huang D 2012 Biomed. Opt. Express 3 3127

    [50]

    van den Berg P, Daoudi K, Steenbergen W 2015 Photoacoustics 3 89

    [51]

    Choi Y S, Lee S J 2009 Appl. Opt. 48 2983

    [52]

    Zakrisson J, Schedin S, Andersson M 2015 Appl. Opt. 54 7442

    [53]

    Rinehart M T, Park H S, Wax A 2015 Biomed. Opt. Express 6 2067

    [54]

    Yi F, Lee C-G, Moon I-K 2012 J. Opt. Soc. Korea 16 115

    [55]

    Fällman E, Schedin S, Jass J, Andersson M, Uhlin B E, Axner O 2004 Biosens. Bioelectron. 19 1429

    [56]

    Andersson M, Fällman E, Uhlin B E, Axner O 2006 Biophys. J. 91 2717

    [57]

    Groen F C, Young I T, Ligthart G 1985 Cytometry A 6 81

    [58]

    Sheng J, Malkiel E, Katz J 2006 Appl. Opt. 45 3893

    [59]

    Langehanenberg P, Kemper B, Dirksen D, Von Bally G 2008 Appl. Opt. 47 D176

    [60]

    Moon I, Daneshpanah M, Javidi B, Stern A 2009 Proc. IEEE Inst. Electr. Electron Eng. 97 990

    [61]

    Moon I, Javidi B, Yi F, Boss D, Marquet P 2012 Opt. Express 20 10295

    [62]

    Gonzalez R, Woods R E 2002 Digital Image Processing (New York: Prentice-Hall) pp588−596

    [63]

    Rencher A C 1998 Multivariate Statistical Inference and Applications (New York: Wiley) pp41−63

    [64]

    Masters B R 2012 J. Biomed. Opt. 17 029901

    [65]

    Mir M, Wang Z, Shen Z, Bednarz M, Bashir R, Golding I, Prasanth S G, Popescu G 2011 P NATL ACAD SCI USA 108 13124

    [66]

    Shaked N T, Finan J D, Guilak F, Wax A P 2010 J Biomed Opt 15 010505

    [67]

    Park Y, Diez-Silva M, Popescu G, Lykotrafitis G, Choi W, Feld M S, Suresh S 2008 Proc. Natl. Acad. Sci. 105 13730

    [68]

    Dubois F, Schockaert C, Callens N, Yourassowsky C 2006 Opt. Express 14 5895

    [69]

    Xu L, Mater M, Ni J 2011 Opt. Express 19 14779

    [70]

    Liebling M, Unser M 2004 J. Opt. Soc. Am. A. Opt. Image. Sci. Vis. 21 2424

    [71]

    Hwang O 2013 Exp. Neurobiol. 22 11

    [72]

    Liu J, Zhu L, Zhang F, Dong M, Qu X 2019 Appl. Opt. 58 4042

    [73]

    Felgner H, Müller O, Schliwa M 1995 Appl. Opt. 34 977

    [74]

    Florin E-L, Pralle A, Stelzer E, Hörber J 1998 Appl. Phys. A 66 S75

    [75]

    Schlenker T, Feranchak A P, Schwake L, Stremmel W, Roman R M, Fitz J G 2000 Gastroenterology 118 395

  • [1] 周宏强, 万玉红, 满天龙. 基于位相变更的非相干数字全息自适应成像. 物理学报, 2018, 67(4): 044202. doi: 10.7498/aps.67.20172202
    [2] 王大勇, 王云新, 郭莎, 戎路, 张亦卓. 基于多角度无透镜傅里叶变换数字全息的散斑噪声抑制成像研究. 物理学报, 2014, 63(15): 154205. doi: 10.7498/aps.63.154205
    [3] 周光照, 佟亚军, 陈灿, 任玉琦, 王玉丹, 肖体乔. 相干X射线衍射成像的数字模拟研究. 物理学报, 2011, 60(2): 028701. doi: 10.7498/aps.60.028701
    [4] 周光照, 王玉丹, 任玉琦, 陈灿, 叶琳琳, 肖体乔. 相干X射线衍射成像三维重建的数字模拟研究. 物理学报, 2012, 61(1): 018701. doi: 10.7498/aps.61.018701
    [5] 陈萍, 唐志列, 王娟, 付晓娣, 陈飞虎. 用Stokes参量法实现数字同轴偏振全息的研究. 物理学报, 2012, 61(10): 104202. doi: 10.7498/aps.61.104202
    [6] 王林, 袁操今, 聂守平, 李重光, 张慧力, 赵应春, 张秀英, 冯少彤. 数字全息术测定涡旋光束拓扑电荷数. 物理学报, 2014, 63(24): 244202. doi: 10.7498/aps.63.244202
    [7] 王淑莹, 章海军, 张冬仙. 基于微球透镜的任选区高分辨光学显微成像新方法研究. 物理学报, 2013, 62(3): 034207. doi: 10.7498/aps.62.034207
    [8] 李俊昌. 数字全息重建图像的焦深研究. 物理学报, 2012, 61(13): 134203. doi: 10.7498/aps.61.134203
    [9] 崔华坤, 王大勇, 王云新, 刘长庚, 赵洁, 李艳. 无透镜傅里叶变换数字全息术中非共面误差的自动补偿算法. 物理学报, 2011, 60(4): 044201. doi: 10.7498/aps.60.044201
    [10] 徐先锋, 韩立立, 袁红光. 两步相移数字全息物光重建误差分析与校正. 物理学报, 2011, 60(8): 084206. doi: 10.7498/aps.60.084206
    [11] 马骏, 袁操今, 冯少彤, 聂守平. 基于数字全息及复用技术的全场偏振态测试方法. 物理学报, 2013, 62(22): 224204. doi: 10.7498/aps.62.224204
    [12] 谷婷婷, 黄素娟, 闫成, 缪庄, 常征, 王廷云. 基于数字全息图的光纤折射率测量研究. 物理学报, 2015, 64(6): 064204. doi: 10.7498/aps.64.064204
    [13] 李俊昌, 楼宇丽, 桂进斌, 彭祖杰, 宋庆和. 数字全息图取样模型的简化研究. 物理学报, 2013, 62(12): 124203. doi: 10.7498/aps.62.124203
    [14] 刘双龙, 刘伟, 陈丹妮, 屈军乐, 牛憨笨. 相干反斯托克斯拉曼散射显微成像技术研究. 物理学报, 2016, 65(6): 064204. doi: 10.7498/aps.65.064204
    [15] 刘诚, 潘兴臣, 朱健强. 基于光栅分光法的相干衍射成像. 物理学报, 2013, 62(18): 184204. doi: 10.7498/aps.62.184204
    [16] 李俊昌, 张亚萍, 许蔚. 高质量数字全息波面重建系统研究. 物理学报, 2009, 58(8): 5385-5391. doi: 10.7498/aps.58.5385
    [17] 卢明峰, 吴坚, 郑明. 数字全息周期像的产生机理及在抑制零级衍射上的应用. 物理学报, 2013, 62(9): 094207. doi: 10.7498/aps.62.094207
    [18] 王建国, 杨松林, 叶永红. 样品表面银膜的粗糙度对钛酸钡微球成像性能的影响. 物理学报, 2018, 67(21): 214209. doi: 10.7498/aps.67.20180823
    [19] 谢静, 张军勇, 岳阳, 张艳丽. 卢卡斯光子筛的聚焦特性研究. 物理学报, 2018, 67(10): 104201. doi: 10.7498/aps.67.20172260
    [20] 王华英, 刘飞飞, 廖薇, 宋修法, 于梦杰, 刘佐强. 优化的数字全息显微成像系统. 物理学报, 2013, 62(5): 054208. doi: 10.7498/aps.62.054208
  • 引用本文:
    Citation:
计量
  • 文章访问数:  1298
  • PDF下载量:  104
  • 被引次数: 0
出版历程
  • 收稿日期:  2020-03-11
  • 修回日期:  2020-04-09
  • 上网日期:  2020-05-20
  • 刊出日期:  2020-08-20

基于数字全息的血红细胞显微成像技术

  • 1. 中国计量大学光学与电子科技学院, 杭州 310018
  • 2. 浙江省现代计量测试技术及仪器重点实验室, 杭州 310018
  • 3. 精密测试技术及仪器国家重点实验室, 清华大学精密仪器系, 北京 100084
  • 通信作者: 金尚忠, jinsz@cjlu.edu.cn ; 曹良才, clc@mail.tsinghua.edu.cn

摘要: 数字全息显微成像有别于传统光学显微成像, 可根据重建全息图获取细胞的生物学参数与形貌信息, 是一种有效的非接触无损三维成像技术. 随着像感器的发展与硬件计算能力的提升, 数字全息显微成像技术在活体生物细胞检测尤其在血红细胞检测领域取得了显著进展和突破. 本文介绍了同轴、离轴以及光镊辅助离轴的数字全息显微技术, 这些技术利用瑞利索末菲反向传播算法、清晰度量化算法、分水岭分割算法、数字重聚焦方法与热涨落方法等来实现血红细胞的形变、空间分布、三维体积信息的高精度提取, 有助于糖尿病、心血管疾病、帕金森氏疾病等病理研究. 数字全息显微成像技术实现了传统三维显微成像技术难以达到的实时性和定量化检测, 由于独有的非接触、无损性特点, 在细胞成像领域应用前景广阔.

English Abstract

    • 数字全息显微技术(digital holographic microscopy, DHM)结合数字全息与显微成像技术, 充分利用两者优势, 不仅能够对待测样本进行放大, 而且能实时记录待测样本的三维信息. 1999年 E. Cuche课题组在待测样本与CCD(charge coupled device)像感器之间放置显微物镜[1], 用于放大待测样本的物光波, 仅基于单幅全息图便可获取有关样品三维结构的精确定量信息, 降低了对CCD分辨率的要求, 实现了高分辨率显微成像. 但由于显微物镜的小孔径限制, 全息图部分高频成分在记录过程中丢失, 重建像质量受到影响, 因此数字全息显微技术仍存在需要改进的地方. 2003年 Ferraro等[2,3]研究了自动跟踪聚焦的数字全息显微镜, 并在随后的研究中提出了运用多波长数字全息技术来校正色差. 2006年Charrière等[4]将数字全息显微技术与电子计算机断层扫描技术(computed tomography, CT)相结合, 通过旋转花粉细胞获取其各方向的相位图分布, 根据滤波反投影算法实现了生物样品的三维成像. 2007年Kemper等[5]将数字全息显微技术与生物医学相结合, 重建了人体肝脏肿瘤细胞与血红细胞(red blood cell, RBC)的三维形貌, 进一步拓展了数字全息显微的实践应用. 2008年Kim等[6]通过数字全息显微术结合生物医学, 揭示了血管的形貌以及视网膜层的光学厚度分布. 国内围绕数字全息成像分辨率提升与生物细胞定量相位测量开展了大量的工作. 在分辨率提升方面, 2010年翟宏琛等[7]利用一种相位掩模板进行多角度照明, 实现了综合孔径数字全息超分辨率成像; 2013年马俊等[8]基于结构光照明的数字全息显微技术, 通过分离和合成傅里叶域中不同频率区域的带宽, 可以获得更高空间分辨率的重建图像; 同年王华英等[9,10]综合分析了数字全息显微成像系统的分辨率影响因素, 对像面数字全息系统展开了讨论与研究; 2015年姚保利等[11]利用散斑照明获得重建物波中的合成数值孔径, 以提高分辨率; 同年王大勇等[12,13]利用动态空间光栅提高了数字全息成像分辨率, 并通过微球体与像面数字全息相结合, 实现了分辨率增强的相衬成像. 在上述研究中, 由于分辨率得到了大幅提升, 使得数字全息显微更有效地分辨生物微观结构如人体细胞与亚细胞结构等. 在数字全息定量相位测量方面, 2007年董可平等[14]成功实现了生物细胞的定量相位测量; 2008年赵建林等[15]通过拍摄两幅全息图并利用相位相减技术实现了活体细胞组织等相位型生物样本的定量测量和有效观察; 2012年马利红等[16]建立了一套预放大式数字全息显微成像系统, 实现了高分辨率的细胞定量相位测量; 2018年肖文等[17]利用离轴马赫泽德实验系统对细胞进行实时成像, 证明了细胞质中肌动蛋白对细胞骨架具有更大的支撑作用. 2019年左超等[18,19]通过数字全息显微结合自适应松弛超像素分辨技术, 在利用相位恢复与数字重聚焦等算法或方法的基础上, 开发了无透镜数字全息显微系统, 集成了定量相位、剖面分析、三维显示和细胞计数等多项功能, 具有无标记、高通量、低成本和微型化的特点, 为面向临床检验应用的显微设备提供了有效的技术支撑. 基于数字全息的定量相位成像技术满足了活体细胞组织等相位型生物样本定量测量和动态观察的需求, 对于生物医学领域的发展具有重要意义.

      相较于当前细胞成像领域广泛应用的激光扫描共聚焦显微(confocal laser scanning microscope, CLSM)技术[20]、近场光学显微(near-field optical microscopy, NSOM)技术[21]和光学相干层析成像(optical coherence tomography, OCT)技术[22]等, 数字全息显微技术具备以下优点: (1)属于面成像, 在成像速度上极具优势, 仅根据单次曝光即可实时记录待测样本完整三维信息[5,23], 利于生物活体组织的动态变化的研究, 传统方法主要基于点扫描成像; (2)属于宽光束照射, 对光功率密度要求较低的同时无需荧光标记, 实现了无扰分析, 而激光扫描共聚焦需在样品上实现激光束聚焦, 对焦点处的光功率密度要求较高, 或需要对样品进行荧光标记, 无法实现无扰分析; (3)可以得到细胞折射率分布信息, 对比当前CT技术具有更高的空间分辨率, 且对比度高, 可实现自动聚焦, 近场光学显微技术仅能实现表面成像, 而光学相干层析成像技术虽能得到生物体组织内部的折射率分布信息, 其空间分辨率约为5 μm, 成像深度在毫米量级, 需要进一步改进才能观察人体细胞和亚细胞结构(1 μm左右). 通过数字全息显微技术能够实现宽视场亚微米级空间分辨率(0.5 μm左右), 不仅可以清晰分辨人体细胞与亚细胞结构, 并且可实现定量、三维和快速追踪的成像, 在细胞成像领域上具有独特的应用价值.

      本文从生物医学细胞成像应用需求出发, 介绍了数字全息显微成像的原理, 综述了同轴、离轴以及光镊辅助离轴数字全息显微系统如何对RBC进行形貌信息提取, 重点分析了瑞利索末菲反向传播算法、清晰度量化算法、分水岭分割算法、数字重聚焦方法与热涨落方法在研究RBC的微形变、三维体积测量与空间分布等方面的应用, 概述了当前数字全息显微成像的发展与应用前景.

    • 数字全息显微光路如图1所示, 其中Object plane代表物平面${x_0} \text{-} {y_0}$, MO plane代表显微物镜平面${x_{\rm{m}}} \text{-} {y_{\rm{m}}}$, Hologram plane代表全息记录平面$x \text{-} y$, Image plane代表透镜成像平面${x_{\rm{i}}} \text{-} {y_{\rm{i}}}$, 样品放置于距离显微物镜平面${z_0}$距离处[24], 显微物镜平面与全息记录面距离为${z_{\rm{m}}}$, 全息记录面与透镜成像平面的距离为$z$, ${{R}}$作为参考光, 与经显微物镜放大后的物光${{O}}$于全息记录面上发生干涉, 形成明暗相间的干涉条纹, 这些条纹记录了物体光波所有振幅与相位信息, 是一个编码和调制的过程.

      图  1  数字全息原理示意图

      Figure 1.  Optical layout of digital holography.

      设物光在物平面的光场分布为${{{O}}_0}\left( {{x_0}, {y_0}} \right)$, 根据菲涅尔衍射近似条件, 物光波经过距离${z_0}$衍射至显微物镜前表面时, 其光场复振幅分布${{{O}}_{\rm{f}}}\left( {{x_{\rm{m}}}, {y_{\rm{m}}}} \right)$

      $\begin{split}\; &{{{O}}_{\rm{f}}}\left( {{x_{\rm{m}}},{y_{\rm{m}}}} \right) = \frac{{\exp \left( {{\rm{j}}k{z_0}} \right)}}{{{\rm{j}}\lambda {z_0}}}\iint {{{{O}}_0}}\left( {{x_0},{y_0}} \right)\\ &~~ \times\exp\left[ {{\rm{j}}k\frac{{{{\left( {{x_{\rm{m}}} - {x_0}} \right)}^2} + {{\left( {{y_{\rm{m}}} - {y_0}} \right)}^2}}}{{2{z_0}}}} \right]{\rm{d}}{x_0}{\rm{d}}{y_0},\end{split}$

      其中, $\lambda $表示物光波长, ${x_{\rm{0}}}$${y_{\rm{0}}}$表示物平面上的空间坐标, ${x_{\rm{m}}}$${y_{\rm{m}}}$表示显微物镜平面上的空间坐标, ${z_{\rm{0}}}$表示物平面与显微物镜平面的距离. 由于显微物镜在全息记录过程中仅作为放大作用的透镜, 依据透镜的相位调制特性, 物光波经过显微物镜传播至其后表面时, 其光场复振幅分布${{{O}}_{\rm{b}}}\left( {{x_m}, {y_m}} \right)$

      ${{{O}}_{\rm{b}}}( {{x_{\rm{m}}},{y_{\rm{m}}}}) \!=\! \exp \left[ { - \frac{{{\rm{j}}k}}{{2f}}\left( {x_{\rm{m}}^2 + y_{\rm{m}}^2} \right)} \right] \cdot {{{O}}_{\rm{f}}}( {{x_{\rm{m}}},{y_{\rm{m}}}} ), $

      其中, $f$为显微物镜的焦距, 物光波由显微物镜后表面经过距离${z_m}$传播至全息记录面时, 其光场复振幅分布${{O}}(x, y)$

      $ \begin{split} {{O}}(x,y) =\; & \frac{{\exp \left( {{\rm{j}}k{z_{\rm{m}}}} \right)}}{{{\rm{j}}\lambda {z_{\rm{m}}}}}\iint {{{{O}}_{\rm{b}}}}\left( {{x_{\rm{m}}},{y_{\rm{m}}}} \right)\\ & \times\exp \left\{\!{\frac{{{\rm{j}}k}}{{2{z_{\rm{m}}}}}\!\left[\!{{{\left( {x \!-\! {x_{\rm{m}}}} \right)}^2} \!+\! {{\left( {y\! -\! {y_{\rm{m}}}} \right)}^2}} \right]{\rm{d}}{x_{\rm{m}}}{\rm{d}}{y_{\rm{m}}}} \right.\\ =\; & - \frac{{\exp \left[ {{\rm{j}}k\left( {{z_0} + {z_{\rm{m}}}} \right)} \right]}}{{{\lambda ^2}{z_0}{z_{\rm{m}}}}}\iint\!\! \iint {{{O}}_{\rm{0}}}\left( {{x_0},{y_0}} \right)\\ & \times\exp \left[ { - \frac{{{\rm{j}}k}}{{2f}}\left( {x_{\rm{m}}^2 + y_{\rm{m}}^2} \right)} \right] \\ & \times\exp \left[ {{\rm{j}}k\frac{{{{\left( {{x_{\rm{m}}} - {x_0}} \right)}^2} + {{\left( {{y_{\rm{m}}} - {y_0}} \right)}^2}}}{{2{z_0}}}} \right] \\ & \times \exp \left[ {{\rm{j}}k\frac{{{{\left( {x - {x_{\rm{m}}}} \right)}^2} + {{\left( {y - {y_{\rm{m}}}} \right)}^2}}}{{2{z_{\rm{m}}}}}} \right]\\ & \times{\rm{d}}{x_{\rm{m}}}{\rm{d}}{y_{\rm{m}}}{\rm{d}}{x_0}{\rm{d}}{y_0}.\\[-10pt] \end{split} $

      根据透镜成像公式$\dfrac{1}{f} = \dfrac{1}{{{z_0}}} + \dfrac{1}{{{z_{\rm{m}}}}}$和显微物镜的放大倍率$M = - {{{z_{\rm{m}}}}}/{{{z_0}}}$, 可得到:

      $\begin{split}{{O}}\left( {x,y} \right) =\; & \frac{{\exp \left[ {{\rm{j}}k\left( {{z_0} + {z_{\rm{m}}}} \right)} \right]}}{M}\exp \left( { - {\rm{j}}k\frac{{{x^2} + {y^2}}}{{2Mf}}} \right) \\ & \cdot {{{O}}_{\rm{0}}}\left( {\frac{x}{M},\frac{y}{M}} \right), \end{split}$

      由于显微物镜的放大作用, 全息记录面记录的物光场相对原物光场存在坐标缩放的关系. 同理, 参考光传播至全息记录面的光学路径与物光类似, 参考光到全息记录面的距离与物光到全息记录面的距离相等, 设参考光在全息记录面的光场分布为${{R}}(x, y)$, 物光和参考光在全息记录面上发生干涉形成全息图, 全息图的光场分布为

      ${{H}}(x,y) = {{O}}(x,y) + {{R}}(x,y),$

      则全息记录面上全息图的强度分布为其光场分布的平方, 其强度分布为

      $\begin{split}{{I}}(x,y) =\;& |{{O}}(x,y){|^2} + |{{R}}(x,y){|^2} \\ &+ {{{R}}^*}(x,y){{O}}(x,y) \!+\! {{R}}(x,y){{{O}}^*}(x,y),\end{split}$

      式中, 第一项表示物光的强度分布, 第二项表示参考光的强度分布, 第一项与第二项合称为零级项或直流项; 第三项表示物光场的复振幅分布, 包含了物光波波前的振幅信息与相位信息, 进行数值重建时可获得待测样本的真实物像[25]; 第四项表示物光场的的共轭复振幅分布, 数值重建时会产生孪生像, 这是全息图重建时获得重建像的最大干扰项, 需要在重建过程中通过算法进行分离.

    • 数字全息重建通过计算机数值模拟重建光波经过数字全息图后的衍射过程, 实现数值重建, 在数字全息显微成像系统中重建的并非物光本身, 而是由显微物镜放大后所成的物光场, 根据近似条件不同可将模拟衍射成像方法分为菲涅尔重建法、卷积重建法、角谱重建法[26]. 三种重建法均基于菲涅尔-基尔霍夫衍射公式进行近似处理或变换后所得. 菲涅尔重建法重建图像的像素尺寸$\delta x$与重建距离z满足关系: $\delta x = {{\lambda z} / {{X_0}}}$, 其中${X_0}$为全息图尺寸. 这意味着重建图像的大小会随着重建距离增加, 而采用卷积重建法和角谱重建法所得图像像素尺寸始终与全息图像素尺寸相同. Kreis等[27]采用菲涅尔重建法与卷积重建法分别对CCD记录下的骰子全息图进行了重建, 其中骰子位于CCD前1.054 m处, 且骰子尺寸大于CCD尺寸. 采用菲涅尔重建法可直接重建完整的骰子图像, 而采用卷积重建法需要平移CCD采集多幅全息图或者补零操作才能完整重建骰子图像, 因此菲涅尔重建法适用于尺寸大于像感器的待测样本. 王大勇等[28]分别利用菲涅尔重建法、卷积重建法与角谱重建法对骰子全息图进行重建, 比较了三种方法的重建速度, 其中菲涅尔重建法的执行时间为0.97 s, 速度最快, 而卷积重建法的执行时间为2.36 s, 速度最慢, 另外角谱重建法的执行时间为1.22 s. 菲涅尔重建法仅通过一次正傅里叶变换即获得重建像; 卷积重建法由于其传播函数在空域中表示, 需要经过两次正傅里叶变换和一次逆傅里叶变换获得重建像; 角谱重建法与卷积重建法原理类似, 其传播函数在频域中表示, 需要经过一次正傅里叶变换和一次逆傅里叶变换获得重建像, 通过比较像质得到角谱重建法的准确性优于卷积重建法. 而根据USAF鉴别率板的重建, 得到在卷积重建与角谱重建中均存在最优重建距离, 当重建距离不同于最优距离时, 重建像的分辨率会下降, 此时角谱重建法的结果优于卷积重建法[29,30]. 卷积重建法仅能在最佳重建距离附近较小的范围内获得高分辨率重建像, 适用性上存在不足, 而角谱重建法重建速度较快, 在最佳重建距离附近较大范围都能获得高分辨率的重建像[31].

      下面以卷积重建法与角谱重建法为例, 介绍生物细胞全息图的重建, 卷积重建过程由记录全息图的光场分布${{H}}(x, y)$以及卷积重建空域内的传播函数$h\left( {{x_i} - x, {y_i} - y} \right)$构成, 通过快速傅里叶变换(fast Fourier transform, FFT)进行运算:

      ${\varGamma} \left( {{x_i},{y_i}} \right) = {{\cal{F}}^{ - 1}}\left[ {{\cal{F}}\{ {{H}}(x,y)\} {\cal{F}}\left\{ {h\left( {{x_i} - x,{y_i} - y} \right)} \right\}} \right]$

      而角谱重建过程由光场分布${{H}}(x, y)$与角谱重建频域内的传播函数$G\left( {{f_{{x_{\rm{i}}}}}, {f_{{y_{\rm{i}}}}}} \right)$构成, 通过FFT进行运算:

      ${{\varGamma}} \left( {{x_{\rm{i}}},{y_{\rm{i}}}} \right) = {{\cal{F}}^{ - 1}}\left[ {{\cal{F}}\{ {{H}}(x,y)\} \cdot G\left( {{f_{{x_{\rm{i}}}}},{f_{{y_{\rm{i}}}}}} \right)} \right]$

      其中, ${\cal{F}}$${{\cal{F}}^{ - 1}}$分别表示快速傅里叶变换与快速傅里叶逆变换, 如图2所示, $x, y$表示全息记录面上的空间坐标, ${x_{\rm{i}}}$${y_{\rm{i}}}$表示重建成像平面上的空间坐标, ${f_{{x_{\rm{i}}}}}$${f_{{y_{\rm{i}}}}}$表示成像面上$({x_{\rm{i}}}, {y_{\rm{i}}})$对应的频率坐标, ${{K}}\left( {{f_x}, {f_y}} \right) = {\cal{F}}\{ {{H}}(x, y)\} $表示在全息记录面上的全息图频谱.

      图  2  数字全息显微重建原理图

      Figure 2.  Optical layout of digital holographic reconstruction.

      重建图像的强度可根据$\left| {{{\varGamma}} \left( {{x_{\rm{i}}}, {y_{\rm{i}}}} \right)} \right|$获得, 由角谱分析方法可对(8)式进行进一步简化[32], 其表达式如下:

      $G\left( {{f_{{x_{\rm{i}}}}},{f_{{y_{\rm{i}}}}}} \right) = \exp \{ {\rm{i}}kz\sqrt {1 - {{\left( {\lambda {f_{{x_{\rm{i}}}}}} \right)}^2} - {{\left( {\lambda {f_{{y_{\rm{i}}}}}} \right)}^2}} \} , $

      其中, $\lambda $表示光源波长, $z$表示全息记录面与成像面的距离. 可通过(7)式对记录面所得全息图进行卷积重建或结合(8)式与(9)式对记录面所得全息图进行角谱重建, 获得经显微物镜放大后待测样本的数值重建像.

      对于同轴全息系统, 由于重建像与孪生像重叠, 需要在数字全息记录过程中引入相移装置, 再根据相移法将所需重建像与孪生像进行分离, 提取出待测样本的物像[33,34]; 对于离轴全息系统, 选择适当的物光与参考光的入射夹角, 使物光与参考光以一定夹角入射至像感器记录面上, 产生干涉全息图, 通过滤波以滤除孪生像, 提取出待测样本的物像. 在离轴全息系统中, 由于物光与参考光存在夹角, 会产生一次项附加相位(呈现斜面分布), 而无论是离轴还是同轴全息系统, 当RBC等弱散射物体的相位信息加载于经显微物镜后形成的球面波上时, 均会引入二次项附加相位(呈现球面分布), 这些附加相位所产生的调制会使待测样本自身的相位分布难以恢复, 引起包裹相位的欠采样, 最终影响到相位解包裹的准确性, 可在参考光路放置与物光光路完全相同的显微物镜进行补偿[35], 或通过数据拟合进行自动相位校正[36], 或采用前后两次拍摄全息图进行相位相减来实现相位校正[15], 此时通过同轴或离轴系统获得的校正相位分布在[–π, π]之间, 还需要进行相位解包裹以实现待测样本相位图像的恢复[37,38].

      为有效测量RBC的微形变, 三维体积及其空间分布, 当前相关RBC测量的数字全息显微系统在重建时主要采用卷积重建法或角谱重建法.

    • RBC是血液中的主要有形成分, 约占血液体积的50%, 生存环境的改变会引起RBC病理性变化, 如糖尿病患者体内的高血糖环境会导致RBC双凹圆盘状丧失、形态拉长、体积减小、变形性降低[39,40], 心血管疾病患者体内RBC体积分布宽度(RDW)升高[41,42], 帕金森氏疾病患者体内存在氧化应激, 氧化应激会引起体内RBC功能性障碍与凋亡[43,44], 因此如何提取RBC的生物学参数与形貌信息将有助于研究人体的糖尿病、心血管疾病、帕金森氏疾病等. 常规的RBC血液测量技术血管造影与超声多普勒分析所提供的血管形状信息与血流速度信息的空间分辨率较差[45,46]; 粒子图像测速(particle image velocimetry, PIV)在血液动力学信息的测量得到广泛应用, 但其仅能在较浅的景深中固有地提供二维平面信息, 在应用范围上存在较大缺陷[47,48]; 光学相干层析成像在血流动力学领域中得到广泛应用, 能实现生物组织微观结构显示, 但RBC的对比度较低[49]; 光声成像(photoacoustic imaging, PAI)能同时实现灌注成像与血流成像, 且RBC的对比度较高, 但其难以完成单个RBC的动态追踪[50]. 本节综述了同轴、离轴以及光镊辅助离轴数字全息显微技术, 这些技术利用瑞利索末菲反向传播算法、清晰度量化算法、分水岭分割算法、数字重聚焦方法等不仅可以获取高分辨率的血管形状与血流速度信息[51], 进一步地还能够检测单个RBC的形状变化[52], 甚至可以完成RBC动态追踪, 实时测量RBC三维体积, 实现无损3D成像[53,54], 对于人体的糖尿病、心血管疾病、帕金森氏疾病等病理研究具有重大意义.

    • 为实现以数字全息显微成像技术为基础, 高效率、低计算强度、低成本地分析细胞的形态变化并跟踪轴向位置提取三维坐标, 瑞典于默奥大学Andersson等[52]设计了一套同轴数字全息系统并提出了一种方法, 通过瑞利索末菲反向传播算法重建复振幅实部信息, 精确确定RBC轴向位置并提供几何形状信息. 其光学系统的分析与实现可根据图3所述结构进行描述, 光学系统围绕Olympus IX71显微镜引入相关光学元件, 采用低成本发光二极管(470 nm, M470 L3-C1)发射LED光[55], 多模光纤安装于压电控制显微镜载物台(PS)上方[56], 载物台上放置单个RBC样品, 其通过人体血液与500 μL的$\rm pH = 7.4$磷酸盐缓冲溶液混合并离心提取制备而成, 多模光纤可以对LED光进行频率滤波, 稳定输出光源, 经过RBC发生衍射的光作为物光, 而未经过待测样本的光作为参考光, 两束光均经过浸水物镜MO(60×/1.40NA)放大, 在CCD像感器上发生相干叠加, 记录下RBC全息图, 其中显微镜的压电控制载物台可进行轴向移动以精确控制物面与记录面的距离$z$, 以10 μm为记录间隔分别记录z = 30—100 μm共8幅RBC全息图像. 根据同轴全息重建算法分别对8幅RBC全息图进行重建, 计算各RBC全息图像在不同重建距离$z'$下的归一化强度, 认定归一化强度最大位置即几何焦点位置, 如图4(a)所示, 插图表示仿真实验使用6.3 μm宽, 折射率$n = 1.40$的RBC形状的Cassini模型, 绿线与红线分别代表实际实验与仿真实验中RBC在不同记录距离z = 30—100 μm下对应的最佳反向重建距离, 可以看出, 在不同的记录距离z下, 通过分别计算RBC全息重建图的归一化强度最大值, 可以精确地得到8个记录距离下对应的几何焦点位置, 作为最佳反向重建距离, 仿真实验与实际实验所得到的反向重建距离曲线一致.

      图  3  同轴数字全息显微成像系统检测血红细胞形变[52]

      Figure 3.  Inline digital holographic microscopic system for detecting micro-deformation of RBC[52].

      图  4  RBC形变与重建的复振幅实部信息之间联系的结果分析 (a)不同记录距离z下最佳反向重建距离z', 绿线代表实际实验, 红线代表模拟实验, 插图代表模拟实验所用RBC; (b)绿线代表实际实验RBC的Re(U), 红线与紫线分别代表模拟实验中初始RBC与变形后RBC的Re(U); (c)所用RBC比(b)所用小约20%, 绿线与红线分别代表实际实验与模拟实验所得Re(U)结果; (b)与(c)所用黑线表示RBC的记录距离[52]

      Figure 4.  The relationship between the RBC deformation and the real part of reconstructed amplitude. (a) The reconstruction distance z' to the focus at different recording distances z of the simulated RBC (red) and the experiment using a real RBC (green). Inset shows the Cassini model of the RBC used in the simulations; (b) the experiment using a real RBC (green), the reconstructed Re(U) for the simulated unaltered RBC and a deformed RBC represented by the red and blue curves, respectively; (c) reconstructed data from a ~20% smaller RBC compared with the one used in (b). Red and green curves represent the Re(U) of simulation and experiment, respectively; Gray vertical line in (b) and (c) indicates position of the RBC[52].

      第一步, 在合适的记录距离下, 将实际实验RBC与仿真实验RBC在不同重建距离z′ = 0—100 μm复振幅实部信息${\rm{Re}} (U)$进行对应比较, 图4(b)中表示绿线与红线分别代表当记录距离z = 60 μm时, 实际实验RBC与仿真实验RBC的复振幅实部信息${\rm{Re}} (U)$随反向重建距离变化而相应的值, 在相变位置z′ = 60 μm周围, 实际实验RBC与仿真实验RBC的${\rm{Re}} (U)$曲线吻合良好; 第二步, 进一步研究RBC微小形态变化对${\rm{Re}} (U)$的影响, 通过改变仿真实验RBC的模型参数, 使初始RBC发生微小形变, 微形变后RBC在不同记录距离$z$下归一化强度峰值分布与初始RBC基本重叠, 但其${\rm{Re}} (U)$变化较大, 如图4(b)中紫线所示, z′ = 60 μm附近波峰值提高, 而波谷位置由z′ = 50 μm处移至z′ = 35 μm附近, 且波谷值出现了下降, 表明RBC微小形变可根据检测${\rm{Re}} (U)$的变化得出; 第三步, 为证明提出实部信息${\rm{Re}} (U)$重建方法的鲁棒性, 作者在实际实验中选取了比初始RBC体积小约20%的新型RBC作为对象, 在仿真实验中更改RBC模型参数, 使其表示为一个5.04 μm宽的RBC形状模型, 其实际实验RBC与仿真实验RBC在不同重建距离z′ = 0—100 μm复振幅实部信息${\rm{Re}} (U)$的曲线如图4(c)所示, 实验结果表明, 在选取不同体积的RBC的情况下, 实际实验与仿真实验的${\rm{Re}} (U)$曲线仍良好吻合, 因此${\rm{Re}} (U)$提供了足够信息以区分细胞不同的几何形状. 综上所述, 所提出瑞利索末菲反向传播算法可以实现细胞形状的辨别, 应用于生物学中数量庞大且深度不同的细胞在特定时间段内形态变化的病理研究.

    • 通过高空间分辨率地提取血流量信息, 进而实现人体循环系统血管疾病的跟踪监测, 这是当前生物医学领域的重要课题之一. 韩国浦项科技大学Lee等[51]提出了一套同轴数字全息显微系统, 通过数字全息显微成像技术对流动的RBC逐一进行追踪, 引入了清晰度量化算法对重建的RBC图像进行评估, 定位RBC的精确深度位置, 可获取单个RBC的运动状态, 实现高空间分辨率的血流量信息测量. 其光学系统的分析与实现可根据图5所述结构进行描述, 注射泵(SP)注射血溶液经过肝素化处理, 防止凝结, 溶液中RBC数量为3.2 × 103 cell/mm3以防止RBC在全息图中出现重叠. He-Ne激光器发出632.8 nm波长红光, 经过水溶液与氟化乙烯丙烯(FEP)微管后进入血溶液, 部分光经过血溶液中的RBC发生衍射作为物光, 而未通过RBC而直接透过的光作为参考光, 两束光由浸水显微物镜MO(20×/0.50$NA$)放大后在CMOS记录面上发生相干叠加, 产生干涉全息图, 由于血溶液折射率$n_{\rm{medium}} = 1.346$, FEP微管折射率$n_{\rm{FEP}}= 1.338$, 水的折射率$n _{\rm{water}}= 1.333$, 三者折射率几乎相同, 避免了折射像畸变问题. 实验初始阶段首先移动CMOS像感器以找到同时满足成像清晰与足够视场两个条件的合适位置, 将CMOS像感器固定, 在测量体积600 μm × 350 μm × 350 μm范围内可观察到185个RBC. 根据同轴全息重建算法对RBC图像进行重建, 重建距离设定为深度方向上以2 μm为间隔建立共计200个重建切片进行RBC图像的重建得到重建图像.

      图  5  同轴数字全息显微成像系统追踪红细胞空间分布[51]

      Figure 5.  Inline digital holographic microscopic system for tracking spatial distribution of RBCS[51].

      为得到流动RBC实际深度位置$z$, 引入图像清晰度量化算法对重建图像进行判断[57]:

      ${\rm{LAP}} (z) = \sum\nolimits_{xy} {{{\left\{ {{\nabla ^2}I(x,y;z)} \right\}}^2}} $

      式中表示重建图像强度梯度的变化状况, 其中参数指标${\rm{LAP}} (z)$可以在忽略恒定强度梯度部分的情况下, 用于测量图像的高频边缘成分进而反映RBC的聚焦值, 在本实验中对比其余图像清晰度量化算法其精度更高[58,59], 参数指标${\rm{LAP}} (z)$的大小随着重建距离$z'$的变化而变化, 可用于描述同一个RBC在不同距离的重建切片下其聚焦值的大小变化, 由于实验在深度方向上以2 μm为间隔建立共计200个重建切片, 则同一个RBC在200个重建距离下存在200个聚焦值, 以重建距离为横坐标, 用光滑曲线表示200个聚焦值, 其聚焦值最大位置表示焦点位置, 定义此位置为RBC的最佳重建距离, 根据计算显微物镜的轴向放大率可进一步得到RBC的深度位置. 图6(a)表示CMOS像感器所采集的RBC全息图, 图中共存在185个RBC, 且RBC的深度互有异同. 首先对RBC进行分割与重建, 根据清晰度量化算法对重建图像进行最佳重建距离判定, 进一步获得RBC的深度位置. 图6(b)6(d)分别表示重建距离z' = 168 μm、z' = 250 μm、z' = 382 μm时RBC重建图像, 其中箭头表示该重建距离下可清晰观察到的单个RBC样品, 依次确定全息图中185个RBC的深度位置, 聚焦RBC中心位置的光强度具有局部最小值, 可通过光强度大小判断每一个RBC的横向位置$(x, y)$.

      图  6  RBC的全息图与重建图像 (a) CMOS像感器拍摄FEP微管内RBC所得全息图; (b), (c), (d)分别表示不同重建深度下RBC重建图像, 箭头表示聚焦的RBC[51]

      Figure 6.  The hologram and the reconstruction images of RBCs. (a) The RBC hologram obtained from CMOS; (b), (c), and (d) represent RBC reconstruction images at different reconstruction depths, respectively. Each focused RBC is shown by an arrow[51].

      综上所述, 可以获取全息图记录的每一个RBC的三维信息, 其相关x, y, z测量精度如表1所示, 由表可知RBC的横向信息与轴向信息的测量精度均较高, 能准确获取RBC的三维信息. 进一步地, 通过CMOS像感器的连续拍摄, 记录RBC的运动轨迹全息图, 通过同轴全息重建算法对RBC进行重建, 分析FEP微流管内RBC的三维运动, 实现RBC动态追踪, 得到RBC的实时空间分布.

      坐标测量精度/μm均方根误差/μm
      $x$ ± 0.30.63
      $y$ ± 0.30.52
      $z$ ± 1.02.05

      表 1  RBC横向信息与轴向信息的测量精度[51]

      Table 1.  The lateral and axial measurement accuracy of RBC[51].

    • 为实现高分辨率和无损地识别不同形状生物细胞[60,61], 韩国朝鲜大学的Moon等[54]提出一套离轴数字全息显微系统用于测量RBC的三维体积, 该系统引入分水岭分割算法, 它是一种基于拓扑理论的数学形态学分割算法, 将图像看作一幅地形图, 图像中每一点像素灰度值表示该点海拔高度, 将每一个局部极小值及其影响区域作为集水盆, 集水盆的边界形成分水岭, 根据图像自身灰度差异进行区域划分[62]. 通过分水岭分割算法对重建RBC相位图像进行区分, 单独提取每一个RBC, 根据单个RBC自身表面积和相位值使用自动算法分别计算每个RBC的体积, 生成全部RBC的三维体积分布. 为验证该方案的可行性, 通过上述方案分别计算生成口腔状RBC与盘状RBC的三维体积分布, 实验表明, 所提出的离轴数字全息显微成像方法可以对各种形状RBC进行三维体积定量分析, 其光学系统的分析与实现可根据图7所述结构进行描述, 激光器发射682 nm波长红光, 由分束器(BS1)将入射光分为两束, 第一束经过待测样本发生衍射, 经显微物镜MO (40×/0.75NA)放大后作为放大物光先后经过反射镜与分束器后入射至CCD像感器, 第二束光作为参考光, 通过分束器(BS2)后与放大物光在CCD像感器上发生干涉, CCD像感器记录生成的全息图.

      图  7  离轴数字全息显微成像系统测量不同形状RBC三维体积[54]

      Figure 7.  Off-axis digital holographic microscopy system for measuring RBCs’ three-dimensional volume of different shapes[54].

      根据离轴全息重建算法对记录的全息图进行重建生成RBC相位图, 单个RBC的体积计算公式如下:

      $ V \cong \frac{{N{p^2}}}{{{M^2}}}\frac{{\lambda \displaystyle\sum \varphi }}{{2{\text{π}}\left( {{n_{{\rm{rbc}}}} - {n_{\rm{m}}}} \right)}}, $

      其中, $N$表示RBC的像素总数, $p$表示像素尺寸大小, $\varphi $表示RBC内每个像素对应的相位值, $M$为显微物镜的放大倍数. RBC折射率${n_{{\rm{rbc}}}} = 1.396$, HEPA介质折射率${n_{\rm{m}}} = 1.3334$.

      为了研究不同形状RBC的三维体积, 在图7中选取了口腔状(stomatcyte)RBC与盘状(discocyte)RBC进行对比. 口腔状RBC与盘状RBC重建相位图像如图8(a)8(b)所示, 通过分水岭分割算法去除相位图像背景与RBC的互相重叠部分, 对RBC进行逐一区分, 如图8(c), 8(d)所示, 可以看出口腔状RBC的重建相位图像与盘状RBC的重建相位图像中相邻细胞间均明显隔开, 能够实现单个口腔状RBC或盘状RBC的提取. 为了对RBC的三维结构进行更合理地分析, 通过标记分水岭分割算法对图8(e)所表示的单个RBC进行标记, 得到标记后的单个RBC的整体C部分, 如图8(h)所示, 经过标记分水岭分割算法计算得到单个RBC的A部分, 如图8(f)所示. 通过减法运算, 将单个RBC的整体C部分减去A部分可得到RBC的B部分, 如图8(g)所示.

      图  8  RBC的相位重建图像 (a)重建的口腔形状RBC相位图像; (b)重建的盘状RBC相位图像; (c)重建后经分水岭算法分割的口腔形状RBC相位图像; (d)重建后经分水岭算法分割的盘状RBC相位图像; (e)分割的单个RBC相位图像; (f)(g)(h)分别经标记分水岭算法进一步分割得到的单个RBC的A、B、C部分[54]

      Figure 8.  The reconstructed phase image for RBCs (a) The reconstructed phase image for RBCs having a stomatocyte shape; (b) the reconstructed RBCs phase image for RBCs having a discocyte shape; (c)the segmented phase image for RBCs having a stomatocyte shape; (d) the segmented phase image for RBCs having a discocyte shape; (e) the segmented phase image for single RBC(f), (g) and (h) represent the A, B and C parts by the marker-controlled watershed algorithm in RBC, respectively[54].

      为了对两种不同形状RBC的三维体积进行定量研究, 首先计算出最佳重建距离下RBC的相位分布, 再根据(11)式分别对每个口腔状或盘状RBC的A与B部分进行计算, 得到平均体积和均方根误差的数值, 结果如表2所示. 由表可知盘状RBC的A部分体积大于口腔状RBC的A部分体积, 而B部分体积小于口腔状RBC的B部分体积, 根据进一步的统计假设检验证明了RBC的三维体积随着其形状的改变会出现较大区别[63]. 综上所述, 根据RBC的体积自动计算算法结合分水岭分割算法可以研究不同形状RBC的三维体积并且实现实时检测RBC的三维体积变化, 对生物医学领域中细胞成像发展研究具有极大价值.

      A部分B部分
      平均体积/μm3均方根
      误差
      平均体积/μm3均方根
      误差
      口腔状7.67.241.514.7
      盘状14.68.032.77.3

      表 2  两种不同形状RBC的A、B部分的三维体积

      Table 2.  The different shapes of RBC’s three-dimensional volume of A and B parts.

    • 数字全息的定量相衬成像在显微生物成像中具备了多种新应用[64], 其中包括微观对象的评估, 如生物细胞生长观测[65], 由渗透压引起的细胞体积变化检测[66], RBC监测[67]. 这些研究中通常使用定性指标即通过重建图像视觉质量来进行焦平面评估, 因此获取结果具有不确定性[6870]. 基于此, 美国杜克大学 Wax等[53]设计了一套离轴数字全息系统并提出了一种数字重聚焦方法, 通过定量指标评估焦平面并描述微观物体的方法, 在尝试多种最佳焦平面的评估方法并对比每种方法的重建图像效果后, 选择振幅方差最小化作为定量指标, 由于其在评估物体透射的相位图像时, 在深度方向聚焦清晰度与聚焦范围均具有较大优势, 且计算量较小, 可找到理想的焦平面位置. 其光学系统的分析与实现可根据图9所述结构进行描述, 低相干激光器发射589 nm黄光, 分束器将黄光分为两束, 由后向反射镜实现两束光的光程匹配, 第一束光通过后向反射镜调整光程入射至RBC发生衍射, RBC置于样品台上, 样品台可移动以控制样品与焦平面距离, 衍射光经过显微物镜MO( 40×/0.75NA)放大后作为放大物光入射至CMOS像感器; 另一束光作为参考光通过后向反射镜调整光程经显微物镜MO(40×/0.75NA)放大后作为放大参考光与放大物光在CMOS像感器上发生相干叠加, 形成全息图, 由CMOS记录生成的全息图.

      图  9  离轴数字全息显微系统研究散焦现象对RBC三维体积测量影响[53]

      Figure 9.  Off-axis digital holographic microscopy system for investigating the effect of defocus on RBC three-dimensional volume measurement[53].

      根据离轴全息重建算法对记录全息图进行重建, 获得RBC重建图像. 为检验该定量指标评估焦平面后重建物体的指标, 定义了微观物体的光学体积, 公式如下

      $\begin{split}OV\;& = \iint_{x,y} \Delta OPL(x,y){\rm{d}}x{\rm{d}}y \\ &= \iint_{x,y} \Delta n(x,y) \cdot h(x,y){\rm{d}}x{\rm{d}}y\end{split}$

      其中, 光程长度的变化$\Delta OPL\left( {x, y} \right) = \Delta \varphi \left( {x, y} \right)/k = \Delta n\left( {x, y} \right) \cdot h\left( {x, y} \right)$用于描述微观物体, $h(x, y)$表示物体的厚度(高度)图, $\Delta n(x, y)$表示物体折射率图, $k = 2{\text{π}}/\lambda $表示波数, 合并$\Delta OPL(x, y)$参数产生物体光学体积$OV$的度量, 若物体的折射率均匀, 即$\Delta n(x, y) = \Delta \bar n$, 说明物体光学体积与物体真实体积线性相关, $OV = \Delta \bar n \cdot V$.

      为进一步比较定性指标与定量指标分析中分别评估得到的焦平面效果, 根据焦平面各自对应的重建物体图像的光学体积对比原物体光学体积真值, 以误差值衡量所得焦平面的效果. 引入直径为8 μm的聚苯乙烯微球(直径变异系数1%)作为样品, 首先根据人工聚焦方法即重建图像视觉质量来进行焦平面评估, 移动样品台通过CMOS像感器记录微球全息图, 在不同重建距离下对微球图像进行重建, 通过人工观察重建微球的图像质量确定最佳焦平面位置, 再计算重建微球的相位分布, 设定此时重建微球球心所在平面位置为z = 0. 根据提出的数字重聚焦方法计算z = ± 50 μm范围内微球球心处的振幅方差, 确定振幅方差最小值位置为z = –2.46 μm, 将其作为最佳焦平面位置, 并计算微球的相位分布. 通过(12)式分别计算人工聚焦方法与数字重聚焦方法所得微球光学体积, 对比微球光学体积真值(1%的直径变异系数对应于3%光学体积变异系数), 结果如表3所示. 表中 ±7.98 fL表示聚苯乙烯微球3%的光学体积变异系数, 人工聚焦方法测得微球OV = 281.7 fL, 以5.90%比率高于微球真值OV = 266.0 fL, 而数字重聚焦方法测得微球光学体积为OV = 266.6 fL, 仅以0.22%比例高于微球光学体积真值, 说明根据数字重聚焦方法所得焦平面更具可靠性, 将其认定为理论最佳焦平面, 此时人工聚焦方法所得$z = 0$处距离理论最佳焦平面距离为2.46 μm, 侧面反映了即使在2.46 μm的轻微散焦距离下, 所测得光学体积结果仍会产生明显误差. 为研究散焦距离变化对微球光学体积测量结果的影响, 通过移动控制台改变微球散焦距离, 实时测量对应微球的光学体积, 对比最佳焦平面处测得微球的光学体积, 可以得到$\dfrac{{{\rm{d}}V}}{{{\rm{d}}z}} \approx 5.75$ fL/μm, 即散焦距离每改变1 μm, 其光学体积测量结果相差2.2%, 当散焦距离达到 ± 10 μm时, 会出现相位包裹伪影进一步降低测量精度, 而根据数字重聚焦方法可以解决散焦的问题, 进一步证明了该方法在待测样本光学体积测量时的高精度优势.

      测量方法光学体积/fL误差比/%
      人工聚焦方法281.75.90
      数字重聚焦方法266.60.22
      真实值266.0 ± 7.980.00

      表 3  人工聚焦方法与数字重聚焦方法测得微球光学体积对比[53]

      Table 3.  Comparison of OV measured by manually-focused and digitally-refocused methods[53].

      为显示生物样品中光学体积测量中散焦的影响, 选取单一RBC作为测试样品, 在约200 μm范围人工选取7个不同焦平面, 分别对RBC进行重建得到相对应的振幅图与相位图, 如图10(a)中首行与中间行(A—G)所示.

      图  10  对单一RBC图像进行数字重聚焦与相应的RBC光学体积测量 (a)通过人工聚焦方法与数字重聚焦方法对单一RBC重建所得振幅图与相位图; (b)A—G的振幅方差分布; (c)RBC在人工聚焦方法所得光学体积(黑线)与数字重聚焦方法所得光学体积(蓝线), 光学体积OV表示为平均值 ± 标准差[53]

      Figure 10.  Digital refocusing of a single red blood cell image and corresponding optical volume measurements. (a) The amplitude and phase images by the manually-focused method and digitally-refocused method from a single RBC; (b) amplitude variance metric of holograms A-G; (c) computed OV of RBC from manually-focused phase images(black) and digitally-refocused phase images(blue). OV reported as mean ± standard deviation[53].

      其中, D图表示重建微球视觉质量最佳时所得微球的振幅图与相位图, D图的焦平面位置即为定性指标中最佳焦平面位置, 分别计算A—G对应的振幅方差最小值可得到7个焦平面相对定量指标中理论最佳焦平面位置($z = 0$)各自的散焦距离, 如图10(b)所示, 人工聚焦方法所得重建RBC图像A—G的光学体积计算结果如图10(c)中黑色线所示, 随着散焦距离的变化, RBC在OV = 5.15—7.03 fL ($ {\rm{\sigma}} = 13.55\%$)范围内波动, 散焦距离较大影响了RBC体积测量数值的可靠性; 而通过数字重聚焦方法获取振幅方差最小值以得到理论最佳焦平面位置, 重建所得RBC相位图如图10(a)末行所示, 通过数字重聚焦方法所得A—G的相位图相似, RBC的光学体积仅在OV = 5.54—5.70 fL ($ {\rm{\sigma}} = 0.99\% $)范围内波动, 进一步证明数字重聚焦方法对光学体积测量的实用性与可靠性. 综上所述, 根据数字全息显微成像方法表征微观物体三维体积时, 重建图像的散焦现象是潜在的重要误差来源, 提出的数字重聚焦方法能够以定量指标评估焦平面, 对重建RBC进行光学体积测量, 实时获得准确RBC的三维体积, 可广泛应用于人体的心血管疾病与糖尿病的病理研究之中.

    • 帕金森氏疾病(PD)作为一种中脑区神经元丢失的神经退行性疾病, 随着PD发病机理研究的不断深入, 在各种病理假设之中, 氧化应激是引起细胞功能障碍与细胞凋亡的主要因素[43,44]. 在随后的研究中, 证实了PD中氧化应激的存在, 根据以往研究发现氧化应激还会诱导细胞变形[71], 因此, 研究氧化应激下细胞的机械性能对PD的研究具有重要意义. 基于此, 祝连庆等[72]设计了一套光镊辅助离轴数字全息显微系统, 用于研究氧化应激下RBC的高度与三维体积变化, 有助于帕金森氏疾病的深入研究. 其光学系统的分析与实现可通过图11所述结构进行描述, 引入含有链霉亲和素涂层的聚苯乙烯微球(直径4 μm)以2∶1的浓度比与人类RBC样本进行混合至磷酸盐缓冲溶液之中(PBS, pH = 7.4), 筛选出RBC水平两侧分别连结一个微球的样品作为实验待测对象, 如图11中待测样本放大图所示. 具有统一参数的水平两侧微球通过光镊控制一侧微球静止, 另一侧恒定慢速移动以实现RBC的拉伸, 考虑到RBC的尺寸并非绝对均匀, 因此, 从满足实验条件细胞中随机选取10个RBC取计算结果平均值, 通过如图11所示光学系统进行全息图的记录, 光学系统由离轴数字全息显微系统与光镊组成, 632.8 nm激光器发射红光通过分束器(BS1)分为物光与参考光, 物光通过油浸物镜MO(60× /1.40NA)放大后与参考光于CMOS像感器记录面上发生相干叠加, 形成干涉全息图, 参考光的反射角度通过反射镜的倾斜度调整, 用于控制干涉条纹的方向与宽度, 1064 nm激光器(5 W, Nd: YAG)作为光镊的照明光源, 发射红外光通过二向色镜(DM)进入倒置显微镜的后孔形成光镊, 其中声光偏转器(AOD)通过声光的相互作用, 改变激光角度从而控制RBC水平两侧微球的位置以实现拉伸.

      图  11  光镊辅助离轴数字全息显微系统测量RBC体积[72]

      Figure 11.  Off-axis digital holographic microscopy system with optical tweezer for measuring RBCs’ three-dimensional volume[72].

      利用离轴全息重建算法对CMOS像感器所记录的全息图进行重建, 并通过重建相位图像的参数计算不同氧化应激浓度下RBC的高度与三维体积, 公式如下:

      $h = \frac{{\lambda \varphi }}{{2{\text{π}}\left( {{n_{{\rm{cell}}}} - {n_{{\rm{medium}}}}} \right)}}, $

      $V = \sum\limits_{i = 1}^M {\sum\limits_{j = 1}^N h } \left( {{x_i},{y_j}} \right), $

      其中, $\lambda $代表全息显微系统中的激光波长(632.8 nm), $\varphi $代表相位值, RBC的折射率$n_{\rm{cell}}=1.37$, 溶液的折射率${n_{{\rm{medium}}}} = 1.34$, $V$表示RBC的体积, $h\left( {{x_i}, {y_j}} \right)$表示RBC上任意点的$\left( {{x_i}, {y_j}} \right)$的高度, $M$$N$表示RBC的横向坐标范围.

      将过氧化氢(H2O2)溶液(0.1 mol/L)在PBS中稀释至50、100、150和200 μmol/L对照组, 置于5个培养皿中, 与符合实验条件附着有微球的RBC样品进行混合. 为具体分析陷阱的拉伸力与细胞变形之间的关系, 需对光阱刚度值进行校准, 校准分为已知力校准方法和热涨落方法[73,74]. 热涨落是系统在平衡状态下相对于其平均状态的随机偏差, 是系统温度的基本体现, 根据Boltzmann统计分析热噪声下随机运动粒子的布朗轨迹, 可以获得粒子位移概率分布图, 接着通过计算粒子的电势获得光阱刚度$k$, 该值被校准为14.23 ± 0.46 pN/μm. 校准后保持激光功率不变以确保光阱强度均匀. 光阱拉伸力$F$由光阱刚度$k$和微球与陷阱之间中心差$\Delta x$共同确定, 控制陷阱拉伸力在0~3 pN之间变化. 将RBC所能达到的最大高度值定义为$H$, 同时为了清楚观察$H$的变化, 在图12中分别显示了$H$大于2 μm的高度部分, 随着氧化应激浓度的增大, RBC的细胞膜出现粗糙度且其边缘呈现锯齿状. 分别计算五项氧化应激浓度下测量所得10个RBC的最大高度与三维体积的平均值, 能够更直观比较RBC在不同浓度氧化应激下高度与三维体积的变化, 绘制的相关折线图如图13所示, H/H0表示拉伸过程中细胞的最大高度$H$与未拉伸前RBC的原始高度${H_0}$的比值, 图13(a)表示0—200 μmol/L氧化应激浓度下RBC高度比$H/{H_0}$与陷阱拉伸力的关系, 反映了RBC的微形变机械性能, 由图可知, 随着氧化应激浓度的增加, RBC的变形能力下降, 因此氧化应激下的RBC其机械性能会受到严重影响. 图13(b)表示0—200 μmol/L氧化应激浓度下RBC三维体积的变化, 由图可知, 对比不加入H2O2的常规PBS中的RBC, 氧化应激使RBC的体积减小, 在50—200 μmol/L的H2O2浓度下, RBC的体积变为初始体积的90%—87.3%, 该结果与过往研究相一致[75]. 综上所述, 在氧化应激下RBC的机械性能会发生变化, 微形变能力下降, 而氧化应激是导致PD的主要因素, 该光镊辅助离轴数字全息显微成像的技术将助于PD的基础研究与临床应用.

      图  12  0−3 pN陷阱拉伸力变化下不同浓度(0−200 μmol/L)氧化应激下重建RBC的高度变化, 颜色深浅代表RBC高度的大小[72]

      Figure 12.  Height change of reconstructed RBCs under different concentrations of oxidative stress (0−200 μmol/L). Four images in each group are corresponding to trap force varying from 0−3 pN. Color bar represents different thickness[72].

      图  13  RBC在不同浓度氧化应激下的性能 (a)不同浓度氧化应激下RBC最大高度H与陷阱拉伸力关系; (b)不同浓度氧化应激下RBC体积[72]

      Figure 13.  Performance of RBC under different oxidative stress. (a) The relationship between the maximum height H of RBC and the trap tensile force under different oxidative stress; (b) the volume of RBC under different oxidative stress[72].

    • 本文简要介绍了数字全息显微成像由全息图记录到全息图重建的原理, 分析讨论了不同数字全息显微系统在RBC成像的应用, 在数字全息显微系统方面, 主要分为同轴、离轴、光镊辅助离轴的数字显微全息成像系统. 同轴数字全息显微成像系统光路简单, 可采用激光作为光源, 但其对光源的相干性以及系统的稳定性要求相对离轴数字全息显微成像系统较低, 可通过弱相干光源(如LED)代替激光作为光源, 采用弱相干光源避免了激光产生的散斑噪声, 可以得到较好的生物细胞图像, 降低成本, 装置结构轻便, 便于集成. 离轴数字全息显微成像系统引入参考光光路, 通过改变参考光入射角度, 实现重建像与孪生像的分离. 光镊辅助离轴的数字全息系统中, 光镊作为辅助检测的手段, 可给予活体细胞不同拉伸力使之产生形变, 并根据离轴数字全息系统对细胞进行实时成像, 进而对活体细胞的性能进行更深一步的研究. 在不同的实验要求下, 可根据具体需求选择合适的全息显微成像系统以实现实验目的.

      在细胞三维信息提取的过程中, 数字重聚焦方法选择图像在不同距离下振幅方差最小值作为定量指标代替图像视觉质量判断的定性指标, 以评估RBC的焦平面位置, 减少由于散焦造成的三维体积测量误差. 或根据图像清晰度量化算法也可以实现RBC焦平面评估, 确定重建距离, 通过重建图像信息进一步获取RBC表面积、相位值, 并根据入射光波长、折射率等参数信息计算RBC高度与三维体积, 也可通过反向重建放大物光复振幅的实部信息来判断RBC几何形状的变化; 对于多个RBC的情况, 引入分水岭分割算法对重建RBC相位图像进行分割, 单独提取每一个RBC后, 再通过上述算法分别计算每个RBC的高度与三维体积并取对应平均值获得RBC整体三维信息. 更进一步地, 通过光镊控制RBC拉伸产生微形变, 根据细胞的三维信息探究RBC在不同浓度氧化应激下的形变能力.

      通过数字全息显微成像技术有效实现细胞三维体积的测量, 并且可以检测出细胞几何形状的微小形变, 对比于当前显微成像领域广泛运用的激光扫描共聚焦显微技术, 近场光学显微技术, 光学相干层析成像技术, 其非接触、无损性、定量化、实时性的优势更有利于细胞生物学参数与细胞性能的研究. 未来数字全息显微成像技术有两个大方向可以进行持续探索, 一方面可进一步发挥先进光学器件优势, 在硬件上实现分辨率的提升, 并设计更优化的重建算法, 通过更少的数据量与更精确的焦平面评估来实现更高分辨率, 更高精度的细胞成像; 另一方面由于系统的分辨率受限于显微物镜数值孔径, 可采用合理的数字全息无透镜成像方法, 待测细胞与像感器之间不存在光学透镜, 因此可避免透镜表面灰尘造成的杂散光影响, 且简化了成像的光路, 具有更大的成像视野, 同时满足细胞显微成像的大视场与高分辨率需求, 便于细胞成像研究方面的集成化运用. 同时, 基于深度学习的图像重建算法的发展, 可根据增加像素数的方式在传统方法的基础上进一步提升细胞成像的分辨率, 提高成像质量.

参考文献 (75)

目录

    /

    返回文章
    返回