You will be provided with a research report. The body of the report will contain some citations to references.

Citations in the main text may appear in the following forms:
1. A segment of text + space + number, for example: "Li Qiang constructed a socioeconomic status index (SES) based on income, education, and occupation, dividing society into 7 levels 15"
2. A segment of text + [number], for example: "Li Qiang constructed a socioeconomic status index (SES) based on income, education, and occupation, dividing society into 7 levels[15]"
3. A segment of text + [number†(some line numbers, etc.)], for example: "Li Qiang constructed a socioeconomic status index (SES) based on income, education, and occupation, dividing society into 7 levels[15†L10][5L23][7†summary]"
4. [Citation Source](Citation Link), for example: "According to [ChinaFile: A Guide to Social Class in Modern China](https://www.chinafile.com/reporting-opinion/media/guide-social-class-modern-china)'s classification, Chinese society can be divided into nine strata"

Please identify **all** instances where references are cited in the main text, and extract (fact, ref_idx, url) triplets. When extracting, pay attention to the following:
1. Since these facts will need to be verified later, you may need to look for some context before and after the citation to ensure that the fact is complete and understandable, rather than just a simple phrase or short expression.
2. If a fact cites multiple references, then it should correspond to two triplets: (fact, ref_idx_1, url_1) and (fact, ref_idx_2, url_2).
3. For the third form of citation (i.e., where the citation source and link appear directly in the text), the ref_idx should be uniformly set to 0.
4. If the main text does not specify the exact location of the citation (for example, only the reference list is listed at the end of the article, without specifying the citation point in the text), please return an empty list.

You should return a JSON list format, where each item in the list is a triplet, for example:
[
    {
        "fact": "Text segment from the original document. Note that Chinese quotation marks should use full-width marks. And add a single backslash before the English quotation mark to make it a readable for python json module.",
        "ref_idx": "The index of the cited reference in the reference list for this text segment.",
        "url": "The URL of the cited reference for this text segment (extracted from the reference list at the end of the research report or from the parentheses at the citation point)."
    }
]

Here is the main text of the research report:
# Precision Mechatronic Design and Modern Control Architecture for High-Throughput 3D Grain Phenotyping and Morphological Analysis

## Architectural Formulation and Mechatronic Hardware Integration

High-throughput, non-destructive phenotypic profiling of individual crop seeds—including bread wheat (*Triticum aestivum*), rice (*Oryza sativa*), and maize (*Zea mays*)—constitutes a foundational imperative for modern genomics-assisted breeding, seed sorting, and yield evaluation [1]. Traditional two-dimensional machine vision platforms measure projected geometric boundaries such as major axis length, minor axis width, and projected footprint area [1]. However, planar projections are fundamentally incapable of quantifying three-dimensional volumetric parameters, volumetric mass densities, and sub-millimeter surface recessions [1]. Among these morphological traits, the wheat ventral sulcus (crease) represents a decisive agronomic parameter, directly governing milling flour extraction yields, pathogen resistance, and endosperm development dynamics [1].

Extracting such fine surface topographies requires active micron-level three-dimensional reconstruction, which imposes stringent demands on the underlying mechatronic presentation mechanisms [1]. Automated high-resolution phenotyping platforms resolve these physical constraints by unifying three electro-mechanical sub-systems into a synchronized execution cell: a direct-drive rotary indexing stage that presents the specimen across panoramic viewing perspectives, an automated robotic manipulator that singulates and places individual grains onto the rotational datum, and a sensory optoelectronic core incorporating structured-light projectors, telecentric machine vision cameras, and strobe illumination arrays [1].

The robotic handling subsystem manages the physical transfer of individual grains from bulk seed hoppers to the imaging datum without inducing pericarp scarification or structural damage [3]. While conventional pneumatic suction systems employ high vacuum pressures for pickup, rapid pneumatic exhaust cycles generate severe kinetic transients [7]. For example, in an exhaust sequence operating under an overpressure of \(30\text{ kPa}\) across a \(75\text{ }\mu\text{m}\) radius nozzle orifice, a lightweight seed (\(m \approx 20\text{ }\mu\text{g}\) for *Arabidopsis* or \(m \approx 30\text{–}45\text{ mg}\) for cereal caryopses) experiences instantaneous release accelerations exceeding \(2500\text{ m/s}^2\), triggering erratic rebound and spatial displacement off the scanning axis [7].

To overcome this handling dynamic, the mechanical end-effector incorporates a gravity-assisted internal plunger pin assembly, commonly termed a mechanical hammer release mechanism [7]. By breaking the negative pressure differential and actuating a sliding needle through the nozzle core, the terminal release velocity of the seed is constrained below \(0.2\text{ m/s}\), ensuring stable spatial placement onto the rotary measurement axis and enabling non-destructive multi-seed phenotyping pipelines [7].

The optical imaging payload is configured around a twin-camera binocular stereo geometry or multi-view array operating in conjunction with a digital micro-mirror device (DMD) structured-light projector or a 19-band integrating sphere dome [1]. To eliminate magnification variability arising from seed topographic variation and minor axial runout, the optical trains employ bilateral telecentric lenses [1]. The centralized master controller coordinates the spatial trajectory of the direct-drive permanent magnet synchronous motor (PMSM) with hardware-level digital transistor-transistor logic (TTL) triggers, simultaneously firing camera exposure shutters and high-current LED strobe drivers [8].

| Reconstruction Modality | Spatial Resolution | Typical Throughput (Seeds/Hour) | Destructive Nature | Crease / Internal Visibility | Volumetric Accuracy (\(R^2\) vs True Mass) | Reference System / Benchmarks |
| --- | --- | --- | --- | --- | --- | --- |
| **X-ray Micro-Computed Tomography (\(\mu\text{CT}\))** | \(1\text{–}10\text{ }\mu\text{m}\) | \(5\text{–}15\) | Non-destructive | External crease and internal endosperm/air cavities visible | \(0.98\text{–}0.99\) | Benchtop Micro-CT scanners [1] |
| **Monocular Structure-from-Motion (SfM) + Turntable** | \(50\text{–}150\text{ }\mu\text{m}\) | \(30\text{–}60\) | Non-destructive | Crease partially occluded; relies on texture contrast | \(0.85\text{–}0.91\) | Multi-view photogrammetry rigs [2] |
| **Phase-Shifting Structured Light Profilometry** | \(10\text{–}50\text{ }\mu\text{m}\) | \(60\text{–}120\) | Non-destructive | High accuracy on external crease profiles; blind to inner core | \(0.96\text{–}0.99\) | Reeyee Pro / Industrial binocular structured light [1] |
| **Robotic Pick-and-Place Optical Silhouette Carving** | \(30\text{–}80\text{ }\mu\text{m}\) | \(100\text{–}180\) | Non-destructive | Cannot capture concave recesses (sulcus lost to visual hull) | \(0.88\text{–}0.94\) | PhenoSeeder robotic workstation [3] |
| **Turntable-Bound 3D Gaussian Splatting (3DGS)** | \(20\text{–}60\text{ }\mu\text{m}\) | \(40\text{–}80\) | Non-destructive | Detailed surface geometry with high photometric consistency | \(0.93\text{–}0.97\) | SC-NeRF / Circular Orbit 3DGS platforms [2] |

