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

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

Info

Publication number
CN111579511B
CN111579511B CN202010545125.8A CN202010545125A CN111579511B CN 111579511 B CN111579511 B CN 111579511B CN 202010545125 A CN202010545125 A CN 202010545125A CN 111579511 B CN111579511 B CN 111579511B
Authority
CN
China
Prior art keywords
seed
seeds
structured light
mechanical arm
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.)
Active
Application number
CN202010545125.8A
Other languages
English (en)
Other versions
CN111579511A (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

Classifications

    • GPHYSICS
    • G01MEASURING; TESTING
    • G01NINVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
    • G01N21/00Investigating or analysing materials by the use of optical means, i.e. using sub-millimetre waves, infrared, visible or ultraviolet light
    • G01N21/17Systems in which incident light is modified in accordance with the properties of the material investigated
    • G01N21/25Colour; Spectral properties, i.e. comparison of effect of material on the light at two or more different wavelengths or wavelength bands
    • G01N21/31Investigating relative effect of material at wavelengths characteristic of specific elements or molecules, e.g. atomic absorption spectrometry
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01BMEASURING LENGTH, THICKNESS OR SIMILAR LINEAR DIMENSIONS; MEASURING ANGLES; MEASURING AREAS; MEASURING IRREGULARITIES OF SURFACES OR CONTOURS
    • G01B11/00Measuring arrangements characterised by the use of optical techniques
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01BMEASURING LENGTH, THICKNESS OR SIMILAR LINEAR DIMENSIONS; MEASURING ANGLES; MEASURING AREAS; MEASURING IRREGULARITIES OF SURFACES OR CONTOURS
    • G01B11/00Measuring arrangements characterised by the use of optical techniques
    • G01B11/02Measuring arrangements characterised by the use of optical techniques for measuring length, width or thickness
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T7/00Image analysis
    • G06T7/80Analysis of captured images to determine intrinsic or extrinsic camera parameters, i.e. camera calibration
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T2207/00Indexing scheme for image analysis or image enhancement
    • G06T2207/30Subject of image; Context of image processing
    • G06T2207/30204Marker
    • G06T2207/30208Marker matrix
    • YGENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
    • Y02TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
    • Y02ATECHNOLOGIES FOR ADAPTATION TO CLIMATE CHANGE
    • Y02A40/00Adaptation technologies in agriculture, forestry, livestock or agroalimentary production
    • Y02A40/10Adaptation technologies in agriculture, forestry, livestock or agroalimentary production in agriculture

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:通过两个向量的向量积可求出Z′轴的方向向量,其中,
其中,上述步骤e中格雷码解码采用以下方法:
设三维空间中的物体在相机和投影仪的成像位置分别为Pc(uc,vc)和Pp(up,vp),且设世界坐标系的原点与相机坐标系的原点重合,则相机的外参矩阵(相机和世界坐标系的转换矩阵)为
投影仪的外参矩阵(投影仪和世界坐标系的转换矩阵),
Rp=R,Tp=T (2)
利用相机模型可得
联立上述方程组,获得Sc,Sp,X,Y,Z,和方程
则该结构光3D成像系统每个像素点的三维解为
其中表示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:通过两个向量的向量积可求出Z′轴的方向向量,其中,
其中,上述步骤e中格雷码解码采用以下方法:
设三维空间中的物体在相机和投影仪的成像位置分别为Pc(uc,vc)和Pp(up,vp),且设世界坐标系的原点与相机坐标系的原点重合,则相机的外参矩阵(相机和世界坐标系的转换矩阵)为
投影仪的外参矩阵(投影仪和世界坐标系的转换矩阵),
Rp=R,Tp=T (2)
利用相机模型可得
联立上述方程组,获得Sc,Sp,X,Y,Z,和方程
则该结构光3D成像系统每个像素点的三维解为
其中表示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 (4)

1.一种基于结构高光谱的种子质量检测方法,其特征为,包括以下步骤:
步骤1:基于结构高光谱的种子质量检测装置标定;
步骤2:种子三维表型信息获取;
步骤3:种子特征光谱信息和二维纹理信息的获取;
步骤4:基于支持向量回归SVR的种子活力预测模型建模;
所述结构高光谱的种子质量检测装置由机械臂基座(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基于结构高光谱的种子质量检测装置标定按照以下步骤标定:
步骤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 calibrationtoolbox中,得到高光谱相机和投影仪的内参矩阵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:通过两个向量的向量积可求出Z'轴的方向向量,其中,
其中,上述步骤e中格雷码解码采用以下方法:
设三维空间中的物体在相机和投影仪的成像位置分别为Pc(uc,vc)和Pp(up,vp),且设世界坐标系的原点与相机坐标系的原点重合,则相机的外参矩阵,即相机和世界坐标系的转换矩阵为
投影仪的外参矩阵,即投影仪和世界坐标系的转换矩阵,
Rp=R,Tp=T (2)
利用相机模型可得
联立上述方程组,获得Sc,Sp,X,Y,Z,和方程
则该结构光3D成像系统每个像素点的三维解为
其中表示Mc和Mp在第i行第j列。
2.根据权利要求1所述的一种基于结构高光谱的种子质量检测方法,其特征为:所述步骤2中种子三维表型信息获取流程为:
步骤a:计算机(8)通过信号线(9)将一组结构光编码信息发送到结构光投影模块,结构光投影模块将编码结构光投射到种子上,计算机(8)再通过信号线(9)和高光谱相机采集种子的图形信息,每张图片存储成jpg格式,然后进行格雷码解码得到第一视角的点云图像;
步骤b:分别将种子旋转90°、180°、270°,同时根据步骤a的方法得到三个视角的点云图;
步骤c:将四个视角的三维点云图根据标定的种子旋转轴Z'变换到统一坐标系,完成点云的粗配准,调用MATLAB的ICP算法和平滑算法实现点云精配准;
步骤d:封装配准后的点云,计算长、宽、高、体积、外表面面积作为种子的三维表型参数。
3.根据权利要求1所述的一种基于结构高光谱的种子质量检测方法,其特征为:所述步骤3种子特征光谱信息和二维纹理信息的获取,按照以下步骤:
步骤a:结构光投影模块向种子投射白光,高光谱相机获取光谱图像立方体;
步骤b:采用MATLAB中OTSU方法进行阈值分割,获取种子的ROI掩模图像,通过ROI掩模图像提取高光谱图像中的种子ROI光谱信息;
步骤c:对光谱调用MATLAB中smooth平滑滤波函数进行平滑滤波;
步骤d:调用MATLAB中PCA算法提取高光谱的前三个特征波长;
步骤e:提取上述三个特征波长对应的灰度图像;
步骤f:用MATLAB调用graycomatrix函数计算上述三个特征波长对应灰度图像的生成灰度共生矩阵,其中包含种子的二维纹理信息。
4.根据权利要求1所述的一种基于结构高光谱的种子质量检测方法,其特征在于所述基于支持向量回归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 CN111579511A (zh) 2020-08-25
CN111579511B true 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)

Families Citing this family (4)

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

Citations (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN213749584U (zh) * 2020-06-15 2021-07-20 南京农业大学 基于结构高光谱的种子质量检测装置

Family Cites Families (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN107133976B (zh) * 2017-04-24 2020-04-14 浙江大学 一种获取三维高光谱信息的方法和装置
CN108444928B (zh) * 2018-03-12 2020-10-09 浙江大学 一种利用种胚谱图特征波段识别谷物种子冻伤情况的方法
MX2020009274A (es) * 2018-03-14 2020-10-01 Monsanto Technology Llc Imaginologia de semillas.
CN109087341B (zh) * 2018-06-07 2022-07-05 华南农业大学 一种近距离高光谱相机与测距传感器的融合方法

Patent Citations (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN213749584U (zh) * 2020-06-15 2021-07-20 南京农业大学 基于结构高光谱的种子质量检测装置

Also Published As

Publication number Publication date
CN111579511A (zh) 2020-08-25

Similar Documents

Publication Publication Date Title
CN111579511B (zh) 基于结构高光谱的种子质量检测方法及装置
CN110276808B (zh) 一种单相机结合二维码测量玻璃板不平度的方法
CN109269430B (zh) 基于深度提取模型的多株立木胸径被动测量方法
CN106595528B (zh) 一种基于数字散斑的远心显微双目立体视觉测量方法
CN109945853B (zh) 一种基于3d点云航拍图像的地理坐标定位系统和方法
CN105046746B (zh) 一种数字散斑人体三维快速扫描方法
CN110378969B (zh) 一种基于3d几何约束的汇聚式双目相机标定方法
CN110728715A (zh) 一种智能巡检机器人像机角度自适应调整方法
CN110672020A (zh) 一种基于单目视觉的立木高度测量方法
CN102221331B (zh) 一种基于不对称双目立体视觉技术的测量方法
CN107084680B (zh) 一种基于机器单目视觉的目标深度测量方法
CN110223355B (zh) 一种基于双重极线约束的特征标志点匹配方法
CN109978940B (zh) 一种sab安全气囊尺寸视觉测量方法
CN109961485A (zh) 一种基于单目视觉进行目标定位的方法
CN109448043A (zh) 平面约束下的立木高度提取方法
CN109272555B (zh) 一种rgb-d相机的外部参数获得及标定方法
CN110608685A (zh) 一种基于树莓派的物体尺寸快速测量方法
CN107144241A (zh) 一种基于景深补偿的双目视觉高精度测量方法
CN111402411A (zh) 一种基于线结构光的散乱物识别抓取方法
CN115375842A (zh) 一种植株三维重建方法、终端及存储介质
CN110060304B (zh) 一种生物体三维信息采集方法
CN116051659A (zh) 一种线阵相机与2d激光扫描仪联合标定方法
CN111462198A (zh) 一种尺度、旋转和辐射不变性的多模态影像配准方法
CN113963067B (zh) 一种采用小靶标对大视场视觉传感器进行标定的标定方法
CN112489141B (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
GR01 Patent grant