时间延迟积分荧光显微成像平场校正技术

时间延迟积分荧光显微成像平场校正技术 Flat field correction technique for time delay integration fluorescence microscopy imaging 1 2 first-author 常松涛（1988-），男，副研究员，2010年于中国科学技术大学获得学士学位，2015年于中国科学院大学获得博士学位，主要从事光机系统设计方面的研究。E-mail： stchang2010@sina.com http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=25416917&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=25416919&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=25416918&type= 22.01333237 32.00399780 常松涛 （1988-），男，副研究员，2010年于中国科学技术大学获得学士学位，2015年于中国科学院大学获得博士学位，主要从事光机系统设计方面的研究。E-mail： stchang2010@sina.com stchang2010@sina.com 1 2 corresp E-mail： hjxia@hfut.edu.cn E-mail： hjxia@hfut.edu.cn 夏豪杰（1979-），男，教授，博士生导师，2006年于合肥工业大学获得博士学位，主要从事光电精密测量技术、微纳测控系统、仪器精度理论和精密仪器设计等领域的研究。E-mail： hjxia@hfut.edu.cn http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=25416920&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=25416922&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=25416921&type= 22.01333618 32.00400162 夏豪杰 （1979-），男，教授，博士生导师，2006年于合肥工业大学获得博士学位，主要从事光电精密测量技术、微纳测控系统、仪器精度理论和精密仪器设计等领域的研究。E-mail： hjxia@hfut.edu.cn hjxia@hfut.edu.cn 合肥工业大学 仪器科学与光电工程学院，安徽 合肥230009 合肥工业大学 仪器科学与光电工程学院，安徽 合肥230009 School of Instrument Science and Opto-electronics Engineering， Hefei University of Technology， Hefei 230009， China School of Instrument Science and Opto-electronics Engineering， Hefei University of Technology， Hefei 230009， China 测量理论与精密仪器安徽省重点实验室，安徽 合肥230009 测量理论与精密仪器安徽省重点实验室，安徽 合肥230009 Anhui Province Key Laboratory of Measuring Theory and Precision Instrument， Hefei 230009， China Anhui Province Key Laboratory of Measuring Theory and Precision Instrument， Hefei 230009， China 时间延迟积分（Time Delay Integration， TDI）图像传感器具有高速、高灵敏度等特点，广泛应用于高通量、大视场的荧光显微成像系统中。显微物镜视场内响应均匀是精确获取荧光能量分布的基础，为提高系统成像质量和测量准确度，研究了适用于TDI荧光显微成像系统的平场校正或响应非均匀性校正方法。根据TDI荧光成像系统的工作原理推导激光诱导荧光成像模型，分析图像均匀性退化机理。提出一种基于微阵列生物芯片的两步式校正方法，将系统响应非均匀性分为高频和低频部分分别校正，高频部分采用传统的两点校正方法，低频部分采用微阵列生物芯片校正。基于高通量TDI荧光显微成像系统开展实验，执行并验证本文的校正方法。实验结果表明：本文方法将TDI荧光成像系统的响应非均匀性由25.21%降低至2.87%，显著提高了系统性能。本文提出的校正方法能够满足TDI荧光显微成像系统的平场校正需求，具有一定的应用价值。 Owing to their high frame rate and sensitivity， time delay integration （TDI） image sensors are widely used in fluorescence microscopy imaging systems with high throughput and large field of view. Uniformity of response in the field of view of the microscope objective is the basis for accurate acquisition of fluorescence energy distribution. To improve the imaging quality and measurement accuracy of the system， this study investigated the response non-uniformity correction method applicable to the TDI fluorescence microscopy imaging system. First， a laser-induced fluorescence imaging model was devised based on the working principle of the TDI fluorescence imaging system， and the mechanism of image uniformity degradation was analyzed. A two-step correction method based on microarray biochips was then proposed， which divided the nonuniformity of system response into high frequency and low-frequency parts to be corrected separately. The former can be corrected by traditional two point correction， whereas the latter is corrected using the proposed microarray biochip-based method. Finally， experiments based on a high-throughput TDI fluorescence microscopy imaging system were conducted to verify the calibration method employed in this study. The experimental results show that the proposed approach reduces the response nonuniformity of the TDI fluorescence imaging system from 25.21% to 2.87%， implying a significant improvement in the system performance. Moreover， it shows that the correction method proposed in this paper can significantly correct the response nonuniformity of TDI fluorescent imaging systems， and is therefore effective and practical. 荧光显微成像 时间延迟积分 非均匀性校正 微阵列生物芯片 fluorescence microscopy imaging time delay integration nonuniformity correction microarray biochip 1 引　言 1 平场校正也称为响应非均匀性校正技术，是光学成像领域的一个重要研究课题，一般而言，红外成像系统和一些对成像质量要求较高的可见光系统需要进行平场校正 ［ 1 1 - 3 3 1-3 ］ 。高通量TDI荧光显微成像系统是一种对像质要求极为苛刻的可见光系统，其应用包括基因测序仪、切片荧光成像系统等 ［ 4 4 - 6 6 4-6 ］ 。高通量的实现途径是大视场和高帧频，要求全视场内响应度均匀，以实现对样本荧光能量的定量测量和分析。然而，图像传感器、光学系统等都会导致视场内的响应非均匀性，使得荧光能量分布产生异常，为了保证后续处理结果中数据的有效率和准确率，必须进行平场校正 ［ 4 4 ， 7 7 - 8 8 7-8 ］ 。 现有的非均匀性校正算法主要针对红外成像领域，可以分为两类：一是基于场景的非均匀性校正技术，二是基于定标的非均匀性校正技术。前者校正精度低，不适用于高精度光学测量系统 ［ 8 8 - 10 10 8-10 ］ 。后者利用均匀校正源（或光源）进行系统成像校正，使得视场内各像素点获得几乎相同的响应度。典型的校正方法有单点校正、两点校正和多点校正 ［ 3 3 ， 11 11 ］ 。目前，显微成像系统的校正方法仍基于传统两点校正原理，采用LED、荧光板、荧光纸和荧光膜等均匀光源对系统响应进行校正 ［ 12 12 - 13 13 12-13 ］ 。高通量时间延迟积分（Time Delay Integration， TDI）荧光显微成像系统基于激光诱导荧光成像原理，并且对响应均匀性的要求较高 ［ 14 14 - 15 15 14-15 ］ 。上述校正方法和校正源存在明显的缺陷，不适用于高通量TDI荧光显微成像系统。例如，基于LED光源的校正方案没有激光诱导荧光的成像过程，因此未考虑激光照明系统的非均匀性，也未考虑杂散光的影响，用于高通量TDI荧光显微成像系统有原理上的缺陷。荧光板、荧光纸和荧光膜等受激荧光光源也作为校正源，其原理是激光照射荧光板等激发出荧光，再由成像系统收集获得图像。其优点是对照明非均匀性进行了校正；但对于高分辨率（亚微米级）荧光成像系统而言，普通的荧光板、荧光纸和荧光膜等校准源颗粒度明显、均匀性不足且荧光易衰减。另外，因缺乏图案类目标，难以将校正源调节到物镜最佳焦面处，校正效果较差。综上，高性能、高通量的显微成像系统对响应均匀性的要求较高，并且系统基于激光诱导荧光成像原理，因此目前缺乏适用的校正源和校正方法。 为解决上述问题，本文将高通量基因测序常用的微阵列生物芯片作为标定样板，提出基于该样板的高通量TDI荧光显微成像系统校正方法。将校正参数分为高频和低频两部分，利用微阵列芯片的宏观均匀性（整体上分布均匀）计算低频非均匀性校正系数，对相机单独进行两点校正获得高频非均匀性校正系数。设计实验逐步验证并修正建立的系统响应均匀性退化数学模型，利用现有的高通量TDI荧光显微成像系统采集图像，验证本文提出的校正方法的效果，形成整套校正方案。本文提出的平场校正方法能够直接采用微阵列实现TDI荧光显微成像系统的准实时两点校正，减小相机固定图形噪声、照明系统非均匀性、成像系统非均匀性对荧光图像的影响。 2 TDI荧光显微成像系统非均匀性来源及两点校正原理 1 2.1　TDI荧光显微成像系统响应非均匀性来源 2 典型的TDI荧光显微成像系统包含APO显微物镜、激光器、滤光片、管镜、TDI相机和载物台。按照来源可将荧光显微成像系统的响应非均匀性分为三类：（1）相机自身响应非均匀性（Photo Response Nonuniformity， PRNU）；（2）激光照明系统导致的照明非均匀性；（3）成像系统导致的非均匀性。相机自身的平场校正主要由相机和图像传感器厂家执行，但对于高性能相机，需要根据应用环境重新校正，一般采用两点非均匀性校正方法。照明非均匀性会受到激光光源、照明光学系统和显微物镜的共同影响，表现为生物芯片上激光照射不均匀，进而受激辐射的荧光强度不均匀，从而影响成像均匀性。成像系统包含显微物镜、管镜、滤光片和分光镜等，导致非均匀性的因素较为复杂，也需要进行校正。相机自身的非均匀性包含高频量和低频量，而激光照明、成像导致的非均匀性绝大部分为低频量。 2.2　TDI成像两点校正原理 2 TDI相机输出为线列，非均匀性一般表现为条纹状，对一维图像进行非均匀性校正即可，计算量及复杂度小于面阵相机。TDI成像系统的输出灰度与成像目标或均匀场景辐射亮度的关系为： <math display="block" id="M1"><msub><mrow><mi>Y</mi></mrow><mrow><mi>i</mi><mo>,</mo><mi>j</mi></mrow></msub><mo>=</mo><msub><mrow><mi>G</mi></mrow><mrow><mi>i</mi><mo>,</mo><mi>j</mi></mrow></msub><mo>×</mo><mi>X</mi><mo>+</mo><msub><mrow><mi>O</mi></mrow><mrow><mi>i</mi><mo>,</mo><mi>j</mi></mrow></msub><mo>+</mo><msub><mrow><mi>σ</mi></mrow><mrow><mi>i</mi><mo>,</mo><mi>j</mi></mrow></msub></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=25416924&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=25416923&type= 41.31733322 3.80999994 ， 其中： <math id="M2"><mi>j</mi></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=25416926&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=25416925&type= 1.10066664 3.64066648 为扫描方向（也称为积分方向、电荷转移方向）的像素序号， <math id="M3"><mi>i</mi></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=25416994&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=25416993&type= 1.10066664 2.79399991 为非扫描方向的像素序号， <math id="M4"><msub><mrow><mi>Y</mi></mrow><mrow><mi>i</mi><mo>,</mo><mi>j</mi></mrow></msub></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=25416930&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=25416929&type= 5.16466665 3.72533321 为成像系统输出图像中 <math id="M5"><mo stretchy="false">(</mo><mi>i</mi><mo>,</mo><mi>j</mi><mo stretchy="false">)</mo></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=25416980&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=25416979&type= 7.11199999 3.97933316 点的灰度， <math id="M6"><msub><mrow><mi>G</mi></mrow><mrow><mi>i</mi><mo>,</mo><mi>j</mi></mrow></msub></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=25416934&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=25416933&type= 5.07999992 3.72533321 ， <math id="M7"><msub><mrow><mi>O</mi></mrow><mrow><mi>i</mi><mo>,</mo><mi>j</mi></mrow></msub></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=25416936&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=25416935&type= 5.16466665 3.72533321 ， <math id="M8"><msub><mrow><mi>σ</mi></mrow><mrow><mi>i</mi><mo>,</mo><mi>j</mi></mrow></msub></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=25416938&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=25416937&type= 3.97933316 3.72533321 分别为响应度、偏置和随机噪声， <math id="M9"><mi>X</mi></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=25416940&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=25416939&type= 2.62466669 2.79399991 为目标或场景的辐射亮度。对扫描方向上 <math id="M10"><mi>H</mi></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=25416954&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=25416953&type= 2.96333337 2.79399991 个像素求均值，以减小随机时间噪声的影响，得到： <math display="block" id="M11"><msub><mrow><mi>Y</mi></mrow><mrow><mi>i</mi></mrow></msub><mo>=</mo><msub><mrow><mi>G</mi></mrow><mrow><mi>i</mi></mrow></msub><mo>×</mo><mi>X</mi><mo>+</mo><msub><mrow><mi>O</mi></mrow><mrow><mi>i</mi></mrow></msub></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=25416944&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=25416943&type= 28.02466965 3.80999994 ， 其中： <math id="M12"><msub><mrow><mi>Y</mi></mrow><mrow><mi>i</mi></mrow></msub><mo>=</mo><mstyle displaystyle="true"><munderover><mo largeop="true">∑</mo><mrow><mi>j</mi><mo>=</mo><mn mathvariant="normal">1</mn></mrow><mrow><mi>H</mi></mrow></munderover></mstyle><mrow><msub><mrow><mi>Y</mi></mrow><mrow><mi>i</mi><mo>,</mo><mi>j</mi></mrow></msub><mo>/</mo><mi>H</mi></mrow></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=25416946&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=25416945&type= 23.28333473 8.97466660 为线列（非扫描方向）第 <math id="M13"><mi>i</mi></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=25416994&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=25416993&type= 1.10066664 2.79399991 个像素在扫描方向的灰度均值， <math id="M14"><msub><mrow><mi>G</mi></mrow><mrow><mi>i</mi></mrow></msub></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=25416950&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=25416949&type= 3.55599999 3.72533321 为响应度均值， <math id="M15"><msub><mrow><mi>O</mi></mrow><mrow><mi>i</mi></mrow></msub></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=25416952&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=25416951&type= 3.72533321 3.72533321 为平均偏置， <math id="M16"><mi>H</mi></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=25416954&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=25416953&type= 2.96333337 2.79399991 为扫描方向上的像素数。校正后每个像素的响应度和偏置相等，得到的理想校正结果为： <math display="block" id="M17"><mover><mrow><mi>Y</mi></mrow><mo>¯</mo></mover><mo>=</mo><mover><mrow><mi>G</mi></mrow><mo>¯</mo></mover><mo>×</mo><mi>X</mi><mo>+</mo><mover><mrow><mi>O</mi></mrow><mo>¯</mo></mover></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=25416956&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=25416955&type= 25.73866653 3.47133350 ， 其中： <math id="M18"><mover><mrow><mi>Y</mi></mrow><mo>¯</mo></mover><mo>=</mo><mstyle displaystyle="true"><munderover><mo largeop="true">∑</mo><mrow><mi>i</mi><mo>=</mo><mn mathvariant="normal">1</mn></mrow><mrow><mi>W</mi></mrow></munderover></mstyle><mrow><msub><mrow><mi>Y</mi></mrow><mrow><mi>i</mi></mrow></msub></mrow><mo>/</mo><mi>W</mi></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=25416958&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=25416957&type= 22.26733208 8.55133343 为输出图像的平均灰度， <math id="M19"><mover><mrow><mi>G</mi></mrow><mo>¯</mo></mover></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=25416960&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=25416959&type= 2.79399991 3.47133350 和 <math id="M20"><mover><mrow><mi>O</mi></mrow><mo>¯</mo></mover></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=25416962&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=25416961&type= 2.96333337 3.47133350 为平均增益与偏置， <math id="M21"><mi>W</mi></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=25417000&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=25416999&type= 3.80999994 2.87866688 为非扫描方向上的像素数。 两点校正的目的是使得各探测元（像素）的响应度一致，满足 式（3） 式（3） 。对于高性能图像传感器而言，像素响应灰度与接收到的能量呈线性关系，即满足线性响应。因此两点校正是较为普遍、效果较好的非均匀性校正方案，需要用到两个亮度不同的均匀辐射场景，通常取亮度接近零的暗场和亮度约使得探测器灰度达到其上限50%的明场，亮度分别设为 <math id="M22"><msup><mrow><mi>X</mi></mrow><mrow><mo stretchy="false">(</mo><mi mathvariant="normal">l</mi><mo stretchy="false">)</mo></mrow></msup></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=25416966&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=25416965&type= 5.16466665 3.04800010 和 <math id="M23"><msup><mrow><mi>X</mi></mrow><mrow><mo stretchy="false">(</mo><mi mathvariant="normal">h</mi><mo stretchy="false">)</mo></mrow></msup></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=25416968&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=25416967&type= 5.67266655 3.04800010 ，暗场通过光学系统封闭、暗室等方式实现，明场通过积分球、LED等光源均匀照射实现。由TDI成像系统响应 公式（2） 公式（2） ，得到： <math display="block" id="M24"><mfenced open="{" close="" separators="|"><mrow><mtable columnalign="left"><mtr><mtd><msubsup><mrow><mi>Y</mi></mrow><mrow><mi>i</mi></mrow><mrow><mo stretchy="false">(</mo><mi mathvariant="normal">l</mi><mo stretchy="false">)</mo></mrow></msubsup><mo>=</mo><msub><mrow><mi>G</mi></mrow><mrow><mi>i</mi></mrow></msub><mo>×</mo><msup><mrow><mi>X</mi></mrow><mrow><mo stretchy="false">(</mo><mi mathvariant="normal">l</mi><mo stretchy="false">)</mo></mrow></msup><mo>+</mo><msub><mrow><mi>O</mi></mrow><mrow><mi>i</mi></mrow></msub></mtd></mtr><mtr><mtd><msubsup><mrow><mi>Y</mi></mrow><mrow><mi>i</mi></mrow><mrow><mo stretchy="false">(</mo><mi mathvariant="normal">h</mi><mo stretchy="false">)</mo></mrow></msubsup><mo>=</mo><msub><mrow><mi>G</mi></mrow><mrow><mi>i</mi></mrow></msub><mo>×</mo><msup><mrow><mi>X</mi></mrow><mrow><mo stretchy="false">(</mo><mi mathvariant="normal">h</mi><mo stretchy="false">)</mo></mrow></msup><mo>+</mo><msub><mrow><mi>O</mi></mrow><mrow><mi>i</mi></mrow></msub></mtd></mtr></mtable></mrow></mfenced></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=25416970&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=25416969&type= 35.05199814 10.75266743 . 进而得到两点校正系数： <math display="block" id="M25"><mfenced open="{" close="" separators="|"><mrow><mtable columnalign="left"><mtr><mtd><msub><mrow><mi>α</mi></mrow><mrow><mi>i</mi></mrow></msub><mo>=</mo><mo stretchy="false">[</mo><mover><mrow><msup><mrow><mi>Y</mi></mrow><mrow><mo stretchy="false">(</mo><mi mathvariant="normal">h</mi><mo stretchy="false">)</mo></mrow></msup></mrow><mo>¯</mo></mover><mo>-</mo><mover><mrow><msup><mrow><mi>Y</mi></mrow><mrow><mo stretchy="false">(</mo><mi mathvariant="normal">l</mi><mo stretchy="false">)</mo></mrow></msup></mrow><mo>¯</mo></mover><mo stretchy="false">]</mo><mo>/</mo><mo stretchy="false">(</mo><msubsup><mrow><mi>Y</mi></mrow><mrow><mi>i</mi></mrow><mrow><mo stretchy="false">(</mo><mi mathvariant="normal">h</mi><mo stretchy="false">)</mo></mrow></msubsup><mo>-</mo><msubsup><mrow><mi>Y</mi></mrow><mrow><mi>i</mi></mrow><mrow><mo stretchy="false">(</mo><mi mathvariant="normal">l</mi><mo stretchy="false">)</mo></mrow></msubsup><mo stretchy="false">)</mo></mtd></mtr><mtr><mtd><msub><mrow><mi>β</mi></mrow><mrow><mi>i</mi></mrow></msub><mo>=</mo><mover><mrow><msup><mrow><mi>Y</mi></mrow><mrow><mo stretchy="false">(</mo><mi mathvariant="normal">h</mi><mo stretchy="false">)</mo></mrow></msup></mrow><mo>¯</mo></mover><mo>-</mo><msub><mrow><mi>α</mi></mrow><mrow><mi>i</mi></mrow></msub><mo>×</mo><msubsup><mrow><mi>Y</mi></mrow><mrow><mi>i</mi></mrow><mrow><mo stretchy="false">(</mo><mi mathvariant="normal">h</mi><mo stretchy="false">)</mo></mrow></msubsup></mtd></mtr></mtable></mrow></mfenced></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=25416972&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=25416971&type= 52.15466690 12.02266598 ， 其中： <math id="M26"><mover><mrow><msup><mrow><mi>Y</mi></mrow><mrow><mo stretchy="false">(</mo><mi mathvariant="normal">h</mi><mo stretchy="false">)</mo></mrow></msup></mrow><mo>¯</mo></mover><mo>=</mo><mstyle displaystyle="true"><munderover><mo largeop="true">∑</mo><mrow><mi>i</mi><mo>=</mo><mn mathvariant="normal">1</mn></mrow><mrow><mi>W</mi></mrow></munderover></mstyle><mrow><msubsup><mrow><mi>Y</mi></mrow><mrow><mi>i</mi></mrow><mrow><mo stretchy="false">(</mo><mi mathvariant="normal">h</mi><mo stretchy="false">)</mo></mrow></msubsup></mrow><mo>/</mo><mi>W</mi></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=25416974&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=25416973&type= 27.85533524 8.55133343 和 <math id="M27"><mover><mrow><msup><mrow><mi>Y</mi></mrow><mrow><mo stretchy="false">(</mo><mi mathvariant="normal">l</mi><mo stretchy="false">)</mo></mrow></msup></mrow><mo>¯</mo></mover><mo>=</mo><mstyle displaystyle="true"><munderover><mo largeop="true">∑</mo><mrow><mi>i</mi><mo>=</mo><mn mathvariant="normal">1</mn></mrow><mrow><mi>W</mi></mrow></munderover></mstyle><mrow><msubsup><mrow><mi>Y</mi></mrow><mrow><mi>i</mi></mrow><mrow><mo stretchy="false">(</mo><mi mathvariant="normal">l</mi><mo stretchy="false">)</mo></mrow></msubsup></mrow><mo>/</mo><mi>W</mi></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=25416976&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=25416975&type= 26.75466537 8.55133343 分别为明场和暗场的灰度均值。 目前，TDI荧光显微成像系统的校正方法仍基于传统两点校正原理，采用LED、荧光板、荧光纸、荧光膜等均匀光源对系统响应进行校正。 3 基于微阵列生物芯片的校正方法 1 3.1　TDI荧光成像系统响应均匀性退化模型 2 TDI荧光成像系统基于激光诱导荧光成像原理，激光经照明光学系统整形后透过物镜均匀投射到生物芯片上，生物芯片上的荧光分子受激发射荧光，荧光经物镜、管镜收集后到达相机，得到图像。根据上述过程建立生物芯片成像及光机系统响应均匀性退化模型。设未校正的暗场图像为 <math id="M28"><msup><mrow><mi>I</mi></mrow><mrow><mo stretchy="false">(</mo><mn mathvariant="normal">0</mn><mo stretchy="false">)</mo></mrow></msup></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=25417153&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=25417151&type= 4.23333359 3.04800010 ， <math id="M29"><mo stretchy="false">(</mo><mi>i</mi><mo>,</mo><mi>j</mi><mo stretchy="false">)</mo></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=25416980&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=25416979&type= 7.11199999 3.97933316 像素点灰度为 <math id="M30"><msubsup><mrow><mi>Y</mi></mrow><mrow><mi>i</mi><mo>,</mo><mi>j</mi></mrow><mrow><mo stretchy="false">(</mo><mn mathvariant="normal">0</mn><mo stretchy="false">)</mo></mrow></msubsup></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=25416982&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=25416981&type= 6.09600019 3.97933316 ，为简化校正过程，暗场通过多层遮光实现，相机接收到的辐亮度几乎为零。相机直接对明场采集图像，得到： <math display="block" id="M31"><msubsup><mrow><mi>Y</mi></mrow><mrow><mi>i</mi></mrow><mrow><mo stretchy="false">(</mo><mn mathvariant="normal">1</mn><mo stretchy="false">)</mo></mrow></msubsup><mo>=</mo><msubsup><mrow><mi>k</mi></mrow><mrow><mi>i</mi></mrow><mrow><mo stretchy="false">(</mo><mn mathvariant="normal">1</mn><mo stretchy="false">)</mo></mrow></msubsup><mo>×</mo><mi>X</mi><mo>+</mo><msubsup><mrow><mi>Y</mi></mrow><mrow><mi>i</mi></mrow><mrow><mo stretchy="false">(</mo><mn mathvariant="normal">0</mn><mo stretchy="false">)</mo></mrow></msubsup></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=25416984&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=25416983&type= 33.52799988 4.06400013 ， 其中：偏置为 <math id="M32"><msubsup><mrow><mi>Y</mi></mrow><mrow><mi>i</mi></mrow><mrow><mo stretchy="false">(</mo><mn mathvariant="normal">0</mn><mo stretchy="false">)</mo></mrow></msubsup><mo>=</mo><mstyle displaystyle="true"><munderover><mo largeop="true">∑</mo><mrow><mi>j</mi><mo>=</mo><mn mathvariant="normal">1</mn></mrow><mrow><mi>H</mi></mrow></munderover></mstyle><mrow><msubsup><mrow><mi>Y</mi></mrow><mrow><mi>i</mi><mo>,</mo><mi>j</mi></mrow><mrow><mo stretchy="false">(</mo><mn mathvariant="normal">0</mn><mo stretchy="false">)</mo></mrow></msubsup><mo>/</mo><mi>H</mi></mrow></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=25416986&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=25416985&type= 26.58533096 8.97466660 ，增益或响应度为 <math id="M33"><msubsup><mrow><mi>k</mi></mrow><mrow><mi>i</mi></mrow><mrow><mo stretchy="false">(</mo><mn mathvariant="normal">1</mn><mo stretchy="false">)</mo></mrow></msubsup></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=25416988&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=25416987&type= 4.48733330 3.97933316 。荧光显微成像系统工作在线性区，即激光能量与激发出的荧光能量呈线性关系。光源和照明光学系统导致入射到生物芯片上的激光能量不均匀，假设生物芯片上的荧光分子是均布的，那么激发出的荧光能量也会呈现非均匀性，即： <math display="block" id="M34"><msubsup><mrow><mi>k</mi></mrow><mrow><mi>i</mi></mrow><mrow><mo stretchy="false">(</mo><mn mathvariant="normal">2</mn><mo stretchy="false">)</mo></mrow></msubsup><mo>=</mo><mfrac><mrow><msub><mrow><mi>E</mi></mrow><mrow><mi>i</mi></mrow></msub><mo>×</mo><msub><mrow><mi>η</mi></mrow><mrow><mi>i</mi></mrow></msub></mrow><mrow><mover accent="true"><mi>E</mi><mo>¯</mo></mover><mo>×</mo><mover accent="true"><mi>η</mi><mo>¯</mo></mover></mrow></mfrac></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=25416990&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=25416989&type= 20.91266632 9.31333351 ， 其中： <math id="M35"><msub><mrow><mi>E</mi></mrow><mrow><mi>i</mi></mrow></msub></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=25416992&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=25416991&type= 3.38666677 3.72533321 为 <math id="M36"><mi>i</mi></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=25416994&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=25416993&type= 1.10066664 2.79399991 像素对应的物方（芯片）区域接收到的激光强度， <math id="M37"><msub><mrow><mi>η</mi></mrow><mrow><mi>i</mi></mrow></msub></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=25416996&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=25416995&type= 2.62466669 4.57200003 为荧光分子的发光效率， <math id="M38"><mover accent="true"><mi>E</mi><mo>¯</mo></mover><mo>=</mo><mstyle displaystyle="true"><munderover><mo largeop="true">∑</mo><mrow><mi>i</mi><mo>=</mo><mn mathvariant="normal">1</mn></mrow><mrow><mi>W</mi></mrow></munderover></mstyle><mrow><msub><mrow><mi>E</mi></mrow><mrow><mi>i</mi></mrow></msub></mrow><mo>/</mo><mi>W</mi></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=25416998&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=25416997&type= 21.67466736 8.55133343 为 <math id="M39"><mi>W</mi></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=25417000&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=25416999&type= 3.80999994 2.87866688 个像素对应物方区域激光强度的均值， <math id="M40"><mover accent="true"><mi>η</mi><mo>¯</mo></mover></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=25417002&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=25417001&type= 1.86266661 3.64066648 为相应位置荧光分子发光效率的均值。加入照明非均匀性后成像系统的响应度为： <math display="block" id="M41"><msubsup><mrow><mi>Y</mi></mrow><mrow><mi>i</mi></mrow><mrow><mo stretchy="false">(</mo><mn mathvariant="normal">2</mn><mo stretchy="false">)</mo></mrow></msubsup><mo>=</mo><msubsup><mrow><mi>k</mi></mrow><mrow><mi>i</mi></mrow><mrow><mo stretchy="false">(</mo><mn mathvariant="normal">2</mn><mo stretchy="false">)</mo></mrow></msubsup><mo>×</mo><msubsup><mrow><mi>k</mi></mrow><mrow><mi>i</mi></mrow><mrow><mo stretchy="false">(</mo><mn mathvariant="normal">1</mn><mo stretchy="false">)</mo></mrow></msubsup><mo>×</mo><mi>X</mi><mo>+</mo><msubsup><mrow><mi>Y</mi></mrow><mrow><mi>i</mi></mrow><mrow><mo stretchy="false">(</mo><mn mathvariant="normal">0</mn><mo stretchy="false">)</mo></mrow></msubsup></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=25417005&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=25417003&type= 43.01066971 4.06400013 . 成像光学系统包括物镜、滤光片、分光镜和管镜，导致非均匀性的因素包括畸变、渐晕、离焦和装调角度误差等，表现为中心视场的响应度比边缘视场略大，其非均匀性系数用 <math id="M42"><msubsup><mrow><mi>k</mi></mrow><mrow><mi>i</mi></mrow><mrow><mo stretchy="false">(</mo><mn mathvariant="normal">3</mn><mo stretchy="false">)</mo></mrow></msubsup></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=25417007&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=25417006&type= 4.57200003 3.97933316 表示，那么响应方程改写为： <math display="block" id="M43"><msubsup><mrow><mi>Y</mi></mrow><mrow><mi>i</mi></mrow><mrow><mo stretchy="false">(</mo><mn mathvariant="normal">3</mn><mo stretchy="false">)</mo></mrow></msubsup><mo>=</mo><msubsup><mrow><mi>k</mi></mrow><mrow><mi>i</mi></mrow><mrow><mo stretchy="false">(</mo><mn mathvariant="normal">3</mn><mo stretchy="false">)</mo></mrow></msubsup><mo>×</mo><msubsup><mrow><mi>k</mi></mrow><mrow><mi>i</mi></mrow><mrow><mo stretchy="false">(</mo><mn mathvariant="normal">2</mn><mo stretchy="false">)</mo></mrow></msubsup><mo>×</mo><msubsup><mrow><mi>k</mi></mrow><mrow><mi>i</mi></mrow><mrow><mo stretchy="false">(</mo><mn mathvariant="normal">1</mn><mo stretchy="false">)</mo></mrow></msubsup><mo>×</mo><mi>X</mi><mo>+</mo><msubsup><mrow><mi>Y</mi></mrow><mrow><mi>i</mi></mrow><mrow><mo stretchy="false">(</mo><mn mathvariant="normal">0</mn><mo stretchy="false">)</mo></mrow></msubsup></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=25417010&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=25417009&type= 52.32400131 4.06400013 . 3.2　新校正方法的理论依据及校正源设计 2 TDI荧光显微成像系统平场校正的难点是设计合适的校正源。原理上，响应度非均匀性包括高频量和低频量，相机引起的非均匀性包含高频和低频，而成像光学系统、激光照明系统引起的非均匀性绝大部分是低频量。本文将低频和高频量分步校正，采用传统的两点校正对相机自身的高频和低频非均匀性进行校正，然后设计一种微阵列生物芯片（一种基因测序使用的非均匀荧光光源）来校正残余的低频非均匀性，校正结果叠加可以实现TDI荧光成像系统的响应非均匀性校正。微阵列芯片如 图1 图1 所示，以光刻方式制备网格线，位置精度优于±10 nm，网格处无荧光样本，在荧光图像中会形成暗线，作为TDI扫描方向调整的参考线。信号点能量在统计上是均匀的，采用多点累加的方式获得各位置每个信号点的灰度。中心与边缘的差异性由像差和渐晕等因素导致，可用于估算这些因素的影响，进而校正响应非均匀性。 利用信号点低频的均匀性，对微阵列生物芯片图像进行低通滤波后获取照明系统和成像系统导致的低频非均匀性，即可实现TDI荧光成像系统的低频非均匀性校正。 3.3　基于微阵列芯片的校正方法 2 在建立的系统响应及均匀性退化数学模型的基础上开展校正工作，实现此类特殊成像系统的平场校正。校正工作实际上是响应均匀性退化模型的逆向推演。校正步骤如下：（1）采集暗场图像 <math id="M44"><msup><mrow><mi>I</mi></mrow><mrow><mo stretchy="false">(</mo><mn mathvariant="normal">0</mn><mo stretchy="false">)</mo></mrow></msup></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=25417015&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=25417014&type= 4.23333359 3.04800010 ，获得相机自身FPN；（2）使用积分球光源或者LED光源照射相机，获得均匀场景或者非均匀平滑场景，采集图像 <math id="M45"><msup><mrow><mi>I</mi></mrow><mrow><mo stretchy="false">(</mo><mn mathvariant="normal">1</mn><mo stretchy="false">)</mo></mrow></msup></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=25417020&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=25417017&type= 4.06400013 3.04800010 ；（3）TDI扫描成像采集微阵列生物芯片图像 <math id="M46"><msup><mrow><mi>I</mi></mrow><mrow><mo stretchy="false">(</mo><mn mathvariant="normal">2</mn><mo stretchy="false">)</mo></mrow></msup></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=25417027&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=25417023&type= 4.23333359 3.04800010 ；（4）经过低通滤波、归一化等处理后，平滑拟合获得低频校正系数，与相机自身校正获得的高频校正系数合并后获得综合性的增益、偏置校正系数。校正工作流程如 图2 图2 所示。 微阵列生物芯片图像减暗背景，得到： <math display="block" id="M47"><mover accent="true"><mrow><msubsup><mrow><mi>I</mi></mrow><mrow><mi>i</mi></mrow><mrow><mi>a</mi></mrow></msubsup></mrow><mo>̂</mo></mover><mo>=</mo><msubsup><mrow><mi>I</mi></mrow><mrow><mi>i</mi></mrow><mrow><mo stretchy="false">(</mo><mn mathvariant="normal">2</mn><mo stretchy="false">)</mo></mrow></msubsup><mo>-</mo><msubsup><mrow><mi>I</mi></mrow><mrow><mi>i</mi></mrow><mrow><mo stretchy="false">(</mo><mn mathvariant="normal">0</mn><mo stretchy="false">)</mo></mrow></msubsup></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=25417043&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=25417039&type= 21.25133324 4.40266657 . 相机自身校正时，去除高频响应度非均匀性。由于此处未引入照明非均匀性和成像光学系统的影响，仅能去除非均匀性中的高频量和少量低频量。高频校正系数为： <math display="block" id="M48"><msubsup><mrow><mi>k</mi></mrow><mrow><mi>i</mi></mrow><mrow><mi>b</mi></mrow></msubsup><mo>=</mo><mfrac><mrow><mover><mrow><msup><mrow><mi>I</mi></mrow><mrow><mo stretchy="false">(</mo><mn mathvariant="normal">1</mn><mo stretchy="false">)</mo></mrow></msup></mrow><mo>¯</mo></mover><mo>-</mo><mover><mrow><msup><mrow><mi>I</mi></mrow><mrow><mo stretchy="false">(</mo><mn mathvariant="normal">0</mn><mo stretchy="false">)</mo></mrow></msup></mrow><mo>¯</mo></mover></mrow><mrow><msubsup><mrow><mi>I</mi></mrow><mrow><mi>i</mi></mrow><mrow><mo stretchy="false">(</mo><mn mathvariant="normal">1</mn><mo stretchy="false">)</mo></mrow></msubsup><mo>-</mo><msubsup><mrow><mi>I</mi></mrow><mrow><mi>i</mi></mrow><mrow><mo stretchy="false">(</mo><mn mathvariant="normal">0</mn><mo stretchy="false">)</mo></mrow></msubsup></mrow></mfrac></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=25417054&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=25417048&type= 22.09799957 9.73666668 . 高频校正后的微阵列图像为： <math display="block" id="M49"><mover accent="true"><mrow><msubsup><mrow><mi>I</mi></mrow><mrow><mi>i</mi></mrow><mrow><mi>b</mi></mrow></msubsup></mrow><mo>̂</mo></mover><mo>=</mo><mfrac><mrow><mover><mrow><msup><mrow><mi>I</mi></mrow><mrow><mo stretchy="false">(</mo><mn mathvariant="normal">1</mn><mo stretchy="false">)</mo></mrow></msup></mrow><mo>¯</mo></mover><mo>-</mo><mover><mrow><msup><mrow><mi>I</mi></mrow><mrow><mo stretchy="false">(</mo><mn mathvariant="normal">0</mn><mo stretchy="false">)</mo></mrow></msup></mrow><mo>¯</mo></mover></mrow><mrow><msubsup><mrow><mi>I</mi></mrow><mrow><mi>i</mi></mrow><mrow><mo stretchy="false">(</mo><mn mathvariant="normal">1</mn><mo stretchy="false">)</mo></mrow></msubsup><mo>-</mo><msubsup><mrow><mi>I</mi></mrow><mrow><mi>i</mi></mrow><mrow><mo stretchy="false">(</mo><mn mathvariant="normal">0</mn><mo stretchy="false">)</mo></mrow></msubsup></mrow></mfrac><mo>×</mo><mover accent="true"><mrow><msubsup><mrow><mi>I</mi></mrow><mrow><mi>i</mi></mrow><mrow><mi>a</mi></mrow></msubsup></mrow><mo>̂</mo></mover></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=25417059&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=25417056&type= 29.54866600 9.73666668 . 采用微阵列生物芯片成像，如 图3 图3 所示。为限制照明非均匀性中的高频信息，可以使得激光照明（临界照明方式）有一定离焦。 图3 图3 在扫描方向取均值，局部放大后如 图4 图4 所示，对每一到两个周期进行低通滤波后，得到灰度分布的骨骼线，即为低频非均匀性分布。 滤波后得到的低频校正系数为： <math display="block" id="M50"><msubsup><mrow><mi>k</mi></mrow><mrow><mi>i</mi></mrow><mrow><mi>c</mi></mrow></msubsup><mo>=</mo><mover><mrow><msup><mrow><mi>I</mi></mrow><mrow><mi>b</mi></mrow></msup></mrow><mo>¯</mo></mover><mo>/</mo><mover accent="true"><mrow><msubsup><mrow><mi>I</mi></mrow><mrow><mi>i</mi></mrow><mrow><mi>b</mi></mrow></msubsup></mrow><mo>̂</mo></mover></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=25417089&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=25417085&type= 16.17133331 4.40266657 ， 其中 <math id="M51"><mover><mrow><msup><mrow><mi>I</mi></mrow><mrow><mi>b</mi></mrow></msup></mrow><mo>¯</mo></mover><mo>=</mo><mstyle displaystyle="true"><munderover><mo largeop="true">∑</mo><mrow><mi>i</mi><mo>=</mo><mn mathvariant="normal">1</mn></mrow><mrow><mi>L</mi></mrow></munderover></mstyle><mrow><mover accent="true"><mrow><msubsup><mrow><mi>I</mi></mrow><mrow><mi>i</mi></mrow><mrow><mi>b</mi></mrow></msubsup></mrow><mo>̂</mo></mover></mrow><mo>/</mo><mi>L</mi></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=25417095&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=25417093&type= 19.81200027 8.55133343 。校正后的图像灰度分布为： <math display="block" id="M52"><mover accent="true"><mrow><msubsup><mrow><mi>I</mi></mrow><mrow><mi>i</mi></mrow><mrow><mi>d</mi></mrow></msubsup></mrow><mo>̂</mo></mover><mo>=</mo><msubsup><mrow><mi>k</mi></mrow><mrow><mi>i</mi></mrow><mrow><mi>c</mi></mrow></msubsup><mo>×</mo><mover accent="true"><mrow><msubsup><mrow><mi>I</mi></mrow><mrow><mi>i</mi></mrow><mrow><mi>b</mi></mrow></msubsup></mrow><mo>̂</mo></mover><mo>=</mo><msubsup><mrow><mi>k</mi></mrow><mrow><mi>i</mi></mrow><mrow><mi>c</mi></mrow></msubsup><mo>×</mo><mfrac><mrow><mover><mrow><msup><mrow><mi>I</mi></mrow><mrow><mo stretchy="false">(</mo><mn mathvariant="normal">1</mn><mo stretchy="false">)</mo></mrow></msup></mrow><mo>¯</mo></mover><mo>-</mo><mover><mrow><msup><mrow><mi>I</mi></mrow><mrow><mo stretchy="false">(</mo><mn mathvariant="normal">0</mn><mo stretchy="false">)</mo></mrow></msup></mrow><mo>¯</mo></mover></mrow><mrow><msubsup><mrow><mi>I</mi></mrow><mrow><mi>i</mi></mrow><mrow><mo stretchy="false">(</mo><mn mathvariant="normal">1</mn><mo stretchy="false">)</mo></mrow></msubsup><mo>-</mo><msubsup><mrow><mi>I</mi></mrow><mrow><mi>i</mi></mrow><mrow><mo stretchy="false">(</mo><mn mathvariant="normal">0</mn><mo stretchy="false">)</mo></mrow></msubsup></mrow></mfrac><mo>×</mo><mo stretchy="false">[</mo><msubsup><mrow><mi>I</mi></mrow><mrow><mi>i</mi></mrow><mrow><mo stretchy="false">(</mo><mn mathvariant="normal">2</mn><mo stretchy="false">)</mo></mrow></msubsup><mo>-</mo><msubsup><mrow><mi>I</mi></mrow><mrow><mi>i</mi></mrow><mrow><mo stretchy="false">(</mo><mn mathvariant="normal">0</mn><mo stretchy="false">)</mo></mrow></msubsup><mo stretchy="false">]</mo></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=25417100&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=25417098&type= 66.37866211 9.73666668 . 因此，校正系数为： <math display="block" id="M53"><mfenced open="{" close="" separators="|"><mrow><mtable columnalign="left"><mtr><mtd><msub><mrow><mi>α</mi></mrow><mrow><mi>i</mi></mrow></msub><mo>=</mo><msubsup><mrow><mi>k</mi></mrow><mrow><mi>i</mi></mrow><mrow><mi>c</mi></mrow></msubsup><mo>×</mo><mfrac><mrow><mover><mrow><msup><mrow><mi>I</mi></mrow><mrow><mo stretchy="false">(</mo><mn mathvariant="normal">1</mn><mo stretchy="false">)</mo></mrow></msup></mrow><mo>¯</mo></mover><mo>-</mo><mover><mrow><msup><mrow><mi>I</mi></mrow><mrow><mo stretchy="false">(</mo><mn mathvariant="normal">0</mn><mo stretchy="false">)</mo></mrow></msup></mrow><mo>¯</mo></mover></mrow><mrow><msubsup><mrow><mi>I</mi></mrow><mrow><mi>i</mi></mrow><mrow><mo stretchy="false">(</mo><mn mathvariant="normal">1</mn><mo stretchy="false">)</mo></mrow></msubsup><mo>-</mo><msubsup><mrow><mi>I</mi></mrow><mrow><mi>i</mi></mrow><mrow><mo stretchy="false">(</mo><mn mathvariant="normal">0</mn><mo stretchy="false">)</mo></mrow></msubsup></mrow></mfrac></mtd></mtr><mtr><mtd><msub><mrow><mi>β</mi></mrow><mrow><mi>i</mi></mrow></msub><mo>=</mo><mo>-</mo><msub><mrow><mi>α</mi></mrow><mrow><mi>i</mi></mrow></msub><mo>×</mo><msubsup><mrow><mi>I</mi></mrow><mrow><mi>i</mi></mrow><mrow><mo stretchy="false">(</mo><mn mathvariant="normal">0</mn><mo stretchy="false">)</mo></mrow></msubsup></mtd></mtr></mtable></mrow></mfenced></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=25417107&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=25417104&type= 31.24200058 16.34066582 . 设 <math id="M54"><msubsup><mrow><mi>I</mi></mrow><mrow><mi>i</mi></mrow><mrow><mo stretchy="false">(</mo><mn mathvariant="normal">2</mn><mo stretchy="false">)</mo></mrow></msubsup></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=25417114&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=25417110&type= 4.40266657 3.97933316 为待校正的微阵列生物芯片图像，则校正方程为： <math display="block" id="M55"><mover accent="true"><mrow><msubsup><mrow><mi>I</mi></mrow><mrow><mi>i</mi></mrow><mrow><mo stretchy="false">(</mo><mn mathvariant="normal">2</mn><mo stretchy="false">)</mo></mrow></msubsup></mrow><mo>̂</mo></mover><mo>=</mo><msub><mrow><mi>α</mi></mrow><mrow><mi>i</mi></mrow></msub><mo>×</mo><msubsup><mrow><mi>I</mi></mrow><mrow><mi>i</mi></mrow><mrow><mo stretchy="false">(</mo><mn mathvariant="normal">2</mn><mo stretchy="false">)</mo></mrow></msubsup><mo>+</mo><msub><mrow><mi>β</mi></mrow><mrow><mi>i</mi></mrow></msub></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=25417123&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=25417118&type= 28.61733437 4.99533367 . 4 实　验 1 为验证本文提出的平场校正方法，基于某TDI荧光显微成像系统（基因测序成像系统）开展实验。实验系统的工作原理如 图5 图5 所示，包括生物芯片（均布间隔为500 nm的携带荧光分子的样本点）、显微物镜（NA为1.0，视场约为2 mm）、双色分光镜（反射激光，透射荧光）、自聚焦滤光片、激发滤光片、发射滤光片（滤除激光，透过荧光）、管镜（接收荧光）和相机。采用DALSA公司的TDI相机HL-HM-16K30H-00R，工作行频为200 kHz，积分级数为128，最大分辨率为16 384 pixel。 采集暗场图像和明场图像，对相机自身的响应非均匀性进行校正，去除非均匀性中的高频量。校正装置如 图6 图6 所示，采用蓝菲积分球作为明场光源，型号为USLR-L20F-NBNL，开口面均匀性优于±1%，积分球直径约为500 mm，开口尺寸约为200 nm。积分球经过照度、色温、光谱辐亮度、均匀性的校准、检测和计量。 得到的暗场图像 <math id="M56"><msup><mrow><mi>I</mi></mrow><mrow><mo stretchy="false">(</mo><mn mathvariant="normal">0</mn><mo stretchy="false">)</mo></mrow></msup></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=25417153&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=25417151&type= 4.23333359 3.04800010 和明场图像 <math id="M57"><msup><mrow><mi>I</mi></mrow><mrow><mo stretchy="false">(</mo><mn mathvariant="normal">1</mn><mo stretchy="false">)</mo></mrow></msup></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=25417161&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=25417157&type= 4.06400013 3.04800010 如 图7 图7 所示，可见暗场存在明显的高频、低频DSNU（或FPN），明场体现出显著的高频响应非均匀性（PRNU）。 1;2; 对微阵列芯片成像，得到图像 <math id="M58"><msup><mrow><mi>I</mi></mrow><mrow><mo stretchy="false">(</mo><mn mathvariant="normal">2</mn><mo stretchy="false">)</mo></mrow></msup></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=25417194&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=25417191&type= 4.23333359 3.04800010 ，代入 式（11） 式（11） 和（13）得到 <math id="M59"><mover accent="true"><mrow><msubsup><mrow><mi>I</mi></mrow><mrow><mi>i</mi></mrow><mrow><mi>b</mi></mrow></msubsup></mrow><mo>̂</mo></mover></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=25417201&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=25417198&type= 2.87866688 4.40266657 ，在扫描方向取灰度均值，得到分布如 图8 图8 所示。较窄的突变暗线为网格线，网格线越窄说明TDI方向调整越准确，边缘和中心的低频非均匀性主要由照明系统和成像系统导致。 去除网格线、低通滤波，采用样条曲线平滑拟合后得到TDI荧光成像系统的低频响应非均匀性，如 图9 图9 所示。对 y 轴归一化后代入 式（14） 式（14） 即可得到 <math id="M60"><msubsup><mrow><mi>k</mi></mrow><mrow><mi>i</mi></mrow><mrow><mi>c</mi></mrow></msubsup></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=25417218&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=25417214&type= 2.87866688 3.97933316 ，进而可以计算校正系数 <math id="M61"><msub><mrow><mi>α</mi></mrow><mrow><mi>i</mi></mrow></msub></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=25417224&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=25417221&type= 2.70933342 3.72533321 和 <math id="M62"><msub><mrow><mi>β</mi></mrow><mrow><mi>i</mi></mrow></msub></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=25417231&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=25417228&type= 2.70933342 4.57200003 。校正系数可直接写入相机配置文件中，实现实时平场校正处理。对HL-HM-16K30H-00R相机进行非均匀性校正测试，不校正时其图像行输出频率为200 kHz，即每行输出时间间隔为5 μs，结果显示加入本文校正方法后仍能保证200 kHz的行输出频率，可见每行图像处理速度已足够快（行运算时间小于5 μs），即能够保证实时输出校正后的图像。 图8 图8 中的微阵列芯片中样本点之间的亮度差异性导致高频量，样本点在整体上是均匀的，因此通常可用 图9 图9 中的非均匀性评估TDI荧光成像系统的响应均匀程度。校正后的图像灰度分布如 图10 图10 所示，校正效果良好。 校正前后的图像如 图11 图11 所示，为便于目视对比，将校正前后图像采用同样的拉伸处理后显示，可见校正后响应均匀性显著提升。 为验证算法的稳定性，获取不同生物芯片、不同强度的图像，采用两种定量评价指标表征本文校正后的明场响应非均匀性、非均匀性（Nonuniformity， NU）和均方根误差（Root-Mean-Square Error， RMSE） ［ 16 16 - 17 17 16-17 ］ ： <math display="block" id="M63"><mi>N</mi><mi>U</mi><mo>=</mo><mfrac><mrow><mi mathvariant="normal">m</mi><mi mathvariant="normal">a</mi><mi mathvariant="normal">x</mi><mo stretchy="false">(</mo><msub><mrow><mi>Y</mi></mrow><mrow><mi>i</mi></mrow></msub><mo stretchy="false">)</mo><mo>-</mo><mi mathvariant="normal">m</mi><mi mathvariant="normal">i</mi><mi mathvariant="normal">n</mi><mo stretchy="false">(</mo><msub><mrow><mi>Y</mi></mrow><mrow><mi>i</mi></mrow></msub><mo stretchy="false">)</mo></mrow><mrow><mstyle displaystyle="true"><munderover><mo largeop="true">∑</mo><mrow><mi>i</mi><mo>=</mo><mn mathvariant="normal">1</mn></mrow><mrow><mi>W</mi></mrow></munderover></mstyle><mrow><msub><mrow><mi>Y</mi></mrow><mrow><mi>i</mi></mrow></msub><mo>/</mo><mi>W</mi></mrow></mrow></mfrac><mo>×</mo><mn mathvariant="normal">100</mn><mi mathvariant="normal">%</mi></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=25417269&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=25417266&type= 56.64199829 13.97000027 ， <math display="block" id="M64"><mtable><mtr><mtd><mi>R</mi><mi>M</mi><mi>S</mi><mi>E</mi><mo>=</mo><mfrac><mrow><mn mathvariant="normal">1</mn></mrow><mrow><mstyle displaystyle="true"><munderover><mo largeop="true">∑</mo><mrow><mi>i</mi><mo>=</mo><mn mathvariant="normal">1</mn></mrow><mrow><mi>W</mi></mrow></munderover></mstyle><mrow><msub><mrow><mi>Y</mi></mrow><mrow><mi>i</mi></mrow></msub><mo>/</mo><mi>W</mi></mrow></mrow></mfrac><mo>×</mo></mtd></mtr><mtr><mtd><mroot><mrow><mfrac><mrow><mn mathvariant="normal">1</mn></mrow><mrow><mi>W</mi></mrow></mfrac><mstyle displaystyle="true"><munderover><mo largeop="true">∑</mo><mrow><mi>i</mi><mo>=</mo><mn mathvariant="normal">1</mn></mrow><mrow><mi>W</mi></mrow></munderover></mstyle><mrow><mo stretchy="false">(</mo><msub><mrow><mi>Y</mi></mrow><mrow><mi>i</mi></mrow></msub><msup><mrow><mo>-</mo><mstyle displaystyle="true"><munderover><mo largeop="true">∑</mo><mrow><mi>i</mi><mo>=</mo><mn mathvariant="normal">1</mn></mrow><mrow><mi>W</mi></mrow></munderover></mstyle><mrow><msub><mrow><mi>Y</mi></mrow><mrow><mi>i</mi></mrow></msub><mo>/</mo><mi>W</mi></mrow><mo stretchy="false">)</mo></mrow><mrow><mn mathvariant="normal">2</mn></mrow></msup></mrow></mrow><mrow/></mroot><mo>×</mo><mn mathvariant="normal">100</mn><mi mathvariant="normal">%</mi></mtd></mtr></mtable></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=25417274&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=25417272&type= 54.35599899 25.31533241 . NU和RMSE都是将图像网格线去除、扫描方向灰度均值处理后得到非扫描方向的响应非均匀性。其中，NU表示响应非均匀性的“峰谷值”评价标准，RMSE表示响应非均匀性的“均方根值”评价标准，均为常规的评价方法。NU和RMSE越小，响应均匀性越好。采集11幅图像，校正前后的非均匀性如 图12 图12 ~ 图13 图13 所示。 由实验结果可知，本文提出的校正方法使得荧光成像系统的响应非均匀性NU从25.21%降低到2.87%，RMSE从3.12%降低到0.48%，校正效果显著。该方法易于实施，广泛适用于各类高性能TDI荧光成像系统。 5 结　论 1 本文根据高性能TDI荧光成像系统的响应均匀性要求，提出一种基于微阵列生物芯片的平场校正方法。建立系统响应非均匀性产生模型，在此理论基础上提出将非均匀性分为高频和低频分别校正，设计了专用的微阵列芯片作为低频非均匀性校正源，将照明系统、成像光学系统和相机自身的响应非均匀性都进行了校正。实验表明，本文方法能够有效地完成TDI荧光成像非均匀性校正，响应非均匀性NU从25.21%降低到2.87%，RMSE从3.12%降低到0.48%，提高了图像质量。此外，由于微阵列芯片是基因测序领域常用的芯片，校正源易于获取，校正方法的实现较为经济、便捷。下一步工作的重点是寻求更为稳定、均匀的校正源。 参考文献 GUTSCHWAGER B ， HOLLANDT J . Nonuniformity correction of imaging systems with a spatially nonhomogeneous radiation source ［J］. Applied Optics ， 2015 ， 54 （ 36 ）： 10599 - 10605 . doi: 10.1364/ao.54.010599 http://dx.doi.org/10.1364/ao.54.010599 Nonuniformity correction of imaging systems with a spatially nonhomogeneous radiation source 程万胜 ， 赵杰 ， 蔡鹤皋 . CCD像素响应非均匀的校正方法 ［J］. 光学 精密工程 ， 2008 ， 16 （ 2 ）： 314 - 318 . doi: 10.3321/j.issn:1004-924X.2008.02.021 http://dx.doi.org/10.3321/j.issn:1004-924X.2008.02.021 CCD像素响应非均匀的校正方法 CHENG W SH ， ZHAO J ， CAI H G . Correction method for pixel response nonuniformity of CCD ［J］. Opt. Precision Eng. ， 2008 ， 16 （ 2 ）： 314 - 318 . （in Chinese） . doi: 10.3321/j.issn:1004-924X.2008.02.021 http://dx.doi.org/10.3321/j.issn:1004-924X.2008.02.021 Correction method for pixel response nonuniformity of CCD JIN Y ， JIANG J ， ZHANG G J . Three-step nonuniformity correction for a highly dynamic intensified charge-coupled device star sensor ［J］. Optics Communications ， 2012 ， 285 （ 7 ）： 1753 - 1758 . doi: 10.1016/j.optcom.2011.12.043 http://dx.doi.org/10.1016/j.optcom.2011.12.043 Three-step nonuniformity correction for a highly dynamic intensified charge-coupled device star sensor ZHANG X ， ZENG F ， LI Y ， et al . Improvement in focusing accuracy of DNA sequencing microscope with multi-position laser differential confocal autofocus method ［J］. Optics Express ， 2018 ， 26 （ 2 ）： 887 - 896 . doi: 10.1364/oe.26.000887 http://dx.doi.org/10.1364/oe.26.000887 Improvement in focusing accuracy of DNA sequencing microscope with multi-position laser differential confocal autofocus method SILVESTRI L ， BRIA A ， SACCONI L ， et al . Confocal light sheet microscopy： micron-scale neuroanatomy of the entire mouse brain ［J］. Optics Express ， 2012 ， 20 （ 18 ）： 20582 . doi: 10.1364/oe.20.020582 http://dx.doi.org/10.1364/oe.20.020582 Confocal light sheet microscopy： micron-scale neuroanatomy of the entire mouse brain GONG H ， XU D ， YUAN J ， et al . High-throughput dual-colour precision imaging for brain-wide connectome with cytoarchitectonic landmarks at the cellular level ［J］. Nature Communications ， 2016 ， 7 ： 12142 . doi: 10.1038/ncomms12142 http://dx.doi.org/10.1038/ncomms12142 High-throughput dual-colour precision imaging for brain-wide connectome with cytoarchitectonic landmarks at the cellular level 张鑫 . 高通量基因测序荧光显微成像光学系统研究 ［D］. 长春 ： 中国科学院长春光学精密机械与物理研究所 ， 2018 . 高通量基因测序荧光显微成像光学系统研究 ZHANG X . Research on High-throughput Gene Sequencing Fluorescence Microscopic Imaging Optical System ［D］. Changchun： Changchun Institute of Optics， Fine Mechanics and Physics ， University of Chinese Academy of Sciences ， 2018 . （in Chinese） Research on High-throughput Gene Sequencing Fluorescence Microscopic Imaging Optical System CAO Y P ， TISSE C L . Single-image-based solution for optics temperature-dependent nonuniformity correction in an uncooled long-wave infrared camera ［J］. Optics Letters ， 2014 ， 39 （ 3 ）： 646 . doi: 10.1364/ol.39.000646 http://dx.doi.org/10.1364/ol.39.000646 Single-image-based solution for optics temperature-dependent nonuniformity correction in an uncooled long-wave infrared camera HE Z ， CAO Y ， DONG Y ， et al . Single-image-based nonuniformity correction of uncooled long-wave infrared detectors： a deep-learning approach ［J］. Applied Optics ， 2018 ， 57 （ 18 ）： D155 - D164 . doi: 10.1364/ao.57.00d155 http://dx.doi.org/10.1364/ao.57.00d155 Single-image-based nonuniformity correction of uncooled long-wave infrared detectors： a deep-learning approach ZHOU D B ， WANG D J ， HUO L J ， et al . Scene-based nonuniformity correction for airborne point target detection systems ［J］. Optics Express ， 2017 ， 25 （ 13 ）： 14210 . doi: 10.1364/oe.25.014210 http://dx.doi.org/10.1364/oe.25.014210 Scene-based nonuniformity correction for airborne point target detection systems MARCOTTE F ， TREMBLAY P ， FARLEY V . Infrared camera NUC calibration： comparison of advanced methods ［J］. SPIE ， 2014 ， 8706 ： 870603 . Infrared camera NUC calibration： comparison of advanced methods 王美昌 ， 于斌 ， 张炜 ， 等 . 基于数字微镜器件的数字线扫描荧光显微成像技术 ［J］. 物理学报 ， 2020 ， 69 （ 23 ）： 238701 . doi: 10.7498/aps.69.20200908 http://dx.doi.org/10.7498/aps.69.20200908 基于数字微镜器件的数字线扫描荧光显微成像技术 WANG M CH ， YU B ， ZHANG W ， et al . Digital line scanning fluorescence microscopy based on digital micromirror device ［J］. Acta Physica Sinica ， 2020 ， 69 （ 23 ）： 238701 . （in Chinese） . doi: 10.7498/aps.69.20200908 http://dx.doi.org/10.7498/aps.69.20200908 Digital line scanning fluorescence microscopy based on digital micromirror device 周豪 . 用于高通量筛选的细胞成像系统研究 ［D］. 杭州 ： 浙江大学 ， 2011 . 用于高通量筛选的细胞成像系统研究 ZHOU H . Research of Cell Image Capture System Used in High Throughput Screening ［D］. Hangzhou ： Zhejiang University ， 2011 . （in Chinese） Research of Cell Image Capture System Used in High Throughput Screening FAN Y ， NGUYEN D T ， AKAY Y ， et al . Engineering a brain cancer chip for high-throughput drug screening ［J］. Scientific Reports ， 2016 ， 6 ： 25062 . doi: 10.1038/srep25062 http://dx.doi.org/10.1038/srep25062 Engineering a brain cancer chip for high-throughput drug screening ZHANG Q ， YANG X ， HU Q ， et al . High axial resolution imaging system for large volume tissues using combination of inclined selective plane illumination and mechanical sectioning ［J］. Biomedical Optics Express ， 2017 ， 8 （ 12 ）： 5767 - 5775 . doi: 10.1364/boe.8.005767 http://dx.doi.org/10.1364/boe.8.005767 High axial resolution imaging system for large volume tissues using combination of inclined selective plane illumination and mechanical sectioning SUI X B ， CHEN Q ， GU G H . A novel non-uniformity evaluation metric of infrared imaging system ［J］. Infrared Physics & Technology ， 2013 ， 60 ： 155 - 160 . doi: 10.1016/j.infrared.2013.04.005 http://dx.doi.org/10.1016/j.infrared.2013.04.005 A novel non-uniformity evaluation metric of infrared imaging system 刘涛 . 红外图像非均匀性校正算法及图像质量评价的研究 ［D］. 杭州 ： 浙江大学 ， 2017 . doi: 10.20965/jdr.2017.p0198 http://dx.doi.org/10.20965/jdr.2017.p0198 红外图像非均匀性校正算法及图像质量评价的研究 LIU T . Non-uniformity Correction Algorithms and Quality Assessment for Infrared Image ［D］. Hangzhou ： Zhejiang University ， 2017 . （in Chinese） . doi: 10.20965/jdr.2017.p0198 http://dx.doi.org/10.20965/jdr.2017.p0198 Non-uniformity Correction Algorithms and Quality Assessment for Infrared Image 10.37188/OPE.20223005.0527 TH742 R318.6 A 1004-924X（2022）05-0527-09 常松涛,夏豪杰.时间延迟积分荧光显微成像平场校正技术[J].光学精密工程,2022,30(05):527-535. 常松涛,夏豪杰.时间延迟积分荧光显微成像平场校正技术[J].光学精密工程,2022,30(05):527-535. CHANG Songtao,XIA Haojie.Flat field correction technique for time delay integration fluorescence microscopy imaging[J].Optics and Precision Engineering,2022,30(05):527-535. CHANG Songtao,XIA Haojie.Flat field correction technique for time delay integration fluorescence microscopy imaging[J].Optics and Precision Engineering,2022,30(05):527-535. 国家自然科学基金面上项目（No.51975179） 2021-08-04 2020-09-22 2022-03-10 2022-03-25 光学精密工程 5 30