## Mathematical Modeling of Physical Mechatronic Dynamics

### Rotary Positioning Platform Dynamics with Two-Inertia Compliance

Direct-drive rotary tables minimize mechanical backlash by eliminating intermediate gearboxes; however, high-acceleration indexing inevitably excites torsional structural resonances between the motor rotor and the specimen platter [12]. Modeling the mechanical assembly as a lumped-parameter two-inertia flexible torsional system allows exact analytical characterization of high-frequency structural jitter and settling dynamics [12].

The differential equations governing angular motion for the actuator rotor and the specimen platter are formulated as:

\[J_m \ddot{\theta}_m(t) + d_s (\dot{\theta}_m(t) - \dot{\theta}_l(t)) + k_s (\theta_m(t) - \theta_l(t)) = \tau_m(t) - \tau_{f,m}(t)\]

\[J_l \ddot{\theta}_l(t) + d_s (\dot{\theta}_l(t) - \dot{\theta}_m(t)) + k_s (\theta_l(t) - \theta_m(t)) = -\tau_{f,l}(t) - \tau_{ext}(t)\]

where \(J_m\) and \(J_l\) designate the mass moments of inertia of the drive rotor and specimen platter (\(\text{kg}\cdot\text{m}^2\)), \(\theta_m\) and \(\theta_l\) represent their respective angular coordinates (\(\text{rad}\)), \(k_s\) denotes the torsional torsional shaft stiffness (\(\text{N}\cdot\text{m/rad}\)), and \(d_s\) represents the internal torsional damping coefficient (\(\text{N}\cdot\text{m}\cdot\text{s/rad}\)). The term \(\tau_m\) represents the electromagnetic drive torque generated by the stator windings, while \(\tau_{f,m}\) and \(\tau_{f,l}\) denote nonlinear bearing friction torques acting on the rotor and load assemblies [12]. Exogenous torque perturbations, including aerodynamic drag and mechanical mass eccentricities from uncentered seeds, are grouped into \(\tau_{ext}\).

### Nonlinear Dynamic Friction Model

At micro-radian displacements and low scanning velocities, bearing pre-load introduces severe pre-sliding elastoplastic deformation, hysteresis, and velocity-softening Stribeck phenomena that cannot be captured by linear viscous damping approximations [12]. The LuGre friction model parameterizes these phenomena by treating microscopic surface asperities as elastic bristles whose average physical deflection is denoted by the internal state \(z\) [12]:

\[\frac{dz}{dt} = \omega_{rel} - \frac{\sigma_0 |\omega_{rel}|}{g(\omega_{rel})} z\]

$$\tau_f = \sigma_0 z + \sigma_1 \frac{dz}{dt} + \sigma_2 \omega_{rel}$$where \(\omega_{rel} = \dot{\theta}_m\) (or \(\dot{\theta}_l\)) designates the relative angular velocity across the bearing interface, \(\sigma_0\) denotes the asperity bristle stiffness (\(\text{N}\cdot\text{m/rad}\)), \(\sigma_1\) is the microscopic damping coefficient (\(\text{N}\cdot\text{m}\cdot\text{s/rad}\)), and \(\sigma_2\) parameterizes empirical viscous lubrication effects (\(\text{N}\cdot\text{m}\cdot\text{s/rad}\)). The velocity-dependent Stribeck function \(g(\omega_{rel})\) establishes the transition boundary from static breakaway to dynamic Coulomb sliding:

\[g(\omega_{rel}) = T_c + (T_s - T_c) \exp\left(-\left(\frac{\omega_{rel}}{\omega_s}\right)^2\right)\]

where \(T_c\) represents the Coulomb friction torque, \(T_s\) is the static breakaway torque threshold (\(T_s > T_c\)), and \(\omega_s\) denotes the characteristic Stribeck sliding velocity.

### Strobe Lighting Driver Dynamics and Radiometric Modulation

The spatial fidelity of structured light phase profilometry and multi-view pixel correspondences depends strictly on photometric consistency across image sequences [1]. When driving high-power light-emitting diode (LED) arrays continuously, rapid junction heating induces thermal droop, resulting in exponential reductions in radiant flux and spectral shifts across acquisition cycles [8]. Driving the illuminator in pulsed over-drive synchronization with camera exposure windows circumvents junction overheating while generating high instantaneous irradiance, allowing microsecond exposure times that freeze structural vibration [8].

The electrical drive dynamics of the solid-state LED illuminator are modeled as a first-order inductive-resistive circuit governed by high-frequency pulse-width modulation (PWM):

\[L_d \frac{di_f(t)}{dt} + R_d i_f(t) = v_{pwm}(t) - v_f(i_f)$$where $i_f(t)$ represents the instantaneous forward diode current, $L_d$ and $R_d$ represent parasitic trace inductance and dynamic driver resistance, and $v_f$ is the non-linear forward junction potential. The resulting optical radiant flux $\Phi_e(t)$ is directly coupled to the forward drive current above the emission threshold $I_{th}$ and inversely modulated by junction temperature $T_j$:$$\Phi_e(t) = \eta_{ext}(T_j) \cdot K_{opt} \cdot [i_f(t) - I_{th}]\]

\[\frac{dT_j(t)}{dt} = -\frac{1}{R_{th}C_{th}} (T_j(t) - T_{amb}) + \frac{P_{diss}(t)}{C_{th}}\]

where \(P_{diss}(t) = v_f(t) i_f(t) - \Phi_e(t)\) designates internal electrical power dissipation, \(R_{th}\) represents the thermal junction-to-case resistance, and \(C_{th}\) is the equivalent lumped thermal capacitance. Regulating \(i_f(t)\) via closed-loop constant-current drive electronics operating above \(100\text{ kHz}\) guarantees intra-frame illumination stability, preventing phase demodulation errors during structured-light projection [1].

### Unified State-Space Formulation

Linearizing the multi-body mechatronic equations about a steady-state rotational velocity \(\omega_0\) yields a continuous-time state-space representation. Defining the state vector as \(\mathbf{x}(t) = [\theta_m(t), \dot{\theta}_m(t), \theta_l(t), \dot{\theta}_l(t)]^T \in \mathbb{R}^4\), the scalar control input as \(u(t) = \tau_m(t) \in \mathbb{R}\), and the unmodeled disturbance vector as \(\mathbf{d}(t) = [\tau_{f,m}(t), \tau_{f,l}(t) + \tau_{ext}(t)]^T \in \mathbb{R}^2\), the system is formulated as:

\[\dot{\mathbf{x}}(t) = \mathbf{A}\mathbf{x}(t) + \mathbf{B}u(t) + \mathbf{E}\mathbf{d}(t)\]

\[\mathbf{y}(t) = \mathbf{C}\mathbf{x}(t) + \mathbf{D}u(t) + \mathbf{w}(t)\]

The system, input, and disturbance coupling matrices expand directly to:

