一种基于图像处理的微米级河鲀表皮体刺轮廓建模方法与流程

本发明属于生物建模技术，具体涉及一种基于图像处理的微米级河鲀表皮体刺轮廓建模方法。

背景技术：

在仿生学中，通常需要提取生物表面微小结构进行特征分析，以建立减阻元模型应用于飞机、船舶、泳衣等表面达到减阻的效果。以往因样件过小建模困难而采取形态轮廓简化，忽略了众多结构特征、模型精度较低。因此，对微米级结构的建模缺少一种操作简单、经济且能够高精度获取轮廓特征的建模方法。

技术实现要素：

发明目的：本发明的主要目的在于提供一种基于图像处理的微米级河鲀表皮体刺轮廓建模方法，以解决现有技术中建模精度低、效率低的问题。

技术方案：一种基于图像处理的微米级河鲀表皮体刺轮廓建模方法，包括如下步骤：

(1)通过显微镜获取轮廓位于同一焦平面的河鲀表皮体刺显微图像；

(2)通过halcon软件获取显微图像中轮廓条线及其轮廓图像；

(3)将轮廓图像导入photoshop中进行去噪处理，并导出ai文件；再通过3dmax软件将ai文件导出为cad文件；

(4)在cad软件中处理转化为cad文件的曲线轮廓；

(5)使用matlab导入所有点的坐标，以x为自变量，y为因变量进行曲线拟合后建立河鲀表皮体刺模型。

进一步的，步骤(1)中样本置于液体，通过液面静置的方式使样本轮廓均位于同一焦平面内，待测的样本中心轴或对称轴为水平线或垂直于水平线。

步骤(2)具体步骤如下：

(21)将原图像转换为rgb单通道图像，选择对比度最高的分量图像进行二值化处理，进而获得完整的图像信息及清晰的目标轮廓线条；

(22)将二值图像中的空洞区域进行填充，并采用腐蚀运算消除边缘毛体刺与干扰噪声；

(23)边缘连接存在小的缝隙，通过膨胀运算实现轮廓的平滑过渡；

(24)通过图像增强中的梯度锐化算子提高图像信噪比，提取出微小结构的真实边缘，获取到完整的河鲀表皮体刺轮廓线。

由于体刺尺寸微小，在显微拍摄过程中无法控制其位置，需对体刺进行坐标系重建与统一；为保证数学模型的精确性，样条点的提取应避免人工操作带来的误差，使用cad软件能够满足要求且操作简便、功能稳定。因此将轮廓图像通过photoshop软件的工作路径功能导出ai文件，经过3dmax软件转换为cad文件。

步骤(4)具体包括如下步骤：

(41)对河鲀表皮体刺轮廓进行坐标系的重置：将表皮体刺两边切线交点为坐标原点，切线所成角的角平分线为y轴，其过原点的垂直方向为x轴，通过坐标变换实现坐标系的统一；

(42)通过cad软件的样条曲线功能对样件轮廓密集取点；

(43)使用list功能批量获取样条曲线上所有点的坐标。

通过所获得的样条点坐标在matlab中进行拟合，得到体刺的轮廓曲线数学模型和构建傅里叶函数表达式。

进一步的，fourier函数表示，阶数取1，表达式如下：

f(x)＝-1.503+1.04×cos(11.28×x)-0.003714×sin(11.28×x)

其数学模型的精度参数为：残差平方和为0.1336，均方差为0.02974，确定系数为0.9954，与生物模型拟合精度高，表达式复杂程度低。

有益效果：与现有技术相比，本发明通过图像获取设备及图像处理技术以获得更高精度和更加贴近生物轮廓的图像数据特征，然后结合制图软件和算法对样本图像做优化处理以获得高精度的样本模型，并且所述方法通过计算机软件即可实现微小样件轮廓的数学模型获取，操作简单、经济，同时可根据结果进行局部特征参数的提取，提高了建模效率及质量，解决了现有技术中微小尺寸样件建模精度较低的问题。

附图说明

图1为本发明所述方法的实施流程图；

图2为本发明所述方法中样本的放置方式示意图；

图3为实施例中样本为河鲀表皮体体刺的提取图像；

图4为实施例中由显微图像精确提取出的河鲀锥状部位轮廓图像；

图5为在cad中对河鲀轮廓进行坐标系重建及与系统坐标系统一处理示意图；

图6为在matlab中对河鲀轮廓样条点进行拟合及函数的选取；

图7为由轮廓数学表达式建立的几何模型函数图像；

图8实施例中为河鲀的建模示意图。

具体实施方式

为了详细地说明本发明所公开的技术方案，下面结合说明书附图及具体实施例做进一步的阐述。

在工业制造方面，仿生学是研究应用的一个重要的技术。本发明所提供的就是一种基于图像处理的微米级河鲀表皮体刺轮廓建模方法。

