基于图像处理技术的薄层色谱定量分析方法技术

本发明专利技术属于数字图像处理和化学分析领域，其特征在于利用数码相机或摄像机在不同波长的光照条件下对薄层板进行拍摄；对所拍摄的图像进行镜头畸变校正、噪声虑除预处理；利用点样时点均匀分布的特点以及纵向距离等分原则确定控制点；根据控制点再利用插值方法构造出图像背景；将构造出的背景图像从镜头畸变校正后的薄层图像中去除，并根据所构造的背景图像对每一像素点进行亮度归一化；对去掉背景图像后的薄层图像进行分割；根据分割结果对每个区域内的灰度进行积分。本方法不仅能够准确的描述所给图像背景信息，分割出目标图像，而且还可以对每一点的光照强度进行修正，图像分割以及样品定量分析结果精度都很高，重复性也非常好，易于推广。

【技术实现步骤摘要】

本专利技术属于计算机数字图像处理领域和化学分析领域。
技术介绍
在化学分析、生物医学、医药研究、临床医学、食品加工等领域中，对某种微量物质成份的检测与鉴定是一项非常必要的工作，其中薄层色谱法或薄层层析法(Thin LayerChromatography，简称TLC)被广为使用。所谓薄层色谱法就是将吸附剂(例如硅胶或氧化铝)铺到平面物体上如玻璃板、塑料片或金属箔形成一薄层，然后用溶剂把样品展开和分离的方法，它是一种简便、经济、微量、适用范围广的色谱方法。薄层色谱法主要由三个步骤组成样品的点样、点样后样品的展开、展开后对样品的测量。样品的点样与展开已经成为非常成熟的技术，而对展开后的样品进行测定，特别是精确的定量测量技术却仍处于不断地探索当中。目前，薄层色谱定量分析基本上都采用了薄层扫描仪法。虽然薄层扫描仪被广泛地使用。但是，薄层扫描仪存在着以下缺点(a)测量速度慢，测量一块薄层板要花费几十分钟，不利于大量样品的检测；(b)精度较低，特别是对于面积较大、不规则的斑点，误差更大；(c)无法直接拍照或录像，降低了可信度；(d)不直观，因无法看到扫描的定位情况，易造成定位不准，从而产生较大误差；(e)操作复杂，由于薄层扫描仪的光学结构比较复杂，因此调整起来非常困难；(f)成本高，由于薄层扫描仪的光学结构比较复杂，因而制作起来工序多，因而成本高；最近几年，国际上有学者提出了利用平板式图像扫描仪法对薄层板上的样品进行定量分析测量。有人利用平板式图像扫描仪将染色后的薄层板进行扫描，得到数字图像后再进行处理。他们用分辨率为100DPI的扫描仪将薄层图像扫描到计算机中，然后利用计算机编程进行定量分析。他们的结论是，此方法不但可以进行定性分析，也可以进行定量分析。对平板式图像扫描仪而言，由于利用线型光源移动扫描，其扫描光线分布比较均匀，因而它对染色后的薄层板进行处理时得到的效果比较好。但是平板式图像扫描仪不能扫描需要特殊光源进行成像的薄层层析板，如荧光照射、紫外线照射等。此外，这种处理方法费时比较长，也无法与其他仪器集成一体。要解决平板式图像扫描仪所带来的问题，人们自然就会想到摄像机。其实，在平板式图像扫描仪被使用之前，利用摄像机成像对薄层色谱进行定量分析就被研究了，这是因为摄像机的成像速度快，安装方便，而且很容易和其它仪器集成使用。然而，在采用了摄像机进行薄层色谱分析后，人们给出的结论却不十分理想。原因是与传统的薄层扫描仪相比，摄像机成像的正确性受到光照均匀性的制约，因而结果不是十分精确，重复性不好。除了他们指出的原因外，摄像机成像时，光源不定、光圈不定、曝光速度不定也给研究带来了困难，也就是重复性不强。正是这些问题，使摄像机成像法在薄层色谱中的定量分析中得不到较好的应用。
技术实现思路
从上面的分析来看，如果能够将摄像机成像过程中光照不均的缺点进行相应的弥补，摄像机成像在薄层色谱中就会有强大的生命力成像速度快、分析精度高、结果直观、易于存储、成本低廉。摄像机成像过程中光照不均的现象是固有的，为此，怎样利用数字图像处理的手段弥补这一现象所带来的缺陷，这就是本专利技术的关键所在。另外，本专利技术的目的还在于提供了一种基于背景估计的一种图像分割方法。本专利技术的基本原理是当被测样品在薄层板上展开并干燥以后，就把该薄层板放置到具有多波长光照的灯箱中，然后利用数码相机或摄像机拍照形成数字图像，最后由计算机估测出薄层板上每一点的光照强度，根据估测出的光照强度把样品分割出来并根据样品灰度的积分来确定样品的实际含量。具体步骤包括利用摄像机在不同波长的光照条件下对薄层板进行拍摄；对所拍摄的图像进行镜头畸变校正、噪声虑除等预处理；利用纵向像素灰度积分投影法确定控制点；根据所得到的控制点，再利用插值方法构造出图像背景；将构造出的背景图像从镜头畸变校正后的薄层图像中去除，同时根据所构造的背景图像对每一像素点进行光照强度归一化；对去掉背景图像后的薄层图像进行分割；根据分割结果对每个样品区域内的灰度进行积分；根据样品灰度积分与样品质量成线性关系的特点，再对比标准品，计算出每个样品的实际含量。本专利技术特征在于，依次按下述步骤步骤1.对载有样品的薄层板利用数码相机或摄像机在所选的波长光照条件下进行拍摄，形成数字图像后存于计算机或DSP处理系统中；步骤2.计算机对所拍摄的图像按以下步骤进行预处理步骤2.1.镜头畸变校正用数码相机或摄像机拍摄标准方格图，再根据网格交点一一对应的关系，把步骤2.1得到的图像校正成标准方格图，记录下各点的校正参数，供以后拍摄薄层板图像时作校正用；步骤2.2.虑除噪声，采用下述滤波方法中的任何一种方法进行平滑滤波、或低通滤波、或维纳滤波，从而得到薄层板的数字图像；步骤3.利用点样时是沿着垂直于展开方向均匀分布的这一特点，从步骤2得到的预处理后图像中找到控制点的横坐标，再按照纵向距离等分的原则，确定控制点的纵坐标，从而得到控制点；步骤4.利用下述各插值方法中的任何一种方法依据步骤3得到的控制点来构造背景图像所述插值方法是拉格朗日插值、或均差的牛顿插值、或分段低次插值、或样条插值；步骤5.从步骤2得到的薄层板数字图像中减去所构造出的背景图像，设f(x，y)是镜头畸变校正后薄层板数字图像中像素(x，y)处的灰度，而g(x，y)是估测出的背景图像中像素(x，y)处的灰度，则相减得到图像的像素(x，y)处的灰度值为(g(x，y)-f(x，y))；再根据所构造的背景图像对每一点进行归一化，得到均匀背景的薄层板数字图像，该图像中对应点(x，y)处的灰度修正后为 或者 其中，Imax为所构造的背景图像的最高灰度，Imin为最小灰度，Ib为所述背景图像中同一点(x，y)处的灰度值，即每一点(x，y)的实际光强与背景光成反比；步骤6.利用全局阈值方法分割出目标图像，所述全局阈值取2或3； 步骤7.对步骤6得到的图像中各个区域内每一象素以灰度值进行累加，即进行灰度积分；步骤8.根据步骤7得到的灰度积分值的大小通过线性比例关系确定各样品点之间定量关系。具体的计算方法为设某一被测样品的灰度积分值为S，而标准品的积分值为T，且标准品的质量为G，则被测样品的质量为 采用本方法进行处理不仅能够准确的描述所给图像背景信息，分割出目标图像，而且还可以对每一点的光照强度进行归一化，分割结果达到很高的准确率，从而使定量分析的结果达到了较高的精度。下面结合附图对本专利技术作进一步说明。附图说明图1为本专利技术的硬件组成框图；图2为本专利技术的整体流程框图；图3是本专利技术的薄层板图像采集装置；图4是本专利技术的镜头畸变校正拍摄图；图5是本专利技术采集到的薄层板图像拍摄原图；图6是本专利技术所采集的薄层板图像的三维立体显示；图7是本专利技术的滤波模板；图8是本专利技术用于选择控制点的示意图，基于图5拍摄图；图9是利用本专利技术所构造图像背景的三维显示；图10是本专利技术得到的去除不均匀光照背景后的薄层板图像的三维显示结果；图11是本专利技术的样品分割结果。具体实施例方式图1给出了本专利技术的硬件组成框图。载有样品的薄层板1被放置到灯箱2中，然后由摄像机或数码相机3对其拍照，形成的图像被送到计算机或DSP处理系统4中，处理的结果被送到打印机6打印或其它设备7中备用。图1中的监视器5用以显示所拍摄的图像和程序运行过程。图2给出本文档来自技高网...