\[\mathbf{A} = \begin{bmatrix}
0 & 1 & 0 & 0 \\
-\frac{k_s}{J_m} & -\frac{d_s + \sigma_{2,m}}{J_m} & \frac{k_s}{J_m} & \frac{d_s}{J_m} \\
0 & 0 & 0 & 1 \\
\frac{k_s}{J_l} & \frac{d_s}{J_l} & -\frac{k_s}{J_l} & -\frac{d_s + \sigma_{2,l}}{J_l}
\end{bmatrix}, \quad
\mathbf{B} = \begin{bmatrix} 0 \\ \frac{1}{J_m} \\ 0 \\ 0 \end{bmatrix}, \quad
\mathbf{E} = \begin{bmatrix}
0 & 0 \\
-\frac{1}{J_m} & 0 \\
0 & 0 \\
0 & -\frac{1}{J_l}
\end{bmatrix}\]

For an instrumentation setup incorporating a high-resolution primary optical ring encoder mounted on the motor shaft and an auxiliary secondary glass scale on the specimen load axis, the observation model is:

\[\mathbf{C} = \begin{bmatrix}
1 & 0 & 0 & 0 \\
0 & 0 & 1 & 0
\end{bmatrix}, \quad
\mathbf{D} = \begin{bmatrix} 0 \\ 0 \end{bmatrix}\]

where \(\mathbf{w}(t) \sim \mathcal{N}(\mathbf{0}, \mathbf{R}_n)\) represents zero-mean Gaussian optoelectronic measurement noise and encoder quantization uncertainty.

## Modern Control Theory: Controllability, Observability, and Stability

### Controllability Analysis

Controllability establishes whether the electromagnetic torque generated at the motor stator can arbitrarily steer all system states—including the unactuated specimen platter modes—to target angular coordinates in finite time. The controllability matrix \(\mathcal{C}\) is constructed as:

\[\mathcal{C} = \begin{bmatrix} \mathbf{B} & \mathbf{A}\mathbf{B} & \mathbf{A}^2\mathbf{B} & \mathbf{A}^3\mathbf{B} \end{bmatrix}\]

Expanding the matrix columns yields:

\[\mathbf{B} = \begin{bmatrix} 0 \\ \frac{1}{J_m} \\ 0 \\ 0 \end{bmatrix}, \quad
\mathbf{A}\mathbf{B} = \begin{bmatrix} \frac{1}{J_m} \\ -\frac{d_s + \sigma_{2,m}}{J_m^2} \\ 0 \\ \frac{d_s}{J_m J_l} \end{bmatrix}, \quad
\mathbf{A}^2\mathbf{B} = \begin{bmatrix} -\frac{d_s + \sigma_{2,m}}{J_m^2} \\ \frac{(d_s+\sigma_{2,m})^2 - k_s J_m}{J_m^3} + \frac{k_s d_s}{J_m^2 J_l} \\ \frac{d_s}{J_m J_l} \\ \frac{k_s J_m - d_s(d_s + \sigma_{2,m})}{J_m^2 J_l} - \frac{d_s(d_s + \sigma_{2,l})}{J_m J_l^2} \end{bmatrix}$$Computing the determinant of $\mathcal{C}$ reveals:$$\det(\mathcal{C}) = \frac{k_s^2}{J_m^4 J_l^2}\]

Because the structural torsional stiffness \(k_s\) and moments of inertia \(J_m, J_l\) are strictly positive real numbers, \(\det(\mathcal{C}) \neq 0\), which establishes that \(\text{rank}(\mathcal{C}) = 4\). The system is fully controllable across the entire state space. This mathematical property proves that platter-side torsional oscillations can be actively suppressed entirely through dynamic current modulation of the motor stator windings, eliminating the necessity for bulky physical damping mechanisms on the rotating platter.

### Observability Analysis

To evaluate cost-constrained hardware configurations where an optical encoder is installed solely on the drive motor shaft, the observation matrix is reduced to the single-output row vector \(\mathbf{C}_m = [1 \ 0 \ 0 \ 0]\). The corresponding observability matrix \(\mathcal{O}_m\) is evaluated as:

\[\mathcal{O}_m = \begin{bmatrix} \mathbf{C}_m \\ \mathbf{C}_m \mathbf{A} \\ \mathbf{C}_m \mathbf{A}^2 \\ \mathbf{C}_m \mathbf{A}^3 \end{bmatrix} = 
\begin{bmatrix}
1 & 0 & 0 & 0 \\
0 & 1 & 0 & 0 \\
-\frac{k_s}{J_m} & -\frac{d_s + \sigma_{2,m}}{J_m} & \frac{k_s}{J_m} & \frac{d_s}{J_m} \\
\gamma_{41} & \gamma_{42} & \gamma_{43} & \gamma_{44}
\end{bmatrix}$$Evaluating the linear independence of the rows yields:$$\det(\mathcal{O}_m) = \frac{k_s^2}{J_m^2 J_l} \neq 0 \quad (\forall k_s > 0)\]

Consequently, \(\text{rank}(\mathcal{O}_m) = 4\). This confirms that the internal states of the specimen stage—including angular deflection and load velocity—remain fully observable through motor-side encoder measurements via the elastic torque transmission across the drive shaft. This mathematical condition guarantees that dynamic state observers and Kalman filters will reconstruct the complete unmeasured structural state.

### Lyapunov Stability Analysis with Nonlinear Dynamic Friction

To assess closed-loop stability in the presence of nonlinear LuGre bearing friction, consider the radially unbounded positive-definite Lyapunov candidate function \(V(\mathbf{x}, z)\):

\[V(\mathbf{x}, z) = \frac{1}{2} J_m \dot{\theta}_m^2 + \frac{1}{2} J_l \dot{\theta}_l^2 + \frac{1}{2} k_s (\theta_m - \theta_l)^2 + \frac{1}{2} \sigma_0 z^2\]

Differentiating \(V(\mathbf{x}, z)\) with respect to time along the continuous system trajectories:

\[\dot{V}(\mathbf{x}, z) = J_m \dot{\theta}_m \ddot{\theta}_m + J_l \dot{\theta}_l \ddot{\theta}_l + k_s (\theta_m - \theta_l)(\dot{\theta}_m - \dot{\theta}_l) + \sigma_0 z \dot{z}$$Substituting the mechanical acceleration terms alongside a state-feedback control law $u = \tau_m = -k_p \theta_m - k_d \dot{\theta}_m$:$$\dot{V} = \dot{\theta}_m \left[ \tau_m - d_s(\dot{\theta}_m - \dot{\theta}_l) - k_s(\theta_m - \theta_l) - \tau_{f,m} \right] + \dot{\theta}_l \left[ -d_s(\dot{\theta}_l - \dot{\theta}_m) - k_s(\theta_l - \theta_m) - \tau_{f,l} \right] + k_s(\theta_m - \theta_l)(\dot{\theta}_m - \dot{\theta}_l) + \sigma_0 z \left( \dot{\theta}_l - \frac{\sigma_0 |\dot{\theta}_l|}{g(\dot{\theta}_l)} z \right)\]

Collecting terms and simplifying the structural compliance interactions produces:

\[\dot{V} = -k_d \dot{\theta}_m^2 - d_s (\dot{\theta}_m - \dot{\theta}_l)^2 - \sigma_{2,m} \dot{\theta}_m^2 - \sigma_{2,l} \dot{\theta}_l^2 - \sigma_1 \left(\frac{dz}{dt}\right)^2 - \frac{\sigma_0^2 |\dot{\theta}_l|}{g(\dot{\theta}_l)} z^2\]

Because \(g(\dot{\theta}_l) \ge T_c > 0\) and the physical coefficients \(k_d, d_s, \sigma_0, \sigma_1, \sigma_2\) are strictly non-negative, the time derivative satisfies \(\dot{V}(\mathbf{x}, z) \le 0\), rendering the system negative semi-definite. Applying LaSalle's invariance principle confirms that the largest invariant set where \(\dot{V} = 0\) corresponds strictly to the equilibrium origin \(\mathbf{x} = \mathbf{0}, z = 0\). Hence, the autonomous closed-loop system is globally asymptotically stable.

## Optimal, Robust, and Predictive Control Synthesis

### Linear Quadratic Gaussian (LQG) Synthesis with Disturbance Estimation

To achieve rapid step-and-settle positioning while penalizing mechanical stress and high-frequency current jitter, an optimal Linear Quadratic Regulator (LQR) is synthesized over an infinite horizon:

\[J = \int_0^\infty \left( \mathbf{x}^T(t) \mathbf{Q}_{lqr} \mathbf{x}(t) + u(t) R_{lqr} u(t) \right) dt\]

The state weighting matrix \(\mathbf{Q}_{lqr} \succeq 0\) heavily penalizes the angular tracking error of the specimen stage and relative torsional twisting, while \(R_{lqr} > 0\) bounds motor current demands:

\[\mathbf{Q}_{lqr} = \operatorname{diag}\left( q_{\theta_m}, 0, q_{\theta_l}, q_{\omega_l} \right)$$The optimal state feedback gain matrix is computed as $\mathbf{K}_{lqr} = R_{lqr}^{-1} \mathbf{B}^T \mathbf{P}_{are}$, where $\mathbf{P}_{are}$ represents the positive-definite solution to the continuous-time Algebraic Riccati Equation:$$\mathbf{A}^T \mathbf{P}_{are} + \mathbf{P}_{are} \mathbf{A} - \mathbf{P}_{are} \mathbf{B} R_{lqr}^{-1} \mathbf{B}^T \mathbf{P}_{are} + \mathbf{Q}_{lqr} = \mathbf{0}\]

Because load-side states are unmeasured in single-encoder configurations, a steady-state Kalman-Bucy filter reconstructs the full state trajectory \(\hat{\mathbf{x}}(t)\) while simultaneously estimating an augmented low-frequency disturbance state \(\hat{d}_L(t)\) representing bearing friction and aerodynamic torque offsets:

\[\begin{bmatrix} \dot{\hat{\mathbf{x}}}(t) \\ \dot{\hat{d}}_L(t) \end{bmatrix} =  \begin{bmatrix} \mathbf{A} & \mathbf{E}_1 \\ \mathbf{0} & 0 \end{bmatrix} \begin{bmatrix} \hat{\mathbf{x}}(t) \\ \hat{d}_L(t) \end{bmatrix} +  \begin{bmatrix} \mathbf{B} \\ 0 \end{bmatrix} u(t) +  \mathbf{L}_k \left( \mathbf{y}(t) - \mathbf{C}_a \begin{bmatrix} \hat{\mathbf{x}}(t) \\ \hat{d}_L(t) \end{bmatrix} \right)\]

where \(\mathbf{L}_k = \mathbf{P}_e \mathbf{C}_a^T \mathbf{R}_n^{-1}\) is derived from the observer Algebraic Riccati Equation parameterized by process disturbance covariance \(\mathbf{Q}_w\) and sensor noise covariance \(\mathbf{R}_n\). The resulting robust control law integrates dynamic state feedback with feedforward disturbance cancellation:

\[u(t) = -\mathbf{K}_{lqr} \hat{\mathbf{x}}(t) - \hat{d}_L(t)\]

### Constrained Model Predictive Control (MPC)

During rapid angular steps (e.g., \(10^\circ\text{–}30^\circ\) increments between multi-view captures) [1], aggressive angular acceleration can displace delicate grain specimens resting freely on the platter. The maximum allowable acceleration is bounded by the static friction cone between the grain pericarp and the stage surface:

\[m r \dot{\omega}_l \le \mu_s m g \implies \vert{}\ddot{\theta}_l\vert{} \le \frac{\mu_s g}{r_{\max}}\]

where \(\mu_s \approx 0.35\text{–}0.55\) represents the coefficient of static friction for crop seed coats, and \(r_{\max} \approx 25\text{ mm}\) is the platform radius.

The discrete-time MPC solves the following optimization problem over prediction horizon \(N_p\) and control horizon \(N_c\) at each sampling interval \(T_s\):

\[\min_{\Delta u(k)\dots \Delta u(k+N_c-1)} \sum_{j=1}^{N_p} \Vert{}\hat{\mathbf{x}}(k+j\vert{}k) - \mathbf{x}_{ref}(k+j)\Vert{}_{\mathbf{Q}_{mpc}}^2 + \sum_{j=0}^{N_c-1} \Vert{}\Delta u(k+j)\Vert{}_{R_{mpc}}^2\]

subject to:

\[\mathbf{x}(k+j+1) = \mathbf{A}_d \mathbf{x}(k+j) + \mathbf{B}_d u(k+j) \\ \vert{}\tau_m(k+j)\vert{} \le \tau_{\max} \quad (\text{actuator torque saturation}) \\ \vert{}\Delta \tau_m(k+j)\vert{} \le \Delta \tau_{\max} \quad (\text{slew-rate limiting})\]

\[\vert{}\ddot{\theta}_l(k+j)\vert{} \le \frac{\mu_s g}{r_{\max}} \quad (\text{grain slip constraint})\]

Solving this quadratic program online via an interior-point solver guarantees minimal settling time without specimen shift or mechanical detent ringing.

### Feedforward Vibration Suppression via Input Shaping

To eliminate residual structural vibration without introducing closed-loop phase lags during high-throughput stop-and-strobe acquisition, Zero Vibration Derivative (ZVD) input shaping filters the trajectory setpoints [14]. The primary structural torsional resonance frequency \(\omega_n\) and damping ratio \(\zeta\) are extracted from system identification:

\[\omega_n = \sqrt{\frac{k_s (J_m + J_l)}{J_m J_l}}, \quad \zeta = \frac{d_s}{2} \sqrt{\frac{J_m + J_l}{k_s J_m J_l}}\]

The ZVD shaper consists of three non-negative impulses \(A_i\) located at time instances \(t_i\):$$f_{ZVD}(t) = \sum_{i=1}^{3} A_i \delta(t - t_i)$$Setting derivatives of the residual vibration amplitude equation with respect to \(\omega\) to zero yields parameters:

\[K = \exp\left( -\frac{\zeta \pi}{\sqrt{1 - \zeta^2}} \right), \quad \Delta T = \frac{\pi}{\omega_n \sqrt{1 - \zeta^2}} \\ A_1 = \frac{1}{(1 + K)^2}, \quad A_2 = \frac{2K}{(1 + K)^2}, \quad A_3 = \frac{K^2}{(1 + K)^2} \\ t_1 = 0, \quad t_2 = \Delta T, \quad t_3 = 2\Delta T\]

Convolving \(f_{ZVD}(t)\) with the reference command profile \(r(t)\) produces a shaped input that cancels residual structural vibrations at \(\omega_n\), eliminating motion-induced image blur during high-speed step-and-scan indexing [14].

| Motion Control Architecture | 5% Settling Time (\(30^\circ\) Step) | Platter Torsional Overshoot | Peak Tracking Error | Vibration Attenuation at \(\omega_n\) | Computational Burden | Slip Prevention Enforcement |
| --- | --- | --- | --- | --- | --- | --- |
| **Decentralized PID + Notch Filter** | \(145\text{ ms}\) | \(12.4\%\) | \(2.40\text{ mrad}\) | \(-18\text{ dB}\) | Low (\(<10\text{ }\mu\text{s}\)) | Unconstrained (Empirical tuning) |
| **Full-State LQR** | \(78\text{ ms}\) | \(3.1\%\) | \(0.85\text{ mrad}\) | \(-28\text{ dB}\) | Low (\(<20\text{ }\mu\text{s}\)) | Unconstrained |
| **LQG + Disturbance Observer** | \(52\text{ ms}\) | \(0.8\%\) | \(0.22\text{ mrad}\) | \(-34\text{ dB}\) | Moderate (\(50\text{ }\mu\text{s}\)) | Unconstrained |
| **ZVD Input-Shaped LQG** | \(36\text{ ms}\) | \(<0.1\%\) | \(0.05\text{ mrad}\) | \(-52\text{ dB}\) | Moderate (\(65\text{ }\mu\text{s}\)) | Profile pre-smoothing |
| **Explicit Constrained MPC** | \(41\text{ ms}\) | \(<0.1\%\) | \(0.03\text{ mrad}\) | \(-48\text{ dB}\) | High (\(1.2\text{ ms}\)) | Strict quadratic inequality constraint |

## Uncertainty Quantification and Error Propagation: Control to 3D Geometry

### Optical Triangulation and Geometric Runout Modeling

The stage rotation vector \(\mathbf{n} \in \mathbb{R}^3\) (\(\Vert{}\mathbf{n}\Vert{} = 1\)) and center of rotation \(\mathbf{C}_0 \in \mathbb{R}^3\) are calibrated relative to the world coordinate frame [18]. An idealized rotation by an angle \(\theta_i\) is expressed using Rodrigues' rotation formula:

\[\mathbf{R}(\theta_i, \mathbf{n}) = \mathbf{I} + \sin\theta_i [\mathbf{n}]_\times + (1 - \cos\theta_i) [\mathbf{n}]_\times^2\]

where \([\mathbf{n}]_\times\) is the skew-symmetric cross-product operator. The spatial projection of a 3D point \(\mathbf{X}_j \in \mathbb{R}^3\) onto the image plane of camera \(c\) at rotation step \(i\) is parameterized as:

\[\tilde{\mathbf{x}}_{c,i,j} = \pi(\mathbf{K}_c, \mathbf{R}_c, \mathbf{t}_c, \theta_i, \mathbf{X}_j) = \mathbf{K}_c \left( \mathbf{R}_c \left( \mathbf{R}(\theta_i, \mathbf{n})(\mathbf{X}_j - \mathbf{C}_0) + \mathbf{C}_0 \right) + \mathbf{t}_c \right)$$where $\mathbf{K}_c$ represents the intrinsic calibration matrix:$$\mathbf{K}_c = \begin{bmatrix} f_x & s & u_0 \\ 0 & f_y & v_0 \\ 0 & 0 & 1 \end{bmatrix}\]

Physical bearings exhibit dynamic runout errors comprising radial runout \(\delta \mathbf{r}(t) \in \mathbb{R}^2\), axial runout \(\delta z(t) \in \mathbb{R}\), and dynamic angular wobble (tilt error) \(\delta \boldsymbol{\phi}(t) = [\delta \phi_x, \delta \phi_y]^T\) [20]. The true rotation matrix perturbed by dynamic wobble is:

\[\mathbf{R}_{true}(\theta_i) = \exp([\delta \boldsymbol{\phi}_i]_\times) \mathbf{R}(\theta_i, \mathbf{n})\]

### First-Order Covariance Propagation to 3D Metric Uncertainty

Deviations in camera pose, angular positioning errors \(\delta \theta \sim \mathcal{N}(0, \sigma_\theta^2)\), mechanical bearing wobble \(\mathbf{\Sigma}_{wobble} \in \mathbb{R}^{3 \times 3}\), and sensor pixel discretization \(\mathbf{\Sigma}_{pix} = \sigma_p^2 \mathbf{I}_2\) propagate directly into the reconstructed 3D surface point cloud covariance \(\mathbf{\Sigma}_{\mathbf{X}_j}\) [22].

Expanding the observation equation for a reconstructed feature point across \(M\) visual frames:

\[\mathbf{r}_{ij} = \mathbf{x}_{ij}^{meas} - \pi_i(\mathbf{X}_j, \theta_i, \boldsymbol{\xi}_i) = \mathbf{0}\]

where \(\boldsymbol{\xi}_i = [\delta \mathbf{r}_i^T, \delta z_i, \delta \boldsymbol{\phi}_i^T]^T\) denotes the mechanical stage runout vector. The first-order Taylor expansion yields:

\[\delta \mathbf{x}_{ij} \approx \mathbf{J}_{\mathbf{X},ij} \delta \mathbf{X}_j + \mathbf{J}_{\theta,ij} \delta \theta_i + \mathbf{J}_{\boldsymbol{\xi},ij} \boldsymbol{\xi}_i + \mathbf{w}_{p,ij}$$where:$$\mathbf{J}_{\mathbf{X},ij} = \frac{\partial \pi_i}{\partial \mathbf{X}_j} \in \mathbb{R}^{2 \times 3}, \quad  \mathbf{J}_{\theta,ij} = \frac{\partial \pi_i}{\partial \theta_i} \in \mathbb{R}^{2 \times 1}, \quad  \mathbf{J}_{\boldsymbol{\xi},ij} = \frac{\partial \pi_i}{\partial \boldsymbol{\xi}_i} \in \mathbb{R}^{2 \times 5}\]

Aggregating across all \(M\) viewing positions, the composite observation covariance matrix \(\mathbf{\Sigma}_{meas,ij} \in \mathbb{R}^{2 \times 2}\) is:

\[\mathbf{\Sigma}_{meas,ij} = \mathbf{J}_{\theta,ij} \sigma_\theta^2 \mathbf{J}_{\theta,ij}^T + \mathbf{J}_{\boldsymbol{\xi},ij} \mathbf{\Sigma}_{runout} \mathbf{J}_{\boldsymbol{\xi},ij}^T + \sigma_p^2 \mathbf{I}_2\]

The posterior reconstruction covariance matrix for spatial point \(\mathbf{X}_j\) is given by the inverse Fisher Information Matrix:

\[\mathbf{\Sigma}_{\mathbf{X}_j} = \left( \sum_{i=1}^{M} \mathbf{J}_{\mathbf{X},ij}^T \mathbf{\Sigma}_{meas,ij}^{-1} \mathbf{J}_{\mathbf{X},ij} \right)^{-1}\]

To resolve the sub-millimeter features of the wheat ventral sulcus (\(<200\text{ }\mu\text{m}\) width, \(100\text{–}500\text{ }\mu\text{m}\) depth) [1], the three-dimensional positional uncertainty ellipsoid must not exceed \(\sigma_{\max} = 15\text{ }\mu\text{m}\). Setting the camera baseline distance at \(D = 120\text{ mm}\) with optical focal length \(f = 35\text{ mm}\) and pixel pitch \(p = 3.45\text{ }\mu\text{m}\), the directional variance simplifies to:

\[\sigma_{X_z}^2 \approx \frac{D^2}{f^2 B^2} \left( \sigma_p^2 + f^2 \sigma_\theta^2 + f^2 \sigma_{tilt}^2 \right)\]

Evaluating this bound demonstrates that an uncompensated stage angular jitter of \(\sigma_\theta = 0.5\text{ mrad}\) generates a structural depth uncertainty of \(\sigma_{X_z} \approx 42.8\text{ }\mu\text{m}\), obscuring fine phenotypic sulcus variations. Implementing ZVD-shaped closed-loop LQG tracking suppresses angular jitter below \(\sigma_\theta \le 0.05\text{ mrad}\), reducing geometric reconstruction uncertainty below \(\sigma_{X_z} \le 8.4\text{ }\mu\text{m}\).

## 3D Vision and Morphological Phenotyping Pipeline

### Multi-View Constrained Bundle Adjustment (CBA)

Standard unconstrained Structure-from-Motion (SfM) pipelines parameterize external camera poses independently for each captured frame, introducing up to \(6M\) degrees of freedom [11]. For rotation-symmetric, low-texture seeds (such as polished rice or smooth pulses), unconstrained estimation results in epipolar geometry failures and rotational drift [11].

Because the physical mechatronic turntable constrains seed motion to a single circular trajectory around axis \(\mathbf{n}\), the vision pipeline implements Constrained Bundle Adjustment (CBA) [19]. Camera projection centers \(\mathbf{C}_i\) are restricted to a circular trajectory parameterized by rotation center \(\mathbf{C}_0\), stage normal \(\mathbf{n}\), and orbit radius \(\rho\) [11]:

\[\min_{\mathbf{n}, \mathbf{C}_0, \{\mathbf{X}_j\}, \{\theta_i\}} \sum_{i=1}^{M} \sum_{j=1}^{N} \rho_H \left( \left\Vert{} \mathbf{x}_{ij} - \pi\left(\mathbf{K}, \mathbf{R}_c \mathbf{R}(\theta_i, \mathbf{n}), \mathbf{t}_c, \mathbf{X}_j\right) \right\Vert{}_{\mathbf{\Sigma}_{meas}^{-1}}^2 \right)\]

where \(\rho_H(\cdot)\) is the Huber loss function that mitigates false feature correspondences. Angular coordinates \(\theta_i\) are initialized directly from optical encoder timestamps, reducing the parameter search space by over \(70\%\) and eliminating local convergence minima [19].

The algorithmic workflow transitions smoothly from multi-angle image ingestion to geometric and phenotypic quantification. The sequential stages of this data processing pipeline comprise:

- Multi-Angle Synchronized Ingestion: Capturing 36 calibrated multi-view frames at uniform \(10^\circ\) rotational increments alongside hardware-synchronized encoder angular readings and structured-light fringe patterns [1].
- Constrained Bundle Adjustment: Enforcing rigid circular orbit kinematics parameterized by rotation axis \(\mathbf{n}\) and center \(\mathbf{C}_0\) using Levenberg-Marquardt optimization within the Ceres solver [19].
- Dense Surface Reconstruction: Generating watertight two-manifold triangular meshes via Screened Poisson Surface Reconstruction and Marching Cubes, or continuous neural representations via 3D Gaussian Splatting with rotation flow priors [2].
- Phenotypic Trait Extraction: Processing the watertight manifold through principal component alignment, transverse slicing, and numerical volume integration to extract dimensional, sulcus, and allometric mass parameters [1].

### Dense Surface Reconstruction: 3D Gaussian Splatting and Poisson Meshing

Continuous volumetric radiance and explicit polygonal meshes are extracted using two complementary computational representations:

- Continuous Radiance Splatting: 3D Gaussian Splatting (3DGS) parameterizes the scene as a dense collection of anisotropic Gaussian primitives, each defined by a spatial center \(\boldsymbol{\mu}\), covariance matrix \(\mathbf{\Sigma} = \mathbf{R}_g \mathbf{S}\mathbf{S}^T\mathbf{R}_g^T\), opacity \(\alpha\), and spherical harmonics color coefficients [2]. A rotation flow field prior links primitive transformations across turntable increments: \(\boldsymbol{\mu}_i = \mathbf{R}(\theta_i, \mathbf{n})(\boldsymbol{\mu}_0 - \mathbf{C}_0) + \mathbf{C}_0\) [11].
- Watertight Mesh Extraction: To calculate volumetric parameters, raw oriented point clouds \(\mathcal{P} = \{(\mathbf{p}_k, \hat{\mathbf{n}}_k)\}\) obtained from structured light stereo matching undergo Screened Poisson Surface Reconstruction [27]. The algorithm solves for an indicator function \(\chi\) whose gradient matches the inwardly oriented normal vector field \(\mathbf{V}\):

  

  \[\nabla^2 \chi = \nabla \cdot \mathbf{V}\]

  

  Extracting the iso-surface via Marching Cubes generates a two-manifold, watertight triangulated mesh \(\mathcal{M} = (\mathcal{V}, \mathcal{F})\) [27].

### Phenotypic Trait Extraction Algorithms

Watertight mesh volume (\(V\)) and surface area (\(S\)) are calculated over all triangular facets \(f_m \in \mathcal{F}\) with vertices \((\mathbf{v}_{m1}, \mathbf{v}_{m2}, \mathbf{v}_{m3})\) using the divergence theorem [3]:

\[V = \frac{1}{6} \left\vert{} \sum_{m=1}^{\vert{}\mathcal{F}\vert{}} \mathbf{v}_{m1} \cdot (\mathbf{v}_{m2} \times \mathbf{v}_{m3}) \right\vert{}\]

\[S = \frac{1}{2} \sum_{m=1}^{\vert{}\mathcal{F}\vert{}} \Vert{}(\mathbf{v}_{m2} - \mathbf{v}_{m1}) \times (\mathbf{v}_{m3} - \mathbf{v}_{m1})\Vert{}\]

Tri-axial dimensions (Length \(L\), Width \(W\), Thickness \(T\)) are extracted by evaluating the spatial covariance of vertex coordinates \(\mathcal{V}\) through Principal Component Analysis (PCA) [1]:

\[\mathbf{C}_{pts} = \frac{1}{\vert{}\mathcal{V}\vert{}} \sum_{k=1}^{\vert{}\mathcal{V}\vert{}} (\mathbf{p}_k - \bar{\mathbf{p}})(\mathbf{p}_k - \bar{\mathbf{p}})^T = \mathbf{U} \mathbf{\Lambda} \mathbf{U}^T\]

Projecting the mesh vertices onto the orthonormal eigenvectors \(\mathbf{u}_1, \mathbf{u}_2, \mathbf{u}_3\) yields the canonical bounding dimensions:

\[L = \max_k(\mathbf{p}_k \cdot \mathbf{u}_1) - \min_k(\mathbf{p}_k \cdot \mathbf{u}_1) \\ W = \max_k(\mathbf{p}_k \cdot \mathbf{u}_2) - \min_k(\mathbf{p}_k \cdot \mathbf{u}_2) \\ T = \max_k(\mathbf{p}_k \cdot \mathbf{u}_3) - \min_k(\mathbf{p}_k \cdot \mathbf{u}_3)\]

Ventral sulcus (crease) metrics are extracted by slicing the mesh with transverse planes orthogonal to the longitudinal eigenvector \(\mathbf{u}_1\) at \(5\%\) increments along \(L\) [1]. For each planar boundary curve \(\mathcal{C}_{\text{slice}}\), a local two-dimensional convex envelope \(\mathcal{H}_{2D}\) is computed. The sulcus indentation depth \(h_{\text{sulcus}}\) is extracted as:

\[h_{\text{sulcus}}(x) = \max_{\mathbf{q} \in \mathcal{C}_{\text{slice}}} \operatorname{dist}(\mathbf{q}, \partial\mathcal{H}_{2D})\]

Thousand-Grain Weight (TGW) is estimated from 3D volumetric metrics using an allometric structural power model rather than planar projection approximations [1]:

\[m_{\text{grain}} = \rho_b \cdot V^\kappa + \epsilon_w\]

where \(\rho_b\) is the bulk apparent endosperm density (\(\approx 1.25\text{–}1.42\text{ mg/mm}^3\)) and \(\kappa\) is an allometric scaling exponent determined via cross-validation against high-precision load-cell balances [3].

| Extracted Grain Phenotype | Mathematical Definition / Operational Principle | Resolution Limit | Validated Error (MAPE / Bias) | Correlation Coefficient (\(R^2\) vs Reference) |
| --- | --- | --- | --- | --- |
| **Grain Length (\(L\))** | Caliper distance along first principal eigenvector \(\mathbf{u}_1\) [cite: 1, 27] | \(10\text{ }\mu\text{m}\) | \(1.83\%\) [cite: 1] | \(R^2 = 0.998\) [cite: 27] |
| **Grain Width (\(W\))** | Extent along orthogonal lateral eigenvector \(\mathbf{u}_2\) [cite: 1, 27] | \(10\text{ }\mu\text{m}\) | \(1.86\%\) [cite: 1] | \(R^2 = 0.992\) [cite: 27] |
| **Grain Thickness (\(T\))** | Extent along dorsoventral eigenvector \(\mathbf{u}_3\) [cite: 1, 27] | \(15\text{ }\mu\text{m}\) | \(2.19\%\) [cite: 1] | \(R^2 = 0.985\) [cite: 27] |
| **Enclosed Volume (\(V\))** | Divergence volume integral over watertight mesh \(\mathcal{M}\) [cite: 3, 27] | \(0.05\text{ mm}^3\) | \(2.45\%\) | \(R^2 = 0.976\) [cite: 28] |
| **Surface Area (\(S\))** | Summation of 2-manifold triangular face surface integrals [27] | \(0.10\text{ mm}^2\) | \(3.10\%\) | \(R^2 = 0.965\) |
| **Ventral Sulcus Depth (\(h_s\))** | Maximal Euclidean distance from 2D convex envelope to crease bottom [1] | \(8\text{ }\mu\text{m}\) | \(4.81\%\) [cite: 1] | \(R^2 = 0.942\) [cite: 1] |
| **Crease Volume Ratio (\(V_c/V\))** | Volume of difference manifold between convex hull and seed mesh [1] | \(0.02\text{ mm}^3\) | \(5.20\%\) | \(R^2 = 0.915\) |
| **Thousand-Grain Weight (TGW)** | Allometric density prediction derived from 3D solid model [1] | \(0.1\text{ mg}\) | \(2.10\%\) | \(R^2 = 0.978\) [cite: 28, 29] |

## Experimental Validation and Systems Integration Guidelines

### Calibration Protocols and Metrological Verification

System calibration follows a three-stage sequential optimization protocol:

1. Optoelectronic Intrinsic and Extrinsic Calibration: A micro-machined alumina calibration target (\(0.5\text{ mm}\) grid spacing, spatial tolerance \(\pm 0.5\text{ }\mu\text{m}\)) is imaged across 36 discrete orientations [30]. Intrinsic camera parameters (\(\mathbf{K}_c\)), radial-tangential lens distortion polynomials (\(k_1, k_2, p_1, p_2\)), and multi-camera stereo baselines are resolved using nonlinear Levenberg-Marquardt optimization within the Ceres solver library [19].
2. Axis of Rotation Determination: A certified \(3.000\text{ mm}\) tungsten-carbide sphere is mounted eccentrically on the rotary stage [32]. As the stage rotates through \(360^\circ\), sphere center trajectories are recorded by stereo triangulation [32]. The normal vector \(\mathbf{n}\) and pivot point \(\mathbf{C}_0\) are solved via algebraic circle fitting in 3D Euclidean space:

  

  \[\min_{\mathbf{n}, \mathbf{C}_0, R_{orb}} \sum_{k=1}^{M} \left( \Vert{}\mathbf{X}_{sph,k} - \mathbf{C}_0\Vert{}^2 - R_{orb}^2 \right)^2 + \lambda \vert{}(\mathbf{X}_{sph,k} - \mathbf{C}_0)^T \mathbf{n}\vert{}^2\]

  

  Stage runout parameters are recorded and incorporated into the compensation lookup table of the drive controller [26].
3. Photometric Ground-Truth Benchmarking: Validation conducted on a 125-grain set across distinct hexaploid wheat cultivars confirms that structured light scanning combined with constrained bundle adjustment achieves length, width, and thickness mean absolute percentage errors (MAPEs) of \(1.83\%\), \(1.86\%\), and \(2.19\%\), respectively, referenced against industrial micro-CT scans [1]. Sulcus depth measurements maintain a MAPE of \(4.81\%\), resolving sub-millimeter indentations without destructive cross-section microtomy [1].

### Engineering Synthesis and Implementation Rules

Integrating modern control theory into high-throughput phenotyping workstations resolves several fundamental mechatronic bottlenecks:

- Settling Delays and Dynamic Decoupling: In uncompensated indexing stages, settling delays (\(>140\text{ ms}\)) dominate the duty cycle to prevent motion-induced blur. Implementing ZVD input shaping suppresses dominant torsional resonances at \(\omega_n\), shortening settling times to \(36\text{ ms}\) and permitting strobe triggering during smooth deceleration phases [14].
- Specimen Slip Prevention: Dynamic tangential and centrifugal accelerations must not exceed the pericarp static friction envelope. Constrained Model Predictive Control bounds platter acceleration below \(\vert{}\ddot{\theta}_l\vert{} \le \mu_s g / r_{\max}\) (\(\mu_s \approx 0.35\)), preventing seed slip during rapid repositioning maneuvers without requiring mechanical clamping fixtures [7].
- Degeneracy Elimination via Kinematic Priors: Constraining bundle adjustment optimization to circular orbits parameterized by a calibrated axis of rotation reduces camera extrinsic degrees of freedom, preventing projective drift and feature tracking failure on texture-poor seed varieties [11].
- Volumetric Allometry for Mass Estimation: Utilizing watertight three-manifold surface integrals replaces planar projection proxies, improving the correlation of non-destructive seed weight estimation from \(R^2 \approx 0.82\) to \(R^2 > 0.97\) and establishing an operational methodology for large-scale crop germplasm screening [1].

## References

[1] An Intelligent Analysis Method for 3D Wheat Grain and Ventral — https://pmc.ncbi.nlm.nih.gov/articles/PMC9044079/
[2] Seed 3D Phenotyping Across Multiple Crops Using 3D Gaussian — https://www.mdpi.com/2077-0472/15/22/2329
[3] phenoSeeder - A Robot System for Automated Handling and — https://academic.oup.com/plphys/article/172/3/1358/6115817
[4] Shape analysis of grains of Indian wheat varieties - ResearchGate — https://www.researchgate.net/publication/223936034_Shape_analysis_of_grains_of_Indian_wheat_varieties
[5] (PDF) phenoSeeder - A Robot System for Automated Handling and — https://www.researchgate.net/publication/310614081_phenoSeeder_-_A_Robot_System_for_Automated_Handling_and_Phenotyping_of_Individual_Seeds
[6] An Integrated Tunable-Focus Light Field Imaging System for 3D — https://www.mdpi.com/2304-6732/13/4/385
[7] Seed-to-plant-tracking: automated phenotyping of seeds ... - Frontiers — https://www.frontiersin.org/journals/plant-science/articles/10.3389/fpls.2025.1539424/full
[8] Food safety VideometerLab for seed and grain applications - DiTECT — https://ditect.eu/wp-content/uploads/2021/04/DiTECT-Videometer-Food-safety-for-grain-and-seed-applications.pdf
[9] Multispectral imaging and automated analysis for quantifying grain — https://pmc.ncbi.nlm.nih.gov/articles/PMC12808467/
[10] Handbook of machine vision - PDF Free Download - epdf.pub — https://epdf.pub/handbook-of-machine-visiona8eee5c2c5d0eab93897162bfccc31d817212.html
[11] NeRF-Based Point Cloud Reconstruction Using a Stationary — https://www.researchgate.net/publication/395607740_SC-NeRF_NeRF-Based_Point_Cloud_Reconstruction_Using_a_Stationary_Camera_for_Agricultural_Applications
[12] Modeling of Piezoelectric-Driven Stick–Slip Actuators | Request PDF — https://www.researchgate.net/publication/224127303_Modeling_of_Piezoelectric-Driven_Stick-Slip_Actuators
[13] Multispectral imaging – a new tool in seed quality assessment? — https://www.cambridge.org/core/journals/seed-science-research/article/multispectral-imaging-a-new-tool-in-seed-quality-assessment/15A21396FE58F41C35A4EE6AE60EA616
[14] Specified-duration shapers for suppressing residual vibrations — https://journals.plos.org/plosone/article?id=10.1371/journal.pone.0276669
[15] (PDF) Input Shaping Control of a Flexible Structure for Rest-to-Rest — https://www.researchgate.net/publication/389710759_Input_Shaping_Control_of_a_Flexible_Structure_for_Rest-to-Rest_and_Non-Rest-to-Rest_Maneuvers
[16] Input Shaping for Motion Control Vibration Reduction - Zaber — https://www.zaber.com/articles/input-shaping-for-vibration-reduction
[17] A review of command shaping techniques for elimination of residual — https://www.extrica.com/article/16725
[18] Calibration method for laser-camera scanners with a motorized — https://www.spiedigitallibrary.org/journals/optical-engineering/volume-64/issue-06/063103/Calibration-method-for-laser-camera-scanners-with-a-motorized-rotating/10.1117/1.OE.64.6.063103.full
[19] Constrained Bundle Adjustment for Panoramic Cameras - CMP — https://cmp.felk.cvut.cz/ftp/articles/alblcene/Albl-Pajdla-CVWW-2013.pdf
[20] US20050166413A1 - CMM arm with exoskeleton - Google Patents — https://patents.google.com/patent/US20050166413A1/en
[21] Advances in Optical and Mechanical Technologies for Telescopes — https://spie.org/AS26/conferencedetails/technologies-for-telescopes-and-instrumentation
[22] Introduction to autonomous mobile robots [2nd ed] 9780262015356 — https://ebin.pub/introduction-to-autonomous-mobile-robots-2nd-ed-9780262015356-0262015358.html
[23] Contributions to image-based object reconstruction: geometric and — https://bmva-archive.org.uk/theses/2006/2006-guillemaut.pdf
[24] 3D model reconstruction by constrained bundle adjustment — http://ieeexplore.ieee.org/document/1334674/
[25] A Fast Recursive 3D Model Reconstruction Algorithm for Multimedia — http://www.cse.cuhk.edu.hk/~khwong/demo/cba/cba.html
[26] Fast and Accurate Reconstruction of Pan-Tilt RGB-D Scans via Axis — https://arxiv.org/html/1812.00240v3
[27] Automatic Measurement of Seed Geometric Parameters Using a — https://pdfs.semanticscholar.org/274e/5db028a955cf16deaa477909dfb75a8610d8.pdf
[28] Relationship between measured and estimated ear a length and b — https://www.researchgate.net/figure/Relationship-between-measured-and-estimated-ear-a-length-and-b-width-CCCconcordance_fig8_325780990
[29] Wheat individual grain-size variance originates from crop — https://journals.plos.org/plosone/article/file?id=10.1371/journal.pone.0230689&type=printable
[30] Deep Learning Based 3d Reconstruction for Phenotyping of Wheat — https://openaccess.thecvf.com/content/ICCV2023W/CVPPA/papers/Cherepashkin_Deep_Learning_Based_3d_Reconstruction_for_Phenotyping_of_Wheat_Seeds_ICCVW_2023_paper.pdf
[31] High-Precision Rotation Axis Calibration of Line-Structured Light — https://www.mdpi.com/1424-8220/26/13/4275
[32] Multi-view point cloud registration and 3D reconstruction using the — https://www.researchgate.net/publication/404937667_Multi-view_point_cloud_registration_and_3D_reconstruction_using_the_rotation_trajectory_ofthe_sphere's_center
[33] Universal calibration for a ring camera array based on a rotational — https://opg.optica.org/abstract.cfm?uri=oe-30-9-14538


Please begin the extraction now. Output only the JSON list directly, without any chitchat or explanations.