CN116612178A - 一种一体化水稻穗部谷粒粒型性状的提取方法 - Google Patents
一种一体化水稻穗部谷粒粒型性状的提取方法 Download PDFInfo
- Publication number
- CN116612178A CN116612178A CN202310900787.6A CN202310900787A CN116612178A CN 116612178 A CN116612178 A CN 116612178A CN 202310900787 A CN202310900787 A CN 202310900787A CN 116612178 A CN116612178 A CN 116612178A
- Authority
- CN
- China
- Prior art keywords
- grain
- image
- rice
- network
- training
- 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.)
- Pending
Links
- 235000013339 cereals Nutrition 0.000 title claims abstract description 357
- 235000007164 Oryza sativa Nutrition 0.000 title claims abstract description 137
- 235000009566 rice Nutrition 0.000 title claims abstract description 137
- 238000000034 method Methods 0.000 title claims description 36
- 240000007594 Oryza sativa Species 0.000 title 1
- 241000209094 Oryza Species 0.000 claims abstract description 137
- 238000012549 training Methods 0.000 claims abstract description 62
- 238000001514 detection method Methods 0.000 claims abstract description 58
- 210000005069 ears Anatomy 0.000 claims abstract description 45
- 238000000605 extraction Methods 0.000 claims abstract description 36
- 238000012360 testing method Methods 0.000 claims abstract description 32
- 230000008439 repair process Effects 0.000 claims abstract description 28
- 230000006870 function Effects 0.000 claims description 16
- 238000010586 diagram Methods 0.000 claims description 14
- 238000003384 imaging method Methods 0.000 claims description 14
- 230000008569 process Effects 0.000 claims description 10
- 238000004364 calculation method Methods 0.000 claims description 9
- 238000002372 labelling Methods 0.000 claims description 9
- 238000012795 verification Methods 0.000 claims description 8
- 238000013528 artificial neural network Methods 0.000 claims description 6
- 230000001629 suppression Effects 0.000 claims description 5
- 238000013135 deep learning Methods 0.000 claims description 4
- 230000004913 activation Effects 0.000 claims description 3
- 238000006243 chemical reaction Methods 0.000 claims description 3
- 239000002131 composite material Substances 0.000 claims description 3
- 230000005684 electric field Effects 0.000 claims description 3
- 238000011156 evaluation Methods 0.000 claims description 3
- 238000003306 harvesting Methods 0.000 claims description 3
- 238000011176 pooling Methods 0.000 claims description 3
- 238000005070 sampling Methods 0.000 claims description 3
- 230000011218 segmentation Effects 0.000 claims description 3
- 238000000638 solvent extraction Methods 0.000 claims description 3
- 230000005764 inhibitory process Effects 0.000 claims description 2
- 238000012545 processing Methods 0.000 abstract description 5
- 238000005516 engineering process Methods 0.000 description 7
- 238000005259 measurement Methods 0.000 description 7
- 239000004464 cereal grain Substances 0.000 description 5
- 239000011159 matrix material Substances 0.000 description 3
- 208000006440 Open Bite Diseases 0.000 description 2
- 238000013527 convolutional neural network Methods 0.000 description 2
- 238000013461 design Methods 0.000 description 2
- 238000012214 genetic breeding Methods 0.000 description 2
- 238000012986 modification Methods 0.000 description 2
- 230000004048 modification Effects 0.000 description 2
- 108090000623 proteins and genes Proteins 0.000 description 2
- 230000008485 antagonism Effects 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 230000000903 blocking effect Effects 0.000 description 1
- 238000010276 construction Methods 0.000 description 1
- 238000011161 development Methods 0.000 description 1
- 230000018109 developmental process Effects 0.000 description 1
- 235000013305 food Nutrition 0.000 description 1
- 230000004927 fusion Effects 0.000 description 1
- 230000014509 gene expression Effects 0.000 description 1
- 238000000691 measurement method Methods 0.000 description 1
- 238000006467 substitution reaction Methods 0.000 description 1
- 238000010200 validation analysis Methods 0.000 description 1
- 235000020985 whole grains Nutrition 0.000 description 1
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T7/00—Image analysis
- G06T7/60—Analysis of geometric attributes
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06N—COMPUTING ARRANGEMENTS BASED ON SPECIFIC COMPUTATIONAL MODELS
- G06N3/00—Computing arrangements based on biological models
- G06N3/02—Neural networks
- G06N3/04—Architecture, e.g. interconnection topology
- G06N3/0464—Convolutional networks [CNN, ConvNet]
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06N—COMPUTING ARRANGEMENTS BASED ON SPECIFIC COMPUTATIONAL MODELS
- G06N3/00—Computing arrangements based on biological models
- G06N3/02—Neural networks
- G06N3/04—Architecture, e.g. interconnection topology
- G06N3/048—Activation functions
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06N—COMPUTING ARRANGEMENTS BASED ON SPECIFIC COMPUTATIONAL MODELS
- G06N3/00—Computing arrangements based on biological models
- G06N3/02—Neural networks
- G06N3/08—Learning methods
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T7/00—Image analysis
- G06T7/60—Analysis of geometric attributes
- G06T7/62—Analysis of geometric attributes of area, perimeter, diameter or volume
-
- 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/20—Special algorithmic details
- G06T2207/20016—Hierarchical, coarse-to-fine, multiscale or multiresolution image processing; Pyramid transform
-
- 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/20—Special algorithmic details
- G06T2207/20081—Training; Learning
-
- 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/20—Special algorithmic details
- G06T2207/20084—Artificial neural networks [ANN]
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- General Physics & Mathematics (AREA)
- Biophysics (AREA)
- Evolutionary Computation (AREA)
- Health & Medical Sciences (AREA)
- Life Sciences & Earth Sciences (AREA)
- Artificial Intelligence (AREA)
- Biomedical Technology (AREA)
- Software Systems (AREA)
- Computational Linguistics (AREA)
- Data Mining & Analysis (AREA)
- Mathematical Physics (AREA)
- General Health & Medical Sciences (AREA)
- Molecular Biology (AREA)
- Computing Systems (AREA)
- General Engineering & Computer Science (AREA)
- Computer Vision & Pattern Recognition (AREA)
- Geometry (AREA)
- Image Analysis (AREA)
Abstract
本发明涉及图像处理技术领域,具体提供了一种一体化水稻穗部谷粒粒型性状的提取方法,包括如下步骤:S1:穗部谷粒检测;S11:采集原始稻穗图像,构建稻穗图像数据集;S12:对稻穗图像数据集进行训练、验证与测试;S2:遮挡谷粒修复;S21:对遮挡谷粒图像进行采集并获得成对谷粒图像训练集;S22:对成对谷粒图像训练集进行训练与测试;S3:谷粒粒型性状提取,得到谷粒长度、宽度、周长以及谷粒投影面积;S4:获取一体化穗部谷粒检测修复模型及一体化穗部谷粒粒型性状提取模型。本方案可对自然形态下受到遮挡的水稻穗部谷粒进行修复而解决穗部谷粒间相互遮挡的问题,提高了水稻穗部谷粒粒型性状提取的效率与精度。
Description
技术领域
本发明涉及图像处理技术领域,尤其涉及一种一体化水稻穗部谷粒粒型性状的提取方法。
背景技术
水稻是重要的粮食作物,已在世界各地被广泛种植。水稻穗部的粒型性状是重要的水稻基因体现,也是水稻基因研究的重要途径,在粒型性状中谷粒长度的增加会使谷粒长宽比提高、谷粒形状由短圆形向细长形变化,但会使谷粒宽度、谷粒厚度下降;谷粒宽度、谷粒厚度协同提高的同时,会使谷粒长宽比下降、谷粒形状由细长形向短圆形改变,谷粒形状与谷粒长宽比具有较高的相关性。
在当前水稻遗传育种的过程中,需要大量表型数据为水稻穗部的粒型基因选择作支撑。传统的表型测量方法主要依赖人工,工作量大、效率低下,且人工测量法受人主观性的影响较大,难以满足当前遗传育种的需求,同时当下劳动力成本逐年增高,水稻穗部粒型性状检测的成本也越来越高。随着成像技术和计算机技术的发展,高通量自动化进行表型测量方法已越来越普遍,但水稻穗部谷粒之间存在大量的相互遮挡,极大地限制了水稻谷粒性状提取的效率。目前,高通量水稻谷粒粒型性状主要分为两种:
1、将穗部谷粒脱粒后铺在平板上,使用振动装置来避免谷粒之间的相互遮挡,随后使用可见光相机进行成像,并使用数字图像处理技术提取谷粒粒型性状。此方法进行脱粒的过程耗时较长,且在脱粒过程中存在谷粒破损的情况,会对性状提取的精度造成一定的影响。
2、通过人工分开水稻枝梗后铺在平板上,使用可见光相机进行成像,使用数字图像处理技术获取穗部未受遮挡的谷粒并提取其粒型性状。该方法需要破坏水稻稻穗的自然形态,且在人工分开枝梗后仍有大量穗部谷粒存在相互遮挡的情况,导致无法使用数字图像处理技术提取其粒型性状,使得粒型性状提取的精度降低。
由上可知,上述方法均不能实现水稻谷粒粒型性状快速准确的提取,水稻穗部谷粒之间的遮挡已成为水稻穗部粒型性状提取的瓶颈。
同时,目前针对水稻穗部谷粒的表型测量主要是通过深度学习目标检测技术实现穗部谷粒的检测与计数。基于自然形态下稻穗上的谷粒间存在严重的相互遮挡,即使通过手动分开水稻枝梗的方法可以在很大程度上减轻谷粒遮挡带来的影响,但人工分开枝梗的的过程中会对稻穗进行破坏,是一个有损的过程,且在执行上述操作后仍无法直接提取部分穗部谷粒的粒型性状,因此目前的检测网络结构难以实现自然形态下水稻穗部谷粒的准确检测与计数。
综上所述,如何设计一种无需对稻穗进行破坏即可解决谷粒之间相互遮挡的问题,可提高谷粒粒型检测效率及检测精度的一体化水稻穗部谷粒粒型性状提取方法,是当下亟需解决的问题。
发明内容
本发明为解决上述问题,提供了一种一体化水稻穗部谷粒粒型性状的提取方法,可对自然形态下水稻穗部谷粒进行准确检测与计数,解决穗部谷粒间相互遮挡的问题,实现一体化穗部谷粒粒型性状提取,可极大提高水稻穗部谷粒粒型性状提取的效率与精度,避免了传统人工测量带来的主观性误差。
为达到上述目的,本发明提出如下技术方案:一种一体化水稻穗部谷粒粒型性状的提取方法,包括如下步骤:
S1:穗部谷粒检测;
S11:对自然形态下的原始稻穗图像进行采集并构建稻穗图像数据集;
S12:通过谷粒检测网络模型对稻穗图像数据集进行训练、验证与测试;
谷粒检测网络模型包括用以提取稻穗图像特征的特征提取网络、用以将特征提取网络在不同深度提取的特征进行融合的特征金字塔网络、用以生成可能存在目标的预测框的区域提案网络、用以对网络检测到的目标进行分类并对检测到的目标边界框位置进行优化的分类回归网络;
S2:遮挡谷粒修复;
S21:对遮挡谷粒图像进行采集并获得包括一张未受遮挡的谷粒图像与一张对应受到遮挡的谷粒图像的成对谷粒图像训练集;
S22:通过谷粒修复网络模型对成对谷粒图像训练集进行训练与测试;
谷粒修复网络模型包括用以将输入的遮挡图像修复为无遮挡图像的生成器、用以判定生成器生成的无遮挡图像真实性的判别器;
S3:谷粒粒型性状提取,得到谷粒长度、宽度、周长以及谷粒投影面积;
S4:获取一体化穗部谷粒粒型性状提取模型;
S41:获取S1中得到的最优谷粒检测网络模型;
S42:获取S2中得到的最优谷粒修复网络模型;
S43:将谷粒修复网络模型集成到谷粒检测网络模型中,置于区域提案网络之后,与分类回归网络并行,得到一体化穗部谷粒检测修复模型。
S44:将性状提取管道集成至一体化穗部谷粒检测修复模型中,得到一体化穗部谷粒粒型性状提取模型。
优选的,S12中的谷粒检测网络模型的目标函数公式如下所示:
其中,表示第i个anchor预测为目标的概率,/>为对应真实标签,如anchor内存在目标则/>为1,反之则为0;/>为一个坐标参数向量,表示预测边界框与真实边界框坐标之间的补偿,/>表示存在目标的anchor所对的真实边界框的坐标;/>表示训练时的批大小;/>表示区域提案网络生成的anchor的数量;
为一个交叉熵损失,用于对anchor中是否包含物体进行分类,具体计算公式如下所示:
;
为一个回归损失,用于使网络得到更精确的边界框,具体计算公式如下所示:
。
优选的,S22中的谷粒修复网络模型的目标函数公式如下所示:
其中,D为判别器,G为生成器,z为噪声,x为真实图像,y为约束条件;
为条件生成式对抗网络的目标函数:
为网络的L1损失函数:
。
优选的,S21包括如下子步骤:
S211:遮挡谷粒图像采集:获取整株水稻稻穗,将稻穗自然放置于扫描仪盖板上,在稻穗上选择一个受遮挡的谷粒,对受遮挡的谷粒进行固定,将固定好的受遮挡的谷粒连同稻穗共同进行成像,获取受遮挡的谷粒图像;
S212:未受遮挡谷粒图像采集:取下S211中的稻穗,仅保留固定在扫描仪盖板上的谷粒,再次进行扫描,获取对应未受遮挡的谷粒图像;
S213:合成遮挡谷粒图像对:对于未受遮挡的谷粒图像,取出图像的红色通道,对红色通道进行二值化,获取图像轮廓并求出谷粒轮廓的外接矩形;再根据外接矩形的坐标在遮挡图像中取出对应受到遮挡的子图,将两个子图在水平方向上进行拼接后得到合成的遮挡谷粒修复数据集;
S214:数据增强:将S213中获得的合成遮挡谷粒图像进行随机水平或垂直翻转,以扩增数据集,训练集构建完成且每一份训练集图像均包括成对的谷粒图像,即包括一张未受遮挡的谷粒图像与一张对应受到遮挡的谷粒图像;
S215:划分数据集:将S214中获得的成对谷粒图像作为数据集,将数据集中的图像按4:1的数量比例划分为训练集和测试集。
优选的,S22包括如下子步骤:
S221:将S215中获得的成对谷粒图像训练集输入至谷粒修复网络模型中进行训练,对神经网络中的参数进行更新,直到达到指定的迭代次数或指定的精度后,停止训练;
S222:将S215中获得的测试集中的合成遮挡谷粒图像输入至谷粒修复网络模型中得到复原后的完整谷粒图像,再将完整谷粒图像与测试集中的真实完整谷粒图像进行对比测试;
选择谷粒的粒长参数、粒宽参数、面积参数、周长参数4个粒型参数进行对比测试,计算4个粒型参数修复值与真实值的平均绝对百分比误差及相关系数R²。
优选的,S3包括如下子步骤:
S31:将单个谷粒图像进行灰度化,提取谷粒与背景之间的对比度最大的r通道的灰度图用于后续的形状提取;
S32:使用阈值分割算法对谷粒灰度图像进行二值化,获取谷粒二值图;
S33:基于谷粒二值图提取谷粒外轮廓,包括轮廓上各点的坐标;
S34:基于S33得到的轮廓坐标计算得到谷粒长度、宽度、周长及谷粒投影面积。
优选的,S34中谷粒长度、宽度、周长及谷粒投影面积的计算方法如下:
谷粒长度:基于谷粒轮廓上的点,求取任意两点之间的距离,取最大距离作为谷粒长度;
谷粒宽度:由求取谷粒长度的两点计算两点之间的直线方程,求取这个直线的法线斜率,随后求取不同截距下法线与轮廓的交点的距离,最大值即为谷粒宽度;
谷粒周长:基于谷粒轮廓上的点,若两点在上下或左右相邻,则两点间的距离定义为1像素大小;若两点在左上、左下、右上或右下方向相邻,则两点之间的距离定义为像素大小;将轮廓上所有相邻像素点之间的距离求和即为谷粒周长;
谷粒面积:基于谷粒轮廓上的点,谷粒面积即为谷粒轮廓内像素的数量。
优选的,S11包括如下子步骤:
S111:稻穗图像采集:在成熟水稻收割后获得水稻主穗,将稻穗随意平铺在可见光扫描成像设备上,对稻穗样本进行成像及存储;
S112:原始图像裁剪:基于S111采集的稻穗图像,将稻穗图像中的冗余部分去除,并裁剪为统一像素大小的图像;
S113:穗部谷粒图像标注:基于S112中得到的裁剪图像,采用边界框标注方式通过软件对所有穗部谷粒进行标注,每个谷粒标注后将边界框的左上角坐标、右下角坐标、边界框面积和框内目标种类记录下来,再将整株稻穗上的谷粒标注完后进行存储;
S114:数据增强:对由S112和S113中获取的稻穗图像及标注文件以一定的概率随机进行水平翻转或垂直翻转,数据集中的每个样本包括裁剪后的稻穗图像及相对应的标注文件;
S115:数据集划分:在经过S114中的数据增强后,将数据集中的样本以2:1:1的比例划分为训练集、验证集和测试集;
S116:数据集格式转换:分别将由S115划分后的数据集样本转换为用于深度学习目标检测方法训练所规定的格式。
优选的,S12包括如下子步骤:
S121:网络训练:将由S116获取的稻穗图像训练集输入到谷粒检测网络模型中进行训练,采用不同的超参数对网络进行训练,训练过程中不断更新神经网络的参数,直到达到指定的迭代次数,停止训练;
S122:网络验证:将由S116获取的稻穗图像验证集输入到谷粒检测网络模型中进行验证,根据目标检测评价指标获取不同超参数下网络的精度,随之获取最优谷粒检测网络;
S123:网络测试:将由S116获取的稻穗图像测试集输入到由S122获取的最优谷粒检测网络中进行测试,得到每个水稻稻穗样本上的谷粒数,计算网络预测值与真实值的平均绝对百分比误差及相关系数R²。
优选的,S12中的特征提取网络包括:输入层,用于输入采集到的稻穗图像;卷积层,用于提取图像中的信息;池化层,用于对由卷积层提取的信息进行挑选,对信息降维;激活层,用于提高网络的非线性拟合能力;输出层,用于输出经过多个卷积池化激活后富含图像高级语义信息的特征图;
特征金字塔网络包括:卷积层,用于对每个阶段输出的特征图进行降维;上采样层:用于对网络得到的特征图进行上采样;
区域提案网络包括:滑动窗口,在由特征金字塔网络得到的特征图上进行滑动生成不同大小的锚框用于目标的预测;分类通道,由若干卷积层构成,用于对锚框中的目标进行分类;回归通道,由若干卷积层构成,用于对锚框的位置进行精调;非极大值抑制,用于消除重复检测的锚框及超出图像边界的锚框;
分类回归网络包括:分类通道,由若干卷积层构成,用于对由区域提案网络生成的提案区域中的目标进行分类;回归通道,由若干卷积层构成,用于对由区域提案网络生成的提案区域的位置进行精调;非极大值抑制,用于消除冗余的检测框,并输出最终目标检测的结果。
与现有技术相比,本发明能够取得如下有益效果:
1、本发明可在不对水稻稻穗进行分穗的情况下高精度实现穗部谷粒的计数并对穗部谷粒进行准确定位,避免了传统人工测量带来的主观性误差,省去了手工分开水稻稻穗枝梗所需的时间,避免了对水稻的自然形态造成破坏,节省了工作量及劳动力成本。
2、本发明可通过基于生成对抗网络的图像修复技术对自然形态下受到遮挡的水稻穗部谷粒进行修复,使受遮挡的谷粒恢复自然形状,极大地避免了由于遮挡导致无法直接进行表型测量的瓶颈,提高了水稻穗部谷粒粒型性状提取的效率与精度。
附图说明
图1是根据本发明实施例提供的一体化穗部谷粒粒型性状提取的整体结构示意图;
图2是根据本发明实施例提供的图像采集设备及采集过程的示意图;
图3是根据本发明实施例提供的穗部谷粒检测网络模型的结构图;
图4是根据本发明实施例提供的谷粒修复网络模型的整体结构图;
图5是根据本发明实施例提供的谷粒粒型性状提取过程示意图。
具体实施方式
在下文中,将参考附图1-5描述本发明的实施例。在下面的描述中,相同的模块使用相同的附图标记表示。在相同的附图标记的情况下,它们的名称和功能也相同。因此,将不重复其详细描述。
为了使本发明的目的、技术方案及优点更加清楚明白,以下结合附图1-5及具体实施例,对本发明进行进一步详细说明。应当理解,此处所描述的具体实施例仅用以解释本发明,而不构成对本发明的限制。
一种一体化水稻穗部谷粒粒型性状的提取方法,包括如下步骤:
S1:穗部谷粒检测。
S11:对自然形态下的原始稻穗图像进行采集并构建稻穗图像数据集,包括如下子步骤:
S111:稻穗图像采集:在成熟水稻收割后获得水稻主穗,将稻穗随意平铺在可见光扫描成像设备上,设置扫描精度为600dpi,对稻穗样本进行成像,并以tiff格式对稻穗样本图像进行存储,单株稻穗的成像时间为60s左右。图像采集过程如图2所示,其中a图为将稻穗随意平铺在扫描成像设备上,b图为对稻穗样本进行成像,c图为扫描得到的原始稻穗图像,d图为裁剪后的稻穗图像。
S112:原始图像裁剪:基于S111采集的稻穗图像,将稻穗图像中的冗余部分去除,并统一裁剪为5700×6800像素大小的图像。
S113:穗部谷粒图像标注:基于S112中得到的裁剪图像,采用边界框标注方式通过labelme软件对所有穗部谷粒进行标注,每个谷粒标注后将边界框的左上角坐标、右下角坐标、边界框面积和框内目标种类记录下来,再将整株稻穗上的谷粒标注完后以json格式进行存储。
S114:数据增强:对由S112和S113中获取的稻穗图像及标注文件以0.2的概率随机进行水平翻转或垂直翻转,数据集中的每个样本包括裁剪后的稻穗图像及相对应的标注文件。
S115:数据集划分:在经过S114中的数据增强后,将数据集中的样本以2:1:1的比例划分为训练集、验证集和测试集。
S116:数据集格式转换:分别将由S115划分后的数据集样本转换为用于深度学习目标检测方法训练所规定的格式。
S12:通过谷粒检测网络模型对稻穗图像数据集进行训练、验证与测试,包括如下子步骤:
S121:网络训练:将由S116获取的稻穗图像训练集输入到谷粒检测网络模型中进行训练,采用不同的超参数对网络进行训练,训练过程中不断更新神经网络的参数,直到达到指定的迭代次数,停止训练;
S122:网络验证:将由S116获取的稻穗图像验证集输入到谷粒检测网络模型中进行验证,根据目标检测评价指标获取不同超参数下网络的精度,随之获取最优谷粒检测网络;
S123:网络测试:将由S116获取的稻穗图像测试集输入到由S122获取的最优谷粒检测网络中进行测试,得到每个水稻稻穗样本上的谷粒数,计算网络预测值与真实值的平均绝对百分比误差及相关系数R²。
本实施例中的谷粒检测网络主要由基于卷积神经网络的Faster R-CNN构成,谷粒检测网络模型包括特征提取网络、特征金字塔网络、区域提案网络及分类回归网络。
特征提取网络用以提取稻穗的图像特征,特征提取网络包括:
输入层,用于输入采集到的稻穗图像;
卷积层,用于提取图像中的信息;
池化层,用于对由卷积层提取的信息进行挑选,对信息降维;
激活层,用于提高网络的非线性拟合能力;
输出层,用于输出经过多个卷积池化激活后富含图像高级语义信息的特征图。
特征金字塔网络用以将特征提取网络在不同深度提取的特征进行融合,特征金字塔网络包括:
卷积层,卷积核大小为1×1,用于对每个阶段输出的特征图进行降维;
上采样层:用于对网络得到的特征图进行上采样,使不同阶段的特征图具有相同的分辨率用于特征图之间的融合;
区域提案网络用以生成可能存在目标的预测框,区域提案网络包括:
滑动窗口,在由特征金字塔网络得到的特征图上进行滑动生成不同大小的锚框用于目标的预测;
分类通道,由若干卷积层构成,用于对锚框中的目标进行分类,分别输出目标属于某个类的概率值;
回归通道,由若干卷积层构成,用于对锚框的位置进行精调,使之与目标真实位置更加匹配;
非极大值抑制,用于消除重复检测的锚框及超出图像边界的锚框。
分类回归网络用以对网络检测到的目标进行分类并对检测到的目标边界框位置进行优化,分类回归网络包括:
分类通道,与区域提案网络中分类通道相同,由若干卷积层构成,用于对由区域提案网络生成的提案区域中的目标进行分类,分别输出目标属于某个类的概率值;
回归通道,与区域提案网络中分类通道相同,由若干卷积层构成,用于对由区域提案网络生成的提案区域的位置进行精调,使之与目标真实位置更加匹配;
非极大值抑制,与区域提案网络中非极大值抑制相同,用于消除冗余的检测框,并输出最终目标检测的结果。
穗部谷粒检测网络的结构如图3所示。谷粒检测网络模型的目标函数公式如下所示:
其中,表示第i个anchor预测为目标的概率,/>为对应真实标签,如anchor内存在目标则/>为1,反之则为0;/>为一个坐标参数向量,表示预测边界框与真实边界框坐标之间的补偿,/>表示存在目标的anchor所对的真实边界框的坐标;/>表示训练时的批大小;/>表示区域提案网络生成的anchor的数量;
为一个交叉熵损失,用于对anchor中是否包含物体进行分类,具体计算公式如下所示:
;
为一个回归损失,用于使网络得到更精确的边界框,具体计算公式如下所示:
。
S2:遮挡谷粒修复。
S21:对遮挡谷粒图像进行采集并获得包括一张未受遮挡的谷粒图像与一张对应受到遮挡的谷粒图像的成对谷粒图像训练集,包括如下子步骤:
S211:遮挡谷粒图像采集:获取整株水稻稻穗,将稻穗自然放置于扫描仪盖板上,在稻穗上选择一个受遮挡的谷粒,对受遮挡的谷粒进行固定,将固定好的受遮挡的谷粒连同稻穗共同进行成像,获取受遮挡的谷粒图像。
S212:未受遮挡谷粒图像采集:取下S211中的稻穗,仅保留固定在扫描仪盖板上的谷粒,再次进行扫描,获取对应未受遮挡的谷粒图像。
S213:合成遮挡谷粒图像对:对于未受遮挡的谷粒图像,取出图像的红色通道,对红色通道进行二值化,获取图像轮廓并求出谷粒轮廓的外接矩形;再根据外接矩形的坐标在遮挡图像中取出对应受到遮挡的子图,将两个子图在水平方向上进行拼接后得到合成的遮挡谷粒修复数据集。
S214:数据增强:将S213中获得的合成遮挡谷粒图像进行随机水平或垂直翻转,以扩增数据集,训练集构建完成且每一份训练集图像均包括成对的谷粒图像,即包括一张未受遮挡的谷粒图像与一张对应受到遮挡的谷粒图像。
S215:划分数据集:将S214中获得的成对谷粒图像作为数据集,将数据集中的图像按4:1的数量比例划分为训练集和测试集。
S22:通过谷粒修复网络模型对成对谷粒图像训练集进行训练与测试,包括如下子步骤:
S221:将S215中获得的成对谷粒图像训练集输入至谷粒修复网络模型中进行训练,对神经网络中的参数进行更新,直到达到指定的迭代次数或指定的精度后,停止训练。
S222:将S215中获得的测试集中的合成遮挡谷粒图像输入至谷粒修复网络模型中得到复原后的完整谷粒图像,再将完整谷粒图像与测试集中的真实完整谷粒图像进行对比测试。
选择谷粒的粒长参数、粒宽参数、面积参数、周长参数4个粒型参数进行对比测试,计算4个粒型参数修复值与真实值的平均绝对百分比误差及相关系数R²。
谷粒修复网络模型包括生成器和判别器,生成器用以将输入的遮挡图像修复为无遮挡图像,生成器包括:
输入层,用于输入人工合成的遮挡谷粒图像或实际的稻穗遮挡谷粒图像子图;
编码层,用于对输入图像编码,将输入图像下采样的同时增加通道数,以提取图像中的高层次语义信息;
解码层,用于对图像进行解码,将编码后的特征图上采样的同时减小通道数,抓取不同层次的语义信息以特征叠加的方式予以整合;编码层与解码层之间采用跳跃连接结构,假设网络总层数为n,则在每第i层和第n-i层之间添加跳跃连接,以此避免了编码-解码结构带来的信息瓶颈问题;
输出层,用于输出修复后的单颗谷粒图像,输出图像的尺寸与输入图像尺寸相同。
判别器用以判定生成器生成的无遮挡图像的真实性,判别器包括:
输入层,用于输入生成器生成的谷粒图像或真实的谷粒图像;
全卷积层,输入图像各区域都独立经过全卷积层,以提取各区域的特征便于网络判定真伪;
输出矩阵,输入图像经过全卷积层后被映射为一个输出矩阵,矩阵的各个值代表着输入图像各区域属于真实图像的概率。
谷粒修复网络模型的目标函数公式如下所示:
其中,D为判别器,G为生成器,z为噪声,x为真实图像,y为约束条件;
为条件生成式对抗网络的目标函数:
为网络的L1损失函数:
。
上述遮挡谷粒修复方法通过前后两次扫描获取的受遮挡谷粒图像和未受遮挡谷粒图像制作的图像对,使修复网络得到有效训练后,即可使用网络的生成器将遮挡谷粒图像修复为未遮挡的完整谷粒图像,对不同品种水稻谷粒复原的综合误差低于4%。利用该方法能够从本质上修复因遮挡造成的图像信息缺失,有效消除谷粒遮挡造成的测量误差,从而为水稻产量性状中最为重要的谷粒粒型测量提供一种高精度的测量途径。
S3:谷粒粒型性状提取,得到谷粒长度、宽度、周长以及谷粒投影面积,包括如下子步骤:
S31:将单个谷粒图像进行灰度化,提取谷粒与背景之间的对比度最大的r通道的灰度图用于后续的形状提取;
S32:使用阈值分割算法对谷粒灰度图像进行二值化,获取谷粒二值图;
S33:基于谷粒二值图提取谷粒外轮廓,包括轮廓上各点的坐标;
S34:基于S33得到的轮廓坐标计算得到谷粒长度、宽度、周长及谷粒投影面积。
S34中谷粒长度、宽度、周长及谷粒投影面积的计算方法如下:
谷粒长度:基于谷粒轮廓上的点,求取任意两点之间的距离,取最大距离作为谷粒长度;
谷粒宽度:由求取谷粒长度的两点计算两点之间的直线方程,求取这个直线的法线斜率,随后求取不同截距下法线与轮廓的交点的距离,最大值即为谷粒宽度;
谷粒周长:基于谷粒轮廓上的点,若两点在上下或左右相邻,则两点间的距离定义为1像素大小;若两点在左上、左下、右上或右下方向相邻,则两点之间的距离定义为像素大小;将轮廓上所有相邻像素点之间的距离求和即为谷粒周长;
谷粒面积:基于谷粒轮廓上的点,谷粒面积即为谷粒轮廓内像素的数量。
S4:获取一体化穗部谷粒粒型性状提取模型,谷粒粒型性状提取细节由图5所示,a图所示为谷粒原始图像,b图所示为通过R通道提取的灰度图像,c图所示为对灰度图像进行二值化后的二值图,d图所示为通过二值图提取出的谷粒轮廓,e图所示为对各种粒型性状进行提取。
S41:获取S1中得到的最优谷粒检测网络模型;
S42:获取S2中得到的最优谷粒修复网络模型。
S43:将谷粒修复网络模型集成到谷粒检测网络模型中,置于区域提案网络之后,与分类回归网络并行,得到一体化穗部谷粒检测修复模型。
S44:将性状提取管道集成至一体化穗部谷粒检测修复模型中,得到一体化穗部谷粒粒型性状提取模型。
应该理解,可以使用上面所示的各种形式的流程,重新排序、增加或删除步骤。例如,本发明公开中记载的各步骤可以并行地执行也可以顺序地执行也可以不同的次序执行,只要能够实现本发明公开的技术方案所期望的结果,本文在此不进行限制。
上述具体实施方式,并不构成对本发明保护范围的限制。本领域技术人员应该明白的是,根据设计要求和其他因素,可以进行各种修改、组合、子组合和替代。任何在本发明的精神和原则之内所作的修改、等同替换和改进等,均应包含在本发明保护范围之内。
Claims (10)
1.一种一体化水稻穗部谷粒粒型性状的提取方法,其特征在于,包括如下步骤:
S1:穗部谷粒检测;
S11:对自然形态下的原始稻穗图像进行采集并构建稻穗图像数据集;
S12:通过谷粒检测网络模型对稻穗图像数据集进行训练、验证与测试;
所述谷粒检测网络模型包括用以提取稻穗图像特征的特征提取网络、用以将特征提取网络在不同深度提取的特征进行融合的特征金字塔网络、用以生成可能存在目标的预测框的区域提案网络、用以对网络检测到的目标进行分类并对检测到的目标边界框位置进行优化的分类回归网络;
S2:遮挡谷粒修复;
S21:对遮挡谷粒图像进行采集并获得包括一张未受遮挡的谷粒图像与一张对应受到遮挡的谷粒图像的成对谷粒图像训练集;
S22:通过谷粒修复网络模型对成对谷粒图像训练集进行训练与测试;
所述谷粒修复网络模型包括用以将输入的遮挡图像修复为无遮挡图像的生成器、用以判定生成器生成的无遮挡图像真实性的判别器;
S3:谷粒粒型性状提取,得到谷粒长度、宽度、周长以及谷粒投影面积;
S4:获取一体化穗部谷粒粒型性状提取模型;
S41:获取S1中得到的最优谷粒检测网络模型;
S42:获取S2中得到的最优谷粒修复网络模型;
S43:将谷粒修复网络模型集成到谷粒检测网络模型中,置于区域提案网络之后,与分类回归网络并行,得到一体化穗部谷粒检测修复模型;
S44:将性状提取管道集成至一体化穗部谷粒检测修复模型中,得到一体化穗部谷粒粒型性状提取模型。
2.根据权利要求1所述的一体化水稻穗部谷粒粒型性状的提取方法,其特征在于,S12中的谷粒检测网络模型的目标函数公式如下所示:
其中,表示第i个anchor预测为目标的概率,/>为对应真实标签,如anchor内存在目标则/>为1,反之则为0;/>为一个坐标参数向量,表示预测边界框与真实边界框坐标之间的补偿,/>表示存在目标的anchor所对的真实边界框的坐标;/>表示训练时的批大小;表示区域提案网络生成的anchor的数量;
为一个交叉熵损失,用于对anchor中是否包含物体进行分类,具体计算公式如下所示:
;
为一个回归损失,用于使网络得到更精确的边界框,具体计算公式如下所示:
。
3.根据权利要求2所述的一体化水稻穗部谷粒粒型性状的提取方法,其特征在于,S22中的谷粒修复网络模型的目标函数公式如下所示:
其中,D为判别器,G为生成器,z为噪声,x为真实图像,y为约束条件;
为条件生成式对抗网络的目标函数:
为网络的L1损失函数:
。
4.根据权利要求1-3中任一项所述的一体化水稻穗部谷粒粒型性状的提取方法,其特征在于,S21包括如下子步骤:
S211:遮挡谷粒图像采集:获取整株水稻稻穗,将稻穗自然放置于扫描仪盖板上,在稻穗上选择一个受遮挡的谷粒,对受遮挡的谷粒进行固定,将固定好的受遮挡的谷粒连同稻穗共同进行成像,获取受遮挡的谷粒图像;
S212:未受遮挡谷粒图像采集:取下S211中的稻穗,仅保留固定在扫描仪盖板上的谷粒,再次进行扫描,获取对应未受遮挡的谷粒图像;
S213:合成遮挡谷粒图像对:对于未受遮挡的谷粒图像,取出图像的红色通道,对红色通道进行二值化,获取图像轮廓并求出谷粒轮廓的外接矩形;再根据外接矩形的坐标在遮挡图像中取出对应受到遮挡的子图,将两个子图在水平方向上进行拼接后得到合成的遮挡谷粒修复数据集;
S214:数据增强:将S213中获得的合成遮挡谷粒图像进行随机水平或垂直翻转,以扩增数据集,训练集构建完成且每一份训练集图像均包括成对的谷粒图像,即包括一张未受遮挡的谷粒图像与一张对应受到遮挡的谷粒图像;
S215:划分数据集:将S214中获得的成对谷粒图像作为数据集,将数据集中的图像按4:1的数量比例划分为训练集和测试集。
5.根据权利要求4所述的一体化水稻穗部谷粒粒型性状的提取方法,其特征在于,S22包括如下子步骤:
S221:将S215中获得的成对谷粒图像训练集输入至谷粒修复网络模型中进行训练,对神经网络中的参数进行更新,直到达到指定的迭代次数或指定的精度后,停止训练;
S222:将S215中获得的测试集中的合成遮挡谷粒图像输入至谷粒修复网络模型中得到复原后的完整谷粒图像,再将完整谷粒图像与测试集中的真实完整谷粒图像进行对比测试;
选择谷粒的粒长参数、粒宽参数、面积参数、周长参数4个粒型参数进行对比测试,计算4个粒型参数修复值与真实值的平均绝对百分比误差及相关系数R²。
6.根据权利要求5所述的一体化水稻穗部谷粒粒型性状的提取方法,其特征在于,S3包括如下子步骤:
S31:将单个谷粒图像进行灰度化,提取谷粒与背景之间的对比度最大的r通道的灰度图用于后续的形状提取;
S32:使用阈值分割算法对谷粒灰度图像进行二值化,获取谷粒二值图;
S33:基于谷粒二值图提取谷粒外轮廓,包括轮廓上各点的坐标;
S34:基于S33得到的轮廓坐标计算得到谷粒长度、宽度、周长及谷粒投影面积。
7.根据权利要求6所述的一体化水稻穗部谷粒粒型性状的提取方法,其特征在于,S34中谷粒长度、宽度、周长及谷粒投影面积的计算方法如下:
谷粒长度:基于谷粒轮廓上的点,求取任意两点之间的距离,取最大距离作为谷粒长度;
谷粒宽度:由求取谷粒长度的两点计算两点之间的直线方程,求取这个直线的法线斜率,随后求取不同截距下法线与轮廓的交点的距离,最大值即为谷粒宽度;
谷粒周长:基于谷粒轮廓上的点,若两点在上下或左右相邻,则两点间的距离定义为1像素大小;若两点在左上、左下、右上或右下方向相邻,则两点之间的距离定义为像素大小;将轮廓上所有相邻像素点之间的距离求和即为谷粒周长;
谷粒面积:基于谷粒轮廓上的点,谷粒面积即为谷粒轮廓内像素的数量。
8.根据权利要求7所述的一体化水稻穗部谷粒粒型性状的提取方法,其特征在于,S11包括如下子步骤:
S111:稻穗图像采集:在成熟水稻收割后获得水稻主穗,将稻穗随意平铺在可见光扫描成像设备上,对稻穗样本进行成像及存储;
S112:原始图像裁剪:基于S111采集的稻穗图像,将稻穗图像中的冗余部分去除,并裁剪为统一像素大小的图像;
S113:穗部谷粒图像标注:基于S112中得到的裁剪图像,采用边界框标注方式通过软件对所有穗部谷粒进行标注,每个谷粒标注后将边界框的左上角坐标、右下角坐标、边界框面积和框内目标种类记录下来,再将整株稻穗上的谷粒标注完后进行存储;
S114:数据增强:对由S112和S113中获取的稻穗图像及标注文件以一定的概率随机进行水平翻转或垂直翻转,数据集中的每个样本包括裁剪后的稻穗图像及相对应的标注文件;
S115:数据集划分:在经过S114中的数据增强后,将数据集中的样本以2:1:1的比例划分为训练集、验证集和测试集;
S116:数据集格式转换:分别将由S115划分后的数据集样本转换为用于深度学习目标检测方法训练所规定的格式。
9.根据权利要求8所述的一体化水稻穗部谷粒粒型性状的提取方法,其特征在于,S12包括如下子步骤:
S121:网络训练:将由S116获取的稻穗图像训练集输入到谷粒检测网络模型中进行训练,采用不同的超参数对网络进行训练,训练过程中不断更新神经网络的参数,直到达到指定的迭代次数,停止训练;
S122:网络验证:将由S116获取的稻穗图像验证集输入到谷粒检测网络模型中进行验证,根据目标检测评价指标获取不同超参数下网络的精度,随之获取最优谷粒检测网络;
S123:网络测试:将由S116获取的稻穗图像测试集输入到由S122获取的最优谷粒检测网络中进行测试,得到每个水稻稻穗样本上的谷粒数,计算网络预测值与真实值的平均绝对百分比误差及相关系数R²。
10.根据权利要求9所述的一体化水稻穗部谷粒粒型性状的提取方法,其特征在于,S12中所述的特征提取网络包括:输入层,用于输入采集到的稻穗图像;卷积层,用于提取图像中的信息;池化层,用于对由卷积层提取的信息进行挑选,对信息降维;激活层,用于提高网络的非线性拟合能力;输出层,用于输出经过多个卷积池化激活后富含图像高级语义信息的特征图;
所述特征金字塔网络包括:卷积层,用于对每个阶段输出的特征图进行降维;上采样层:用于对网络得到的特征图进行上采样;
所述区域提案网络包括:滑动窗口,在由特征金字塔网络得到的特征图上进行滑动生成不同大小的锚框用于目标的预测;分类通道,由若干卷积层构成,用于对锚框中的目标进行分类;回归通道,由若干卷积层构成,用于对锚框的位置进行精调;非极大值抑制,用于消除重复检测的锚框及超出图像边界的锚框;
所述分类回归网络包括:分类通道,由若干卷积层构成,用于对由区域提案网络生成的提案区域中的目标进行分类;回归通道,由若干卷积层构成,用于对由区域提案网络生成的提案区域的位置进行精调;非极大值抑制,用于消除冗余的检测框,并输出最终目标检测的结果。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202310900787.6A CN116612178A (zh) | 2023-07-21 | 2023-07-21 | 一种一体化水稻穗部谷粒粒型性状的提取方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202310900787.6A CN116612178A (zh) | 2023-07-21 | 2023-07-21 | 一种一体化水稻穗部谷粒粒型性状的提取方法 |
Publications (1)
Publication Number | Publication Date |
---|---|
CN116612178A true CN116612178A (zh) | 2023-08-18 |
Family
ID=87682293
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202310900787.6A Pending CN116612178A (zh) | 2023-07-21 | 2023-07-21 | 一种一体化水稻穗部谷粒粒型性状的提取方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN116612178A (zh) |
Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US8779233B1 (en) * | 2010-07-14 | 2014-07-15 | Iowa State University Research Foundation, Inc. | QTL regulating ear productivity traits in maize |
CN207408272U (zh) * | 2017-08-10 | 2018-05-25 | 武汉谷丰光电科技有限公司 | 基于线阵列相机和x射线双模式成像的水稻粒型参数测量装置 |
CN110969654A (zh) * | 2018-09-29 | 2020-04-07 | 北京瑞智稷数科技有限公司 | 基于收割机的玉米高通量表型测量的方法及装置、收割机 |
US20220121919A1 (en) * | 2020-10-16 | 2022-04-21 | X Development Llc | Determining Cereal Grain Crop Yield Based On Cereal Grain Trait Value(s) |
-
2023
- 2023-07-21 CN CN202310900787.6A patent/CN116612178A/zh active Pending
Patent Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US8779233B1 (en) * | 2010-07-14 | 2014-07-15 | Iowa State University Research Foundation, Inc. | QTL regulating ear productivity traits in maize |
CN207408272U (zh) * | 2017-08-10 | 2018-05-25 | 武汉谷丰光电科技有限公司 | 基于线阵列相机和x射线双模式成像的水稻粒型参数测量装置 |
CN110969654A (zh) * | 2018-09-29 | 2020-04-07 | 北京瑞智稷数科技有限公司 | 基于收割机的玉米高通量表型测量的方法及装置、收割机 |
US20220121919A1 (en) * | 2020-10-16 | 2022-04-21 | X Development Llc | Determining Cereal Grain Crop Yield Based On Cereal Grain Trait Value(s) |
Non-Patent Citations (2)
Title |
---|
LEJUN YU ET AL.: "An integrated rice panicle phenotyping method based on X-ray and RGB scanning and deep learning", 《THE CROP JOURNAL》, pages 42 - 56 * |
余乐俊: "水稻稻穗性状提取关键技术研究", 《中国博士学位论文全文数据库农业科技辑》, pages 047 - 17 * |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN111986099B (zh) | 基于融合残差修正的卷积神经网络的耕地监测方法及*** | |
CN113392775B (zh) | 一种基于深度神经网络的甘蔗幼苗自动识别与计数方法 | |
CN106610969A (zh) | 基于多模态信息的视频内容审查***及方法 | |
Zharkova et al. | Solar feature catalogues in EGSO | |
CN111259925B (zh) | 基于k均值聚类和宽度突变算法的田间麦穗计数方法 | |
CN111160096A (zh) | 禽蛋异常的识别方法、装置及***、存储介质、电子装置 | |
CN113269224B (zh) | 一种场景图像分类方法、***及存储介质 | |
CN112435214B (zh) | 基于先验框线性缩放的花粉检测方法、装置和电子设备 | |
Gaillard et al. | Voxel carving‐based 3D reconstruction of sorghum identifies genetic determinants of light interception efficiency | |
CN111079999A (zh) | 一种基于cnn和svm的洪水灾害易发性预测方法 | |
CN115272828A (zh) | 一种基于注意力机制的密集目标检测模型训练方法 | |
CN115439654B (zh) | 动态约束下的弱监督农田地块精细化分割方法及*** | |
CN116092179A (zh) | 一种改进的yolox跌倒检测*** | |
CN115797216B (zh) | 一种基于自编码网络的碑文文字修复模型及修复方法 | |
CN115861823A (zh) | 一种基于自监督深度学习的遥感变化检测方法和装置 | |
CN115239672A (zh) | 缺陷检测方法及装置、设备、存储介质 | |
CN116758539A (zh) | 一种基于数据增强的胚胎图像卵裂球识别方法 | |
Xia et al. | A Deep Learning Application for Building Damage Assessment Using Ultra-High-Resolution Remote Sensing Imagery in Turkey Earthquake | |
CN117727027A (zh) | 一种豆荚分类方法 | |
CN116805415A (zh) | 一种基于轻量化改进YOLOv5的笼养肉鸡健康状态识别方法 | |
CN116612178A (zh) | 一种一体化水稻穗部谷粒粒型性状的提取方法 | |
CN114626445B (zh) | 基于光流网络与高斯背景建模的大坝白蚁视频识别方法 | |
CN112287787B (zh) | 一种基于梯度直方图特征的作物倒伏分级方法 | |
CN114387517A (zh) | 一种基于高分辨率遥感影像的温室大棚智能提取方法 | |
CN114821165A (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 | ||
RJ01 | Rejection of invention patent application after publication |
Application publication date: 20230818 |