SPECTROSCOPY METHOD AND SYSTEM

Abstract

We describe a method of making a prediction of when a tuber will sprout. In some embodiments the method comprises making a measurement of an optical reflectance of the tuber in an eye region of the tuber; making a reference measurement of optical reflectance of the or another tuber, processing the tuber eye measurement in combination with the reference measurement; and making a prediction of when the tuber will sprout in response.

Claims

1. A method of making a prediction of when a tuber will sprout, the method comprising: making a tuber eye measurement of an optical reflectance of said tuber in an eye region of the tuber; making a reference measurement of optical reflectance of the said or another tuber; processing said tuber eye measurement in combination with said reference measurement; and making a prediction of when said tuber will sprout responsive to a result of said processing.

2. A method as claimed in claim 1 wherein said reference measurement comprises a measurement of a non-eye region of the said or another tuber and wherein said processing includes determining a ratio of said tuber eye measurement and said reference measurement.

3. A method as claimed in claim 1 wherein said tuber eye measurement comprises a measurement at a wavelength in the range 600 nm to 750 nm, in particular in the range 650 nm to 700 nm, more particularly at a wavelength of approximately 690 nm.

4. A method as claimed in claim 2 using a wavelength-sensitive optical sensing system to measure features of an optical reflectance spectrum of said eye region of said tuber to make said measurements of said eye region of the same said tuber at different respective wavelengths.

5. A method as claimed in claim 1 comprising making a further optical reflectance measurement on the said or another tuber and compensating for a reference background signal level in said optical reflectance measurements when processing said tuber eye measurement.

6. A method as claimed in claim 1 wherein said processing comprises providing said tuber eye measurement in combination with said reference measurement to a mathematical model providing a prediction of when said tuber will sprout responsive to said tuber eye measurement and said reference measurement, wherein said mathematical model has been trained using a training data set defining measured time to sprouting for a set of said tuber eye measurements and reference measurements.

7. A method as claimed in claim 6 further comprising providing a plurality of said mathematical models, for a plurality of different tuber cultivars, and selecting a said model for making said prediction responsive to a cultivar of said tuber.

8. A method as claimed in claim 6 wherein said model comprises a partial least squares model.

9. A method of identifying a potato or batch of potatoes using the method of claim 1 to make a prediction of when the potato or one or more potatoes of the batch will sprout, and making a decision on sale or use of said potato or said batch responsive to said prediction.

10. A method of processing potatoes using the method of claim 1, the method comprising: making a tuber eye measurement of an optical reflectance of a potato tuber in an eye region of the tuber; making a reference measurement of optical reflectance of the said or another tuber; making a prediction of when said tuber will sprout responsive to said tuber eye measurement in combination with said reference measurement; and processing said potatoes responsive to said prediction, wherein said processing includes at least separating potatoes predicted to sprout sooner from potatoes predicted to sprout relatively later.

11. A method as claimed in claim 10 further comprising one or both of: i) storing potatoes selected by said processing in a selected storage area; and ii) dispatching potatoes selected by said processing.

12. A method as claimed in claim 1, for making a prediction of when a tuber will sprout, the method comprising: using an optical reflectance measurement system for making a tuber eye measurement of an optical reflectance of an eye region of the tuber, and for making a reference measurement of optical reflectance of the said or another tuber; and using a data processor, coupled to said optical reflectance measurement system, to process said tuber eye measurement in combination with said reference measurement and to make a prediction of when said tuber will sprout responsive to a result of said processing.

13. A method as claimed in claim 1, wherein making said prediction of when said tuber will sprout comprises identifying a change in gradient of said optical reflectance over time, and/or identifying a reduction in a value of said optical reflectance to below a threshold level or by a threshold amount, and making a prediction of a point of intervention to inhibit said sprouting.

14. A method for predicting sprouting of one or more root vegetables having at least one bud, using the method of claim 1, the method comprising: measuring a reflectance of at least one of said one or more root vegetables to obtain a reflectance spectrum; and applying a multivariate model to said reflectance spectrum to predict said sprouting of said one or more root vegetables.

