CN108169173B - 一种大视场高分辨三维衍射层析显微成像方法 - Google Patents

一种大视场高分辨三维衍射层析显微成像方法 Download PDF

Info

Publication number
CN108169173B
CN108169173B CN201711484784.XA CN201711484784A CN108169173B CN 108169173 B CN108169173 B CN 108169173B CN 201711484784 A CN201711484784 A CN 201711484784A CN 108169173 B CN108169173 B CN 108169173B
Authority
CN
China
Prior art keywords
dimensional
field
resolution
spectrum
intensity
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Active
Application number
CN201711484784.XA
Other languages
English (en)
Other versions
CN108169173A (zh
Inventor
左超
李加基
陈钱
孙佳嵩
冯世杰
张玉珍
顾国华
张佳琳
范瑶
丁君义
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Nanjing University of Science and Technology
Original Assignee
Nanjing University of Science and Technology
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by Nanjing University of Science and Technology filed Critical Nanjing University of Science and Technology
Priority to CN201711484784.XA priority Critical patent/CN108169173B/zh
Publication of CN108169173A publication Critical patent/CN108169173A/zh
Application granted granted Critical
Publication of CN108169173B publication Critical patent/CN108169173B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
    • G01N21/00—Investigating or analysing materials by the use of optical means, i.e. using sub-millimetre waves, infrared, visible or ultraviolet light
    • G01N21/17—Systems in which incident light is modified in accordance with the properties of the material investigated
    • G01N21/41—Refractivity; Phase-affecting properties, e.g. optical path length

Landscapes

  • Physics & Mathematics (AREA)
  • Health & Medical Sciences (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • Chemical & Material Sciences (AREA)
  • Analytical Chemistry (AREA)
  • Biochemistry (AREA)
  • General Health & Medical Sciences (AREA)
  • General Physics & Mathematics (AREA)
  • Immunology (AREA)
  • Pathology (AREA)
  • Investigating Or Analysing Materials By Optical Means (AREA)

Abstract

本发明公开了一种基于LED阵列照明的大视场高分辨三维傅里叶叠层衍射层析显微成像方法，首先进行原始强度图像采集，根据LED阵列中每颗LED在空间中的坐标位置计算出在照明系统中每个LED灯对应的入射光的空间频率，然后初始化被测物体的大视场高分辨三维频谱，且该初始化矩阵满足每个方向上的最小采样数与最终的成像分辨率要求，并将每个照明角度下所拍摄的强度图像迭代至初始化的三维频谱中，并进行多轮迭代，迭代得出被测物体的三维频谱，并将该三维频谱变换至空域，最终得到被测的三维物体大视场高分辨的折射率信息分布。本发明无须采用高放大倍率的物镜，在保证较大的成像视场前提下也可达到较高的重构分辨率。

Description

一种大视场高分辨三维衍射层析显微成像方法
技术领域
本发明属于光学显微测量、三维折射率成像技术，特别是一种基于LED阵列的大视场高分辨三维傅里叶叠层衍射层析显微成像方法。
背景技术
大部分活细胞和未染色的生物标本都是无色透明的，这是因为细胞内各部份细微结构的折射率和厚度的不同，当光波通过时，波长和振幅并不发生变化，仅相位发生变化，这种相位差人眼无法观察。这就需要通过一些化学或者生物手段来对细胞进行染色，从而使其在显微镜下可见。由于生物细胞内部不同物质对有色染料吸收的吸收率不同，故最终在显微镜下就可以看见生物细胞内部的结构和表现出不同的生物特性(龚志锦,詹熔洲.病理组织制片和染色技术[M].上海科学技术出版社,1994)。但由于这种染色手段会对细胞和组织的正常生理活动造成影响，更有可能会杀死细胞。所以通过光学调控对生物样品的无标记观测手段应运而生，其中，泽尼克相衬成像和微分干涉相衬成像作为传统的无标记成像手段，通过将生物组织的相位延迟转化为强度信息，从而使未染色的细胞和生物组织在镜下可见。此外，待测物体的定量相位信息或者由折射率分布不均引起的光程差可利用干涉的方法来获取，利用干涉成像手段来获取定量相位信息的的典型方法为数字全息显微镜(Digital Holography Microscopy,DHM)，但这类成像方法需要激光光源和复杂的干涉光路，且对环境的振动较为敏感，不利于成像质量的提升(M.K.Kim,Digital HolographicMicroscopy,Springer,New York 2011)。
计算层析成像(Computed Tomography，CT)作为生物医学领域常用的成像手段之一，利用X射线和不同入射角度下的生物组织对X射线的不同吸收率来重构出物体的三维结构信息(Beerlink A,Mell M,Tolkiehn M,et al.Hard x-ray phase contrast imagingof black lipid membranes[J].Applied Physics Letters,2009,95(20):203703-203703-3)。当入射光波长变短至可见光波段的时候，被探测物体与入射光波会发生衍射现象，所以必须利用另一种成像模型来重构生物组织的三维分布。传统的光学衍射层析成像是通过改变干涉系统中物光的入射光波与参考光之间的夹角，从而可以得到每个入射光角度下物体相位分布，最终对生物物体三维折射率分布进行重构。Choi.W等人在2007年提出了衍射层析相位显微镜(Choi W,Fang-Yen C,Badizadegan K,et al.Tomographic phasemicroscopy[J].Nature Methods,2007,4(9):717-719)，该成像方法是将振镜引入马赫曾德离轴全息的干涉系统中，入射光从各个方向扫描样片，重构出了未染色生物细胞和线虫(C.elegans)的三维折射率分布。2009年Sung.Y等人又将基于全息干涉和振镜旋转的方法用于活细胞的动态观察中(Sung,Y.,Choi,W.,Fang-Yen,C.,Badizadegan,K.,Dasari,R.R.,&Feld,M.S.Optical diffraction tomography for high resolution live cellimaging[J].Optics express,2009,17(1),266-277.)。传统三维折射率层析成像方法利用干涉的方法需先获取在每个入射方向上的散射复振幅分布，该方法不仅需要复杂的干涉系统，而且振镜的机械转动也会影响最终的三维重构效果。南京理工大学智能计算成像实验室在2017年提出了一种利用LED照明和光强传输方程的非干涉的三维衍射层析成像方法(Li J,Chen Q,Zhang J,et al.Optical diffraction tomography microscopy withtransport of intensity equation using a light-emitting diode array[J].Optics&Lasers in Engineering,2017,95:26-34.)，该方法通过拍摄在两个或多个不同焦面上的不同角度强度图像堆栈，再利用基于光强传输方程的相位解法直接求解出在不同入射光方向下的相位分布，最终重构出待测样品的三维折射分布。
但是，上述所有三维衍射层析成像方法拍摄的强度图像均为明场图像，即照明数值孔径小与或等于物镜数值孔径NAill≤NAobj，而一些物体高频信息的散射光且没有被充分利用。所以如何将傅里叶衍射层析定理和傅里叶叠层成像技术有机结合，利用低倍物镜的大视场成像和由增大照明角度提升分辨率的方法来同时兼顾大视场和高分辨的三维折射率衍射层析成像是一个亟待解决的一个技术难题。
发明内容
本发明的目的在于提供一种基于LED阵列照明的大视场高分辨三维傅里叶叠层衍射层析显微成像方法，不仅无须干涉，简化了成像系统光路，而且可高度兼容传统明场显微镜，通过在低倍物镜下增大照明角度就解决了传统三维折射率衍射层析成像高分辨率和大视场不能同时兼顾的技术难题。
实现本发明目的的技术方案为：一种基于LED阵列照明的大视场高分辨三维傅里叶叠层衍射层析显微成像方法，步骤如下：
步骤一，原始强度图像采集，将LED光源作为传统显微镜的照明光源，LED面板上的每颗LED依次被点亮，入射光经过聚光镜后从各个不同的入射角下照射样品，并且在被测的厚样品为聚焦情况下，由相机在同步触发信号下采集一系列不同照明角度下的强度图像
Figure BDA0001534566910000021
其中某一个位置下采集的光强为
Figure BDA0001534566910000022
步骤二，根据LED阵列中每颗LED在空间中的坐标位置计算出在照明系统中每个LED灯对应的入射光的空间频率kin＝(kx,ky)；
步骤三，初始化被测物体的大视场高分辨三维频谱
Figure BDA0001534566910000031
三维频谱矩阵像素个数分别为Nx,Ny,Nz，且该初始化矩阵满足每个方向上的最小采样数与最终的成像分辨率要求；
步骤四，三维频谱迭代重构，将每个照明角度下所拍摄的强度图像迭代至初始化的三维频谱中，并进行多轮迭代；
步骤五，迭代得出被测物体的三维频谱，并将该三维频谱变换至空域，最终得到被测的三维物体大视场高分辨的折射率信息分布n(r)。
本发明与现有技术相比，其显著优点：(1)采用LED阵列作为该显微成像方法的照明光源，且与聚光镜有机组合，这样既保证了照明方向的编程可控，同时也保证了照明数值孔径最高可达到传统阿贝聚光镜的最大值1.20。(2)图像采集过程中不需要机械移动或振镜转动，仅由LED阵列的编程来实现照明光角度的改变。(3)所采集的强度图像中不仅有明场强度图像，还有包含高频衍射信息的暗场强度图像，将暗场图像迭代至三维频谱中，提高了三维成像的轴向分辨率与横向分辨率。(4)无须采用高放大倍率的物镜，在保证较大的成像视场前提下也可达到较高的重构分辨率。
下面结合附图对本发明作进一步详细描述。
附图说明
图1是单颗LED产生入射光线与实际装置示意图。
图2是本发明基于LED阵列照明的大视场高分辨三维傅里叶叠层衍射层析显微成像方法的流程图。
图3是三维傅里叶衍射层析定理示意图。
图4是明场和暗场情况下不同方向的入射光在三维频谱上对应的位置。
图5是最终重构之后的三维频谱的形状示意图。
图6是利用该方法重构出的实球藻在不同位置上的二维折射率分布和三维折射率分布渲染图。
具体实施方式
如图1所示，传统的明场显微镜的照明光源被替换为高密度LED阵列。LED阵列放置在聚光镜的前焦面位置，后端成像系统与传统明场显微镜一致。其中f为聚光镜的焦距，一般在10-20mm之间，并且LED阵列的中心处于成像系统的光轴上。LED阵列中包括若干个(至少261个)LED单元，它们等间隔排布形成一个二维阵列。其中每个LED单元均为红绿蓝三色LED单元，其典型波长为红光633nm、绿光525nm和蓝光465nm。每个LED单元之间中心间距d典型值1-4mm。LED阵列并不需要进行单独加工，一般在市场上可直接购置。表1给出了一个市面上可购置的LED阵列的产品参数。在此LED阵列中，LED单元共有32行、32列，一共1024个，每个LED单元的亮度在2000cd/m2以上。
表1 LED阵列的物理参数
Figure BDA0001534566910000041
LED阵列中每个LED单元均可通过单独点亮，点亮LED单元的具体方法为现有常规技术，实现电路可以采用(但不限于)单片机、ARM、或者可编程逻辑器件等现有技术即可实现，具体实现方法可参考相关文献(如郭宝增,邓淳苗：基于FPGA的LED显示屏控制系统设计[J].液晶与显示,2010,25(3):424-428)。
结合图2，本发明利用上述装置实现的大视场高分辨三维傅里叶叠层衍射层析显微成像方法包括以下步骤：
步骤一，原始强度图像采集。在传统的明场显微镜中，照明光源为高密度LED阵列，每颗单色LED发出的光线均为中心波长为λ0的准单色光，且在不同介质中的归一化波长为λ。该LED阵列放置于聚光镜的前焦面位置上，且该阵列中心的LED与成像系统的光轴重合。LED面板上的每颗LED依次被点亮，入射光经过聚光镜之后从各个不同的入射角下照射样品，并对被测的厚样品进行聚焦。通过产生同步触发信号控制相机采集光强图像堆栈
Figure BDA0001534566910000051
其中某一个位置下采集的光强为
Figure BDA0001534566910000052
区别于传统的三维衍射层析成像方法中只有明场图像，在拍摄的强度图像序列不仅有包含透射光线的明场图像也有只存在衍射光的暗场图像。
步骤二，计算入射光线的空间频率。根据LED阵列中每颗LED在空间中的坐标位置计算出在照明系统中每个LED灯对应的入射光的空间频率kin＝(kx,ky),如图1(b)所示。根据LED面板上被点亮的LED空间位置计算出所对应的入射光空间频率kin,当kin小于所用显微物镜可探测空间频率NAobj/λ时，此时所拍摄的图像为明场图像，而当其中kin大于NAobj/λ时，
Figure BDA0001534566910000059
为暗场光线，其中NAobj为所用物镜的数值孔径。
步骤三，初始化被测物体的三维频谱
Figure BDA0001534566910000053
被测物体的初始化三维频谱
Figure BDA0001534566910000054
的像素尺寸必须满足最终的成像分辨率，且三维矩阵像素个数Nx,Ny,Nz满足每个方向上的最小采样数，其中k为三维频域变量。当最大照明数值孔径为NAill，且NAill≥NAobj时，该成像方法最终可达的横向分辨率的空间频率为：
Figure BDA0001534566910000055
轴向最高可达空间频率为：
Figure BDA0001534566910000056
故最终初始化的三维频谱的像素分辨率必须大于系统的横向分辨率：
Figure BDA0001534566910000057
代入Δux,y后可得关系式：
Figure BDA0001534566910000058
其中ΔXpixel为初始化之后的三维物函数的像素尺寸。
步骤四，三维频谱的强度迭代重构。当入射光线的空间频率为kin＝(kx,ky)时，利用傅里叶衍射层析定理公式：
Figure BDA0001534566910000061
即可得到在某确定入射光空间频率下，对应的埃瓦尔德球壳上的物体三维频谱分量，如图3(a)所示。其中
Figure BDA0001534566910000062
为聚焦位置处的被测物体的二维傅里叶频谱，
Figure BDA0001534566910000063
即为埃瓦尔德球壳的定义函数。经过物镜光瞳限制之后的埃瓦尔德球壳上的频谱信息可以被物镜接收到，而后将该频谱信息在反投影至三维对应位置的埃瓦尔德球壳上，如图3(b)和(c)所示。将三维球壳上的频谱成分映射至二维频谱空间中，即可得到在对应入射光角度下的聚焦位置二维复振幅频谱：
Figure BDA0001534566910000064
其中
Figure BDA0001534566910000065
为二维散射场的傅里叶变换，k2D为二维频域变量，
Figure BDA0001534566910000066
为三维频谱中与入射光线对应的埃瓦尔德球壳上的频谱分量，kz定义为投影系数
Figure BDA0001534566910000067
对于明场图像来说，所拍摄的总场包含着入射分量与散射分量，而暗场图像却只有被物体所散射之后的分量，所以拍摄的强度图与散射场和入射光线的关系为：
Figure BDA0001534566910000068
其中I(xT)为实际所拍摄的强度图像，U(xT)为包含散射场Us(xT)和入射场Uin(xT)的总和场，xT为二维空间变量。
在三维衍射层析成像中，Born近似与Rytov近似是较为常用的两种使散射模型简化的近似方法。其中Born近似可描述为在透射分量与散射分量的总和中，散射成分所占比重极小，即弱散射近似。而Rytov近似则将总场与一阶散射场之间用对数联系起来，这两种近似方法可简化为：
Figure BDA0001534566910000069
其中Ub(xT)表示一阶散射场，对上式进行进一步的化简，对各个变量分别除以入射分量Uin(xT)，化简后的公式可得：
Figure BDA0001534566910000071
其中UBN(xT)，UNS(xT)和UN(xT)分别表示经由入射光场归一化之后的一阶散射场，散射场和测量的总场，即可表示为
Figure BDA0001534566910000072
和
Figure BDA0001534566910000073
故可以分别得到在Born近似和Rytov近似下所测量的明场和暗场的情况下，总场与一阶散射场之间关系的表达式：
对于Born近似
Figure BDA0001534566910000074
对于Rytov近似
Figure BDA0001534566910000075
利用上述关系可以分别得到在Born近似和Rytov近似下将所拍摄的强度图像
Figure BDA00015345669100000711
迭代至二维一阶散射场中的公式为：
对于Born近似
Figure BDA0001534566910000076
对于Rytov近似
Figure BDA0001534566910000077
其中
Figure BDA0001534566910000078
为原始估计的一阶散射场，
Figure BDA0001534566910000079
为迭代更新光强信息之后的散射场，
Figure BDA00015345669100000710
为实际拍摄的强度图像，系数α作为调节项可加快收敛速度减小噪声影响。
再对
Figure BDA0001534566910000081
进行二维傅里叶变换得到
Figure BDA0001534566910000082
将由强度收敛之后的一阶散射势的傅里叶变换反投影到三维频谱中对应入射光线下的埃瓦尔德球壳上，即可实现一个照明角度下的三维频谱迭代重构，即表示为：
Figure BDA0001534566910000083
经过上述公式即可完成一次对三维傅里叶频谱的强度迭代收敛过程。其中图4(a)表示在相干照明情况下在三维频域空间中埃瓦尔德球所覆盖的范围，而图4(b)是明场和暗场情况下不同方向的入射光线在三维频谱上对应的埃瓦尔德球的位置。
步骤五，进行被测物体的三维频谱完整迭代，得到被测物体的三维折射率分布，图5为最终进行强度迭代重构之后的三维频谱的形状示意图，其三维形状为图4(b)的三维渲染图形。根据步骤四进行单次三维频谱光强迭代，将每个照明角度下的强度图
Figure BDA0001534566910000084
依次进行迭代，通过数次的光强迭代之后，所拍摄的强度信息逐渐收敛至三维散射势频谱中，最终即可得到所测量样片的三维频谱信息。最后进行三维逆傅里叶变换得到三维物体的物函数V(r)分布：
Figure BDA0001534566910000085
其中IFT3表示三维傅里叶变换，再利用物函数散射势计算公式即可得到待测样品的三维折射率分布：
Figure BDA0001534566910000086
其中k0为真空中对应波长的波数，nm(r)为物体所处介质的折射率。图6(a)-(d)为利用该方法所重构出的实球藻在不同轴向平面位置上的二维折射率分布，而图6(e)则是待测的实球藻三维折射率分布渲染图。

Claims (3)

1.一种基于LED阵列照明的大视场高分辨三维傅里叶叠层衍射层析显微成像方法，其特征在于步骤如下：
步骤一，原始强度图像采集，将LED光源作为传统显微镜的照明光源，LED面板上的每颗LED依次被点亮，入射光经过聚光镜后从各个不同的入射角下照射样品，并且在被测的厚样品为聚焦情况下，由相机在同步触发信号下采集一系列不同照明角度下的强度图像
Figure FDA0002530654570000011
其中某一个位置下采集的光强为
Figure FDA0002530654570000012
步骤二，根据LED阵列中每颗LED在空间中的坐标位置计算出在照明系统中每个LED灯对应的入射光的空间频率kin＝(kx,ky)，其中kx和ky分别表示空间频率在x方向和y方向的分量；
步骤三，初始化被测物体的大视场高分辨三维频谱
Figure FDA0002530654570000013
三维频谱矩阵像素个数分别为Nx,Ny,Nz，且该三维频谱矩阵满足每个方向上的最小采样数与最终的成像分辨率要求，其中k是指三维频域变量；
步骤四，三维频谱迭代重构，将每个照明角度下所拍摄的强度图像迭代至初始化的三维频谱中，并进行多轮迭代；
步骤五，迭代得出被测物体的三维频谱，并将该三维频谱变换至空间域，最终得到被测的三维物体大视场高分辨的折射率信息分布n(r)，其中r是指三维空间变量；
在步骤四中：当入射光线的空间频率为kin＝(kx,ky)时，利用傅里叶衍射层析定理公式：
Figure FDA0002530654570000014
即可得到在某确定入射光空间频率下，对应的埃瓦尔德球壳上的物体三维频谱分量；其中
Figure FDA0002530654570000015
为聚焦位置处的被测物体的二维傅里叶频谱，
Figure FDA0002530654570000016
即为埃瓦尔德球壳的定义函数，
Figure FDA0002530654570000017
是指三维散射势频谱分布，λ为不同介质中的归一化波长；将三维球壳上的频谱成分映射至二维频谱空间中，即可得到在对应入射光角度下的聚焦位置二维复振幅频谱：
Figure FDA0002530654570000021
其中
Figure FDA0002530654570000022
为二维散射场的傅里叶变换，j是指虚数单位值，k2D为二维频域变量，
Figure FDA0002530654570000023
为三维频谱中与入射光线对应的埃瓦尔德球壳上的频谱分量，kz定义为投影系数
Figure FDA0002530654570000024
对于明场图像来说，所拍摄的总场包含着入射分量与散射分量，而暗场图像却只有被物体所散射之后的分量，所以拍摄的强度图与散射场和入射光线的关系为：
Figure FDA0002530654570000025
其中I(xT)为实际所拍摄的强度图像，U(xT)为包含散射场Us(xT)和入射场Uin(xT)的总和场，xT为二维空间变量；
在三维衍射层析成像中，Born近似与Rytov近似方法简化为：
Figure FDA0002530654570000026
其中Ub(xT)表示一阶散射场，对上式进行进一步的化简，对各个变量分别除以入射分量Uin(xT)，化简后的公式可得：
Figure FDA0002530654570000027
其中UBN(xT)，USN (xT)和UN(xT)分别表示经由入射光场归一化之后的一阶散射场，散射场和测量的总场，即可表示为
Figure FDA0002530654570000028
和
Figure FDA0002530654570000029
故可以分别得到在Born近似和Rytov近似下所测量的明场和暗场的情况下，总场与一阶散射场之间关系的表达式：
对于Born近似
Figure FDA0002530654570000031
对于Rytov近似
Figure FDA0002530654570000032
利用上述关系可以分别得到在Born近似和Rytov近似下将所拍摄的强度图像
Figure FDA00025306545700000312
迭代至二维一阶散射场中的公式为：
对于Born近似
Figure FDA0002530654570000033
对于Rytov近似
Figure FDA0002530654570000034
其中
Figure FDA0002530654570000035
为原始估计的一阶散射场，
Figure FDA0002530654570000036
为迭代更新光强信息之后的散射场，
Figure FDA0002530654570000037
为实际拍摄的强度图像，系数α作为调节项可加快收敛速度减小噪声影响；
再对
Figure FDA0002530654570000038
进行二维傅里叶变换得到
Figure FDA0002530654570000039
将由强度收敛之后的一阶散射势的傅里叶变换反投影到三维频谱中对应入射光线下的埃瓦尔德球壳上，即可实现一个照明角度下的三维频谱迭代重构，即表示为：
Figure FDA00025306545700000310
经过上述公式即可完成一次对三维傅里叶频谱的强度迭代收敛过程。
2.根据权利要求1所述的基于LED阵列照明的大视场高分辨三维傅里叶叠层衍射层析显微成像方法，其特征在于在步骤三中：被测物体的初始化三维频谱
Figure FDA00025306545700000311
的像素尺寸必须满足最终的成像分辨率，且三维矩阵像素个数分别Nx,Ny,Nz满足每个方向上的最小采样数，当最大照明数值孔径为NAill，且NAill≥NAobj时，该成像方法最终可达的横向分辨率的空间频率为：
Figure FDA0002530654570000041
轴向最高可达空间频率为：
Figure FDA0002530654570000042
故最终初始化的三维频谱的像素分辨率必须大于系统的横向分辨率：
Figure FDA0002530654570000043
代入Δux,y后可得关系式：
Figure FDA0002530654570000044
其中ΔXpixel为初始化之后的三维物函数的像素尺寸，NAobj为显微物镜的数值孔径。
3.根据权利要求1所述的基于LED阵列照明的大视场高分辨三维傅里叶叠层衍射层析显微成像方法，其特征在于在步骤五中：将每个照明角度下的强度图
Figure FDA0002530654570000045
依次进行迭代，通过数次的光强迭代之后，所拍摄的强度信息逐渐收敛至三维散射势频谱中，最终即可得到所测量样片的三维频谱信息，最后进行三维逆傅里叶变换得到三维物体的物函数V(r)分布为：
Figure FDA0002530654570000046
其中IFT3表示三维逆傅里叶变换，再利用物函数散射势计算公式即可得到待测样品的三维折射率分布：
Figure FDA0002530654570000047
其中，k0为真空中对应波长的波数，nm(r)为物体所处介质的折射率。
CN201711484784.XA 2017-12-29 2017-12-29 一种大视场高分辨三维衍射层析显微成像方法 Active CN108169173B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201711484784.XA CN108169173B (zh) 2017-12-29 2017-12-29 一种大视场高分辨三维衍射层析显微成像方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201711484784.XA CN108169173B (zh) 2017-12-29 2017-12-29 一种大视场高分辨三维衍射层析显微成像方法

Publications (2)

Publication Number Publication Date
CN108169173A CN108169173A (zh) 2018-06-15
CN108169173B true CN108169173B (zh) 2020-08-11

Family

ID=62516311

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201711484784.XA Active CN108169173B (zh) 2017-12-29 2017-12-29 一种大视场高分辨三维衍射层析显微成像方法

Country Status (1)

Country Link
CN (1) CN108169173B (zh)

Families Citing this family (12)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN108919487A (zh) * 2018-08-22 2018-11-30 中国科学院重庆绿色智能技术研究院 成像曲面非均匀采样方法及doe元件设计方法和曲面成像方法
CN109580457B (zh) 2018-11-01 2021-05-04 南京理工大学 基于led阵列编码照明的三维衍射层析显微成像方法
CN111105415A (zh) * 2019-12-31 2020-05-05 北京理工大学重庆创新中心 一种基于深度学习的白细胞大视场图像检测系统及方法
CN111610621B (zh) * 2020-01-19 2022-04-08 北京大学 一种双模态显微成像系统和方法
CN111258046A (zh) * 2020-02-26 2020-06-09 清华大学 基于前置微透镜阵列的光场显微系统及方法
CN113495043B (zh) * 2020-03-18 2023-11-10 江苏长电科技股份有限公司 微尺度衍射图像检测装置
CN111445554B (zh) * 2020-05-09 2022-09-23 南京大学 一种多会聚角混合扫描电子显微成像方法
CN112580808B (zh) * 2020-12-30 2024-04-26 杭州电子科技大学 一种基于改进型led矩阵无透镜成像方法
CN114413783A (zh) * 2021-12-16 2022-04-29 中国科学院上海光学精密机械研究所 基于仪器三维传递函数的相干扫描干涉仪表面重构方法
CN114967397B (zh) * 2022-04-25 2023-04-25 上海交通大学 一种无透镜全息三维成像构建方法及装置
WO2023221741A1 (zh) * 2022-05-19 2023-11-23 南京理工大学 一种基于非干涉合成孔径的光强传输衍射层析显微成像方法
CN115144371A (zh) * 2022-07-19 2022-10-04 南京理工大学 基于波长扫描的无透镜傅里叶叠层衍射层析显微成像方法

Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN105158894A (zh) * 2015-09-29 2015-12-16 南京理工大学 基于彩色led阵列照明的无透镜相位显微层析装置及图像重构方法
CN105911692A (zh) * 2016-06-24 2016-08-31 南京理工大学 一种傅立叶叠层显微成像系统的最优系统参数选择方法
CN106204434A (zh) * 2016-06-24 2016-12-07 南京理工大学 一种面向大视场高分辨率显微成像的图像迭代重构方法
CN107065159A (zh) * 2017-03-24 2017-08-18 南京理工大学 一种基于大照明数值孔径的大视场高分辨率显微成像装置及迭代重构方法

Family Cites Families (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US10652444B2 (en) * 2012-10-30 2020-05-12 California Institute Of Technology Multiplexed Fourier ptychography imaging systems and methods
AU2016209275A1 (en) * 2015-01-21 2017-06-29 California Institute Of Technology Fourier ptychographic tomography
US10168525B2 (en) * 2015-01-26 2019-01-01 California Institute Of Technology Multi-well fourier ptychographic and fluorescence imaging

Patent Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN105158894A (zh) * 2015-09-29 2015-12-16 南京理工大学 基于彩色led阵列照明的无透镜相位显微层析装置及图像重构方法
CN105911692A (zh) * 2016-06-24 2016-08-31 南京理工大学 一种傅立叶叠层显微成像系统的最优系统参数选择方法
CN106204434A (zh) * 2016-06-24 2016-12-07 南京理工大学 一种面向大视场高分辨率显微成像的图像迭代重构方法
CN107065159A (zh) * 2017-03-24 2017-08-18 南京理工大学 一种基于大照明数值孔径的大视场高分辨率显微成像装置及迭代重构方法

Non-Patent Citations (6)

* Cited by examiner, † Cited by third party
Title
Diffraction tomography with Fourier ptychography;ROARKE HORSTMEYER等;《Optica》;20160727;第3卷(第8期);第828-831页 *
Efficient quantitative phase microscopy using programmable annular LED illumination;JIAJI LI等;《BIOMEDICAL OPTICS EXPRESS》;20170926;第8卷(第10期);第5687-4705页 *
Lensless phase microscopy and diffraction tomography with multi-angle and multi-wavelength illuminations using a LED matrix;Chao Zuo等;《OPTICS EXPRESS》;20150522;第23卷(第11期);第14314-14328页 *
Sampling criteria for Fourier ptychographic microscopy in object space and frequency space;JIASONG SUN等;《OPTICS EXPRESS》;20170705;第24卷(第14期);第15765-15781页 *
傅里叶叠层显微成像技术：理论、发展和应用;孙佳嵩等;《光学学报》;20161031;第36卷(第10期);第1011005-1至1011005-9页 *
高精度数字全息显微衍射层析成像;马利红等;《光子学报》;20141031;第43卷(第10期);第1011005-1至1011005-5页 *

Also Published As

Publication number Publication date
CN108169173A (zh) 2018-06-15

Similar Documents

Publication Publication Date Title
CN108169173B (zh) 一种大视场高分辨三维衍射层析显微成像方法
Zuo et al. Wide-field high-resolution 3D microscopy with Fourier ptychographic diffraction tomography
CN109580457B (zh) 基于led阵列编码照明的三维衍射层析显微成像方法
Kuś et al. Holographic tomography: hardware and software solutions for 3D quantitative biomedical imaging
CA2703102A1 (en) Depth of field extension for optical tomography
CN108761752A (zh) 傅立叶重叠关联成像系统、设备和方法
CN105973845A (zh) 光学测量装置和光学测量方法
US9404857B2 (en) White light diffraction tomography of unlabeled live cells
Xue et al. Deep-learning-augmented computational miniature mesoscope
US10848683B2 (en) Method for observing a sample by lensless imaging, with a spatial dispersion in the sample taken into account
US20160048010A1 (en) Structured standing wave microscope
Lee et al. Color-coded LED microscopy for quantitative phase imaging: Implementation and application to sperm motility analysis
Wang et al. Three-dimensional refractive index reconstruction of red blood cells with one-dimensional moving based on local plane wave approximation
Huang et al. Multiplane digital holography based on extrapolation iterations
Sung Snapshot three-dimensional absorption imaging of microscopic specimens
Coe et al. Computational modeling of optical projection tomographic microscopy using the finite difference time domain method
Juntunen et al. Hyperspectral three-dimensional absorption imaging using snapshot optical tomography
Sokol et al. The Use of Digital Interferometry Devices to Analyze the State of Red Blood Cell Membranes
Bai et al. Absorption and phase decoupling in transport of intensity diffraction tomography
CN114324245A (zh) 基于部分相干结构光照明的定量相位显微装置和方法
Bergoënd et al. Diffraction tomography for biological cells imaging using digital holographic microscopy
Li et al. Three-dimensional optical diffraction tomographic microscopy with optimal frequency combination with partially coherent illuminations
Wang et al. Imaging and positioning through scattering media noninvasively by bi-directional exposure
Chamgoulov et al. Optical computed-tomography micoscope using digital spatial light modulation
Li et al. Three-dimensional tomographic microscopy technique with multi-frequency combination with partially coherent illuminations

Legal Events

Date Code Title Description
PB01 Publication
SE01 Entry into force of request for substantive examination
SE01 Entry into force of request for substantive examination
GR01 Patent grant
GR01 Patent grant