CN109285196B - 一种花毛茛花色变化特征确定方法及系统 - Google Patents

一种花毛茛花色变化特征确定方法及系统 Download PDF

Info

Publication number
CN109285196B
CN109285196B CN201811294109.5A CN201811294109A CN109285196B CN 109285196 B CN109285196 B CN 109285196B CN 201811294109 A CN201811294109 A CN 201811294109A CN 109285196 B CN109285196 B CN 109285196B
Authority
CN
China
Prior art keywords
value
flower color
change
color
day
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
CN201811294109.5A
Other languages
English (en)
Other versions
CN109285196A (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.)
INSTITUTE OF QUALITY STANDARD AND DETECTION TECHNOLOGY YUNNAN ACADEMY OF AGRICULTURAL SCIENCES
Original Assignee
INSTITUTE OF QUALITY STANDARD AND DETECTION TECHNOLOGY YUNNAN ACADEMY OF AGRICULTURAL SCIENCES
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 INSTITUTE OF QUALITY STANDARD AND DETECTION TECHNOLOGY YUNNAN ACADEMY OF AGRICULTURAL SCIENCES filed Critical INSTITUTE OF QUALITY STANDARD AND DETECTION TECHNOLOGY YUNNAN ACADEMY OF AGRICULTURAL SCIENCES
Priority to CN201811294109.5A priority Critical patent/CN109285196B/zh
Publication of CN109285196A publication Critical patent/CN109285196A/zh
Application granted granted Critical
Publication of CN109285196B publication Critical patent/CN109285196B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T7/00—Image analysis
    • G06T7/90—Determination of colour characteristics
    • 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
    • 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/30181—Earth observation
    • G06T2207/30188—Vegetation; Agriculture