首先是样本提取过程：取成年健康菊黄东方鲀(河鲀)为实验对象，剥取完整鱼皮浸泡一夜至黏液与表皮分离，清洗干净，使用生物解剖工具对河鲀表皮进行解剖并标记。配置5％浓度氢氧化钾溶液浸泡表皮，置于kh-45a恒温箱内6-8小时加速鱼皮溶解。鱼皮溶解后，倒出浑浊的氢氧化钾溶液，避免长时间浸泡对体体刺造成损伤。滴入清水反复清洗，获得干净的乳白色河鲀表皮体体刺。配置1％浓度氢氧化钾溶液及1％浓度茜红素溶液浸没体体刺，静置2-3分钟进行体体刺染色，获得体体刺试验样本，用于显微镜观察。

对于河鲀表皮体刺的轮廓建模包括如下过程：

(1)通过显微镜及软件cellsensstandard获取样件的显微图像，其样件观测位置及方式根据所需测量轮廓决定。其中，河鲀表皮体刺所需测量轮廓为锥状部位曲线，使用液面静置方式置于显微镜下观测；体刺尺寸微小，根部呈爪状，无法使用支架等方式对体刺锥状部位进行水平放置。可以采用增加液体密度的方式以支撑河鲀体刺，保持样本短时间静止状态。步骤(1)中样本轻置于平静水面保持水平状态，使体刺轮廓在显微镜中位于同一焦平面，获得清晰完整准确的体刺锥状部位图像。

(2)使用halcon软件和相应的算法精确获取显微图像中所需轮廓样条线；

2.1、将原图像转换为rgb单通道图像，选择对比度最高的分量图像进行二值化处理，进而获得完整的图像信息及清晰的目标轮廓线条；

2.2、将二值图像中的空洞区域进行填充，并采用腐蚀运算消除边缘毛体刺与干扰噪声；

2.3、边缘连接存在小的缝隙，可通过膨胀运算实现轮廓的平滑过渡；

2.4、使用图像增强中的梯度锐化算子提高图像信噪比，使图像轮廓更为清晰，准确地提取出微小结构的真实边缘

(3)将轮廓图像导入photoshop中进行去噪处理，并导出ai文件；再通过3dmax软件将ai文件导出为cad文件；

(4)在cad软件中处理转化为cad文件的曲线轮廓；

4.1针对样件进行坐标系的重置；

4.2将新建立的坐标系与cad系统坐标系重合；

4.3使用样条曲线功能对样件轮廓密集取点，曲率变化较小处取点可较为稀疏，曲率变化越大处取点约密集，以保证样条曲线与样件轮廓线高度重合；

4.4使用list功能批量获取样条曲线上所有点的坐标；

(5)使用matlab导入所有点的坐标，以x为自变量，y为因变量进行曲线拟合。选取不同的函数及阶数拟合轮廓曲线，根据残差平方和(sse)、均方差(rmse)值及确定系数(r-square)值判断拟合精度，根据需求获得合适的数学模型；

当微小结构为旋成体，根据所得轮廓数学模型可继而在三维软件nx10.0中建立几何模型；然后根据所得轮廓数学模型，在matlab中可继而获得轮廓曲线的局部曲率、坐标等曲线特征数值。

详细地说，本发明所述的方法具体建模过程如图1所示。对于河鲀表皮体刺这一类比较小的生物轮廓，需要获取其轮廓的函数表达式，则提取时确保体刺锥状部分是完整的，由于体刺尺寸微小，难以实现拆分，因此该专利提供的提取方式不需要对样件进行拆分并能够保证体刺锥状部分的完整性。

通过显微镜及软件cellsensstandard获取样件的显微图像，其样件观测位置及方式根据所需测量轮廓决定。其中，河鲀表皮体刺所需测量轮廓为锥状部位曲线，使用液面静置方式置于显微镜下观测，如图2-图8所示。然后使用halcon软件编写算法精确获取显微图像中所需轮廓样条线；接着将轮廓图像导入photoshop中进行去噪处理，并导出ai文件；再通过3dmax软件将ai文件导出为cad文件。由于体刺尺寸微小，在显微拍摄过程中无法控制其位置，需对体刺进行坐标系重建与统一；为保证数学模型的精确性，样条点的提取应避免人工操作带来的误差，使用cad软件能够满足要求且操作简便、功能稳定。因此将轮廓图像通过photoshop软件的工作路径功能导出ai文件，经过3dmax软件转换为cad文件。对体刺轮廓进行坐标系的重置：体刺两边切线交点为坐标原点，切线所成角的角平分线为y轴，其过原点的垂直方向为x轴，通过坐标变换实现坐标系的统一。

本实施例中构建的河鲀轮廓线表达式具体如下：

f(x)＝a0+a1×cos(x×ω)+b1×sin(x×ω)

表1.河鲀表皮体体刺轮廓曲线的方程系数

即：f(x)＝-1.503+1.04×cos(11.28×x)-0.003714×sin(11.28×x)。

本发明所提供的方法仅需通过计算机软件即可实现微小样件轮廓的数学模型获取，操作简单、经济，同时可根据结果进行局部特征参数的提取，提高了建模效率及质量，解决了现有技术中微小尺寸样件建模精度较低的问题。