15. A method as claimed in claim 14, wherein said measuring comprises performing a plurality of reflectance measurements at two or more locations on each of said at least one of said one or more root vegetables.

16. A method as claimed in claim 15, wherein said two or more locations are at least one apical bud and at least one location away from an apical bud on said at least one of said one or more root vegetables.

17. A method as claimed in claim 16, wherein said measurement of said reflectance at said at least one location away from an apical bud is used to obtain a reference background reflectance spectrum.

18. A method as claimed in claim 14, wherein said multivariate model is a linear mathematical model.

19. A method as claimed in claim 14, wherein said multivariate model is a Partial Least Squares model.

20. A method as claimed in claim 19, wherein said applying said multivariate model comprises predicting said sprouting using a calculated regression coefficient for a previous said measured reflectance.

21. A method as claimed in claim 14, wherein said measuring comprises measuring said reflectance in a non-contact, non-destructive mode.

22. A method as claimed in claim 14, wherein said reflectance spectrum is obtained for optical wavelengths between 500 nm and 1200 nm.

23. A method as claimed in claim 14, wherein said root vegetables are potato tubers.

24. A method as claimed in claim 14 for predicting sprouting of one or more root vegetables having at least one bud, the method comprising: using a spectrometer to measure a reflectance of at least one of said one or more root vegetables to obtain a reflectance spectrum; and using a processor configured to apply a multivariate model to said reflectance spectrum to predict said sprouting of said one or more root vegetables.

25. A method for predicting sprouting of one or more root vegetables, as claimed in claim 14, the method comprising: measuring a reflectance of at least one of said one or more root vegetables at a plurality of points in time to obtain a plurality of reflectance spectra; fitting a curve to each of said reflectance spectra; and comparing said fitted curves to predict said sprouting of said one or more root vegetables.

26. A method as claimed in claim 25, wherein said fitted curves are polynomial curves.

27. A method as claimed in claim 25 further comprising subtracting each of said curves from an initial reflectance spectrum of said plurality of reflectance spectra prior to said comparison.

28. A method as claimed in claim 25 further comprising subtracting each of said measured reflectance spectra from an initial reflectance spectrum of said plurality of reflectance spectra to obtain a plurality of normalised reflectance spectra, and wherein said fitting comprises fitting a curve to each of said normalised reflectance spectra.

29. A method as claimed in claim 25, wherein said measuring comprises measuring a reflectance at an apical bud of a said root vegetable.

30. A method as claimed in claim 1 for monitoring tuber sprouting to predict a point of intervention to inhibit the sprouting, the method comprising: making a time series of optical measurements on said tuber at at least one wavelength to provide time series spatial data; processing said time series optical data to monitor the evolution over time of a spectral feature in said optical measurements; and predicting a point of intervention to inhibit the sprouting, from said evolution over time.

31. A method as claimed in claim 30 comprising making said optical measurements at a plurality of wavelengths or over a wavelength range; wherein said spectral feature is an integrated spectral feature comprising a sum or integration over said plurality of wavelengths or wavelength range; and wherein said processing comprises monitoring the evolution of said integrated spectral feature over time.

32. A method as claimed in claim 30 wherein predicting comprises identifying a change in gradient of said evolution over time.

33. A method as claimed in claim 32 wherein said identifying of said change in gradient comprise identifying a reduction in gradient.

34. A method as claimed in claim 32 wherein said identifying of said change in gradient comprises identifying when $\frac{I}{t}>C$ where C is a threshold value less than 0, l represents a value of said spectral feature, and t represents time.

35. A method as claimed in claim 32 wherein predicting comprises identifying a reduction in a value of said monitored spectral feature to below a threshold level or by a threshold amount.

36. A method as claimed in claim 30 wherein said tuber is a potato, wherein said at least one wavelength is a wavelength in the range of 600 nm to 750 nm, and wherein said optical measurements comprise measurements of an eye of said potato.

37. A method as claimed in claim 30 further comprising inhibiting said sprouting by applying a sprouting suppressant to said tuber in response to said predicting.