【技术保护点】
基于图像处理技术的薄层色普定量分析法，其特征在于：步骤１．对载有样品的薄层板利用数码相机或摄像机在所选的波长光照条件下进行拍摄，形成数字图像后存于计算机或ＤＳＰ处理系统中；步骤２．计算机对所拍摄的图像按以下步骤进行预处理：　 　步骤２．１．镜头畸变校正：用数码相机或摄像机拍摄标准方格图，再根据网格交点一一对应的关系，把步骤２．１得到的图像校正成标准方格图，记录下各点的校正参数，供以后拍摄薄层板图像时作校正用；步骤２．２．虑除噪声，采用下述滤波方法中的 任何一种方法进行：平滑滤波、或低通滤波、或维纳滤波，从而得到薄层板的数字图像；步骤３．利用点样时是沿着垂直于展开方向均匀分布的这一特点，从步骤２得到的预处理后图像中找到控制点的横坐标，再按照纵向距离等分的原则，确定控制点的纵坐标，从 而得到控制点；步骤４．利用下述各插值方法中的任何一种方法依据步骤３得到的控制点来构造背景图像：所述插值方法是拉格朗日插值、或均差的牛顿插值、或分段低次插值、或样条插值；步骤５．从步骤２得到的薄层板数字图像中减去所构造出的背景 图像，设ｆ（ｘ，ｙ）是镜头畸变校正后薄层板数字图像中像素（ｘ，ｙ）处的灰度，而ｇ（ｘ，ｙ）是估测出的背景图像中像素（ｘ，ｙ）处的灰度，则相减得到图像的像素（ｘ，ｙ）处的灰度值为（ｇ（ｘ，ｙ）－ｆ（ｘ，ｙ））；再根据所构造的背景图像对每一点进行归一化，得到均匀背景的薄层板数字图像，该图像中对应点（ｘ，ｙ）处的灰度修正后为Ｉ↓［ｍｉｎ］／Ｉ↓［ｂ］.Ｉ，或者Ｉ↓［ｍａｘ］／Ｉ↓［ｂ］.Ｉ，其中，Ｉ↓［ｍａｘ］为所构造的背景图像的最高灰度，Ｉ↓［ｍｉｎ］为最小灰度，Ｉ↓［ｂ］为所述背景图像中同一点（ｘ，ｙ）处的灰度值，即每一点（ｘ，ｙ）的实际光强与背景光成反比；步骤６．利用全局阈值方法分割出目标图像，所述全局阈值取２或３；步骤７．对步骤６得到的图像中各个区域内每一象素以灰度值进行累加，即进行灰度积分； 步骤８．根据步骤７得到的灰度积分值的大小通过线性比例关系确定各样品点之间定量关系。具体的计算方法为：设某一被测样品的灰度积分值为Ｓ，而标准品的积分值为Ｔ，且标准品的质量为Ｇ，则被测样品的质量为Ｓ／ＴＧ。...

【技术特征摘要】

【专利技术属性】
技术研发人员：张利，
申请(专利权)人：清华大学，
类型：发明
国别省市：11[中国|北京]