基于水稻三维模型的表型参数提取及生物量估测.pdf

河南农业科学,2 0 2 3,52(7)：144-153Journal of Henan Agricultural Sciencesdoi:10.15933/ki.1004-3268.2023.07.015基于水稻三维模型的表型参数提取及生物量估测程志强，方圣辉（武汉大学遥感信息工程学院，湖北武汉430 0 7 2）摘要：基于二维图像的水稻表型检测方法缺乏空间维度信息，难以提取全面的水稻表型参数，因此建立水稻三维模型进而提取全面的水稻表型参数具有重要意义。体积是指示水稻生长状况的重要参数，如何在不破坏水稻植株的前提下获取水稻体积参数是当下研究中仍有待解决的问题。基于上述考虑，提出一种基于水稻三维模型的表型参数提取及生物量估测方法，实现在不破坏水稻植株的前提下提取水稻体积参数并估测水稻生物量。以盆栽水稻为对象，首先利用Alpha-shape算法重建其三维模型，然后提取并评估水稻株高、茎粗、植被覆盖度、体积参数，并基于体积参数估测水稻生物量。结果表明，空间雕刻方法可以重建高精度的水稻三维模型，并准确测算水稻表型参数及生物量，其中水稻株高、茎粗、单株植被覆盖度均方根误差（RMSE）和平均绝对百分比误差（MAPE）分别为6 3.2 7 mm、4.0 1mm、5.0 4%和7.15%、14.91%、12.59%，水稻生物量RMSE和MAPE分别为6 2.44g、19.2 5%，均优于现有三维重建软件的测算结果。关键词：水稻；三维重建；Alpha-shape；表型参数提取；体积；生物量中图分类号：S126Rice Phenotypic Parameters Extraction and Biomass Estimation(School of Remote Sensing and Information Engineering,Wuhan University,Wuhan 430072,China)Abstract:The phenotype detection methods based on two-dimensional images lack spatial dimensioninformation and have dfficulty in extracting comprehensive rice phenotype parameters.Therefore,it is ofgreat significance to establish a three-dimensional model of rice to extract comprehensive rice phenotypeparameters.The volume is an important parameter indicating the growth status of rice.How to obtain ricevolume parameters without damaging the rice plant is still a problem to be solved in current research.Based on the above considerations,a method of rice phenotypic parameters extraction and biomassestimation based on three-dimensional model is proposed in this paper,which can extract rice volumeparameters and estimate rice biomass without damaging the rice plant.This study focused on potted rice,first reconstructed its three-dimensional model using the Alpha-shape algorithm,and then extracted andevaluated rice height,stem thickness,vegetation coverage,volume parameters,and estimated ricebiomass based on the volume parameter.The experimental results showed that space carving could文献标志码：ABased on Three-Dimensional ModelCHENG Zhiqiang,FANG Shenghui文章编号：10 0 4-32 6 8（2 0 2 3)0 7-0 144-10收稿日期：2 0 2 3-0 3-10基金项目：湖北省重点研发计划项目（2 0 2 0 BBB058）作者简介：程志强（1998-），男，四川井研人，在读硕士研究生，研究方向：水稻植株三维重建和表型分析。E-mail;通信作者：方圣辉（196 4-），男，湖北武汉人，教授，主要从事定量遥感等研究。E-mail：s h f a n g w h u.e d u.c n第7 期reconstruct a high-precision three-dimensional model of rice and accurately measure rice phenotypeparameters and biomass,the RMSE and MAPE of rice height,stem thickness,and single plant vegetationcoverage were 63.27 mm,4.01 mm,5.04%and 7.15%,14.91%,12.59%,and the RMSE and MAPE ofrice biomass were 62.44 g and 19.25%,which were better than the results obtained by existingthree-dimensional reconstruction software.Key words:Rice;3D reconstruction;Alpha-shape;Phenotypic parameters extraction;Volume;Biomass水稻一直是我国重要的粮食作物之一，对维系我国粮食安全和社会稳定发挥着重要作用叫。为了实现水稻作物的遗传改良和产量提升，当下许多研究关注于水稻作物表型2-。水稻作物表型提取方法已经经历了人工测量、二维图像分析和专业数字化仪器测量等阶段4。其中，人工测量、专业数字化仪器测量的方法5 劳动量大且具有一定的破坏性。相较于二维图像分析技术，三维重建技术可以提取二维图像不具备的空间信息，成为近年来研究的热点。当前，植株三维重建方法主要分为两类。一类是利用传感器直接获得植物点云的主动视觉法。其中，PAULUS等7 使用三维激光扫描仪实现了大麦器官的自动扫描建模，TANG等18 使用时差法和可见光成像提取玉米苗期植株性状。然而，这类方法因受到设备昂贵、环境要求苛刻、重建效率低下等限制无法投入大规模应用。另一类是基于多视角图像的摄影测量方法，如运动恢复重建（Structurefrommotion，SFM)19、空间雕刻（Structure from spacecarving,SSC)i等。这类方法能够捕获植株高密度空间点云信息和色彩纹理信息，从而实现更优的三维重建效果。胡鹏程等利用SFM方法实现了盆栽黄瓜、茄子、青椒的三维重建，但该方法依赖于特征点匹配结果，难以完成水稻等纹理稀疏、叶片细窄植物的三维重建。SSC算法利用物体的多视角图像和二值轮廓，采用视角交集的方法计算植株的可视化外壳，最终保留体素构成的目标三维模型，可以高效稳定地完成叶片细窄植物的三维重建。基于该方法，吴丹等112 通过旋转拍摄盆栽水稻图像，成功实现了对不同时期和品种的水稻样本进行三维重建。魏天翔等在空间雕刻的基础上加入了深度学习背景分割技术，实现了对田间单株水稻的重建。水稻体积(Volume of rice,VR)是一个与水稻生物量信息密切相关、指示水稻生长状况的重要参数指标，但在当下研究中却并未得到广泛关注。传统程志强等：基于水稻三维模型的表型参数提取及生物量估测1.1水稻体素三维重建1.1.1序列图像采集图像采集地点位于武汉大学花山水稻试验基地，采集设备为CanonEOS6DMark单反相机，拍摄分辨率为16 0 0 像素2 40 0像素，固定焦距f为2 4mm，保持拍摄光圈、ISO和曝光时间等参数一致，以最大程度地还原水稻真实色彩。145的体积测量方法往往会对水稻植株造成破环，而最近提出的基于植株三维模型的体积测量方法如点云包围盒方法3、凸包法14,则由于水稻结构复杂、内部存在大量遮挡，难以获取准确的水稻体积。因此，在不破坏水稻植株的前提下，如何精确获取水稻的体积参数，仍是当下研究中函待解决的问题。Alpha-shape算法最初用于构建点集轮廓，目前已广泛应用于建筑物轮廓提取和三维重建15。该算法可以构建反映物体在空间中表面积、体积、疏密程度等特征的模型，但如何利用它提取水稻点云轮廓并进行三维曲面重建仍缺乏深入研究。鉴于此，利用Alpha-shape算法生成水稻植株三维曲面模型，精确提取株高（Plantheight，PH）、茎粗（Stemthickness，ST）、植被覆盖度（Vegetation coverage,VC)、体积等多类重要表型参数，并使用水稻体积估算水稻生物量（Biomassofrice,BR），为水稻表型研究提供一种无损可靠的参数测算方法。1材料和方法图1展示了基于水稻三维模型提取表型参数和估测生物量的过程。首先，采集了水稻的多视角图像，利用张氏标定法确定相机内外参数16 并通过色彩信息分离水稻与背景。接下来，通过空间雕刻方法生成水稻体素模型。随后，实现水稻表型参数的提取和生物量的估测。具体而言，基于水稻体素模型构建曲面模型，提取株高、茎粗、植被覆盖度和水稻体积等表型参数，并利用水稻体积来估测水稻生物量。最后，将获得的表型参数提取结果与人工测量值进行比对分析，以验证模型的精度。146河南农业科学第52 卷水稻参数采集Rice parameter acquisition体素模型重建Voxel model reconstruction序列图像采集Sequence image acquisition相机标定Camera calibrationSpace carving工模型着色Model coloring水稻曲面模型重构Reconstruction ofricesurfacemodel模型精度检测Model accuracy test植被覆盖度水稻株高Vegetation coveragePlant height of rice7俯视图轮廓提取TopviewSilhouette extraction重投影Reprojection空间雕刻水稻茎粗Stem thickness of rice最小外接矩形茎部聚类去噪Minimumbounding rectangleStem cluster denoisingAABB包围盒模型分割Axis-alignedbounding boxModel segmentation水稻表型参数分析Analysisofricephenotypic parameters水稻生物量Biomass of rice线性回归计算Linear regression calculation水稻体积Volume ofrice图1水稻三维重建流程图Fig.1Diagram of the 3D reconstruction flow of rice为了对水稻进行全面观察和研究，选用了1穴2株种植的盆栽水稻作为试验对象。在合适的水稻观察点放置棋盘格标定板并进行多角度拍摄，拍摄距离保持在0.6 1.0 m，以保证完整地捕捉水稻在各个方位视角的表型信息。本研究共包含4个时期（幼苗期、分期、拔节期和抽穗期)的水稻样本，每株水稻采集了2 0 张照片，包含水稻俯视图、正反视图和左右视图。1.1.2水稻轮廊提取水稻轮廓提取是水稻体素三维重建的必要步骤，旨在去除投影在水稻轮廓外的体素块，从而分离出精确的水稻体素模型。在可见光范围内，抽穗期以前的水稻呈现为绿色，而抽穗期至成熟期的水稻颜色则逐渐变黄。根据上述特性，基于色彩信息提取水稻轮廓。HSV空间 7 能够直观地表达色彩的明暗、色调以及鲜艳程度，常用于分割指定颜色的物体。因此，采用HSV分量阈值对水稻轮廓进行提取。HSV分量阈值设置如表1所示，水稻原始图像和水稻轮廓提取结果如图2 所示。表1HSV阈值设置Tab.1HSV threshold settingsHSV阅值区间Threshold interval of HSV绿色Green黄色Yellowba.水稻图像；b.背景分割；c.轮廓提取a.Rice image;b.Background segmentation;c.Silhouette extraction图2 水稻原始图像和水稻轮廓提取结果Fig.2Original rice image and silhouette extraction result1.1.3相机标定相机标定旨在建立包围盒中每个体素的三维空间坐标与对应的二维像素坐标之间的映射关系，该映射关系如公式(1)所示。X.X.uPYZS=KRtZ1L1其中，X,Y,z,1表示世界坐标系下的坐标点，u,u,1则表示该点在某视图图像中的像素点坐标。世界坐标系是以棋盘格固定内角点为原点，棋盘格所在平面为oxy平面的坐标系。P是将空间坐标转换为像素坐标的投影矩阵，由相机内参数矩阵HS0.160.400.161.000.100.130.171.00YL1VK和相机外参数矩阵Rt相乘所得。R、t 分别是相0.181.00机的旋转矩阵和平移向量，它们通过张氏标定法获0.181.00得，用于表示相机在世界坐标系中的位置，如图3所(1)第7 期示。s为尺度因子。相机位姿Camerapose-1 000-800wu/z-600-400-2000-6000-400-200Y/mm图3相机在世界坐标系中的位置Fig.3Position of the camera in world coordinate system1.1.4空间雕刻在水稻体素三维重建的过程中，需要获取相机标定得到的世界坐标范围，在范围内建立包围盒并进行体素化。然而，如果体素块过粗略计算缩小包围盒Rough calculation ofshrinking boundingbox标定范围Calibration range程志强等：基于水稻三维模型的表型参数提取及生物量估测标定板Calibration targetX/mm0200147大，会导致水稻建模结果过于粗糙甚至出现错误；而过小则会极大地增加程序的时间复杂度和空间复杂度。为平衡三维重建精度和程序运行速度，首先采用边长为3mm的体素块初始化网格，然后以边长为1mm的体素块重新建立更加准确的包围盒，执行空间雕刻算法。空间雕刻流程如图4所示，首先根据公式（1)计算包围盒中每个体素块角点在某一视图图像的像素坐标，去除坐标在水稻轮廓外的体素块。接下来，遍历所有视图图像，重复上述步骤进行水稻外500壳的细化，最终生成水稻体素模型。随后，将水稻体素模型转换为点云模型。最后遍历点云，使用投-500影变换公式计算每个点云在不同视图的投影坐标，以点云在所有视图中的RGB均值C(w)作为该点的色彩值，C(w)计算过程如公式(2)所示。C(nu)=c(P,u)n其中，n代表视图数量，P代表第i个视图对应的投影矩阵，W表示水稻三维点云中的任意点。1001张图像雕刻结果500Result after1carving(2)5张图像雕刻结果9X6标定板9x6 chessboardResultafter5carvings重投影Reprojection色彩重投影与高度渣染Color reprojection andheight renderingFig.4 Process of space carving1.2水稻表型参数提取及生物量估测1.2.1水稻曲面模型重构及表型参数提取Alpha-shape是一种边界点的提取算法，适用于无序点云的边界提取和将点云模型转换为曲面模型。该算法基于空间中的3个点构建半径为的球，以判定点云模型的边界。参数直接决定了算法的精度和模型的重建效果，通过调整值，利用该算法建所有图像雕刻结果Result after allcarvings图4空间雕刻流程立不同精度的水稻曲面模型。图5展示了利用Alpha-shape算法重建分期水稻曲面模型的效果。当=100时，水稻内部结构被高度包围，曲面模型接近于水稻外部点云构建的凸包，无法展现水稻的具体结构。随着的减小，重建出的曲面模型更接近于水稻原始图像和点云的结果。当=10时，曲面模型的重建结果较好，水稻的几何结构表达效果148较为精细，但存在部分叶片粘连的问题。当减小到1时，水稻的曲面模型与点云模型基本一致，此时的模型接近于真实水稻。当继续减小时，由于小于水稻点云之间的距离，无法建立有效的曲面模型。通过调整参数值，构建出不同精度的水稻曲面模型，反映出水稻在空间伸展时的不同形态和疏河南农业科学密特征。当较大时，可以通过计算凸包体积获取水稻在生长过程中占据的实际空间，当较小时，水稻模型更接近真实形态，此时可以利用Alpha-shape算法将曲面模型分割为多个三角体区域，计算各三角体区域的体积并相加得到整个曲面模型的体积，从而实现水稻体积的提取。第52 卷水稻原图像Original image of rice=100=75=50点云Point cloudFig.5 Curved surface models of rice with different a parameters将水稻株高定义为植株最高点到盆沿所在平面的垂直距离。在本研究中，标定板水平放置在盆沿所在平面上，因此将标定板所在平面定义为世界坐标系中的Z=0平面，水稻三维模型的高度即为水稻株高参数。水稻株高的提取方式如图6 所示。株高参数计算Calculation of plantheight盆沿所在高度Height of theedgeofthepot单株水稻模型Ricemodel图6 水稻株高参数提取示意图Fig.6 Schematic diagram of extraction of rice plant height在水稻进入分期后，由于分数量增加且根部遮挡较少，在此时期测量水稻茎粗可以获得更加准确的结果。本研究将水稻茎粗定义为整株水稻盆沿平面往上1cm根茎横截面的最长直径，在点云模型中表现为根茎横截面最大坐标值与最小坐标值的差值。在提取根茎横截面时，可能存在叶片遮挡引起的噪点，因此在计算坐标差值以前需要进行=10图5不同参数下的水稻曲面模型去噪处理。水稻茎粗提取的具体过程如图7 所示。根茎截取Stem segmentation水稻模型Ricemodel株高Plant height株高提取Extraction of plant height=5Fig.7 Schematic diagram of extraction of rice stem thickness水稻植被覆盖度被定义为水稻植株在地面的垂直投影面积占最小外接矩形面积的百分比，人工测量值则为水稻拍摄俯视图植株部分与其最小外接矩形面积的百分比。为了提取该参数，需要将水稻曲面模型平面投影至Z=0平面，然后计算在二维图像中的占比。植被覆盖度具体提取方式如图8所示。1.2.2水稻生物量估测通常，植物生物量指单位面积内有机物质干质量的总量18,然而水稻作为一种对水分极为敏感的粮食作物，其鲜质量仍然具有重要的研究意义。因此,将水稻生物量定义为水稻=1茎秆提取Stem extraction图7 水稻茎粗参数提取示意图点云去噪Point clouddenoising茎粗Stem thickness茎粗计算Calculation of stemthickness第7 期植株地上部分的鲜质量。在数据采集阶段，以单株水稻地上部分的鲜质量作为水稻生物量的人工测量值，以=1时水稻曲面模型的体积作为水稻体积计算值，建立两者的线性回归预测模型，并基于该模型估测水稻生物量。图8 水稻植被覆盖度参数提取示意图Fig.8 Schematic diagram of schematic diagram ofextraction of rice vegetation coverage程志强等：基于水稻三维模型的表型参数提取及生物量估测水稻体素三维重建结果采集不同时期的水稻多视角图像，使用VisualSFMu19和ContextCapture20软件生成了相应的点云模型。与此同时，使用空间雕刻方法实现了水稻三维重建，并与2 种软件的结果进行对比分析。表2 展示了各种重建方法生成的平均水稻点云最小外接矩形数量。Minimumbounding rectangle1492结果与分析2.1表2 不同重建方法的平均水稻点云数量Tab.2Average number of rice point clouds with differentreconstruction methodVisualSFMContext Capture1.771031.47106通过对比结果可以发现，VisualSFM重建点云数量极少，这是因为VisualSFM没有重建出有效的水稻模型，而ContextCapture重建的水稻点云数量庞大，其中存在大量的背景余。由于VisualSFM无法重建水稻模型，不再对其进行讨论。图9展示空间雕刻SSC6.43105aCd植株IPlant Ia.ContextCapture重建的水稻模型；b.ContextCapture重建的水稻俯视深度图；c.空间雕刻重建的水稻模型；d.空间雕刻重建的水稻俯视深度图a.Reconstructed rice model of Context Capture;b.Reconstructed top view depth maps of Context Capture;c.Reconstructed rice model of SSC;d.Reconstructed top view depth maps of SSCFig.9Reconstruction results of Context Capture and SsC植株Plant II图9Context Capture和空间雕刻重建结果植株Plant II植株IVPlant IV150了使用Context Capture软件和空间雕刻方法得到不同时期的水稻三维模型。重建结果显示，ContextCapture 的重建结果存在严重的水稻叶片和茎秆扭曲失真，而相比之下，空间雕刻方法能够重建更完整的水稻模型，还原出水稻真实形态。2.2水稻表型参数提取结果2.2.1水稻体积由于水稻体积真值无法直接获取，将Alpha-shape提取的水稻体积与凸包法提取的水稻体积进行对比，用于估算本研究水稻体积提取的准确性。具体而言，提取了不同值下的水稻体积，并计算了它们与凸包法提取的体积之间的决定系数（R)和平均绝对百分比误差（MAPE）。表3不同对应体积参数计算结果与凸包法所计算水稻体积之间的决定系数及平均绝对百分比误差Tab.3 R and MAPE of rice volume calculated by Alpha-shape with different and convex hull指标IndicatorR2MAPE/%2.2.2水稻其他表型参数利用Context Capture和空间雕刻重建水稻模型，分别提取水稻俯视图的最小外接矩形和包围盒，并对三维模型进行点云分割和去噪处理，从而获得了水稻的株高、植被覆盖度河南农业科学由表3可知，当=1000时，2 种方法得到的水稻体积趋于相等，决定系数高达0.99，MAPE仅为3.97%。随着的减小，水稻曲面模型逐渐变得复杂和真实，求得的水稻体积逐渐减少，与凸包法得到的体积差值变大，相关性也逐渐减小。当=1时，2 种方法的决定系数仅为0.42，而MAPE则高达97.14%，说明由于水稻的茎叶延展特性不同，相同的水稻外包模型内部结构差异明显，也进一步说明随着的减小，Alpha-shape算法构建的水稻模型能更加真实地反映水稻生长情况，本研究在=1的情况下求得的水稻体积最接近真实值。10005000.990.993.9710.35第52 卷100750.510.5626.1139.72和茎粗等参数，并计算了参数计算值与人工测量值之间的MAPE和均方根误差（RMSE)指标。图10 为水稻表型参数计算值与人工测量值散点图，参数评估结果显示，空间雕刻方法提取的水500.5555.29100.4784.5750.4691.2810.4297.14a1800160014001200100080060040020000200株高人工测量值/mmManualmeasurement of PHd120010008006004002000020040060080010001200株高人工测量值/mmManual measurement of PHa.ContextCapture株高散点图；b.ContextCapture茎粗散点图；c.ContextCapture植被覆盖度散点图；d.空间雕刻株高散点图；e.空间雕刻茎粗散点图；f.空间雕刻植被覆盖度散点图a.Scatter plots of PH of Context Capture;b.Scatter plots of ST of Context Capture;c.Scatter plots of VC of Context Captured.Scatter plots of PH of space carving;e.Scatter plots of ST of space carving;f.Scatter plots of VC of space carvingFig.10 Scatter plot of measured and calculated rice phenotypic parameter valuesR=0.50RMSE=214.40 mmMAPE=24.41%40060080010001200R=0.97RMSE=63.27mmMAPE=7.15%b9080RMSE=11.79mm70MAPE=29.55%60504030201000e7060504030201000图10 水稻不同表型参数计算值与人工测量值散点图R=0.651530茎粗人工测量值/mmManual measurement ofSTR=0.96RMSE=4.01mmMAPE=14.91%1530茎粗人工测量值/mmManual measurement of TC10080604020045604560R=0.53RMSE-14.13%MAPE=30.03%750f90%8070605040302010075020植被覆盖度人工测量值/%Manual measurementofvcR=0.92RMSE=5.04%MAPE-12.59%2040植被覆盖度人工测量值/%Manual measurement ofvc4060608080100100第7 期稻表型参数的R、RM SE和MAPE均优于专业三维建模软件Context Capture，验证了空间雕刻方法在水稻表型参数提取方面的有效性。2.3水稻生物量估测结果图11展示了人工测量的水稻生物量与水稻体积参数的散点图。结果显示，体积参数与水稻生物量信息的相关性较弱，R仅为0.34。考虑到水稻孕穗以后，相同体积下稻穗质量明显更大，因此植株是否含穗会对水稻生物量信息产生影响，从而影响最终的提取结果。为了分析水稻是否含穗对水稻生物量测量的影响，将数据样本分为了含穗和不含穗2 类，并将它们作为独立数据样本进行回归分析。70060050040030020010000500100015002000250030003500水稻体积计算值/mm3Calculated value of VR图11水稻生物量测量值和体积计算值散点图Fig.11 Scatter plot of measured BR and calculated VR图12 展示了分类后的水稻生物量人工测量值和体积参数提取结果的散点图，其中蓝色点代表不含穗的水稻植株样本，橙色点代表含穗的植株样本。与图11相比，图12 中2 类样本各自计算的回归结果更具相关性，其中含穗水稻样本R达到了0.77,不含穗的水稻样本R达到了0.8 3，相较于分类前有了明显的提高。此外，含穗水稻植株和不含穗水稻植株的分布态势也符合含穗水稻植株密度更大的现实规律。70060050040030020010000500100015002000250030003500水稻体积计算值/mm3Calculated values of VR图12分类水稻生物量测量值和体积计算值散点图Fig.12 Scatter plot of categorical measured BR andcalculatedVR程志强等：基于水稻三维模型的表型参数提取及生物量估测RMSE=62.44 gMAPE=19.25%50040030020010000500100015002000250030003500水稻生物量计算值/mm3Calculated valueof BR图13水稻生物量测量值和计算值散点图R2=0.34y=0.13x+161.91R=0.77y=0.68x-49.60不含穗植株Plantswithoutspikes含穗植株Plantswithspikes151图13展示了将植株分为两类后水稻生物量的人工测量值与系统计算值的散点图及y=x参照。可以看到，分为2 类回归计算的水稻生物量和人工检测值之间存在较强的相关性，RMSE和MAPE均可控制在较小范围内，证明了本研究利用水稻体积参数估测水稻生物量的有效性及可靠性。700R2=0.93600不含穗植株Plantswithoutspikes含穗植株PlantswithspikesFig.13Scatter plot of measured BR and calculated BR2类数据分别计算的结果显示，不含穗植株的RMSE和MAPE分别为53.38 g、17.52%，含穗植株的RMSE和MAPE分别为8 0.40 g、19.95%。可以看出，在不含穗情况下水稻植株体积与生物量的相关程度大于含穗情况，这是因为含穗水稻在不同时期稻穗生物量体积比例不同，从而影响最终的提取结果。3结论与讨论本研究结果表明，使用空间雕刻方法重建水稻三维模型比SFM方法更为优越，在点云余、表型参数估测等方面均表现更好，其中水稻株高、茎粗、单株植被覆盖度的MAPE分别为7.15%、14.9 1%和12.59%，均低于使用ContextCapture软件重建提取的水稻表型误差结果。这是因为空间雕刻算法不需要进行特征匹配，所以其细化和分割体素的过程不会扭曲和穴余水稻模型，从而提取的相关表型更R2=0.83加准确2 1。尽管如此，SFM方法仍然是目前多视角y=0.15x+82.05图像三维重建的主流方法，这是因为其操作简单、具有高自动化和鲁棒性的优点2 。相比之下，虽然空间雕刻算法能取得更加精确的重建结果，但是它需要提前进行相机标定，实现过程更为复杂2 3。因此,在未来的研究中，可以考虑结合SFM方法和空间雕刻方法，在不借助标定板的情况下实现水稻的三维重建。J=X152此外，植物冠层的体积特征对于指示植株生物量信息和生长状况具有重要意义。针对植株体积参数计算的需求，李庆等2 4 证明Alpha-shape算法可用于计算冠层结构参数，杨征鹤等2 5 基于该算法成功实现了对番茄植株的体积计算，与凸包法相比平均误差降低了15.3%，越来越多的研究证明了该算法在植株三维模型体积计算方面的有效性。或许是因为水稻生长情况复杂，虽然吴丹等12 成功重建了水稻三维模型，但是并未考虑基于水稻模型实现水稻体积和生物量的估测。本研究使用Alpha-shape将体素模型进一步改进获得了曲面模型，从而提取了水稻模型的体积并用于估测水稻生物量，按照抽穗期将水稻分为2 类，讨论是否含穗对估测的影响。研究结果显示，针对抽穗期以前的水稻，预测的生物量较为准确，水稻体积与生物量人工测量值的R达到了0.8 3。然而，对于抽穗期以后水稻，水稻体积与生物量人工测量值的R?为0.7 7,预测误差RMSE和MAPE也更高。如前所述，这是由于不同的水稻植株穗质量比例和生长情况各不相同导致。本研究存在一些局限性：首先，由于计算机内存和处理器性能的限制，本研究在仅使用CPU进行建模时耗时较长，难以投入大规模应用；其次，本研究在使用体积参数估测水稻生物量时，仅考虑了水稻是否含穗对结果的影响，而由于水稻生长情况的复杂性，植株之间存在较大差异，因此在综合考虑其他因素之后或许会得到更准确的估测值；最后，本研究仅提取了水稻株高、茎粗、植被覆盖度和体积等表型参数，而如叶倾角、展开度、植株骨架等更多的表型参数的研究还有待进一步开展。在后续研究中，将考虑采用GPU加速和优化体素存储结构等方式来提高算法的效率，同时还将探索提取和研究更多水稻表型参数，以便更全面地了解水稻的生长形态和特性。综上所述，本研究通过空间雕刻方法重建了水稻高精度的三维模型，并基于此提取出了更全面的水稻表型参数，通过体积参数成功预测了水稻生物量，从三维层面为水稻表型研究提供了参考。此外，本研究发现，空间雕刻算法对重建物体没有限制，理论上适用于任何对象的重建，具有广泛的应用价值，因此对如小麦、高梁等缺乏特征纹理信息的禾本科农作物重建均具有参考价值。河南农业科学参考文献：1孙艳华.水稻高产栽培技术与病虫害防控J.江西农业,2 0 2 2(2 4):2-5.SUN Y H.Cultivation techniques for high-yielding riceand prevention and control of diseases and pests J.Jiangxi Agriculture,2022(24):2-5.2朱业宝，陈立喆，张丹，等.福建省水稻地方品种表型性状遗传多样性分析J.福建农业学报，2 0 2 1，36（10)：1119-1125.ZHU Y B,CHEN L Z,ZHANG D,et al.Genetic diversityon phenotypic traits of rice landraces in Fujian J.FujianJournal of Agricultural Sciences，2 0 2 1,36（10):1119-1125.3万亮，杜晓月，陈硕博，等.基于无人机多源图谱融合的水稻稻穗表型监测J.农业工程学报，2 0 2 2，38（9）：162-170.WANL,DUXY,CHENSB,et al.Rice paniclephenotyping using UAV-based multi-source spectralimage data fusion J.Transactions of the Chinese Societyof Agricultural Engineering,2022,38(9):162-170.4魏天翔，汪小邑，施印炎，等.基于三维重建的秸秆地单株水稻生长形态检测研究.南京农业大学学报，2021,44(6):1197-1208.WEI T X,WANG X C,SHI Y Y,et al.Detection of singlerice early growth morphol