38. Apparatus for determining the condition of a tuber, the system comprising: a light source to stimulate the tuber to promote chlorophyll production and/or sprouting, wherein said light source is configured to provide substantially continuous optical stimulation to the tuber; an optical instrument to measure an optical response of the tuber; a controller to control said optical instrument to make a time series of optical measurements on said tuber at intervals during said substantially continuous optical stimulation to determine a time evolution of an optical response of the tuber, wherein said time series evolution is determined over a period of less than twenty four hours; a data processor to analyse said time evolution of said optical response to determine a condition of the tuber.

39. Apparatus as claimed in claim 38 wherein said condition of said tuber comprises a prediction of when said tuber will sprout or defines a sprouting propensity of said tuber.

40. Apparatus as claimed in claim 39 wherein said data processor is configured to analyse said optical response to determine a condition of the tuber by determining one or more of i) a time interval until a threshold change in said optical response; ii) a rate of change of said optical response; and iii) a curve fit to said optical response.

41. Apparatus as claimed in claim 39 wherein said optical response comprises an optical reflectance or absorption response in the range 600 nm to 750 nm and/or an integrated optical response over a wavelength band.

42. Apparatus as claimed in claim 39 wherein said optical response comprises an optical reflectance or absorption response within 50 nm of an absorption band of chlorophyll.

43. Apparatus as claimed in claim 38 wherein said time series evolution is determined over a period of less than one hour.

44. Apparatus as claimed in claim 38 wherein said controller is configured to control said optical instrument to make measurements at least every 5 minutes.

45. A method of determining the sprouting propensity of a tuber, the method comprising: applying substantially continuous stimulation to an eye region of the tuber to drive said eye region of the tuber to develop at a faster than natural rate; making a time series of optical measurements on said eye region of said tuber at at least one wavelength to provide time series optical data to determine a time evolution of an optical response of said eye region of the tuber during said stimulation, wherein said time series evolution is determined over a period of less than twenty four hours; and determining a sprouting propensity of the tuber from said time evolution of said optical response.

46. A method as claimed in claim 45 wherein said substantially continuous stimulation comprises light.

47. A method as claimed in claim 46 wherein said development of said eye region of the tuber comprises chlorophyll production in said eye region of the tuber; wherein said sprouting propensity comprises a prediction of when said tuber will sprout, and wherein said determining of said sprouting propensity comprises determining a speed of response of said eye region of said tuber to said stimulation.

48. A method as claimed in claim 45 wherein said determining of said sprouting propensity comprises determining one or more of i) a time interval until a threshold change in said optical response; ii) a rate of change of said optical response; and iii) a curve fit to said optical response.

49. A method as claimed in claim 45 wherein said time series evolution is determined over a period of less than one hour.

Description

BRIEF DESCRIPTION OF THE DRAWINGS

[0065] These and other aspects of the invention will now be further described, by way of example only, with reference to the accompanying figures in which:

[0066] FIG. 1 shows an example experimental setup of an embodiment of the invention;

[0067] FIG. 2 shows a series of reflectance spectra of an apical bud;

[0068] FIG. 3 shows a set of visible/near-infrared data;

[0069] FIGS. 4a and 4b show cross-validation predictions and regression coefficients;

[0070] FIG. 5 shows further shows cross-validated predictions;

[0071] FIG. 6 shows regression coefficients for various types of potatoes;

[0072] FIG. 7 shows leave-tuber-out cross-validations;

[0073] FIG. 8 shows a boxplot summary of predictions;

[0074] FIG. 9 shows an example experimental setup of an embodiment of the invention;

[0075] FIGS. 10a to 10d show predicted sprouting age versus actual sprouting age for different types of tubers;

[0076] FIG. 11 shows regression coefficients related to PLS modelling for the Mozart tubers analysed as shown in FIG. 10a;

[0077] FIGS. 12a and 12b show visible IR spectra and time analysis of the spectra, respectively;

[0078] FIG. 13 shows calculated area versus sprouting age for Maris Piper tubers;

[0079] FIG. 14 shows calculated area versus sprouting age for different tubers;