Landscapes

  • Physics & Mathematics (AREA)
  • Engineering & Computer Science (AREA)
  • General Physics & Mathematics (AREA)
  • Chemical & Material Sciences (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • Analytical Chemistry (AREA)
  • Biochemistry (AREA)
  • General Health & Medical Sciences (AREA)
  • Health & Medical Sciences (AREA)
  • Immunology (AREA)
  • Pathology (AREA)
  • Spectroscopy & Molecular Physics (AREA)
  • Computer Vision & Pattern Recognition (AREA)
  • Theoretical Computer Science (AREA)
  • Spectrometry And Color Measurement (AREA)
  • Measurement Of The Respiration, Hearing Ability, Form, And Blood Characteristics Of Living Organisms (AREA)

Abstract

本发明公开一种花毛茛花色变化特征确定方法及系统。该方法包括：获取连续多天测量的花毛茛的花色测量值；根据花色测量值计算每种花色测量值对应的花色变化速度；根据花色变化速度计算每种花色测量值对应的花色变化加速度绝对值；根据花色变化速度计算每种花色测量值在每连续两天内的速度变化状态值；根据花色变化加速度绝对值计算每种花色测量值在每天的变化速度趋势值；根据速度变化状态值和速度变化趋势值计算每种花色测量值在每天的变化特征值；针对每种花色测量值，将每天的变化特征值求和得到每种花色测量值的变化动态值。本发明的方法及系统，能够动态地说明花毛茛花色整体的变化状态和趋势，从而科学准确的确定花毛茛花色变化特征。

Description

一种花毛茛花色变化特征确定方法及系统
技术领域
本发明涉及生物技术领域，特别是涉及一种花毛茛花色变化特征确定方法及系统。
背景技术
花毛茛隶属于毛茛科毛茛属，是近几年发展迅猛的盆栽、切花和花坛使用花卉。我国花毛茛年销售切花量达2500-4000万枝，利润在45％以上。在欧洲、日本等发达地区，花毛茛冬春季节产量及产值已达到切花市场的第二位，占据国际花卉拍卖市场重要地位，具很大研究和应用价值。花毛茛花色丰富，除常见色系外，还具有绿色、棕色、紫蓝色等罕见色系；同一色系存在浅、中、深过渡颜色。绝大多数花毛茛品种在成花不同阶段花瓣颜色上存在明度、彩度、甚至色调的变化。因此，如何科学准确的描述花毛茛花色变化特征，是获得有效的品种描述，评估花色市场经济价值的突破口，是应用具观赏优势的花色变化，控制不利的花色变化，花色育种改良，花色测试评价等的数据支撑。
目前在对花色表型变化研究方面，主要通过色差仪测量值L*(亮度)、a*(红绿度)、b*(黄蓝度)反映颜色状态，通过在测量值基础上计算得到的C*(彩度：
Figure BDA0001850706420000011
)反映颜色饱和度，通过色差值
Figure BDA0001850706420000012
反映颜色变化。然而，由ΔE*反映的颜色变化，并不能评价颜色变化状态、趋势以及变化速度等。通过各时期花色L*、a*、b*、C*测量值增减反映的花色变化也不能动态地说明整体的变化状态、趋势和速度特征，无法应用于花毛茛花色变化特征的科学准确描述。
发明内容
本发明的目的是提供一种花毛茛花色变化特征确定方法及系统，能够动态地说明毛茛花色整体的变化状态、趋势和速度特征，从而科学准确的确定花毛茛花色变化特征。
为实现上述目的，本发明提供了如下方案：
一种花毛茛花色变化特征确定方法，包括：
获取色差仪在连续多天测量的花毛茛的花色测量值；
根据所述花色测量值计算每种花色测量值对应的花色变化速度；
根据每种花色测量值对应的花色变化速度计算每种花色测量值对应的花色变化加速度绝对值；
根据每种花色测量值对应的花色变化速度计算每种花色测量值在每连续两天内的速度变化状态值；
根据每种花色测量值对应的花色变化加速度绝对值计算每种花色测量值在每天的变化速度趋势值；
根据每连续两天内的速度变化状态值和每天的速度变化趋势值计算每种花色测量值在每天的变化特征值；
针对每种花色测量值，将每天的变化特征值求和得到每种花色测量值的变化动态值；所述变化动态值用于确定花色变化的整体变化特征。
可选的，所述获取色差仪在连续连续多天测量的花毛茛的花色测量值，具体包括：
获取色差仪在连续多天测量的花毛茛的花色亮度、花色红绿度和花色黄蓝度；
根据所述花色红绿度和所述花色黄蓝度计算彩度；
根据所述花色亮度、所述花色红绿度和所述花色黄蓝度计算色差值。
可选的，所述根据所述花色测量值计算每种花色测量值对应的花色变化速度，具体包括：
利用公式
Figure BDA0001850706420000021
计算花色变化速度；其中j表示天数，vj表示第j天的花色变化速度；tj表示第j天，tj+1表示第j+1天，lj表示第j天的花色测量值，lj+1表示第j+1天的花色测量值。
可选的，所述根据每种花色测量值对应的花色变化速度计算每种花色测量值对应的花色变化加速度绝对值，具体包括：
利用公式
Figure BDA0001850706420000031
计算花色变化加速度绝对值；其中vj+1表示第j+1天的花色变化速度。
可选的，所述根据每种花色测量值对应的花色变化速度计算每种花色测量值在每连续两天内的速度变化状态值，具体包括：
利用公式
Figure BDA0001850706420000032
计算每种花色测量值在每连续两天内的速度变化状态值；其中sv(tj,tj+1)表示第j天到第j+1天期间的速度变化状态值，t为时间。
可选的，所述根据每种花色测量值对应的花色变化加速度绝对值计算每种花色测量值在每天的变化速度趋势值，具体包括：
利用公式
Figure BDA0001850706420000033
计算每种花色测量值在每天的变化速度趋势值；其中λ(aj)表示第j天的变化速度趋势值，ε为常数。
可选的，所述根据每连续两天内的速度变化状态值和每天的速度变化趋势值计算每种花色测量值在每天的变化特征值，具体包括：
利用公式Yjv(tj,tj+1)＝sv(tj,tj+1)·λ(aj)计算每种花色测量值在每天的变化特征值；其中Yjv(tj,tj+1)表示第j天的变化特征值。
可选的，所述针对每种花色测量值，将每天的变化特征值求和得到每种花色测量值的变化动态值，具体包括：
利用公式
Figure BDA0001850706420000034
将每天的变化特征值求和得到每种花色测量值的变化动态值；其中
Figure BDA0001850706420000035
表示连续多天的变化动态值。
可选的，在所述针对每种花色测量值，将每天的变化特征值求和得到每种花色测量值的变化动态值之后，还包括：
求取每种花色测量值对应的变化动态值的绝对值；
将各个变化动态值的绝对值进行对比，确定绝对值最大的变化动态值；
将绝对值最大的变化动态值对应的花色测量值类型作为花色变化主导因子。
本发明还公开一种花毛茛花色变化特征确定系统，包括：
花色测量值获取模块，用于获取色差仪在连续多天测量的花毛茛的花色测量值；
变化速度计算模块，用于根据所述花色测量值计算每种花色测量值对应的花色变化速度；
加速度绝对值计算模块，用于根据每种花色测量值对应的花色变化速度计算每种花色测量值对应的花色变化加速度绝对值；
速度变化状态值计算模块，用于根据每种花色测量值对应的花色变化速度计算每种花色测量值在每连续两天内的速度变化状态值；
变化速度趋势值计算模块，用于根据每种花色测量值对应的花色变化加速度绝对值计算每种花色测量值在每天的变化速度趋势值；
变化特征值计算模块，用于根据每连续两天内的速度变化状态值和每天的速度变化趋势值计算每种花色测量值在每天的变化特征值；
变化动态值计算模块，用于针对每种花色测量值，将每天的变化特征值求和得到每种花色测量值的变化动态值；所述变化动态值用于确定花色变化的整体变化特征。
根据本发明提供的具体实施例，本发明公开了以下技术效果：本发明所公开的花毛茛花色变化特征确定方法及系统，根据每种花色测量值对应的花色变化速度计算每种花色测量值对应的花色变化加速度绝对值；根据每种花色测量值对应的花色变化速度计算每种花色测量值在每连续两天内的速度变化状态值；根据每种花色测量值对应的花色变化加速度绝对值计算每种花色测量值在每天的变化速度趋势值；根据每连续两天内的速度变化状态值和每天的速度变化趋势值计算每种花色测量值在每天的变化特征值；针对每种花色测量值，将每天的变化特征值求和得到每种花色测量值的变化动态值；所述变化动态值用于确定花色变化的整体变化特征。通过计算连续多个时间段的花色变化速度以及花色变化的加速度绝对值确定花色变化的整体变化特征，能够克服传统花色变化确定方法缺少整体性的的缺点，能够动态地说明毛茛花色整体的变化状态、趋势和速度特征，从而科学准确的确定花毛茛花色变化特征。
附图说明
为了更清楚地说明本发明实施例或现有技术中的技术方案，下面将对实施例中所需要使用的附图作简单地介绍，显而易见地，下面描述中的附图仅仅是本发明的一些实施例，对于本领域普通技术人员来讲，在不付出创造性劳动性的前提下，还可以根据这些附图获得其他的附图。
图1为本发明实施例一的花毛茛花色变化特征确定方法的方法流程图；
图2为本发明实施例四的花毛茛花色变化特征确定系统的系统结构图。
具体实施方式
下面将结合本发明实施例中的附图，对本发明实施例中的技术方案进行清楚、完整地描述，显然，所描述的实施例仅仅是本发明一部分实施例，而不是全部的实施例。基于本发明中的实施例，本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其他实施例，都属于本发明保护的范围。
本发明的目的是提供一种花毛茛花色变化特征确定方法及系统，能够动态地说明毛茛花色整体的变化状态、趋势和速度特征，从而科学准确的确定花毛茛花色变化特征。
为使本发明的上述目的、特征和优点能够更加明显易懂，下面结合附图和具体实施方式对本发明作进一步详细的说明。
实施例一：
图1为本发明实施例一的花毛茛花色变化特征确定方法的方法流程图
参见图1，该花毛茛花色变化特征确定方法，包括：
一种花毛茛花色变化特征确定方法，包括：
步骤101：获取色差仪在连续多天测量的花毛茛的花色测量值。
该步骤101具体包括：
获取色差仪在连续多天测量的花毛茛的花色亮度L*、花色红绿度a*和花色黄蓝度b*；根据所述花色红绿度a*和所述花色黄蓝度b*计算彩度
Figure BDA0001850706420000051
根据所述花色亮度L*、所述花色红绿度a*和所述花色黄蓝度b*计算色差值
Figure BDA0001850706420000061
步骤102：根据所述花色测量值计算每种花色测量值对应的花色变化速度。具体为：
利用公式
Figure BDA0001850706420000062
计算花色变化速度；其中j表示天数，vj表示第j天的花色变化速度；tj表示第j天，tj+1表示第j+1天，lj表示第j天的花色测量值，lj+1表示第j+1天的花色测量值。
步骤103：根据每种花色测量值对应的花色变化速度计算每种花色测量值对应的花色变化加速度绝对值。具体为：
利用公式
Figure BDA0001850706420000063
计算花色变化加速度绝对值；其中vj+1表示第j+1天的花色变化速度。
本发明采用加速度绝对值而非原始值的原因在于：加速度值的正负与花色好坏无关，但绝对值大小与花色变化加速或减速有关。
步骤104：根据每种花色测量值对应的花色变化速度计算每种花色测量值在每连续两天内的速度变化状态值。具体为：
利用公式
Figure BDA0001850706420000064
计算每种花色测量值在每连续两天内的速度变化状态值；其中sv(tj,tj+1)表示第j天到第j+1天期间的速度变化状态值，t为时间。
步骤105：根据每种花色测量值对应的花色变化加速度绝对值计算每种花色测量值在每天的变化速度趋势值。具体为：
利用公式
Figure BDA0001850706420000065
计算每种花色测量值在每天的变化速度趋势值；其中λ(aj)表示第j天的变化速度趋势值，ε为常数。
步骤106：根据每连续两天内的速度变化状态值和每天的速度变化趋势值计算每种花色测量值在每天的变化特征值。具体为：
利用公式Yjv(tj,tj+1)＝sv(tj,tj+1)·λ(aj)计算每种花色测量值在每天的变化特征值；其中Yjv(tj,tj+1)表示第j天的变化特征值。
步骤107：针对每种花色测量值，将每天的变化特征值求和得到每种花色测量值的变化动态值；所述变化动态值用于确定花色变化的整体变化特征。具体为：
利用公式
Figure BDA0001850706420000071
将每天的变化特征值求和得到每种花色测量值的变化动态值；其中
Figure BDA0001850706420000072
表示连续多天的变化动态值。
步骤108：确定花色变化的主导因子。具体为：
求取每种花色测量值对应的变化动态值的绝对值；将各个变化动态值的绝对值进行对比，确定绝对值最大的变化动态值；将绝对值最大的变化动态值对应的花色测量值类型作为花色变化主导因子。
作为一种可选的实施方式，分别计算花色亮度L*、花色红绿度a*和花色黄蓝度b*三个花色参数对应的变化动态值的绝对值，选择绝对值较大的变化动态值对应的花色参数为花色变化主导因子。
作为一种可选的实施方式，分别计算花色亮度L*、花色红绿度a*和花色黄蓝度b*、彩度C*和色差值E*五个花色参数对应的变化动态值的绝对值，选择绝对值较大的变化动态值对应的花色参数为花色变化主导因子。
主导因子的确定可以明确花色变化的主要影响因素，从而明确花色变化方向，对花卉育种，销售过程中充分利用具有优势的花色变化，减小不利变化具有重要作用。
实施例二：
选择花毛茛品种“娇艳洋牡丹砖红”主茎顶花，每天在固定时间，使用色差仪对待测区域颜色进行测定，记录L*，a*，b*值，测定3次，取平均值。本实例中，测定时间为9:30am-10:00am，测定时长为外轮花瓣展开至外轮花瓣松动，测定区域为外轮花瓣上表面上三分之一处中心位置。在L*，a*，b*值基础上，根据
Figure BDA0001850706420000073
和
Figure BDA0001850706420000074
计算获得彩度C*和色差值E*；根据
Figure BDA0001850706420000075
其中Yjv(tj,tj+1)＝sv(tj,tj+1)·λ(aj)，
Figure BDA0001850706420000076
Figure BDA0001850706420000077
计算获得花色的变化动态值。表1为实施例二的花色变化速度动态值。
表1
Figure BDA0001850706420000081
参见表1，结果表明，花色在整个测定时期总体亮度变大(L*值为正，1.2)，红度减小(a*值为负，-3.1)，黄度变化不大，微减小(b*值为负，-0.2)，彩度减小(C*值为负，-2.5)，色差值减小(E*值为负，-1.0)；红绿度变化速度|-3.1|>亮度变化速度|1.2|>黄蓝度变化速度|-0.2|，即红绿度快速减小是测定时期内影响花色变化的主要因子。以此类推，各时期变化特征同样可通过数据量化说明。此外，应用英国皇家园艺学会比色卡(RHSCC)测定颜色，结果表明，连续九天色号对应颜色值饱和度减小，亮度变大，证明本发明方法科学有效。
实施例三：
选择花毛茛品种“娇艳洋牡丹洋红”主茎顶花，每天在固定时间，使用色差仪对待测区域颜色进行测定，记录L*，a*，b*值，测定3次，取平均值。本实例中，测定时间为9:30am-10:00am，测定时长为外轮花瓣展开至外轮花瓣松动，测定区域为外轮花瓣上表面上三分之一处中心位置。在L*，a*，b*值基础上，根据
Figure BDA0001850706420000082
和
Figure BDA0001850706420000083
计算获得彩度C*和色差值E*；根据
Figure BDA0001850706420000091
其中Yjv(tj,tj+1)＝sv(tj,tj+1)·λ(aj)，
Figure BDA0001850706420000092
Figure BDA0001850706420000093
计算获得花色的变化动态值。表2为实施例三的花色变化速度动态值。
表2
Figure BDA0001850706420000094
参见表2，结果表明，花色在整个测定时期总体亮度基本没变，微减小(L*值为负，-0.4)，红度减小(a*值为负，-1.3)，黄度变化不大，微减小(b*值为负，-0.3)，彩度减小(C*值为负，-1.3)，颜色值微减小(E*值为负，-0.7)，表明整体颜色变化不大；红绿度变化速度|-1.3|，即红度减小是测定时期内影响花色变化的主要因子。以此类推，各时期变化特征同样可通过数据量化说明。此外，应用英国皇家园艺学会比色卡(RHSCC)测定颜色，结果表明，连续九天色号对应颜色值饱和度减小，亮度变大，证明本发明方法科学有效。
实施例四：
图2为本发明花毛茛花色变化特征确定系统实施例的系统结构图。
参见图2，该花毛茛花色变化特征确定系统，包括：
本发明还公开一种花毛茛花色变化特征确定系统，包括：
花色测量值获取模块201，用于获取色差仪在连续多天测量的花毛茛的花色测量值。具体为：
获取色差仪在连续多天测量的花毛茛的花色亮度L*、花色红绿度a*和花色黄蓝度b*；根据所述花色红绿度a*和所述花色黄蓝度b*计算彩度
Figure BDA0001850706420000101
根据所述花色亮度L*、所述花色红绿度a*和所述花色黄蓝度b*计算色差值
Figure BDA0001850706420000102
变化速度计算模块202，用于根据所述花色测量值计算每种花色测量值对应的花色变化速度；具体为：
利用公式
Figure BDA0001850706420000103
计算花色变化速度；其中j表示天数，vj表示第j天的花色变化速度；tj表示第j天，tj+1表示第j+1天，lj表示第j天的花色测量值，lj+1表示第j+1天的花色测量值。
加速度绝对值计算模块203，用于根据每种花色测量值对应的花色变化速度计算每种花色测量值对应的花色变化加速度绝对值；具体为：
利用公式
Figure BDA0001850706420000104
计算花色变化加速度绝对值；其中vj+1表示第j+1天的花色变化速度。
本发明采用加速度绝对值而非原始值的原因在于：加速度值的正负与花色好坏无关，但绝对值大小与花色变化加速或减速有关。
速度变化状态值计算模块204，用于根据每种花色测量值对应的花色变化速度计算每种花色测量值在每连续两天内的速度变化状态值；具体为：
利用公式
Figure BDA0001850706420000105
计算每种花色测量值在每连续两天内的速度变化状态值；其中sv(tj,tj+1)表示第j天到第j+1天期间的速度变化状态值，t为时间。
变化速度趋势值计算模块205，用于根据每种花色测量值对应的花色变化加速度绝对值计算每种花色测量值在每天的变化速度趋势值；具体为：
利用公式
Figure BDA0001850706420000111
计算每种花色测量值在每天的变化速度趋势值；其中λ(aj)表示第j天的变化速度趋势值，ε为常数。
变化特征值计算模块206，用于根据每连续两天内的速度变化状态值和每天的速度变化趋势值计算每种花色测量值在每天的变化特征值；具体为：
利用公式Yjv(tj,tj+1)＝sv(tj,tj+1)·λ(aj)计算每种花色测量值在每天的变化特征值；其中Yjv(tj,tj+1)表示第j天的变化特征值。
变化动态值计算模块207，用于针对每种花色测量值，将每天的变化特征值求和得到每种花色测量值的变化动态值；所述变化动态值用于确定花色变化的整体变化特征。具体为：
利用公式
Figure BDA0001850706420000112
将每天的变化特征值求和得到每种花色测量值的变化动态值；其中
Figure BDA0001850706420000113
表示连续多天的变化动态值。
主导因子确定模块208，用于确定花色变化的主导因子。具体为：
求取每种花色测量值对应的变化动态值的绝对值；将各个变化动态值的绝对值进行对比，确定绝对值最大的变化动态值；将绝对值最大的变化动态值对应的花色测量值类型作为花色变化主导因子。
根据本发明提供的具体实施例，本发明公开了以下技术效果：本发明所公开的花毛茛花色变化特征确定方法及系统，根据每种花色测量值对应的花色变化速度计算每种花色测量值对应的花色变化加速度绝对值；根据每种花色测量值对应的花色变化速度计算每种花色测量值在每连续两天内的速度变化状态值；根据每种花色测量值对应的花色变化加速度绝对值计算每种花色测量值在每天的变化速度趋势值；根据每连续两天内的速度变化状态值和每天的速度变化趋势值计算每种花色测量值在每天的变化特征值；针对每种花色测量值，将每天的变化特征值求和得到每种花色测量值的变化动态值；所述变化动态值用于确定花色变化的整体变化特征。通过计算连续多个时间段的花色变化速度以及花色变化的加速度绝对值确定花色变化的整体变化特征，能够动态地说明毛茛花色整体的变化状态、趋势和速度特征，从而科学准确的确定花毛茛花色变化特征。
对于实施例公开的系统而言，由于其与实施例公开的方法相对应，所以描述的比较简单，相关之处参见方法部分说明即可。
本文中应用了具体个例对本发明的原理及实施方式进行了阐述，以上实施例的说明只是用于帮助理解本发明的方法及其核心思想；同时，对于本领域的一般技术人员，依据本发明的思想，在具体实施方式及应用范围上均会有改变之处。综上所述，本说明书内容不应理解为对本发明的限制。

Claims (4)

1.一种花毛茛花色变化特征确定方法，其特征在于，包括：
获取色差仪在连续多天测量的花毛茛的花色测量值；
根据所述花色测量值计算每种花色测量值对应的花色变化速度；
根据每种花色测量值对应的花色变化速度计算每种花色测量值对应的花色变化加速度绝对值；
根据每种花色测量值对应的花色变化速度计算每种花色测量值在每连续两天内的速度变化状态值；
根据每种花色测量值对应的花色变化加速度绝对值计算每种花色测量值在每天的变化速度趋势值；
根据每连续两天内的速度变化状态值和每天的速度变化趋势值计算每种花色测量值在每天的变化特征值；
针对每种花色测量值，将每天的变化特征值求和得到每种花色测量值的变化动态值；所述变化动态值用于确定花色变化的整体变化特征；
所述根据所述花色测量值计算每种花色测量值对应的花色变化速度，具体包括：
利用公式
Figure FDA0002573395200000011
计算花色变化速度；其中j表示天数，vj表示第j天的花色变化速度；tj表示第j天，tj+1表示第j+1天，lj表示第j天的花色测量值，lj+1表示第j+1天的花色测量值；
所述根据每种花色测量值对应的花色变化速度计算每种花色测量值对应的花色变化加速度绝对值，具体包括：
利用公式
Figure FDA0002573395200000012
计算花色变化加速度绝对值；其中vj+1表示第j+1天的花色变化速度；
所述根据每种花色测量值对应的花色变化速度计算每种花色测量值在每连续两天内的速度变化状态值，具体包括：
利用公式
Figure FDA0002573395200000013
计算每种花色测量值在每连续两天内的速度变化状态值；其中sv(tj,tj+1)表示第j天到第j+1天期间的速度变化状态值，t为时间；
所述根据每种花色测量值对应的花色变化加速度绝对值计算每种花色测量值在每天的变化速度趋势值，具体包括：
利用公式
Figure FDA0002573395200000021
计算每种花色测量值在每天的变化速度趋势值；其中λ(aj)表示第j天的变化速度趋势值，ε为常数；
所述根据每连续两天内的速度变化状态值和每天的速度变化趋势值计算每种花色测量值在每天的变化特征值，具体包括：
利用公式Yjv(tj,tj+1)＝sv(tj,tj+1)·λ(aj)计算每种花色测量值在每天的变化特征值；其中Yjv(tj,tj+1)表示第j天的变化特征值；
所述针对每种花色测量值，将每天的变化特征值求和得到每种花色测量值的变化动态值，具体包括：
利用公式
Figure FDA0002573395200000022
将每天的变化特征值求和得到每种花色测量值的变化动态值；其中
Figure FDA0002573395200000023
表示连续多天的变化动态值。
2.根据权利要求1所述的一种花毛茛花色变化特征确定方法，其特征在于，所述获取色差仪在连续连续多天测量的花毛茛的花色测量值，具体包括：
获取色差仪在连续多天测量的花毛茛的花色亮度、花色红绿度和花色黄蓝度；
根据所述花色红绿度和所述花色黄蓝度计算彩度；
根据所述花色亮度、所述花色红绿度和所述花色黄蓝度计算色差值。
3.根据权利要求1所述的一种花毛茛花色变化特征确定方法，其特征在于，在所述针对每种花色测量值，将每天的变化特征值求和得到每种花色测量值的变化动态值之后，还包括：
求取每种花色测量值对应的变化动态值的绝对值；
将各个变化动态值的绝对值进行对比，确定绝对值最大的变化动态值；
将绝对值最大的变化动态值对应的花色测量值类型作为花色变化主导因子。
4.一种花毛茛花色变化特征确定系统，其特征在于，包括：
花色测量值获取模块，用于获取色差仪在连续多天测量的花毛茛的花色测量值，具体包括：
获取色差仪在连续多天测量的花毛茛的花色亮度L*、花色红绿度a*和花色黄蓝度b*；根据所述花色红绿度a*和所述花色黄蓝度b*计算彩度
Figure FDA0002573395200000031
根据所述花色亮度L*、所述花色红绿度a*和所述花色黄蓝度b*计算色差值
Figure FDA0002573395200000032
变化速度计算模块，用于根据所述花色测量值计算每种花色测量值对应的花色变化速度，具体包括：
利用公式
Figure FDA0002573395200000033
计算花色变化速度；其中j表示天数，vj表示第j天的花色变化速度；tj表示第j天，tj+1表示第j+1天，lj表示第j天的花色测量值，lj+1表示第j+1天的花色测量值；
加速度绝对值计算模块，用于根据每种花色测量值对应的花色变化速度计算每种花色测量值对应的花色变化加速度绝对值，具体包括：
利用公式
Figure FDA0002573395200000034
计算花色变化加速度绝对值；其中vj+1表示第j+1天的花色变化速度；
速度变化状态值计算模块，用于根据每种花色测量值对应的花色变化速度计算每种花色测量值在每连续两天内的速度变化状态值，具体包括：
利用公式
Figure FDA0002573395200000035
计算每种花色测量值在每连续两天内的速度变化状态值；其中sv(tj,tj+1)表示第j天到第j+1天期间的速度变化状态值，t为时间；
变化速度趋势值计算模块，用于根据每种花色测量值对应的花色变化加速度绝对值计算每种花色测量值在每天的变化速度趋势值，具体包括：
利用公式
Figure FDA0002573395200000036
计算每种花色测量值在每天的变化速度趋势值；其中λ(aj)表示第j天的变化速度趋势值，ε为常数；
变化特征值计算模块，用于根据每连续两天内的速度变化状态值和每天的速度变化趋势值计算每种花色测量值在每天的变化特征值，具体包括：
利用公式Yjv(tj,tj+1)＝sv(tj,tj+1)·λ(aj)计算每种花色测量值在每天的变化特征值；其中Yjv(tj,tj+1)表示第j天的变化特征值；
变化动态值计算模块，用于针对每种花色测量值，将每天的变化特征值求和得到每种花色测量值的变化动态值；所述变化动态值用于确定花色变化的整体变化特征，具体包括：
利用公式
Figure FDA0002573395200000041
将每天的变化特征值求和得到每种花色测量值的变化动态值；其中
Figure FDA0002573395200000042
表示连续多天的变化动态值。
CN201811294109.5A 2018-11-01 2018-11-01 一种花毛茛花色变化特征确定方法及系统 Active CN109285196B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201811294109.5A CN109285196B (zh) 2018-11-01 2018-11-01 一种花毛茛花色变化特征确定方法及系统

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201811294109.5A CN109285196B (zh) 2018-11-01 2018-11-01 一种花毛茛花色变化特征确定方法及系统

Publications (2)

Publication Number Publication Date
CN109285196A CN109285196A (zh) 2019-01-29
CN109285196B true CN109285196B (zh) 2020-11-13

Family

ID=65174742

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201811294109.5A Active CN109285196B (zh) 2018-11-01 2018-11-01 一种花毛茛花色变化特征确定方法及系统

Country Status (1)

Country Link
CN (1) CN109285196B (zh)

Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
WO2009016710A1 (ja) * 2007-07-30 2009-02-05 M & M Holdings Co., Ltd. 置物
CN103971011A (zh) * 2014-05-20 2014-08-06 中国农业科学院农业资源与农业区划研究所 全球陆地植被生物量变化监测方法
CN105890764A (zh) * 2016-03-31 2016-08-24 杭州市农业科学研究院 基于计算机图像处理的分辨植物花色深浅的方法
CN107067016A (zh) * 2016-11-10 2017-08-18 浙江农林大学 基于叶片颜色的枇杷属内植物分类方法

Patent Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
WO2009016710A1 (ja) * 2007-07-30 2009-02-05 M & M Holdings Co., Ltd. 置物
CN103971011A (zh) * 2014-05-20 2014-08-06 中国农业科学院农业资源与农业区划研究所 全球陆地植被生物量变化监测方法
CN105890764A (zh) * 2016-03-31 2016-08-24 杭州市农业科学研究院 基于计算机图像处理的分辨植物花色深浅的方法
CN107067016A (zh) * 2016-11-10 2017-08-18 浙江农林大学 基于叶片颜色的枇杷属内植物分类方法

Non-Patent Citations (5)

* Cited by examiner, † Cited by third party
Title
FLORAL COLOR CHANGE AND INSECT POLLINATION: A DYNAMIC RELATIONSHIP;Martha R. Weiss et al;《IsraelJournal of Plant Sciences》;20130419;第45卷;第185-199页 *
Flower color variation: A model for the experimental;Michael T. Clegg et al;《Proc Natl Acad Sci USA 》;20000620;第97卷(第13期);第7016-7023页 *
两个牡丹品种开花过程中花色变化的研究;杨琴等;《园艺学报》;20150531;第930-938页 *
月季花发育过程中花色变化的生理生化研究;程怡等;《西北植物学报》;20140430;第34卷(第4期);第733-739页 *
绚丽海棠实生后代叶色及花色变化研究;田涛等;《北方园艺》;20170731;第69-73页 *

Also Published As

Publication number Publication date
CN109285196A (zh) 2019-01-29

Similar Documents

Publication Publication Date Title
CN103565419B (zh) 测色装置、方法和适宜粉底的选择方法
EP1541960A3 (en) Parametric optimization of optical metrology model
CN104159053B (zh) 一种大色域显示设备三原色优化设置方法
CN109285196B (zh) 一种花毛茛花色变化特征确定方法及系统
CN104062010B (zh) 一种优化定标算法的分光光源颜色照度测量仪器
CN108734074A (zh) 指纹识别方法以及指纹识别装置
CN106017868B (zh) 基于亮度参数的多通道led照明系统的光谱匹配方法
CN109472210A (zh) 一种植物的识别方法、装置、储存介质和设备
CN103559712A (zh) 黑瓜子色选方法
CN109580926A (zh) 一种润州凤头白鸭日龄综合鉴定方法
CN109283145B (zh) 一种花毛茛花色测试时期确定方法及系统
ATE322521T1 (de) Herstellung von farbstoffmischungen zur gewinnung bestimmter zielfarben
CN115099992B (zh) 蘑菇种植房的数据处理方法及处理系统
CN106248591B (zh) 一种烟叶成熟度区分用器具及其区分方法
CN103575673A (zh) 可测量色彩情感的分光光度计
CN106061000B (zh) 基于光谱峰值辐亮度参数的多通道照明系统光谱匹配方法
CN102945555B (zh) 一种天然绿色蚕茧颜色的计算机分类方法
CN108956486B (zh) 一种染料评估方法
Pridmore Model of saturation and brightness: Relations with luminance
KR101435994B1 (ko) 후지 및 홍로 사과나무의 수피 발색에 의한 질소수준 판정방법
CN101551274A (zh) 一种新型棉花色泽仪及棉花色泽测量方法
CN104766301B (zh) 一种基于图像的单色采集算法
CN208333480U (zh) 一种基于田间精准定位仪
Qingfang Game Analysis System of Landscape Garden Plant Configuration Optimization based on Multi-View Image Comprehensive Evaluation Index Matrix Algorithm
CN110907036B (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