CN111579511A - 基于结构高光谱的种子质量检测方法及装置 - Google Patents

基于结构高光谱的种子质量检测方法及装置 Download PDF

Info

Publication number
CN111579511A
CN111579511A CN202010545125.8A CN202010545125A CN111579511A CN 111579511 A CN111579511 A CN 111579511A CN 202010545125 A CN202010545125 A CN 202010545125A CN 111579511 A CN111579511 A CN 111579511A
Authority
CN
China
Prior art keywords
seed
structured light
mechanical arm
seeds
information
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.)
Granted
Application number
CN202010545125.8A
Other languages
English (en)
Other versions
CN111579511B (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 Agricultural University
Original Assignee
Nanjing Agricultural University
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 Agricultural University filed Critical Nanjing Agricultural University
Priority to CN202010545125.8A priority Critical patent/CN111579511B/zh
Publication of CN111579511A publication Critical patent/CN111579511A/zh
Application granted granted Critical
Publication of CN111579511B publication Critical patent/CN111579511B/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/25—Colour; Spectral properties, i.e. comparison of effect of material on the light at two or more different wavelengths or wavelength bands
    • G01N21/31—Investigating relative effect of material at wavelengths characteristic of specific elements or molecules, e.g. atomic absorption spectrometry
    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01B—MEASURING LENGTH, THICKNESS OR SIMILAR LINEAR DIMENSIONS; MEASURING ANGLES; MEASURING AREAS; MEASURING IRREGULARITIES OF SURFACES OR CONTOURS
    • G01B11/00—Measuring arrangements characterised by the use of optical techniques
    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01B—MEASURING LENGTH, THICKNESS OR SIMILAR LINEAR DIMENSIONS; MEASURING ANGLES; MEASURING AREAS; MEASURING IRREGULARITIES OF SURFACES OR CONTOURS
    • G01B11/00—Measuring arrangements characterised by the use of optical techniques
    • G01B11/02—Measuring arrangements characterised by the use of optical techniques for measuring length, width or thickness
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T7/00—Image analysis
    • G06T7/80—Analysis of captured images to determine intrinsic or extrinsic camera parameters, i.e. camera calibration
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T2207/00—Indexing scheme for image analysis or image enhancement
    • G06T2207/30—Subject of image; Context of image processing
    • G06T2207/30204—Marker
    • G06T2207/30208—Marker matrix

Landscapes

  • Physics & Mathematics (AREA)
  • General Physics & Mathematics (AREA)
  • Spectroscopy & Molecular Physics (AREA)
  • Engineering & Computer Science (AREA)
  • Health & Medical Sciences (AREA)
  • Theoretical Computer Science (AREA)
  • Computer Vision & Pattern Recognition (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • Chemical & Material Sciences (AREA)
  • Analytical Chemistry (AREA)
  • Biochemistry (AREA)
  • General Health & Medical Sciences (AREA)
  • Immunology (AREA)
  • Pathology (AREA)
  • Image Processing (AREA)

Abstract

本发明提出一种基于结构高光谱的种子质量检测装置，由机械臂基座、机械臂、真空吸种模块、结构光投影模块、白光光源、光纤、高光谱相机、计算机、信号线、种子组成；一种基于结构高光谱的种子质量检测方法，包括以下步骤：基于结构高光谱的种子质量检测装置标定、种子三维表型信息获取、种子特征光谱信息和二维纹理信息的获取、基于支持向量回归(SVR)的种子活力预测模型建模；本发明采用结构光三维成像技术，可精确测量种子的长度、宽度、体积等三维尺寸信息，通过结构光调制技术可检测种子不同深度的高光谱信息，拓展信息维度，提高种子活力指数的检测精度。

Description

基于结构高光谱的种子质量检测方法及装置
技术领域
本发明属于种子质量自动化检测的技术领域，尤其是一种基于结构高光谱的种子质量检测方法及装置，属于农业工程领域。
背景技术
农作物产量是是保证粮食安全的重要指标，是国家安全的基石。种子质量的优劣很大程度上决定了粮食的产量，进行种子表型检测是增加农作物产量的重要因素，进行种子质量检测具有重要意义。衡量种子质量的标准主要包括发芽率、种子净度、纯度、活力、损伤裂纹等。其中，种子物理尺寸、活力指数等直接关系到作物的产量和质量。
传统的种子物理尺寸表型检测通过提取相机拍照后的二维图像进行种子直径、长度的测量，但无法实现种子体积的精确测量。此外，传统采用发芽法、四唑法等进行种子活力检测为有损检测，而基于近红外、高光谱的检测方法因无法提取不同深度的种子光谱信息，检测精度难以提升。
发明内容
针对现有种子表型参数和活力指数检测精度低的不足，本发明提出一种基于结构高光谱的种子质量检测方法及装置。
为实现上述目的，本发明采用如下的技术方案：
一种基于结构高光谱的种子质量检测装置，其特征为，由机械臂基座(1)、机械臂(2)、真空吸种模块(3)、结构光投影模块(4)、白光光源(5)、光纤(6)、高光谱相机(7)、计算机(8)、信号线(9)、种子(10)组成，其中，机械臂(2)由机械臂升降台(21)、机械臂大臂(22)、机械臂小臂(23)、机械臂连接法兰(24)、手臂连接轴(25)组成，且机械臂升降台(21)的底部固装于机械臂基座(1)上，机械臂升降台(21)的滑动端与机械臂大臂(22)一端固装，机械臂大臂(22)另一端与机械臂小臂(23)一端通过手臂连接轴(25)滑动连接，机械臂小臂(23)另一端通过机械臂连接法兰(24)与真空吸种模块(3)的法兰连接件(31)固装；白光光源(5)通过光纤(6)为结构光投影模块(4)提供光源，结构光投影模块(4)投射结构光到种子(10)上，高光谱相机(7)采集种子(10)上投射的结构光图像并通过信号线(9)传输到计算机(8)。
一种基于结构高光谱的种子质量检测装置，所述的真空吸种模块(3)，其特征为：由法兰连接件(31)、真空管转接头(32)、气吸空心针(33)、真空发生器(34)和软气管(35)组成，其中，真空发生器(34)输出端和软气管(35)一端相连，软气管(35)另一端通过真空管转接头(32)连接至法兰连接件(31)一端，法兰连接件(31)另一端与气吸空心针(33)一端相连，气吸空心针(33)另一端可吸附种子(10)，且法兰连接件(31)，气吸空心针(33)，真空管转接头(32)，软气管(35)均为真空结构，并共同构成内部空心连通体，可吸附种子(10)。
一种基于结构高光谱的种子质量检测方法，其特征为，包括以下步骤：
步骤1：基于结构高光谱的种子质量检测装置标定；
步骤2：种子三维表型信息获取；
步骤3：种子特征光谱信息和二维纹理信息的获取；
步骤4：基于支持向量回归(SVR)的种子活力预测模型建模。
所述的一种基于结构高光谱的种子质量检测方法，其特征为：所述步骤1基于结构高光谱的种子质量检测装置标定按照以下步骤标定：
步骤a：将标准棋盘格标定板垂直置于气吸空心针下方，标定板垂直轴线与气吸空心针同轴，标定板平面中心与结构光投影模块出光口中心的连线与标定板平面成90°，通过电脑8生成条纹形格雷码；
步骤b：电脑8通过信号线9将一组条纹形格雷码顺次发送到结构光投影模块，结构光投影模块顺次将编码结构光投射到标准棋盘格标定板上，电脑8再通过信号线9和高光谱相机7顺次采集标准棋盘格标定板上的结构光图像信息，并顺次存储保存为以***数字命名的jpg格式文件；
步骤c：电脑8将步骤b中存储的一组以***数字命名的jpg图片格式文件保存于angle1文件夹中，再分别将标准棋盘格标定板旋转2种不同角度，并根据步骤b的方法分别把2组保存为以***数字命名的jpg格式文件分别存储在angle2和angle3文件夹中；
步骤d：采用张正友棋盘格标定法，将步骤c所述angle1、angle2、angel3文件夹导入到MATALB calibration toolbox中，得到高光谱相机和投影仪的内参矩阵KL、KR，以及投影仪坐标系转换到高光谱相机坐标系的旋转R矩阵和平移T矩阵，其中内参矩阵包含每个像素的物理尺寸、焦距、图像物理坐标的扭曲因子、图像原点相对于光心成像点的纵横偏移量；
步骤e：根据步骤d得到的高光谱相机和投影仪的内参矩阵KL，KR，以及投影仪坐标系转换到高光谱相机坐标系的旋转R矩阵和平移T矩阵进行格雷码解码，得到三幅点云图像，分别命名为point1，point2point3；
步骤f：在point1点云图像中任取一个点作为参考点PR(x，y，z)，其在point2、point3中对应的点为PR1(x，y，z)，PR2(x，y，z)；
步骤g：基于3个点PR(x，y，z)、PR1(x，y，z)、PR2(x，y，z)的坐标信息，拟合出圆并计算出圆心坐标作为种子旋转坐标系的坐标原点P0(x，y，z)；
步骤h：设P0(x，y，z)为坐标原点，PR(x，y，z)为机械臂旋转台工作表面上的参考点，PR1(x，y，z)和PR2(x，y，z)为其在两个不同视场中所对应的坐标；
步骤i：通过
Figure BSA0000211498140000021
两个向量的向量积可求出Z′轴的方向向量，其中，
Figure BSA0000211498140000022
Figure BSA0000211498140000023
其中，上述步骤e中格雷码解码采用以下方法：
设三维空间中的物体在相机和投影仪的成像位置分别为Pc(uc，vc)和Pp(up，vp)，且设世界坐标系的原点与相机坐标系的原点重合，则相机的外参矩阵(相机和世界坐标系的转换矩阵)为
Figure BSA0000211498140000024
投影仪的外参矩阵(投影仪和世界坐标系的转换矩阵)，
Rp＝R，Tp＝T (2)
利用相机模型可得
Figure BSA0000211498140000025
Figure BSA0000211498140000031
联立上述方程组，获得Sc，Sp，X，Y，Z，和方程
Figure BSA0000211498140000032
Figure BSA0000211498140000033
则该结构光3D成像***每个像素点的三维解为
Figure BSA0000211498140000034
其中
Figure BSA0000211498140000035
和
Figure BSA0000211498140000036
表示Mc和Mp在第i行第j列。
所述的一种基于结构高光谱的种子质量检测方法，其特征为：所述步骤2中种子三维表型信息获取流程为：
步骤a：电脑8通过信号线9将一组结构光编码信息发送到结构光投影模块，结构光投影模块将编码结构光投射到种子上，电脑8再通过信号9和高光谱相机采集种子的图形信息，每张图片存储成jpg格式，然后进行格雷码解码得到第一视角的点云图像；
步骤b：分别将种子旋转90°、180°、270°，同时根据步骤a的方法得到三个视角的点云图；
步骤c：将四个视角的三维点云图根据上述标定的种子旋转轴Z′变换到统一坐标系，完成点云的粗配准，调用MATLAB的ICP算法和平滑算法实现点云精配准；
步骤d：封装配准后的点云，计算长、宽、高、体积、外表面面积作为种子的三维表型参数。
所述的一种基于结构高光谱的种子质量检测方法，其特征为：所述步骤3种子特征光谱信息和二维纹理信息的获取，按照以下步骤：
步骤a：结构光投影模块向种子投射白光，高光谱相机获取光谱图像立方体；
步骤b：采用MATLAB中OTSU方法进行阈值分割，获取种子的ROI掩模图像，通过ROI掩模图像提取高光谱图像中的种子ROI光谱信息；
步骤c：对光谱调用MATLAB中smooth平滑滤波函数进行平滑滤波；
步骤d：调用MATLAB中PCA算法提取高光谱的前三个特征波长；
步骤e：提取上述三个特征波长对应的灰度图像；
步骤f：用MATLAB调用graycomatrix函数计算上述三个特征波长对应灰度图像的生成灰度共生矩阵，其中包含种子的二维纹理信息。
所述的基于支持向量回归(SVR)的种子质量预测模型建模，按照以下方法实施，调用MATLAB中Lib-SVM工具箱中的SVR算法，以上述前三个特征波长、灰度共生矩阵作为输入，以种子活力指数的数值为输出，建立种子活力预测模型，可用于种子活力预测。
与现有技术相比，本发明的有益效果为：
(1)本发明采用结构光三维成像技术，可精确测量种子的长度、宽度、体积等三维尺寸信息。
(2)本发明通过结构光调制技术可检测种子不同深度的高光谱信息，拓展信息维度，提高种子活力指数的检测精度。
附图说明
图1为基于结构高光谱种子质量检测装置结构示意图；
图2为基于结构高光谱种子质量检测方法流程图；
具体实施方式
下面结合附图及具体实施例对本发明作进一步说明。
实施例1
如图1所示，一种基于结构高光谱的种子质量检测装置，其特征为，由机械臂基座(1)、机械臂(2)、真空吸种模块(3)、结构光投影模块(4)、白光光源(5)、光纤(6)、高光谱相机(7)、计算机(8)、信号线(9)、种子(10)组成，其中，机械臂(2)由机械臂升降台(21)、机械臂大臂(22)、机械臂小臂(23)、机械臂连接法兰(24)、手臂连接轴(25)组成，且机械臂升降台(21)的底部固装于机械臂基座(1)上，机械臂升降台(21)的滑动端与机械臂大臂(22)一端固装，机械臂大臂(22)另一端与机械臂小臂(23)一端通过手臂连接轴(25)滑动连接，机械臂小臂(23)另一端通过机械臂连接法兰(24)与真空吸种模块(3)的法兰连接件(31)固装；白光光源(5)通过光纤(6)为结构光投影模块(4)提供光源，结构光投影模块(4)投射结构光到种子(10)上，高光谱相机(7)采集种子(10)上投射的结构光图像并通过信号线(9)传输到计算机(8)。
实施例2
一种基于结构高光谱的种子质量检测装置，所述的真空吸种模块(3)，其特征为：由法兰连接件(31)、真空管转接头(32)、气吸空心针(33)、真空发生器(34)和软气管(35)组成，其中，真空发生器(34)输出端和软气管(35)一端相连，软气管(35)另一端通过真空管转接头(32)连接至法兰连接件(31)一端，法兰连接件(31)另一端与气吸空心针(33)一端相连，气吸空心针(33)另一端可吸附种子(10)，且法兰连接件(31)，气吸空心针(33)，真空管转接头(32)，软气管(35)均为真空结构，并共同构成内部空心连通体，可吸附种子(10)。
实施例3
如图2所示，一种基于结构高光谱的种子质量检测方法，其特征为，包括以下步骤：
步骤1：基于结构高光谱的种子质量检测装置标定；
步骤2：种子三维表型信息获取；
步骤3：种子特征光谱信息和二维纹理信息的获取；
步骤4：基于支持向量回归(SVR)的种子活力预测模型建模。
实施例4
权力要求2所述的一种基于结构高光谱的种子质量检测方法，其特征为：所述步骤1基于结构高光谱的种子质量检测装置标定按照以下步骤标定：
步骤a：将标准棋盘格标定板垂直置于气吸空心针下方，标定板垂直轴线与气吸空心针同轴，标定板平面中心与结构光投影模块出光口中心的连线与标定板平面成90°，通过电脑8生成条纹形格雷码；
步骤b：电脑8通过信号线9将一组条纹形格雷码顺次发送到结构光投影模块，结构光投影模块顺次将编码结构光投射到标准棋盘格标定板上，电脑8再通过信号线9和高光谱相机7顺次采集标准棋盘格标定板上的结构光图像信息，并顺次存储保存为以***数字命名的jpg格式文件；
步骤c：电脑8将步骤b中存储的一组以***数字命名的jpg图片格式文件保存于angle1文件夹中，再分别将标准棋盘格标定板旋转2种不同角度，并根据步骤b的方法分别把2组保存为以***数字命名的jpg格式文件分别存储在angle2和angle3文件夹中；
步骤d：采用张正友棋盘格标定法，将步骤c所述angle1、angle2、angel3文件夹导入到MATALB calibration toolbox中，得到高光谱相机和投影仪的内参矩阵KL、KR，以及投影仪坐标系转换到高光谱相机坐标系的旋转R矩阵和平移T矩阵，其中内参矩阵包含每个像素的物理尺寸、焦距、图像物理坐标的扭曲因子、图像原点相对于光心成像点的纵横偏移量；
步骤e：根据步骤d得到的高光谱相机和投影仪的内参矩阵KL，KR，以及投影仪坐标系转换到高光谱相机坐标系的旋转R矩阵和平移T矩阵进行格雷码解码，得到三幅点云图像，分别命名为point1，point2point3；
步骤f：在point1点云图像中任取一个点作为参考点PR(x，y，z)，其在point2、point3中对应的点为PR1(x，y，z)，PR2(x，y，z)；
步骤g：基于3个点PR(x，y，z)、PR1(x，y，z)、PR2(x，y，z)的坐标信息，拟合出圆并计算出圆心坐标作为种子旋转坐标系的坐标原点P0(x，y，z)；
步骤h：设P0(x，y，z)为坐标原点，PR(x，y，z)为机械臂旋转台工作表面上的参考点，PR1(x，y，z)和PR2(x，y，z)为其在两个不同视场中所对应的坐标；
步骤i：通过
Figure BSA0000211498140000051
两个向量的向量积可求出Z′轴的方向向量，其中，
Figure BSA0000211498140000052
Figure BSA0000211498140000053
其中，上述步骤e中格雷码解码采用以下方法：
设三维空间中的物体在相机和投影仪的成像位置分别为Pc(uc，vc)和Pp(up，vp)，且设世界坐标系的原点与相机坐标系的原点重合，则相机的外参矩阵(相机和世界坐标系的转换矩阵)为
Figure BSA0000211498140000054
投影仪的外参矩阵(投影仪和世界坐标系的转换矩阵)，
Rp＝R，Tp＝T (2)
利用相机模型可得
Figure BSA0000211498140000055
Figure BSA0000211498140000056
联立上述方程组，获得Sc，Sp，X，Y，Z，和方程
Figure BSA0000211498140000061
Figure BSA0000211498140000062
则该结构光3D成像***每个像素点的三维解为
Figure BSA0000211498140000063
其中
Figure BSA0000211498140000064
和
Figure BSA0000211498140000065
表示Mc和Mp在第i行第j列。
实施例5
所述的一种基于结构高光谱的种子质量检测方法，其特征为：所述步骤2中种子三维表型信息获取流程为：
步骤a：电脑8通过信号线9将一组结构光编码信息发送到结构光投影模块，结构光投影模块将编码结构光投射到种子上，电脑8再通过信号9和高光谱相机采集种子的图形信息，每张图片存储成jpg格式，然后进行格雷码解码得到第一视角的点云图像；
步骤b：分别将种子旋转90°、180°、270°，同时根据步骤a的方法得到三个视角的点云图；
步骤c：将四个视角的三维点云图根据上述标定的种子旋转轴Z′变换到统一坐标系，完成点云的粗配准，调用MATLAB的ICP算法和平滑算法实现点云精配准；
步骤d：封装配准后的点云，计算长、宽、高、体积、外表面面积作为种子的三维表型参数。
实施例6
所述的一种基于结构高光谱的种子质量检测方法，其特征为：所述步骤3种子特征光谱信息和二维纹理信息的获取，按照以下步骤：
步骤a：结构光投影模块向种子投射白光，高光谱相机获取光谱图像立方体；
步骤b：采用MATLAB中OTSU方法进行阈值分割，获取种子的ROI掩模图像，通过ROI掩模图像提取高光谱图像中的种子ROI光谱信息；
步骤c：对光谱调用MATLAB中smooth平滑滤波函数进行平滑滤波；
步骤d：调用MATLAB中PCA算法提取高光谱的前三个特征波长；
步骤e：提取上述三个特征波长对应的灰度图像；
步骤f：用MATLAB调用graycomatrix函数计算上述三个特征波长对应灰度图像的生成灰度共生矩阵，其中包含种子的二维纹理信息。
实施例7
所述的基于支持向量回归(SVR)的种子质量预测模型建模，按照以下方法实施，调用MATLAB中Lib-SVM工具箱中的SVR算法，以上述前三个特征波长、灰度共生矩阵作为输入，以种子活力指数的数值为输出，建立种子活力预测模型，可用于种子活力预测。

Claims (7)

1.一种基于结构高光谱的种子质量检测装置，其特征为，由机械臂基座(1)、机械臂(2)、真空吸种模块(3)、结构光投影模块(4)、白光光源(5)、光纤(6)、高光谱相机(7)、计算机(8)、信号线(9)、种子(10)组成，其中，机械臂(2)由机械臂升降台(21)、机械臂大臂(22)、机械臂小臂(23)、机械臂连接法兰(24)、手臂连接轴(25)组成，且机械臂升降台(21)的底部固装于机械臂基座(1)上，机械臂升降台(21)的滑动端与机械臂大臂(22)一端固装，机械臂大臂(22)另一端与机械臂小臂(23)一端通过手臂连接轴(25)滑动连接，机械臂小臂(23)另一端通过机械臂连接法兰(24)与真空吸种模块(3)的法兰连接件(31)固装；白光光源(5)通过光纤(6)为结构光投影模块(4)提供光源，结构光投影模块(4)投射结构光到种子(10)上，高光谱相机(7)采集种子(10)上投射的结构光图像并通过信号线(9)传输到计算机(8)。
2.权利要求1一种基于结构高光谱的种子质量检测装置，所述的真空吸种模块(3)，其特征为：由法兰连接件(31)、真空管转接头(32)、气吸空心针(33)、真空发生器(34)和软气管(35)组成，其中，真空发生器(34)输出端和软气管(35)一端相连，软气管(35)另一端通过真空管转接头(32)连接至法兰连接件(31)一端，法兰连接件(31)另一端与气吸空心针(33)一端相连，气吸空心针(33)另一端可吸附种子(10)，且法兰连接件(31)，气吸空心针(33)，真空管转接头(32)，软气管(35)均为真空结构，并共同构成内部空心连通体，可吸附种子(10)。
3.使用权利要求1一种基于结构高光谱的种子质量检测装置的一种基于结构高光谱的种子质量检测方法，其特征为，包括以下步骤：
步骤1：基于结构高光谱的种子质量检测装置标定；
步骤2：种子三维表型信息获取；
步骤3：种子特征光谱信息和二维纹理信息的获取；
步骤4：基于支持向量回归(SVR)的种子活力预测模型建模。
4.权利要求2所述的一种基于结构高光谱的种子质量检测方法，其特征为：所述步骤1基于结构高光谱的种子质量检测装置标定按照以下步骤标定：
步骤a：将标准棋盘格标定板垂直置于气吸空心针下方，标定板垂直轴线与气吸空心针同轴，标定板平面中心与结构光投影模块出光口中心的连线与标定板平面成90°，通过电脑8生成条纹形格雷码；
步骤b：电脑8通过信号线9将一组条纹形格雷码顺次发送到结构光投影模块，结构光投影模块顺次将编码结构光投射到标准棋盘格标定板上，电脑8再通过信号线9和高光谱相机7顺次采集标准棋盘格标定板上的结构光图像信息，并顺次存储保存为以***数字命名的jpg格式文件；
步骤c：电脑8将步骤b中存储的一组以***数字命名的jpg图片格式文件保存于angle1文件夹中，再分别将标准棋盘格标定板旋转2种不同角度，并根据步骤b的方法分别把2组保存为以***数字命名的jpg格式文件分别存储在angle2和angle3文件夹中；
步骤d：采用张正友棋盘格标定法，将步骤c所述angle1、angle2、angel3文件夹导入到MATALB calibration toolbox中，得到高光谱相机和投影仪的内参矩阵KL、KR，以及投影仪坐标系转换到高光谱相机坐标系的旋转R矩阵和平移T矩阵，其中内参矩阵包含每个像素的物理尺寸、焦距、图像物理坐标的扭曲因子、图像原点相对于光心成像点的纵横偏移量；
步骤e：根据步骤d得到的高光谱相机和投影仪的内参矩阵KL，KR，以及投影仪坐标系转换到高光谱相机坐标系的旋转R矩阵和平移T矩阵进行格雷码解码，得到三幅点云图像，分别命名为point1，point2point3；
步骤f：在point1点云图像中任取一个点作为参考点PR(x，y，z)，其在point2、point3中对应的点为PR1(x，y，z)，PR2(x，y，z)；
步骤g：基于3个点PR(x，y，z)、PR1(x，y，z)、PR2(x，y，z)的坐标信息，拟合出圆并计算出圆心坐标作为种子旋转坐标系的坐标原点P0(x，y，z)；
步骤h：设P0(x，y，z)为坐标原点，PR(x，y，z)为机械臂旋转台工作表面上的参考点，PR1(x，y，z)和PR2(x，y，z)为其在两个不同视场中所对应的坐标；
步骤i：通过
Figure FSA0000211498130000021
两个向量的向量积可求出Z′轴的方向向量，其中，
Figure FSA0000211498130000022
Figure FSA0000211498130000023
其中，上述步骤e中格雷码解码采用以下方法：
设三维空间中的物体在相机和投影仪的成像位置分别为Pc(uc，vc)和Pp(up，vp)，且设世界坐标系的原点与相机坐标系的原点重合，则相机的外参矩阵(相机和世界坐标系的转换矩阵)为
Figure FSA0000211498130000024
投影仪的外参矩阵(投影仪和世界坐标系的转换矩阵)，
Rp＝R，Tp＝T (2)
利用相机模型可得
Figure FSA0000211498130000025
Figure FSA0000211498130000026
联立上述方程组，获得Sc，Sp，X，Y，Z，和方程
Figure FSA0000211498130000027
Figure FSA0000211498130000028
则该结构光3D成像***每个像素点的三维解为
Figure FSA0000211498130000029
其中
Figure FSA00002114981300000210
和
Figure FSA00002114981300000211
表示Mc和Mp在第i行第j列。
5.权利要求2所述的一种基于结构高光谱的种子质量检测方法，其特征为：所述步骤2中种子三维表型信息获取流程为：
步骤a：电脑8通过信号线9将一组结构光编码信息发送到结构光投影模块，结构光投影模块将编码结构光投射到种子上，电脑8再通过信号9和高光谱相机采集种子的图形信息，每张图片存储成jpg格式，然后进行格雷码解码得到第一视角的点云图像；
步骤b：分别将种子旋转90°、180°、270°，同时根据步骤a的方法得到三个视角的点云图；
步骤c：将四个视角的三维点云图根据上述标定的种子旋转轴Z′变换到统一坐标系，完成点云的粗配准，调用MATLAB的ICP算法和平滑算法实现点云精配准；
步骤d：封装配准后的点云，计算长、宽、高、体积、外表面面积作为种子的三维表型参数。
6.权利要求2所述的一种基于结构高光谱的种子质量检测方法，其特征为：所述步骤3种子特征光谱信息和二维纹理信息的获取，按照以下步骤：
步骤a：结构光投影模块向种子投射白光，高光谱相机获取光谱图像立方体；
步骤b：采用MATLAB中OTSU方法进行阈值分割，获取种子的ROI掩模图像，通过ROI掩模图像提取高光谱图像中的种子ROI光谱信息；
步骤c：对光谱调用MATLAB中smooth平滑滤波函数进行平滑滤波；
步骤d：调用MATLAB中PCA算法提取高光谱的前三个特征波长；
步骤e：提取上述三个特征波长对应的灰度图像；
步骤f：用MATLAB调用graycomatrix函数计算上述三个特征波长对应灰度图像的生成灰度共生矩阵，其中包含种子的二维纹理信息。
7.权利要求2所述的基于支持向量回归(SVR)的种子质量预测模型建模，按照以下方法实施，调用MATLAB中Lib-SVM工具箱中的SVR算法，以上述前三个特征波长、灰度共生矩阵作为输入，以种子活力指数的数值为输出，建立种子活力预测模型，可用于种子活力预测。
CN202010545125.8A 2020-06-15 2020-06-15 基于结构高光谱的种子质量检测方法及装置 Active CN111579511B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202010545125.8A CN111579511B (zh) 2020-06-15 2020-06-15 基于结构高光谱的种子质量检测方法及装置

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202010545125.8A CN111579511B (zh) 2020-06-15 2020-06-15 基于结构高光谱的种子质量检测方法及装置

Publications (2)

Publication Number Publication Date
CN111579511A true CN111579511A (zh) 2020-08-25
CN111579511B CN111579511B (zh) 2024-07-23

Family

ID=72112819

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202010545125.8A Active CN111579511B (zh) 2020-06-15 2020-06-15 基于结构高光谱的种子质量检测方法及装置

Country Status (1)

Country Link
CN (1) CN111579511B (zh)

Cited By (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN113379824A (zh) * 2021-06-10 2021-09-10 浙江大学 一种基于双视角点云配准的类圆形水果纵横径测量方法
CN115205319A (zh) * 2022-09-19 2022-10-18 山东省金春雨种业科技有限公司 用于选种过程中的种子特征提取及分类方法
CN117405607A (zh) * 2023-12-15 2024-01-16 慧诺云谱(海南)科技有限公司 种子表型测量***、方法及相关设备

Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN107133976A (zh) * 2017-04-24 2017-09-05 浙江大学 一种获取三维高光谱信息的方法和装置
CN109087341A (zh) * 2018-06-07 2018-12-25 华南农业大学 一种近距离高光谱相机与测距传感器的融合方法
WO2019174572A1 (zh) * 2018-03-12 2019-09-19 浙江大学 一种利用种胚谱图特征波段识别谷物种子冻伤情况的方法
US20190281781A1 (en) * 2018-03-14 2019-09-19 Monsanto Technology Llc Seed Imaging
CN213749584U (zh) * 2020-06-15 2021-07-20 南京农业大学 基于结构高光谱的种子质量检测装置

Patent Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN107133976A (zh) * 2017-04-24 2017-09-05 浙江大学 一种获取三维高光谱信息的方法和装置
WO2019174572A1 (zh) * 2018-03-12 2019-09-19 浙江大学 一种利用种胚谱图特征波段识别谷物种子冻伤情况的方法
US20190281781A1 (en) * 2018-03-14 2019-09-19 Monsanto Technology Llc Seed Imaging
CN109087341A (zh) * 2018-06-07 2018-12-25 华南农业大学 一种近距离高光谱相机与测距传感器的融合方法
CN213749584U (zh) * 2020-06-15 2021-07-20 南京农业大学 基于结构高光谱的种子质量检测装置

Cited By (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN113379824A (zh) * 2021-06-10 2021-09-10 浙江大学 一种基于双视角点云配准的类圆形水果纵横径测量方法
CN113379824B (zh) * 2021-06-10 2022-04-01 浙江大学 一种基于双视角点云配准的类圆形水果纵横径测量方法
CN115205319A (zh) * 2022-09-19 2022-10-18 山东省金春雨种业科技有限公司 用于选种过程中的种子特征提取及分类方法
CN115205319B (zh) * 2022-09-19 2023-01-10 山东省金春雨种业科技有限公司 用于选种过程中的种子特征提取及分类方法
CN117405607A (zh) * 2023-12-15 2024-01-16 慧诺云谱(海南)科技有限公司 种子表型测量***、方法及相关设备
CN117405607B (zh) * 2023-12-15 2024-02-23 慧诺云谱(海南)科技有限公司 种子表型测量***、方法及相关设备

Also Published As

Publication number Publication date
CN111579511B (zh) 2024-07-23

Similar Documents

Publication Publication Date Title
CN111579511B (zh) 基于结构高光谱的种子质量检测方法及装置
CN110118528B (zh) 一种基于棋盘靶标的线结构光标定方法
CN109269430B (zh) 基于深度提取模型的多株立木胸径被动测量方法
CN105046746B (zh) 一种数字散斑人体三维快速扫描方法
CN109029284A (zh) 一种基于几何约束的三维激光扫描仪与相机标定方法
CN110378969B (zh) 一种基于3d几何约束的汇聚式双目相机标定方法
CN109579695B (zh) 一种基于异构立体视觉的零件测量方法
CN110672020A (zh) 一种基于单目视觉的立木高度测量方法
CN111336948B (zh) 基于成像平面转换的非标定手持式廓形检测方法及装置
CN110443879B (zh) 一种基于神经网络的透视误差补偿方法
CN113137920A (zh) 一种水下测量的设备与水下测量方法
CN109978940B (zh) 一种sab安全气囊尺寸视觉测量方法
CN104224212A (zh) Ct***、其扫描定位方法及校准方法
CN107167073A (zh) 一种线阵结构光三维快速测量装置及其测量方法
CN110223355B (zh) 一种基于双重极线约束的特征标志点匹配方法
CN114283203B (zh) 一种多相机***的标定方法及***
CN111402411A (zh) 一种基于线结构光的散乱物识别抓取方法
CN110608685A (zh) 一种基于树莓派的物体尺寸快速测量方法
CN109272555B (zh) 一种rgb-d相机的外部参数获得及标定方法
CN107869954A (zh) 一种双目视觉体积重量测量***及其实现方法
CN115375842A (zh) 一种植株三维重建方法、终端及存储介质
CN116051659A (zh) 一种线阵相机与2d激光扫描仪联合标定方法
CN109461183A (zh) 一种空间坐标点与点云位置点三维对比的方法
CN109493418B (zh) 一种基于LabVIEW的三维点云获取方法
CN108614277A (zh) 双激光单摄像头三维成像扫描台及扫描、成像方法

Legal Events

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