[0080] FIGS. 15a and 15b show predicted sprouting age versus actual sprouting age for different Maris Piper tubers;

[0081] FIG. 16 shows Vis/NIR spectra of a Mozart tuber; and

[0082] FIGS. 17a and 17b show integrated feature intensity versus time for different tubers.

DETAILED DESCRIPTION OF PREFERRED EMBODIMENTS

[0083] Sprouting is a major cause of losses, and therefore commercial losses, in stored potatoes. There is therefore a need for predicting the onset of sprouting in potato tubers.

[0084] In order to predict the onset of sprouting in potato tubers, visible/near-infrared spectra were collected using a portable, low-cost StellarNet spectrometer, equipped with a fibre-optic probe operating in the wavelength range of 500 nm to 1200 nm. The spectra were collected through a non-contact, non-destructive reflectance measurement, taken from the potato skin at two different locations: the apical buds (eyes) and a portion of the smooth skin located away from the eyes (background). In some experiments, pseudo-absorbance spectra were calculated for each eye, by ratioing to the previously collected background spectrum.

[0085] It was discovered that potatoes tend to age in a consistent manner which may be tracked by the visible/near-infrared spectra. A model based on a standard chemometric analysis method allowed predicting the length of time until sprouting of individual potatoes and potato batches.

[0086] Embodiments of the present invention were used to predict sprouting in various types of potatoes under different ageing conditions. In particular, studies related to monitoring of a Maris Piper potato forced into sprouting by exposure to light over a three-day period, longer-term monitoring of Maris Piper potatoes stored in the dark between analyses, and longer-term monitoring of Desiree and King Edward potatoes stored in the dark between analyses.

Maris Piper Potato Exposed to Light Over a Three-Day Period

[0087] In a first attempt at detecting changes occurring before sprouting, an eye on a single Maris Piper tuber was monitored over a period of three days. Spectra were acquired hourly. This study took place at ambient laboratory temperature (21 C.) and the tuber was illuminated throughout, promoting quick sprouting. FIG. 1 shows an example experimental setup used in this investigation.

[0088] FIG. 2 shows a series of spectra obtained from an apical bud, collected over three days.

[0089] These data clearly show a profound difference between the initial (before sprouting) and final (after sprouting) spectra, and a gradual progression between the two states. This allowed predicting the onset of potato sprouting by pattern recognition methods applied to the pre-sprouted spectra.

Longer-Term Monitoring of Maris Piper Potatoes Stored in the Dark Between Analyses

[0090] In a second series of experiments, 4 tubers were monitored every few days over a period of weeks. The tubers were from the 2012 harvest, and had been stored and treated with sprout suppressant as per standard commercial practice. Each tuber had a number of readily identifiable apical buds.

[0091] The tubers were stored at 4 C. in the dark between analyses. On each tuber, different sites were identified and labelled as follows: site 0 corresponded to a background, and sites 1 to 3 were non-sprouted apical buds.

[0092] Replicate analyses (at least two-fold) were made of each site on each day of analysis, repositioning the fibre probe for each acquisition. This may be particularly useful as there may be a variance in the data arising from repositioning of the collection optics when using a handheld device.

[0093] FIG. 3 shows a set of visible/near-infrared data from two eye sites on two different Maris Piper tubers which ultimately sprouted at around 4 weeks (here called analysis day 12). The data shown in these graphs relates to the ratio of tuber eye measurements and reference measurements taken on the respective tubers. In other words, background ambient light is taken into consideration. Note that in all cases the x-axes are in local data points, not wavelength units.

[0094] Left hand plots show the complete spectral range, ratioed to the site 0 backgrounds. Right hand plots show an expansion of the marked feature, which rapidly becomes prominent as the tubers age, in advance of the onset of sprouting.

[0095] For the majority of the sites monitored, a clear and systematic trend was seen in the spectra as the tubers aged, in advance of the visible signs of sprouting. As well as an overall spectral shape change, which represents an increasing difference in the visible/near-infrared reflectance of the eyes and the reference background skin, changes also occurred which led to the emergence of individual spectral bands in the spectra, specifically an increasingly negative feature at around 690 nm (data points 300-400 in the left hand figures). Background locations did not change quite so much, which is in agreement with literature studies which have shown that the dormancy-breaking chemical processes occur predominantly at the eyes.

