Automatic extraction device and experiment of rice seed investigation parameters
-
摘要:目的
传统考种方法测量精度和效率难以满足现代水稻育种研究的需求，设计一种水稻谷粒图像与质量信息同步采集装置，实现水稻谷粒考种参数的自动提取。
方法采用掩膜法自动提取稻谷区域图像，根据稻谷投影面积、数量规律获取稻谷总粒数；根据空粒、实粒颖壳轮廓差异识别空粒；基于角点间距均值标定法，结合轮廓最小外接矩形法获取粒长、粒宽，结合链码法获取粒周长；采用正方形面积均值标定法结合像素累加法获取粒面积。分析摄像头高度、谷粒数量、谷粒种类、规则图形类型对谷粒性状参数提取精度的影响。
结果摄像头高度对稻谷总数、空粒数、长、宽测量精度有明显影响，稻谷种类对宽度测量精度有明显影响，规则图形类型对周长和面积测量精度有明显影响。采用本文提出方法测量总粒数、空粒数、粒长、粒宽、粒周长、粒面积的决定系数(R2)分别为0.99830、0.98780、0.99610、0.78290、0.99510和0.99998，测量的平均精度分别为99.47%、87.17%、96.55%、96.36%、98.00%和95.86%，测量效率为16.52粒/s。
结论本文所采用的稻谷谷粒考种参数自动提取方法可行，可为全自动考种机的研发提供技术参考。
Abstract:ObjectiveThe measurement accuracy and efficiency of traditional seed investigation methods could not meet the needs of modern rice breeding research. A synchronous collection device of rice grain image and quality information was designed to automatically extract rice seed investigation parameters.
MethodThe image of grain region was automatically extracted by mask method, and the total number of rice grains was obtained according to the law between rice projection area and rice quantity. Empty grains were identified according to the difference of the hull contour between empty grains and full grains. Based on the mean value calibration method of corner spacing, the grain length and width were obtained by combining the minimum circumscribed rectangle method of contour, and the grain perimeter was obtained by combining the chain code method. The square area mean calibration method and pixel accumulation method were used to obtain the grain area. The effects of camera height, grain quantity, grain type and regular graph type on the extraction accuracy of grain character parameters were analyzed.
ResultThe camera height had a obvious impact on the measurement accuracies of total number, empty number, length and width of rice grains, the grain type had a obvious impact on the measurement accuracy of width, and the regular graph type had a obvious impact on the measurement accuracies of grain area and perimeter. The determination coefficients (R2) of total grain number, empty grain number, grain length, grain width, grain perimeter and grain area measured by the proposed method were 0.99830, 0.98780, 0.99610, 0.78290, 0.99510 and 0.99998 respectively, the average measurement accuracies were 99.47%, 87.17%, 96.55%, 96.36%, 98.00% and 95.86% respectively, and the measurement efficiency was 16.52 grains per second.
ConclusionThe automatic extraction method of rice seed investigation parameters used in this paper is feasible, and can provide technical references for the development of automatic seed investigation machine.
-
Keywords:
- Rice /
- Seed investigation /
- Parameter extraction /
- Image processing
-
桉树Eucalyptus spp.在中国南方的种植已经超过367万hm2[1]，因其显著的速生性和丰产性而被列为我国短周期工业原料林造林的首选树种之一，它用不到3%的林地面积解决了18%的工业用材需求，其良好的短期经济效益为其大面积的发展起到了积极的作用[2]。然而，由于桉树引种和推广种植的地区大多为台风频发的东南沿海地区，又因杂交无性系的推出使其品种繁多，其适生性和抗风性又各不同，造成明显产量差异，因此有必要对广泛推广种植的品种作系统的多点对比，为其造林品种的选择提供科学的依据。目前关于桉树适生性和抗性方面的研究多集中于桉树的抗病性[3-5]和抗寒性[6-7]，虽然前人亦有相关桉树的适生性研究[8]和抗风性对比试验[9-13]，但多为单一地区的研究，缺乏对多个立地、多个桉树无性系的系统比较，且研究的品种和立地情况也各有不同。本文对10个目前主要推广种植的桉树无性系，在广东和广西的6个不同立地的生长情况进行对比研究，以期对其抗性和适生性作出科学评价，为“适地适树”的选择提供科学参考。
1. 材料与方法
1.1 试验地概况
广东3个试验地位于雷州林业局的唐家林场（20°50′01″ N，109°51′76″ E）和遂溪林场（21°21′81″ N，110°01′66″ E），及湛江林业局属下的吴川林场（21°18′37″ N，110°31′13″ E）。唐家林场林地为玄武岩发育的厚有机质层赤红壤，土壤黏重，肥力较好；吴川林场林地为花岗岩发育的赤红壤，中壤土，肥力中等；遂溪林场林地为浅海沉积物发育的砂壤土，松散而少黏粒，肥力中下。广西3个试验地：钦廉林场（21°45′46″ N，109°03′08″E），浅海沉积物发育的中壤土，肥力中上；防城林地（21°48′02″N，108°25′16″E），沙页岩发育的赤红壤，土层较薄，中壤土，肥力中下；东门林场（22°28′N，107°51′E），石灰岩发育的赤红壤，土层深厚黏重，肥力较好，相对干旱。各试验地土壤分析结果见表 1。除广西东门林场试验地属于南亚热带湿润季风气候外，其余5个试验地均属于典型热带北缘海洋性季风气候，干湿季明显，夏季盛行东南风、南风及西南风，冬季多北风和东北风，每年6—11月为热带风暴和台风季节，这为各无性系的抗风性比较提供了条件。
表 1 各试验地土壤分析结果Table 1. The analysis results of soil in six experimental plots试验地 取样深度/cm pH w/(g·kg-1) w/(mg·kg-1) 土壤质地 土层深度/m 有机质 全氮 全磷 全钾 速效氮 速效磷 速效钾 有效硼 吴川 0~20 4.22 6.14 0.31 0.13 1.04 26.68 3.27 14.13 0.76 中壤 >2.0 20~40 4.27 4.84 0.26 0.12 1.11 19.47 1.85 11.79 0.66 唐家 0~20 4.31 7.15 0.43 0.11 0.75 48.60 4.29 116.11 1.12 黏土 >2.0 20~40 4.40 5.63 0.34 0.11 0.84 42.32 1.30 145.79 1.09 遂溪 0~20 4.19 5.41 0.20 0.17 0.75 22.51 2.75 16.51 0.39 轻壤 >2.0 20~40 4.26 3.99 0.14 0.17 0.77 16.13 1.99 12.52 0.31 钦廉 0~20 3.90 24.28 0.98 0.15 10.60 82.86 1.79 37.27 1.58 中壤 >1.5 20~40 4.08 19.57 0.77 0.15 11.49 66.11 1.07 80.43 0.91 防城 0~20 4.10 14.67 0.75 0.15 15.20 53.48 0.53 47.14 1.67 中壤 >0.8 20~40 4.15 9.67 0.63 0.15 17.28 39.60 0.41 29.60 1.54 东门 0~20 4.27 17.02 0.80 0.40 3.65 66.14 0.85 50.41 0.79 重壤 >1.5 20~40 4.32 13.97 0.69 0.35 3.81 54.37 0.53 35.45 0.66 1.2 供试材料
供试桉树无性系有来自广西东门林场的尾巨桉E. urophylla×E. grandis系列DH3226、DH3213、DH3327、DH3222和巨赤桉E. grandis× camaldulensis DH2012，尾细桉E.urophylla×E.tereticornis系列中有来自雷州林业局林业科学研究所的M1和来自热带林业研究所的LH9211及LH9224，天然杂交种有雷州林业局选育的石岭1号（SH1）和来自湛江的尾叶桉E. urophylla自然杂交种U6，合计10个无性系，全部使用扦插苗，苗高25~30 cm，挑选生长基本均匀一致的苗木，统一于2002年6月雨后造林。
1.3 试验设计和造林
6个试验地，全部采用随机区组设计，广东各试验地重复4次，广西各试验地重复3次。6月间定植，每小区10行每行20株，株行距2 m×3 m，小区面积1 200 m2。3种桉树基肥（S-芬兰的Kemira“星王”复合肥N15P15K15，R-热林所Ritf桉树专用配方肥N5P18K8B0.1，L-雷林1号桉树专用配方肥N10P10K10B0.1）处理10个无性系(各肥料名称中元素N、P、K、B及下标分别表示N、P2O5、K2O、B的质量百分比)，等质量施加3种桉树基肥，造林前统一每穴施桉树基肥500 g，造林后2个半月结合抚育除草，每株分别追施尿素50 g和氯化钾30 g。第2和第3年4—5月间抚育追施桉树专用肥N16P4K8B0.02 250 g·株-1，广东各试验地机械开沟小穴种植，广西防城试验地人工垦带后挖穴种植，其余2个试验地采用机耕全垦小穴造林。
1.4 调查指标及方法
试验的调查指标包括成活率、树高、胸径、发病率和风害率等，每小区从第3行第3株开始，每木进行检测，观测3行每行10株，即每小区有效观测30株，3次重复。成活率调查在植后1个月内，各试验地各处理调查3个重复，每重复观测30株。风害调查在台风“伊布都”和“黑格比”登陆广东后进行，其中心风力分别为11级和14级，台风过后4 d内完成。
风害率=受害株数/调查总株数×100%；风折率=折断株数/调查总株数×100%[9-11]。文中单株材积和蓄积量均取自当地采伐数据，由于年均出材量数据有所缺失，本文中年均出材量（MAI）=单株材积×种植密度（每公顷株数）/材积年份。
1.5 模糊判别的隶属函数
单一指标的排序结果无法恰当反映多指标的综合效应，所以采用模糊数学中隶属函数的方法对各桉树无性系的适生性进行综合评价[14]。用于分析的隶属函数R(Xi)计算方程为：
$$R(X_i)=(X_i-X_{\text{min}})/(X_{\text{max}}-X_{\text{min}})，$$ (1) $$R(X_i)=1-(X_i-X_{\text{min}})/(X_{\text{max}}-X_{\text{min}})，$$ (2) 式中，Xi为指标测定值，Xmin和Xmax分别为所有处理的某一指标的最小值和最大值。如果某一指标与生长表现为正相关，则采用式(1) 计算隶属值，反之则用式(2)。累加具体各指标的隶属值，并求出平均值进行比较。
1.6 数据统计与分析
所有试验数据采用Microsoft Excel 2010输入并整理绘图，使用SPSS19.0软件进行方差分析、相关分析和Duncan’s多重比较。
2. 结果与分析
2.1 成活率分析
在同一生长环境条件及相同栽培管理的情况下，造林成活率、抗逆性（如抗风、旱、寒、病虫害等）能力的大小，基本可以反映品系（树种/种源/无性系）对环境的适应程度及遗传差异[12]。造林成活率的高低是树种适生性的第一反应，调查结果（表 2）表明，10个桉树无性系的成活率平均为96.0%，LH9224的成活率(平均93.7%)相对较差，显著低于其余品系（95.2%~97.4%）；试验点间，黏重土壤（唐家）和降雨量相对较少的立地（东门）有相对较低的成活率，分别为93.9%和92.2%，低于平均值。
表 2 各试验地各桉树无性系成活率1）Table 2. The survival rates of ten Eucalyptus clones at six sites% 试验地 DH3226 DH3213 DH3327 DH3222 DH2012 SH1 M1 LH9211 LH9224 U6 平均值 唐家 92.8±0.9 94.2±1.1 91.7±1.2 93.6±1.4 93.9±1.5 93.9±1.2 97.8±1.2 91.4±1.2 93.6±1.6 96.7±1.1 93.9b 吴川 98.4±1.6 96.1±1.2 98.2±0.7 95.6±1.2 99.7±2.1 98.2±1.1 99.5±1.4 97.7±1.4 87.8±2.3 98.4±0.9 96.9a 遂溪 97.9±1.3 99.1±1.5 98.5±0.9 99.1±1.1 97.3±1.2 98.5±1.2 98.2±1.6 95.8±1.3 97.3±0.9 98.8±2.3 98.1a 钦廉 96.3±1.1 97.8±1.6 91.9±1.2 98.1±1.2 99.3±0.9 96.3±1.2 94.8±0.8 96.3±1.2 98.1±2.2 98.5±1.3 96.7a 防城 98.5±1.2 97.8±1.6 99.3±1.1 98.1±0.9 98.5±2.1 97.8±1.4 97.4±1.2 98.1±1.1 98.9±0.9 97.8±2.1 98.2a 东门 93.0±0.9 86.7±1.3 90.0±1.1 93.7±2.3 95.2±1.2 94.1±2.1 96.7±1.5 95.9±1.7 86.3±2.1 90.7±1.8 92.2b 平均值 96.1a 95.2a 94.9a 96.3a 97.3a 96.4a 97.4a 95.8a 93.7b 96.8a 96.0 1）表中同行或同列数据后，凡是有一个相同小写字母者，表示无性系或试验地间差异不显著（Duncan’s法，P>0.05）。 由表 2可知，桉树成活率的平均值在无性系间的差异小于4%，立地间的差异小于6%，说明绝大部分杂交桉树无性系仍然保持桉树适生性强的优良特性，不论在黏重的立地（唐家）还是沙质土地（遂溪），不论在降雨量大（防城年降雨量>2 000 mm）或小（东门年降雨<1 200 mm）的地方均有理想的成活率。
2.2 抗病性分析
对各试验地各无性系桉树的抗病性调查发现，当年生桉树以LH9224的发病率最高，主要为青枯病造成的死亡或黄株落叶，在唐家、吴川和钦廉其发病率均超过10%，6个试验地其平均发病率为6.7%；其次为M1无性系，其平均发病率为4.08%；以LH9211抗病性最强，在6个试验地1 620株有效调查株中，没有1株病例，其次为DH2012和SH1，发病率只为0.12%（图 1）。从立地类型上分析，桉树在黏重的土壤上发病率较高，同在广东雷州半岛地区的3个试验地，唐家土壤黏重，其平均发病率为3.63%，而中壤土的吴川林地发病率为1.62%，沙质土的遂溪林地发病率为0；东门林场发病率较低则可能是其降雨量明显减少所致（图 2）。
![图 1 10个桉树无性系发病率]()
图 1 10个桉树无性系发病率柱子上方，凡是有一个相同大写字母者表示无性系间差异不显著（Duncan’s法，P>0.01）。Figure 1. The disease incidences of ten Eucalyptus clones![图 2 6个试验地桉树无性系发病率]()
图 2 6个试验地桉树无性系发病率柱子上方，凡是有一个相同大写字母者表示无性系或试验地间差异不显著（Duncan’s法，P>0.01）。Figure 2. The disease incidences of Eucalyptus clones at six sites2.3 抗风性分析
由于桉树大多种植在台风多发的东南沿海地区，所以桉树抗风性强弱是其推广种植的重要指标[6-7]。试验林营建后第2年受台风“伊布都”的影响，广东3个试验地的桉树都受到较大程度的台风危害，经过试验地的风力8~10级。对桉树试验林的抗风情况调查结果列于表 3，认为抗风无性系（风害率1%~10%）有LH9211、LH9224和SH1，较抗风无性系（风害率10%~20%）有U6和M1，抗风性较差的无性系（风害率20%~40%）有DH3213、DH2012和DH3222，抗风性差的无性系（风害率>50%）为DH3327和DH3226。调查还发现DH2012无性系易风折，其次为DH3226，而M1则易折顶（树顶1.0~1.5 m折断），在吴川受风害严重的小区DH2012风折率高达56%，主要在树干基部1~2 m处折断，且发生在生长最好的IIR5小区，在唐家IIL1小区DH3226风折率也高达52%，在吴川M1风折率高达46%。
表 3中结果同时表明，几乎任何无性系其风害率都随风力的增大而显著增大，吴川风力8级，其平均风害率为16.1%，遂溪风力9级，其平均风害率为25.8%，唐家风力10级，其平均风害率为38.2%。
表 3 广东各试验地1.5年生桉树无性系风害率1）Table 3. The wind damage rates of 1.5-year-old Eucalyptus clones at three Guangdong sites% 试验地 DH3226 DH3213 DH3327 DH3222 DH2012 SH1 M1 LH9211 LH9224 U6 平均值 唐家 70.5±3.2E 48.6±3.2C 80.5±3.4EF 58.3±3.2CD 40.0±2.8C 16.9±1.3B 29.0±1.4BC 3.7±0.4A 18.3±0.9B 16.3±0.8B 38.2 遂溪 65.8±4.3C 22.9±2.3B 65.4±2.1C 27.8±1.3B 15.0±1.2A 11.1±0.6A 19.8±1.2AB 9.2±0.7A 6.2±0.3A 14.5±0.7A 25.8 吴川 55.8±2.1D 7.4±1.2A 37.9±0.9C 13.5±0.7B 25.0±1.7BC 1.1±0.1A 3.2±0.2A 2.1±0.1A 3.4±0.2A 11.8±1.1B 16.1 平均值 64.0 26.3 61.3 33.2 26.7 9.7 17.3 5.0 9.3 14.2 26.7 1）表中同行数据后，凡是有一个相同大写字母者，表示无性系间差异不显著（Duncan’s法，P>0.01）; 试验地唐家、遂溪和吴川的风力依次为10、9和8级。 桉树试验林6年生时，台风“黑格比”在广东阳江地区登陆，中心风力14级，主要对其邻近的吴川和遂溪2个试验地试验林造成了危害，吴川风力12级，遂溪风力11级。对所有观察株（每无性系360株）的调查结果（表 4）显示，6年生桉树试验林以M1的抗风性最差，平均风折率达到17.5%，其次为DH2012，平均风折率12.4%，以LH9211和DH3222最为抗风，风折率为1.8%，风倒情况少有发生。
表 4 6年生桉树无性系受台风危害的风折率1）Table 4. The breakage rates of 6-year-old Eucalyptus clones under typhoon% 试验地 DH3226 DH3213 DH3327 DH3222 DH2012 SH1 M1 LH9211 LH9224 U6 吴川 10.5±1.3C 3.9±0.6A 2.1±0.1A 2.1±0.2A 13.2±1.3CD 2.3±0.4A 13.2±0.9CD 2.1±0.3A 4.5±0.2AB 5.8±0.4B 遂溪 3.8±0.5A 5.6±0.3AB 8.9±0.5B 1.4±0.3A 11.6±1.0C 5.0±0.6AB 21.7±2.0D 1.4±0.3A 8.0±0.7B 1.4±0.1A 平均值 7.2 4.8 5.5 1.8 12.4 3.7 17.5 1.8 6.2 3.6 1）表中同行数据后，凡是有一个相同大写字母者，表示无性系间差异不显著（Duncan’s法，P>0.01）；吴川和遂溪的风力分别为12和11级。 表 3、表 4还表明，不同林龄桉树受台风危害的情况有所不同，幼林易受风害，当风力8级以上时，幼林平均风害率达到16.1%，而当风力11级时成林平均风害率仅为6.9%。幼林主要以风倒或倾斜为主，因其枝叶特别茂盛，有头重脚轻的情况，但无性系间亦有差别，DH2012风折情况同样严重就是例证；成林则主要以风折为主，少见风倒情况，不同林龄不同无性系其抗风性表现亦不相同。总体而言，以来自广东的LH9211和LH9224抗风性最强，其次为SH1和U6，以来自广西东门的DH系列无性系抗风性相对较差。
2.4 速生性分析
短轮伐期桉树的砍伐周期一般为6年。对6.5年生桉树无性系各立地的生长调查结果列于表 5和表 6，从树木的速生性来看，除受台风影响最为严重的唐家试验地外，其余各试验地无论树高和胸径均以DH3213表现最为良好，以U6表现为差，其他无性系在各试验地生长表现各异。
表 5 6.5年生桉树无性系在广东不同试验地的生长表现1）Table 5. The growth traits of 6.5-year-old Eucalyptus clones at different sites in Guangdong Province无性系 吴川 遂溪 唐家 树高/m 胸径/cm 树高/m 胸径/cm 树高/m 胸径/cm DH3226 13.62±2.5E 10.29±0.7D 13.67±1.2BC 11.60±0.6A 14.94±0.6C 11.18±0.9DE DH3213 17.22±1.4A 13.09±1.3A 14.94±2.3A 12.26±0.7A 16.46±1.4B 12.84±0.8ABC DH3327 14.41±1.3D 10.56±1.2BD 11.56±2.2D 9.93±0.7BCD 12.28±1.1D 8.88±1.0F DH3222 15.25±1.2C 11.27±2.1BC 14.87±2.1A 11.81±0.3A 14.43±0.8C 10.97±0.7E DH2012 12.90±0.7F 9.45±1.0E 13.12±1.2C 9.73±1.4BCD 16.51±0.6B 12.15±0.6BCDE SH1 16.41±0.9B 10.60±0.9CD 14.21±2.1AB 9.52±2.2CDE 18.52±2.2A 12.10±0.8CDE M1 15.48±1.1C 11.27±0.6BC 10.99±0.9D 8.74±1.1E 16.72±2.4B 13.49±1.1A LH9211 17.03±1.3AB 11.77±0.7B 13.52±1.3BC 10.23±2.1BC 19.79±1.6A 13.34±0.9AB LH9224 15.68±1.5C 11.46±0.7B 14.00±2.1ABC 10.57±0.3B 17.14±1.4B 12.55±1.2ABC U6 12.44±1.1F 10.22±0.8D 10.92±0.8D 9.08±0.8DE 14.97±1.3C 12.23±1.1BCD 平均值 15.04 11.01 13.47 10.42 16.32 11.99 F值 75.34** 28.56** 38.93** 23.81** 18.14** 73.17** 1）表中同列数据后，凡是有一个相同大写字母者，表示无性系间差异不显著（Duncan’s法，P>0.01）。 表 6 6.5年生桉树无性系在广西不同试验地的生长表现1）Table 6. The growth traits of 6.5-year-old Eucalyptus at different sites in Guangxi Province无性系 钦廉 东门 防城 树高/m 胸径/cm 树高/m 胸径/cm 树高/m 胸径/cm DH3226 17.68±1.2CD 12.95±0.8BCD 17.72±0.8B 13.26±0.7B 13.48±1.1CD 10.66±0.6BC DH3213 20.34±1.3A 15.10±1.0A 18.86±1.2A 15.01±0.6A 15.91±1.3A 12.75±1.1A DH3327 18.69±1.4BC 13.55±0.7B 16.77±1.4CD 12.39±0.9CD 13.45±1.2CD 10.00±0.6CD DH3222 18.84±2.1B 13.15±1.1BC 17.70±1.6B 12.87±0.9BC 14.45±1.1B 10.96±0.9B DH2012 18.20±0.9BCD 12.05±0.4DE 17.36±1.7BC 11.99±2.1DE 13.82±1.2BC 10.98±0.8B SH1 17.69±0.7CD 11.71±0.6E 15.79±0.9EF 11.08±1.2F 13.25±1.0CD 9.68±0.7D M1 17.63±0.8D 12.33±2.0CDE 16.91±1.3BCD 11.87±1.6DEF 14.46±0.6B 10.76±1.3B LH9211 18.42±0.9BCD 12.85±1.2BCD 15.59±1.3F 11.43±0.8EF 13.46±0.9CD 11.02±0.9B LH9224 18.93±0.9B 13.17±0.8BC 16.45±1.2DE 12.51±1.1BCD 13.57±0.5CD 11.02±1.0B U6 14.91±0.6E 11.64±1.1E 13.95±1.1G 11.28±0.9EF 12.94±0.7D 10.39±0.6BC 平均值 18.13 12.85 19.07 13.50 13.88 10.82 F值 29.59** 18.61** 40.93** 27.86** 16.90** 16.25** 1）表中同列数据后，凡是有一个相同大写字母者，表示无性系间差异不显著（Duncan’s法，P>0.01）。 8.5年生时，保存最为完好的钦廉试验地为代表的单株材积见图 3，速生性最好的为DH3213，其单株材积达到0.222 m3·株-1，年均出材量（MAI）为38.87 m3·hm-2·a-1，其次为DH3327和DH3222，单株材积分别为0.167和0.146 m3·株-1，MAI为29.12和25.46 m3·hm-2·a-1，再次为抗风性强的LH9224，单株材积为0.140 m3·株-1，MAI为24.52 m3·hm-2·a-1，表现最差的U6单株材积为0.08 m3·株-1，MAI为17.19 m3·hm-2·a-1。
![图 3 钦廉8.5年生桉树无性系单株材积]()
图 3 钦廉8.5年生桉树无性系单株材积柱子上方，凡是有一个相同大写字母者，表示无性系间差异不显著（Duncan’s法，P>0.01）。Figure 3. The individual volumes of 8.5-year-old Eucalyptus clones in Qinlian图 4结果则表明，立地对桉树产量的影响极为明显，好的立地比差的立地增产1倍以上，6.5年生时，遂溪的平均出材蓄积量为79.6 m3·hm-2，MAI为12.24 m3·hm-2·a-1，而东门的平均出材蓄积量为175.3 m3·hm-2，MAI为26.97 m3·hm-2·a-1。
![图 4 6.5年生桉树无性系各试验地平均出材蓄积量]()
图 4 6.5年生桉树无性系各试验地平均出材蓄积量柱子上方，凡是有一个相同大写字母者，表示试验地间差异不显著（Duncan’s法，P>0.01）。Figure 4. The average volume production of 6.5-year-old Eucalyptus clones at different sites2.5 综合评价
采用模糊数学中隶属函数的方法对各桉树无性系的适生性进行综合评价，隶属值越大，说明桉树的生长表现越好。由表 7结果可知，在广东和广西6个试验地条件下，无性系LH9211、DH3213、SH1、DH3222和U6的隶属值均不小于3.5，高于其他5个无性系，其中，最大隶属值与最小隶属值相差近1倍。表 8表明，由于较低的成活率和最高发病率，种植于唐家的桉树生长综合表现较差，隶属值最低，其他5个立地的隶属值相差不大。就具体某个地区而言：从表 9可以看出，唐家最适合M1和LH9211这2个无性系，隶属值相对较高；从表 10可以看出，吴川适宜种植DH3213、LH9211和SH1这3个无性系，这3个无性系的隶属值排在前三；从表 11可以看出，适宜在遂溪种植的桉树无性系为DH3213和DH3222。
表 7 不同桉树无性系适生性综合评价Table 7. The adaptability evaluation of different Eucalyptus clones无性系 得分 排名 成活率 发病率 风害率 风折率 单株材积 总和 DH3226 0.65 0.94 0.00 0.65 0.40 2.64 7 DH3213 0.41 0.96 0.64 0.81 1.00 3.82 2 DH3327 0.32 0.80 0.05 0.76 0.60 2.53 8 DH3222 0.70 0.91 0.52 1.00 0.46 3.59 4 DH2012 0.97 0.98 0.63 0.32 0.19 3.09 6 SH1 0.73 0.98 0.92 0.88 0.10 3.61 3 M1 1.00 0.39 0.79 0.00 0.22 2.40 9 LH9211 0.57 1.00 1.00 1.00 0.27 3.84 1 LH9224 0.00 0.00 0.93 0.72 0.41 2.06 10 U6 0.84 0.94 0.84 0.88 0.00 3.50 5 表 8 不同试验地桉树的适生性综合评价Table 8. The comprehensive evaluation of Eucalyptus plantation at different sites试验地 得分 排名 成活率 发病率 岀材蓄积量 总和 唐家 0.28 0.00 0.50 0.78 6 吴川 0.78 0.55 0.20 1.53 5 遂溪 0.98 1.00 0.00 1.98 2 钦廉 0.75 0.47 0.87 2.09 1 防城 1.00 0.71 0.07 1.78 4 东门 0.00 0.95 1.00 1.95 3 表 9 唐家10个桉树无性系适生性综合评价Table 9. The adaptability evaluation of ten Eucalyptus clones in Tangjia无性系 得分 排名 成活率 风害率 树高 胸径 总和 DH3226 0.22 0.13 0.35 0.50 1.20 8 DH3213 0.44 0.42 0.56 0.86 2.28 5 DH3327 0.05 0.00 0.00 0.00 0.05 9 DH3222 0.34 0.29 0.29 0.45 1.37 7 DH2012 0.39 0.53 0.56 0.71 2.19 6 SH1 0.39 0.83 0.83 0.70 2.75 3 M1 1.00 0.67 0.59 1.00 3.26 1 LH9211 0.00 1.00 1.00 0.97 2.97 2 LH9224 0.34 0.81 0.65 0.80 2.60 4 U6 0.83 0.84 0.36 0.73 2.76 3 表 10 吴川10个桉树无性系适生性综合评价Table 10. The adaptability evaluation of ten Eucalyptus clones in Wuchuan无性系 得分 排名 成活率 风害率 风折率 树高 胸径 总和 DH3226 0.89 0.00 0.24 0.25 0.23 1.61 10 DH3213 0.70 0.88 0.84 1.00 1.00 4.42 1 DH3327 0.87 0.33 1.00 0.41 0.30 2.91 7 DH3222 0.66 0.77 1.00 0.59 0.50 3.52 4 DH2012 1.00 0.56 0.00 0.10 0.00 1.66 9 SH1 0.87 1.00 0.98 0.83 0.32 4.00 3 M1 0.98 0.96 0.00 0.64 0.50 3.08 5 LH9211 0.83 0.98 1.00 0.96 0.64 4.41 2 LH9224 0.00 0.96 0.78 0.68 0.55 2.97 6 U6 0.89 0.80 0.67 0.00 0.21 2.57 8 表 11 遂溪10个桉树无性系适生性综合评价Table 11. The adaptability evaluation of ten Eucalyptus clones in Suixi无性系 得分 排名 成活率 风害率 风折率 树高 胸径 总和 DH3226 0.64 0.00 0.88 0.68 0.81 3.01 6 DH3213 1.00 0.72 0.79 1.00 1.00 4.51 1 DH3327 0.82 0.01 0.63 0.16 0.34 1.96 9 DH3222 1.00 0.64 1.00 0.98 0.87 4.49 2 DH2012 0.45 0.85 0.50 0.55 0.28 2.63 8 SH1 0.82 0.92 0.82 0.82 0.22 3.60 3 M1 0.73 0.77 0.00 0.02 0.00 1.52 10 LH9211 0.00 0.95 1.00 0.65 0.42 3.02 5 LH9224 0.45 1.00 0.67 0.77 0.52 3.41 4 U6 0.91 0.86 1.00 0.00 0.10 2.87 7 3. 讨论与结论
杂交桉树无性系保留了桉树适应性强的特性，各无性系造林成活率均能达到93%以上，平均为96.0%，对比桉树实生苗造林[13]没有太大的差异，保留了桉树强大生命力的优点，在桉树的适宜种植地区，其成活率均不成为问题。
桉树无性系的选择一要考虑其速生性，二要考虑其抗风性，特别是台风频发的沿海地区。试验中DH系列有相对理想的速生性，但其抗风能力严重不足，在风力9级时即可造成65%以上的风害。因此，建议台风频发地区可多种植抗风性强的尾细桉系列的LH9224无性系，以减少因风倒伏以致重新种植的麻烦，产量亦能达中等以上水平，其在福建闽南地区亦有良好表现[12]。在内陆少有风害的地区，DH3213和DH3327是较为理想的种植品系，其能使产量提高30%以上，并取得较为理想的经济效益，DH系列的速生性也被前人的试验所肯定[13, 15]。
10个桉树无性系，在唐家和吴川的综合评价隶属值较差，因此，对于成活率差和发病率最高的唐家，只能选择成活率最高的M1和具抗病性的LH9211，而对于出材率较差的吴川和遂溪，适宜选择速生性最佳的DH3213无性系。
立地条件的好坏显著地影响桉树单位面积的产量，良好立地条件下的桉树产量可以比较差立地条件下桉树产量多1倍甚至3倍以上，因此，林地的选择极为重要，虽然桉树的适生性很强，但不同无性系产值有着极大的差别，因此，选好地、种好桉树具有极为现实的意义。
桉树各无性系的抗病性虽然有差别，但目前桉树种植更多地使用组培苗，并在造林地附近直接培育进行造林，一定程度上减少了病害的传播。桉树成林的病害也可能在桉树与相思Acacia confusa、桉树与西南桦Betula spp、桉树与红椎Castanopsis hystrix等的混交种植模式中得以减少，改变桉树大面积纯林种植的方式是减少桉树病虫害发生的有效途径。
-
![]()
图 1 水稻考种信息采集装置
1：安装支架；2：图像采集设备；3：光源；4：透明标定膜；5：质量计量装置；6：计算机
Figure 1. Acquisition device of rice seed investigation information
1: Mounting bracket; 2: Image acquisition equipment; 3: Light source; 4: Transparent calibration film; 5: Mass metering device; 6: Computer
![]()
图 3 谷粒区域提取及总粒数和空粒数计数流程图
Figure 3. Flowchart of regional extraction of grains and counting of total and empty grains
![]()
图 6 各因素各水平对稻谷总计数精度的影响
Figure 6. Effects of various factors and levels on counting accuracy of total grains
![]()
图 7 稻谷总数计数值与实际值的关系
Figure 7. Relationship between count value and actual value of total grains
![]()
图 8 各因素各水平对稻谷空粒计数精度的影响
Figure 8. Effects of various factors and levels on counting accuracy of empty grains
![]()
图 9 稻谷空粒计数值与实际值的关系
Figure 9. Relationship between count value and actual value of empty grains
![]()
图 10 各因素各水平对稻谷长度、宽度测量精度的影响
Figure 10. Effects of various factors and levels on the measurement accuracy of grain length and width
![]()
图 11 稻谷长度、宽度测量值与实际值的关系
Figure 11. Relationship between measured value and actual value of grain length and width
![]()
图 12 各因素各水平对规则图形周长、面积测量精度的影响
Figure 12. Effects of various factors and levels on the perimeter and area measurement accuracy of regular patterns
![]()
图 13 规则图形周长、面积测量值与实际值的关系
Figure 13. Relationship between measured value and actual value of perimeter and area of regular patterns
表 1 测量效率测试试验结果
Table 1 Results of measurement efficiency test
测试组号
Group
number稻谷粒数
Rice grain
quantity每组所需时间/s
Time required
for each group每秒测量谷粒数
Measured grain
quantity per second1 152 9.12 16.67 2 143 9.20 15.54 3 174 9.61 18.11 4 168 9.46 17.76 5 141 9.19 15.34 6 127 9.00 14.11 7 181 9.65 18.76 8 130 8.97 14.49 9 158 9.11 17.34 10 162 9.47 17.11 
下载: 导出CSV
-
[1] 邓若玲, 潘威杰, 王志琪, 等. 农作物表型技术及其智能装备研究进展与展望[J]. 现代农业装备, 2021, 42(1): 2-9. [2] FURBANK R T, TESTER M. Phenomics-technologies to relieve the phenotyping bottleneck[J]. Trends in Plant Science, 2011, 16(12): 635-644. doi: 10.1016/j.tplants.2011.09.005
[3] AFONNIKOV D A, GENAEV M A, DOROSHKOV A V, et al. Methods of high-throughput plant phenotyping for large-scale breeding and genetic experiments[J]. Russian Journal of Genetics, 2016, 52(7): 688-701. doi: 10.1134/S1022795416070024
[4] 周洪垒. 基于图像处理的水稻考种系统的设计与实现[D]. 成都: 电子科技大学, 2019. [5] FIORANI F, SCHURR U. Future scenarios for plant phenotyping[J]. Annual Review of Plant Biology, 2013, 64: 267-291. doi: 10.1146/annurev-arplant-050312-120137
[6] DENG R L, TAO M, HUANG X N, et al. Automated counting grains on the rice panicle based on deep learning method[J]. Sensors, 2021, 21(1): 281. doi: doi.org/10.3390/s21010281.
[7] VIRLET N, SABERMANESH K, SADEGHI-TEHRAN P, et al. Field scanalyzer: An automated robotic field phenotyping platform for detailed crop monitoring[J]. Functional Plant Biology, 2016, 44(1): 143-153.
[8] GUO Q H, WU F F, PANG S X, et al. Crop 3D: A LiDAR based platform for 3D high-throughput crop phenotyping[J]. Science China-Life Sciences, 2018, 61(3): 328-339. doi: 10.1007/s11427-017-9056-0
[9] 宋鹏, 张晗, 王成, 等. 玉米籽粒考种信息获取装置设计与试验[J]. 农业机械学报, 2017, 48(12): 19-25. doi: 10.6041/j.issn.1000-1298.2017.12.002 [10] 胡正方, 向阳, 熊瑛, 等. 基于机器视觉的千粒质量测量仪的设计与试验[J]. 湖南农业大学学报(自然科学版), 2021, 47(4): 476-481. doi: 10.13331/j.cnki.jhau.2021.04.017 [11] 崔亮. 基于机器视觉的农作物种子计数检测系统[D]. 太原: 中北大学, 2016. [12] 杨万能. 水稻产量相关性状参数自动提取的数字化技术研究[D]. 武汉: 华中科技大学, 2011. [13] 宋初苏, 蓝景针, 夏世峰, 等. ZPXG-18型转盘斜刮式自动光电数粒仪和千粒重仪的设计[J]. 浙江农业学报, 2011, 23(5): 1023-1028. doi: 10.3969/j.issn.1004-1524.2011.05.030 [14] ARVIDSSON S, PÉREZ-RODRÍGUEZ P, MUELLER-ROEBER B. A growth phenotyping pipeline for Arabidopsis thaliana integrating image analysis and rosette area modeling for robust quantification of genotype effects[J]. The New Phytologist, 2011, 191(3): 895-907. doi: 10.1111/j.1469-8137.2011.03756.x
[15] YANG W N, GUO Z L, HUANG C L, et al. Combining high-throughput phenotyping and genome-wide association studies to reveal natural genetic variation in rice[J]. Nature Communications, 2014, 5(1): 1-9. doi: 10.1038/ncomms6087
[16] 段凌凤, 杨万能. 水稻表型组学研究概况和展望[J]. 生命科学, 2016, 28(10): 1129-1137. doi: 10.13376/j.cbls/2016151 [17] HAIN R, KÄHLER C J, TROPEA C. Comparison of CCD, CMOS and intensified cameras[J]. Experiments in Fluids, 2007, 42(3): 403-411. doi: 10.1007/s00348-006-0247-1
[18] 袁洪波, 赵努东, 程曼. 基于图像处理的田间杂草识别研究进展与展望[J]. 农业机械学报, 2020, 51(S2): 323-334. doi: 10.6041/j.issn.1000-1298.2020.S2.038 [19] 夏庆凯, 于大泳. 基于图像处理的工业生产手机屏幕裂纹检测[J]. 农业装备与车辆工程, 2021, 59(1): 112-115. doi: 10.3969/j.issn.1673-3142.2021.01.024 [20] 项荣, 应义斌, 蒋焕煜. 田间环境下果蔬采摘快速识别与定位方法研究进展[J]. 农业机械学报, 2013, 44(11): 208-223. doi: 10.6041/j.issn.1000-1298.2013.11.037 [21] 孔明, 徐志玲, 徐勇, 等. 实际边缘轮廓下的螺栓综合尺寸质量检测[J]. 计量学报, 2021, 42(6): 724-730. doi: 10.3969/j.issn.1000-1158.2021.06.06 [22] 贾建楠, 吉海彦. 基于病斑形状和神经网络的黄瓜病害识别[J]. 农业工程学报, 2013, 29(S1): 115-121. [23] 李擎, 唐欢, 迟健男, 等. 基于改进最大类间方差法的手势分割方法研究[J]. 自动化学报, 2017, 43(4): 528-537. doi: 10.16383/j.aas.2017.c150862 [24] 伍艳莲, 赵力, 姜海燕, 等. 基于改进均值漂移算法的绿色作物图像分割方法[J]. 农业工程学报, 2014, 30(24): 161-167. doi: 10.3969/j.issn.1002-6819.2014.24.019 [25] 胡炼, 刘海龙, 何杰, 等. 智能除草机器人研究现状与展望[J]. 华南农业大学学报, 2023, 44(1): 34-42. doi: 10.7671/j.issn.1001-411X.202203060 [26] 贺锋, 钟宏民, 胡友旺. 基于图像处理的焊缝跟踪检测方法研究[J]. 应用激光, 2020, 40(5): 847-854. doi: 10.14128/j.cnki.al.20204005.847 [27] 胡园园. 施肥与灌溉对水稻颖壳不闭合影响研究[D]. 合肥: 安徽农业大学, 2011. [28] 赵凯旋, 刘晓航, 姬江涛. 基于EfficientNet与点云凸包特征的奶牛体况自动评分[J]. 农业机械学报, 2021, 52(5): 192-201. doi: 10.6041/j.issn.1000-1298.2021.05.021 [29] 赵强, 黄远, 赵春海. 基于凸包的电池片边缘缺陷检测算法[J]. 计算机测量与控制, 2020, 28(6): 129-134. doi: 10.16526/j.cnki.11-4762/tp.2020.06.027 [30] 李鹏阳, 杨文辉, 曹利平, 等. 圆盘类零件几何量机器视觉测量技术研究[J]. 西安工业大学学报, 2020, 40(6): 598-604. doi: 10.16185/j.jxatu.edu.cn.2020.06.002 [31] 张文霞, 王春光, 王海超, 等. 无芒隐子草叶片卷曲度和厚度测量方法[J]. 农业机械学报, 2021, 52(5): 184-191. doi: 10.6041/j.issn.1000-1298.2021.05.020 [32] 古亭. 基于机器视觉的蚕茧计数与识别研究[D]. 杭州: 中国计量大学, 2020. [33] 汪晓艳, 丁义明. 一种去除图像中Cauchy噪声的滤波算法[J]. 数学物理学报, 2018, 38(4): 823-832. doi: 10.3969/j.issn.1003-3998.2018.04.018 -
期刊类型引用(5)
1. 许文英. 西南桦的生长特征及栽培技术探讨. 农村科学实验. 2024(07): 82-84 . 
百度学术
2. 连晓倩，田云龙，黄健，钟兆全，周垂帆，黄伟，吴鹏飞. 杉木林分密度对套种闽楠树型、光合能力及生物量分配的影响. 林业科学研究. 2024(05): 13-22 . 百度学术
3. 余清海. 不同基质·遮阴强度对盐肤木种子萌发的影响. 安徽农业科学. 2023(23): 90-92+103 . 百度学术
4. 代嵩华，黄印冉，冯树香，闫淑芳，刘易超，陈丽英，樊彦聪. 不同遮阴强度对白桦幼苗生长情况的影响. 河北林业科技. 2022(02): 25-27+36 . 百度学术
5. 朱秀征，靳超，房丽莎，刘震，王艳梅，蔡齐飞，李志，耿晓东. 不同密度下山桐子自然整枝差异性分析. 河南科学. 2022(08): 1257-1263 . 百度学术
其他类型引用(3)
