一种拟南芥轮廓表型的检测方法.pdf

(19)中华人民共和国国家知识产权局*CN103471523A*(12)发明专利申请(10)申请公布号(10)申请公布号CNCN103471523103471523A(43)申请公布日2013.12.25(21)申请号201310456227.2(22)申请日2013.09.30(71)申请人北京林业大学地址100083北京市海淀区清华东路35号(72)发明人田野张俊梅柯秋红聂凤梅(51)Int.Cl.G01B11/24(2006.01)权权利要求书2页利要求书2页说明书10页说明书10页附图3页附图3页(54)发明名称一种拟南芥轮廓表型的检测方法(57)摘要一种拟南芥轮廓表型的检测方法，该方法具体包括：在拟南芥的种植盆中，放置标定板，利用相机采集拟南芥的RGB图像；对采集后的图像进行预处理，实现图像的自动校正和标定，其中图像校正是为了校正图像的畸变，图像标定是为了获得单位像素的真实尺寸；对预处理后的图像进行分割，将拟南芥与背景分割，从图像中提取出来；分割出拟南芥图像之后，提取拟南芥的轮廓表型参数，用椭圆傅里叶描述子定量地描述其整体轮廓。通过对整体轮廓表型参数的分析，可以描述不同基因的拟南芥在整体轮廓形态和生长取向上面的差异，从而可以推断出不同基因的功能及对拟南芥植物的影响。CN103471523ACN1034752ACN103471523A权利要求书1/2页1.一种拟南芥轮廓表型的检测方法，包括以下步骤：S1.在拟南芥的种植盆中，放置标定板，利用相机采集拟南芥的RGB图像；S2.对采集后的图像进行预处理，实现图像的自动校正和标定，其中图像校正是为了校正图像的畸变，图像标定是为了获得单位像素的真实尺寸；步骤S2包括定位黑白棋盘格、角点检测、图像校正和图像标定；所述定位黑白棋盘格包括以下步骤：根据蓝色的RGB特征，提取蓝色边框；对图像的内部进行孔洞填充，再减去原来的蓝色边框图像，得到新图像；对新图像进行开运算，去除噪点；再进行闭运算，连接断点，得到黑白棋盘格的区域；所述角点检测包括以下步骤：计算黑白棋盘格区域中的每一个点的水平方向和垂直方向上的一阶导数，得到三幅新图像：水平一阶导数的平方，垂直一阶导数的平方和两个一阶导数的乘积；用高斯滤波对三幅图像进行滤波，去除噪声；由上述三幅图像组成相关矩阵，计算准则函数，判断其中的像素点是否为角点；所述图像校正和图像标定包括以下步骤：通过角点检测得到棋盘格的每个正方形的顶点在图像中坐标，并根据其在真实世界中的空间位置关系，得到二者的变换矩阵；求取变换矩阵的逆，作用于图像，实现图像校正；通过角点坐标得到黑白棋盘格的总像素个数，并根据它的真实尺寸，得到单位像素的真实尺寸；S3.对预处理后的图像进行分割，将拟南芥与背景分割，从图像中提取出来；步骤S3包括图像初分割和去除噪点；所述图像初分割包括以下步骤：对每个像素的RGB值进行归一化获得rgb，提取3g-2.4r-b的色差图，用0作为阈值对图像进行二值化；将结果减去原来检测出的蓝色边框区域，判断所获前景区域像素的G的灰度值是否大于50，若是，则保留，否则去除；提取前景连通区域中的具有最多像素的区域，即为植物区域；所述去除噪点包括以下步骤：对所得图像进行开运算，得到新的图像，其将只保留大叶片区域和植物的中心区域，而去除图像中的细节部分，包括植物的茎和叶片周围的噪点；对去除的细节部分的各个连通区域计算其像素个数，像素个数小于12的直接去除，像素个数大于等于12的连通区域，分别单独与开运算得到的新图像叠加，再计算叠加后图像中的区域个数，如果区域个数减少，说明该连通区域是茎，必须保留，否则，如果叠加后图像中的区域个数增加或者不变，则说明该连通区域是叶片周围的噪点，必须去除；判断完所有的细节部分后，将保留的所有连通区域全部与开运算得到的新图像叠加，得到最终的植物区域；S4.分割出拟南芥图像之后，提取拟南芥的轮廓表型参数，用椭圆傅里叶描述子定量地描述其整体轮廓；其中步骤S4包括如下步骤：S4.1对植物区域的孔洞进行填充；S4.2通过轮廓查找算法，提取填充后的植物区域的外轮廓点；S4.3将轮廓点按照x和y坐标分别存储，从而将2维(2D)图像变成2个1维(1D)信号；S4.4对S4.3中所述的两个1D信号分别进行数据点镜像填充，使其扩展到2的整数次幂倍，得到两个新的1D信号；S4.5对S4.4得到两个1D信号分别进行4层Haar小波分解，去除噪声，取第4层的近2CN103471523A权利要求书2/2页似层系数作为新的信号进行处理；S4.6对S4.5得到的两个新的1D信号进行傅里叶变换，获得10层的傅里叶系数，每层系数包含4个描述子；S4.7对每层的四个系数求归一化的模数，从而用10个傅里叶描述子描述植物的整体轮廓。2.根据权利要求1所述的拟南芥轮廓表型的检测方法