[0096] The findings were consistent across different tubers, and persisted despite the large sampling variability that arose from the positioning of the fibre-optic probe relative to the tuber surface (replicate measurements were made to mitigate this source of variance, since a single measurement may be performed within a few seconds).

[0097] By using a non-eye location on the same tuber to acquire a reference background, other confounding factors such as ambient lighting may be accounted for.

[0098] In order to develope a predictive model for time until sprouting, a multivariate modelling may be applied to the visible/near-infrared spectra, in this example a Partial Least Squares modelling.

[0099] In this example, a continuum regression model is used to allow for making a prediction of when a tuber will sprout. The continuum regression model is of the form y=Kx, where y is a prediction of time until a tuber will sprout, K is a scalar or vector constant, and x represents a scalar or vector wavelength ratio or spectrum ratio.

[0100] The skilled person will appreciate that the continuum regression model used here is just one example of a model which may be exploited to make a prediction of when a tuber will sprout. The skilled person will immediately understand that any alternative model may be used to make this prediction.

[0101] FIG. 4a shows cross-validated predictions of tuber age (data from tubers 1 and 2 predicting age of tubers 5 and 8).

[0102] FIG. 4b shows regression coefficients as a function of wavelength for the data in FIG. 4a.

[0103] The modelling work showed that it was possible to apply Partial Least Squares modelling to the visible/near-infrared spectra from one collection of tubers, and use this model to predict the age of other tubers, as in this example Maris Piper potato tubers, from their spectra.

Longer-Term Monitoring of Desiree and King Edward Potatoes Stored in the Dark Between Analyses

[0104] This study examined two further commercially important cultivars: King Edward and Desiree.

[0105] Three batches of potatoes were examined. The first batch was delivered direct from the commercial store soon after harvest. Six tubers were selected for study, three of each cultivar. The tubers were monitored using the same experimental protocol as described above in relation to the Maris Piper potatoes stored in the dark between analyses. In this instance, however, all potatoes were analysed until sprouting, and this was found to occur at somewhat different times for each tuber.

[0106] The visible/near-infrared spectra were analysed using Partial Least Squares regression, in this example using calculated time until sprouting as the dependent variate.

[0107] FIG. 5 shows cross-validated predictions obtained by treating each potato separately. It is to be noted that potato Kind Edward 2 failed to sprout before becoming spoiled due to storage.

[0108] As can be seen, the error in prediction for all potatoes was approximately +/seven days, which represents a commercial useful outcome in practice.

[0109] The regression coefficients were examined for each of the Desiree and King Edward potatoes of FIG. 5, and are shown in FIG. 6.

[0110] The obtained regression coefficients were very similar, in particular across the Desiree tubers. Furthermore, they were also highly similar to the equivalent regression vector obtained from Maris Piper tubers discussed previously. Specifically, there was a large feature at approximately 675 nm that was negatively associated with increasing age. A predominant feature occurred at 640 nm for the King Edward tubers, wherein the feature was shifted slightly with respect to the one observed for Maris Piper tubers. The skilled person will appreciate that the precise location of the feature is cultivar dependent.

[0111] It is important to note that the data from each tuber were treated entirely separately here. The fact that the coefficients were so similar was strong evidence for a common effect across all potatoes.

[0112] In a further modelling work, leave-tuber-out cross-validation was employed (as for Maris Piper potatoes stored in the dark between analysis), as this more closely mimics the required real-world situation, in which different tubers are used in the model from those that will be subsequently tested. The results are shown in FIG. 7.

[0113] The most effective calibrations were obtained from the Desiree tubers, which were able to predict the left-out tuber's time-until-sprouting with a precision of approximately +/seven days.

[0114] Two further batches of potatoes were obtained from the same commercial store, at intervals of two months. Unfortunately, the Batch 2 tubers sprouted within a day of arrival at the laboratory, so no proper measurements could be taken in advance of sprouting. Batch 3 tubers also sprouted within a few days of arrival at the laboratory, but one set of measurements were able to be taken in that time.

[0115] Where measurements were performed after sprouting, the data was less reliable due to the physical changes in the apical bud area and the consequent reflectance changes. However, this data was used as independent test data with which to challenge the models described for the first batch (Batch 1 above).

[0116] FIG. 8 shows a boxplot summary of the predictions for Batch 3 Desiree tubers using a model developed from Batch 1 data.

[0117] It can be seen that there was a general trend implying that Batch 1 data could make credible predictions about Batch 3, as seen for the 95% confidence interval around the median prediction for time point 4 days. At 4 days before sprouting, the data predicts 4 days to sprout.

[0118] As described above, Batch 2 and Batch 3 tubers sprouted soon after arrival at the laboratory. Hence, the Batch 2 and Batch 3 data in FIG. 8 can be disregarded, since no proper measurements could be taken in advance of sprouting. Hence, it was clear that more tubers exhibiting a longer dormancy period would be needed to carry out proper model validation.

[0119] It will be appreciated that not all potatoes always behave in a predictable manner. This is because sampling handling and storage is difficult: any kind of shock (light, temperature change, humidity change, etc.) seems to be able to cause the potatoes to break dormancy suddenly and be driven to sprouting quickly. It may therefore be preferable to analyse potatoes in situ in the potato storage facility.

[0120] FIG. 9 shows an example experimental setup 900 used in embodiments of this invention. A system 902, which may comprise a light source, a spectrometer, optical fibres and a reflectance probe, is used to measure one or more reflectance spectra. One or more data sets 904 are obtained using the system 902 and may then be further processed and stored in computer 908 comprising a processor and memory.

[0121] It will be understood that controlling of system 902 and processing and/or storing data sets 904 may be performed using the same single computer 908, or different computers.

[0122] A model 906, as outlined above, which may be stored in memory of computer 908 or a different memory, may then be applied to the one or more data sets 904. The skilled person will appreciate that the model may be applied to the data sets 904 first before the obtained data is input in computer 908 for further processing. Alternatively or additionally, the data sets 904 may first be input and stored in computer 908, and model 906 may be applied to the data sets 904 stored in computer 908 for further processing. In the latter case, it will be appreciated that model 906 may be in direct communication with computer 908, or alternatively, model 906 is stored in memory of computer 908.

[0123] It will be understood that one model 906 may be applied per cultivar, i.e. model 906 may be different for, e.g. Kind Edwards tubers, Maris Piper tubers or Desiree tubers.

[0124] As a result of the processing in the processor of computer 908, a training algorithm 910 may further be developed. This is an optional feature in use merely for training purposes, and the skilled person will appreciate that this functionality is not essential for the method and system for making a prediction of when a tuber will sprout. Hence, training algorithm 910 may be stored separately from the system comprising computer 908. A separate storing of the training algorithm 910 allows for a remote updating of model 906 in a network.

[0125] It will be understood that alternatively, training algorithm 910 may be stored in memory of computer 908 where it is part of the analysis system.

[0126] Training algorithm 910 may then be used to train model 906 in order to improve the system and method for making a prediction of when a tuber will sprout.

[0127] It will be appreciated that training algorithm 910 may further comprise information, such as, but not limited to, measurement parameters, for example a preferred wavelength of light used which may depend on the cultivar investigated, duration of one or more measurements, or number or intervals of measurements to be taken. The training data may be processed in computer 908 and applied to future measurements.

[0128] Data sets 904 and/or training algorithm 910 and/or model 906 may be stored in memory for measurements being performed at a later stage. It will be understood that the memory of computer 908 may be used, or one or more memories outside computer 908 may be exploited to store data sets 904 and/or training algorithm 910 and/or model 906.

Harvest 2014 Monitoring

[0129] Further experiments were conducted on additional tubers obtained from the 2014 harvest. In order to examine the longer NIR wavelengths, two StellarNet spectrometers (with separate Vis and NIR channels) were exploited to increase the wavelength range to 500-2300 nm.

[0130] Four different samples of tubers were collected on Jun. 10, 2014 from various locations: King Edward and Mozart tubers were obtained from a first location, and Maris Piper tubers from second and third locations.

[0131] The tubers were prepared and analysed using the methods as outlined above. Spectra were collected from four locations on each tuber, a background location that was an area of skin free from any apical buds and three separate eye locations of three apical buds. For the present harvest, six tubers were studied from each sample to provide a larger data set. Intending to mimic commercial stores, these tubers were kept in a cold, ventilated unit to ensure a stable temperature (5.7 C.) and relative humidity (80%).

[0132] The selected potato tubers were monitored over an 18 week period. The spectral data was firstly analysed using Partial Least Squares regression (PLS). This method was used to demonstrate that Vis/NIR spectra collected from a sample of potato tubers may be used to predict the sprouting age of another.

[0133] FIG. 10a shows actual sprouting ages (calculated by defining the day on which sprout growth was first visible as day zero) against the predicted sprouting ages calculated by the optimal PLS model for six Mozart tubers. Corresponding graphs are shown in FIG. 10b for six King Edwards tubers, in FIG. 10c for six Maris Piper tubers from a first location, and in FIG. 10d for six Maris Piper tubers from a second location.

[0134] It can be seen that very similar features are observed in the figures from analysing King Edwards and Maris Piper tubers in the same way as Mozart tubers. Correlation coefficients between the actual and predicted sprouting ages were found to be as high as 0.92.

[0135] FIG. 11 shows regression coefficients related to the PLS modelling created for the six Mozart tubers analysed as shown in FIG. 10a. In each case, the pair of plots represents the regression coefficients from the two spectrometers (i.e. the separate Vis and NIR channels).

[0136] These regression coefficients reiterate the importance of the change in the spectra that occurs between 600-750 nm, evident as the greatest feature in all plots of FIG. 11. The similarity between the six plots of the six separate Mozart tubers stresses that these individual tubers all behaved in a similar way to one another, illustrating that the change in this feature during ageing is common for all tubers.

[0137] Due to the extended wavelength range of 500-2300 nm, the experiments specified above, which were conducted over a period of three days, were repeated using a tuber from the 2014 harvest.

[0138] The result is shown in FIG. 12a, which only shows a part of the measured wavelength range (500-2300 nm), in this example from about 620 nm to 720 nm. The curves show the initial spectra (red) taken from an unsprouted tuber at time zero through to the final spectra (cyan) taken from the same tuber, which sprouted a few days later.

[0139] As can be seen in FIG. 12a, the largest change in intensity was observed in the 600-700 nm region.

[0140] In order to evaluate the change from the initial to the final spectra, a method was applied to calculate the size of the increasingly negative feature between 660 nm and 680 nm. Polynomial curves (from 600 nm to 750 nm) were hereby fitted to each spectrum and then subtracted from the original time zero spectrum.

[0141] The result is shown in FIG. 12b, in which the area of these spectra section is calculated and plotted against time. A sudden change occurring in the tuber after approximately 400 minutes was observed.

[0142] In order to assess whether the sudden change at around 400 minutes as shown in FIG. 12b may be linked to dormancy breaking, a similar method of data analysis was also performed on the long-term monitoring data of the King Edward, Mozart and Maris Piper tubers.

[0143] FIG. 13 shows the calculated area difference between each spectra and a related fitted polynomial curve from 600 nm to 750 nm for Maris Piper tubers. In contrast to the single-beam spectra outlined above, the polynomial fitting has here been applied to the eye:background ratio spectra.

[0144] The 690 nm feature only begins to appear at the point at which the tubers had been recorded to start showing (barely) visible signs of sprouting. From this point onwards, the summed area steadily decreases over time.

[0145] The same analysis as the one shown in FIG. 13 has been applied to various tubers. FIG. 14 shows the calculated area difference versus sprouting age for King Edwards tubers, Mozart tubers, and Maris Piper tubers from two different sources. For this plot, the ratio (eye/background spectra), for each set of tuber samples was used to calculate the area difference between a each spectra and fitted to a polynomial curve from 600 nm to 750 nm.

[0146] The analysis displayed in FIG. 14 shows that the results were similar for the Mozart tubers and each of the two Maris Piper tubers. This method analysis may therefore be used as a simple alternative to the PLS regression approaches to predict when tubers are about to break dormancy.

[0147] FIGS. 15a and 15b shows predicted sprouting age versus actual sprouting age for different Maris Piper tubers.

[0148] For the plot shown in FIG. 15a, PLS predictions were calculated using B & C Farming Maris Piper tubers to create the training set that was used to predict the G's Fresh Maris Piper tubers sprouting ages. For the plot shown in FIG. 15b, G's Fresh tubers were used to obtain a training set to predict sprouting ages of the B & C Farming tubers.

[0149] FIG. 16 shows a series of Vis/NIR spectra collected from a single eye and tuber, in this example a Mozart tuber. For clarity, the wavelength labels on the x-axis have been omitted. It can be seen that the 690 nm feature changes as a function of date, and generally becomes more pronounced over time.

[0150] The above examples show that for separate harvest years and different tuber cultivars, dormancy breaking can be predicted using Vis/NIR spectroscopy. The observed change in the spectral data which occurs between 600 nm and 750 nm (with the central position of the spectral feature concerned varying slightly for different cultivars) allows for predicting sprouting in varies types of tubers.

[0151] The PLS modelling applied to the 2014 harvest data set has been shown to be effective for all types of cultivars monitored. As shown, in this example, for Maris Piper tubers from different locations (which may further encompass different farming practices and/or a different local climate), the method allows for predicting one another's sprouting ages, as shown in FIG. 16, reinforcing that the dominant feature between 600 nm and 750 nm is independent of growing conditions.

[0152] A further method of analysis may be implemented to investigate more specifically the spectral data between 600 nm and 750 nm, where the change in a baseline-corrected summed area for this section of the spectrum is plotted against time.

[0153] The results in this example are very similar for the Mozart and Maris Piper tubers. Once the summed area started to fall below a value of around 0.25, the first signs of tuber sprouting were observed.

[0154] FIG. 17a shows integrated feature intensity of the 690 nm feature versus time (days from the start of the study) for different batches/cultivars of potatoes, in this example from the 2014 harvest. The intensity is, in this example, integrated over a wavelength range, for example from 600 nm to 750 nm. The optical measurements are taken on a tuber eye, and in this example additionally at a non-eye region of the tuber which allows for background correction of the tuber eye measurements.

[0155] FIG. 17b shows integrated feature intensity of the 690 nm feature versus time for tubers which were forced to sprout.

[0156] As can be seen, the rate at which the change in integrated feature intensity occurs over time, as well as the starting point of the curves, may vary significantly between types of tubers. The gradient is largest, in this example, for King Edwards tubers, and a much larger gradient may be observed for tubers which are forced to sprout compared to naturally aging tubers (note the different x-axis scale of FIGS. 17a and 17b, respectively).

[0157] The King Edwards tubers exhibit a lower integrated feature intensity at the start of the measurements (day zero). This shows that the King Edwards tubers already started sprouting prior to day zero of the measurements.

[0158] As outlined above, monitoring the feature intensity (or integrated feature intensity where integration is employed) allows for determining a point in time at which intervention of sprouting may be desired. This point in time may be determined by the change in gradient of the (integrated) feature intensity over time, particularly if the gradient is above a threshold. Additionally or alternatively, the point in time at which intervention may be desired may be determined by the (integrated) feature intensity dropping by a threshold value and/or being below a (integrated) feature intensity level.

[0159] Monitoring the progression of a batch in this manner may be particularly useful for identifying a point along a time course at which the sprouting may be intervened, e.g. by spraying a suppressant known to those skilled in the art onto the tubers, and/or using temperature control methods, and/or atmospheric control methods, and/or in-field treatments as outlined above.

[0160] No doubt many other effective alternatives will occur to the skilled person. It will be understood that the invention is not limited to the described embodiments and encompasses modifications apparent to those skilled in the art and lying within the spirit and scope of the claims appended hereto.