WO2020115316A1 - Procédé d'identification d'un paramètre interne d'un œuf - Google Patents

Procédé d'identification d'un paramètre interne d'un œuf Download PDF

Info

Publication number
WO2020115316A1
WO2020115316A1 PCT/EP2019/084077 EP2019084077W WO2020115316A1 WO 2020115316 A1 WO2020115316 A1 WO 2020115316A1 EP 2019084077 W EP2019084077 W EP 2019084077W WO 2020115316 A1 WO2020115316 A1 WO 2020115316A1
Authority
WO
WIPO (PCT)
Prior art keywords
egg
excitation
thermograms
thermogram
parameter
Prior art date
Application number
PCT/EP2019/084077
Other languages
English (en)
Inventor
Bart De Ketelaere
Karlien D'huys
Niels WOUTERS
Original Assignee
Katholieke Universiteit Leuven
Moba B.V.
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by Katholieke Universiteit Leuven, Moba B.V. filed Critical Katholieke Universiteit Leuven
Publication of WO2020115316A1 publication Critical patent/WO2020115316A1/fr

Links

Classifications

    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
    • G01N33/00—Investigating or analysing materials by specific methods not covered by groups G01N1/00 - G01N31/00
    • G01N33/02—Food
    • G01N33/08—Eggs, e.g. by candling
    • A—HUMAN NECESSITIES
    • A01—AGRICULTURE; FORESTRY; ANIMAL HUSBANDRY; HUNTING; TRAPPING; FISHING
    • A01K—ANIMAL HUSBANDRY; AVICULTURE; APICULTURE; PISCICULTURE; FISHING; REARING OR BREEDING ANIMALS, NOT OTHERWISE PROVIDED FOR; NEW BREEDS OF ANIMALS
    • A01K43/00—Testing, sorting or cleaning eggs ; Conveying devices ; Pick-up devices
    • A01K43/04—Grading eggs
    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
    • G01N25/00—Investigating or analyzing materials by the use of thermal means
    • G01N25/72—Investigating presence of flaws
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T7/00—Image analysis
    • G06T7/0002—Inspection of images, e.g. flaw detection
    • G06T7/0004—Industrial image inspection
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T2207/00—Indexing scheme for image analysis or image enhancement
    • G06T2207/10—Image acquisition modality
    • G06T2207/10016—Video; Image sequence
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T2207/00—Indexing scheme for image analysis or image enhancement
    • G06T2207/10—Image acquisition modality
    • G06T2207/10048—Infrared image
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T2207/00—Indexing scheme for image analysis or image enhancement
    • G06T2207/20—Special algorithmic details
    • G06T2207/20048—Transform domain processing
    • G06T2207/20056—Discrete and fast Fourier transform, [DFT, FFT]
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T2207/00—Indexing scheme for image analysis or image enhancement
    • G06T2207/30—Subject of image; Context of image processing
    • G06T2207/30108—Industrial image inspection
    • G06T2207/30128—Food products

Definitions

  • the present invention is related to a method for identifying an internal parameter of an egg.
  • JP5967766 describes an egg inspection apparatus comprising a support for supporting an egg; a temperature measurement unit for contactlessly measuring the temperature in a predetermined region of the egg supported by the support while moving relative to said egg, said temperature measurement unit comprising at least one temperature sensor for outputting temperature data relating to the temperature; and an inspection unit for inspecting the state of the egg in accordance with a change occurring, as a result of the abovementioned relative movement, in the temperature data output from the temperature measurement unit.
  • JP5967766 a single temperature data point per location on an egg is used in the analysis.
  • a single thermal image may not provide the required level of contrast between egg shell adjacent to the air cell and egg shell adjacent to the albumen to allow reliable separation and determination of the size of the air cell.
  • WO2017/137837 discloses an apparatus to identify upside-down eggs of a batch of eggs based on the presence of a heated zone in an air cell of each egg, for use in the poultry industry and especially in hatcheries.
  • US2008/0149033 concerns a method and apparatus for candling incubated eggs, wherein eggs are designated as live, inverted, or having a side air cell. Summary
  • an active thermography method for identifying at least one internal or subsurface parameter of at least one egg, the method comprising: thermally exciting at least a surface of the at least one egg with a means for excitation such that the egg is provided in a non-equilibrium thermal condition for a non-equilibrium period having duration transient / obtaining a thermogram at a plurality of time intervals (N) of the at least one egg in during the non-equilibrium period, resulting in N thermograms; determining at least one parameter in dependence upon the N thermograms; identifying the at least one internal or subsurface quality parameter of the at least one egg based on the at least one parameter.
  • the obtaining the thermogram at a plurality of time intervals of the at least one egg can result in a temperature data cube comprising the N thermograms.
  • Determining at least one parameter in dependence upon the N thermograms may comprise transforming a temperature-time signal of each spatial pixel of the datacube into a corresponding phase-frequency signal, resulting in a phase data cube comprising N phase images and identifying the at least one internal or subsurface quality parameter of the at least one egg based on the at least one parameter may comprise identifying the at least one internal or subsurface quality parameter of the at least one egg based on at least one of the phase images.
  • the method may further comprise a segmentation step, said segmentation step comprising segmenting at least one phase image based on at least one thermogram.
  • the segmentation step may comprise: extracting a mask of the at least one egg, said mask extraction comprising using at least a first thermogram and a second thermogram, wherein the first thermogram is taken immediately before thermal excitation and the second thermogram is taken immediately after thermal excitation; and applying said egg mask to a selected phase image.
  • the method may further comprise determining a time region of interest in the phase data cube and selecting the selected phase image as a phase image within the time region of interest.
  • the means for excitation may be hot air.
  • the means for excitation may be at least one source configured to thermally heat or cool the at least one egg.
  • the means for excitation may be an optical source, a mechanical source, or an inductive excitation means.
  • the method may further comprise a transforming step, said transforming step configured to transform the at least one thermogram in phase and/or amplitude data.
  • thermograms may be obtained at a time interval during a period t.
  • the plurality of thermograms may be provided in a 3D data cube.
  • the plurality of thermograms may be obtained at regular or irregular time intervals.
  • the plurality of thermograms may be obtained at a combination of regular and irregular time intervals.
  • the means for excitation may be configured to provide an excitation pulse for a period (tpulse).
  • the means for excitation may be configured to provide a periodic signal.
  • Se may be maximal 5 second and preferably 3 seconds and more preferably less than 2 seconds.
  • An air cell size of the at least one egg may be extracted as internal quality parameter troue may be 1 second.
  • the first of the plurality of thermograms may be obtained 3 seconds after thermal excitation is ended.
  • the plurality of thermograms may be obtained at a sampling frequency of 60 Hz or higher.
  • thermography to detect or monitor internal or subsurface quality parameter of at least one egg.
  • the active thermography may be pulsed active thermography. Also, the active thermography can be lock-in thermography using a periodical signal.
  • the internal or subsurface quality parameter of the at least one egg may be an air chamber.
  • an active thermography assembly for identifying at least one internal or subsurface quality parameter of at least one egg, said assembly comprising:
  • processing means configured to identify the at least one internal or subsurface quality parameter of the at least one egg based on the at least one thermogram.
  • the assembly may further comprise at least one optical filter configured to cancel out at least a part of the reflected radiation coming from the at least one egg.
  • the at least one egg is in particular an unfertilized egg, in particular a dead egg, free from any embryos.
  • the egg can be a poultry egg.
  • the at least one egg is preferably produced for direct or indirect human consumption.
  • the at least one egg is preferably an uncooked egg, i.e. a raw egg.
  • An egg shell of the egg is preferably substantially intact.
  • Figure 1 is a schematic cross-section of an egg
  • Figure 2 is a schematic plan view of an apparatus for performing active thermography
  • Figure 3 is a flow chart of a method of identifying the size of an air cell in an egg
  • Figure 4 is a schematic representation of a data cube
  • Figure 5 is an example phase image of an egg
  • FIGS. 6a - 6d illustrate stages in the segmentation of a thermogram
  • Figures 7a - 7d illustrate stages in the segmentation of a phase image
  • Figure 8 is an example phase profile for an air cell pixel and an egg pixel
  • Figure 9 is a flow chart of a method of determining optimal parameters for determining an air cell size
  • Figure 10 illustrates boxplots of F values for each of 10 Fourier frequencies for a pulse duration of 1 second and a cool down time of 7 seconds;
  • Figure 11 illustrates boxplots of F values for pulse duration Is, for six cool down time values, each cool down time value being associated with a Fourier frequency;
  • Figure 12 illustrates boxplots of F values for five pulse duration values, each pulse duration value being associated with a Fourier frequency and a cool down time value;
  • Figure 13 illustrates a sphere and a spherical cap thereof
  • Figure 14 shows an egg and a circle which fits the contour of the egg.
  • the egg 1 may be, for example, a chicken egg or a duck egg.
  • the egg 1 includes a yolk 2 which is surrounded by an albumen 3.
  • the albumen 3 is enclosed by an inner shell membrane 4.
  • the egg includes an external shell layer 5 (or “shell layer”) and an outer shell membrane 6 between the external shell layer 5 and the inner shell membrane 4.
  • the egg 1 includes an air cell 7 between the inner shell membrane 4 and the outer shell membrane 6.
  • the air cell 7 is formed when the egg 1 cools down after being laid, as a result of contraction of the egg contents. This causes air from outside the egg to be drawn inside the egg, forming the air cell 7.
  • the components of the egg 1 change in structure and composition following the laying of the egg.
  • the air cell 7 tends to increase in size due to moisture loss, and loss of carbon dioxide through pores of the egg shell 5.
  • the carbon dioxide can be replaced by air drawn into the egg 1 from outside the egg, which may cause the volume (or size) of the air cell 7 to increase.
  • Thermal properties of the air cell 7 may be different to thermal properties of the albumen 3. If the egg 1 is heated, the heating rate of the air cell 7 may be different to that of the albumen 3.
  • a thermal contrast may be observed between the air cell 7 and the albumen. In this way, the size of the air cell 7 may be determined.
  • FIG 2 a schematic plan view of a pulsed thermography apparatus 8 for investigating an egg 1 is shown.
  • the apparatus includes first and second excitation means in the form of first and second halogen lamps 9, 10 respectively.
  • the first and second halogen lamps 9, 10 may be, for example, Hedler H25 available from Hedler Systemlicht GmbH, Germany.
  • the first and second lamps 9, 10 are positioned so as to illuminate the egg 1 with excitation radiation 11 from first and second directions 12, 13 respectively.
  • the first and second lamps 9, 10 may include reflectors (not shown) and frosted glass diffusers (not shown) which can help to provide uniform illumination of the egg 1.
  • the apparatus 8 includes detection means in the form of an infrared camera 14 arranged to receive radiation 15 which may be reflected, transmitted, and/or radiated by the egg 1.
  • the detection means is capable of detecting and optionally recording thermal images, or thermograms.
  • One suitable detection means is the FUR SC7600 available from FUR Systems Inc., USA.
  • the apparatus 8 includes a filter 16 arranged between the egg 1 and the infrared camera 14.
  • the filter 16 is configured to block direct reflections of excitation radiation 11 from first and second lamps 9, 10.
  • the filter 16 may be, for example, a high pass optical filter.
  • One suitable filter is available from Spectrogon US Inc., having a cut-off frequency of 4080 nm. However, the choice of suitable filter depends on the heat sources and camera specifications.
  • the camera 14 and the egg 1 maintain the same mutual orientation during the obtaining of a respective N thermograms of that egg 1.
  • the camera is not substantially moved (remains stationary) during the obtaining of the respective N thermograms of that egg 1; similarly, the egg is not substantially moved (remains stationary) during the obtaining of the respective N thermograms of that egg 1.
  • the apparatus 8 includes a control module 17.
  • the first and second lamps 9, 10 and the infrared camera 14 are coupled to the control module 17.
  • the control module 17 is configured to provide activation signals to first and second excitation means 9, 10 and to receive data from detection means 14.
  • the control module 17 may also be configured to provide activation signals to detection means 14, for example, if it is preferred that the detection means 14 is only active (for example, detecting a thermal image) at specified times.
  • the control module 17 is further configured to process data received from detection means 14 as will be described in further detail in the following.
  • the control module 17 may include, for example, an optical excitation module 18 (OTvis 4000, edevis GmbH, Germany) configured to synchronise and control the first and second lamps 9, 10 and infrared camera 14.
  • the control module 17 may include memory 19 configured to store data received from infrared camera 14 and processor 20 configured to process data received from infrared camera 14.
  • the memory 19 may also store processed data received from processor 20.
  • the control module 17 provides an activation pulses to the first and second halogen lamps 9, 10, causing the lamps to emit an excitation pulse in the form of radiation 11 which has the same duration as the activation pulse.
  • the activation pulse has an activation pulse duration (or on-time) of t pU
  • the activation pulse duration may be, for example, in the range of 0.5 seconds to 5 seconds. However, the activation pulse duration may be longer, for example up to 10 seconds, up to 20 seconds, up to 30 seconds, or longer than 30 seconds.
  • Energy provided by an excitation pulse 11 from first and second lamps 9, 10 is absorbed at the surface (or shell 5) of the egg 1, and subsequently propagates by conduction under the surface, or shell 5, of the egg 1.
  • the presence of an air cell 7 reduces the thermal diffusion rate, so that surface or shell temperature above the air cell 7 reaches a higher temperature than the surrounding surface or shell and remains at a higher temperature during a cooling phase which follows the heating phase initiated by the radiation pulse 11. This non-equilibrum thermal condition persists for a time t transient .
  • time evolution of surface temperature of the egg 1 may be monitored after thermal excitation by radiation pulse 11.
  • An activation pulse is applied to first and second halogen lamps 9, 10.
  • the activation pulse causes the first and second halogen lamps 9, 10 to each emit an excitation pulse of radiation of the same duration as the activation pulse.
  • the excitation pulses are incident on the egg 1 and cause an internal surface of the egg 1 to be thermally excited (step S301).
  • a plurality of thermograms are obtained by the infrared camera 14 at times f, spaced apart by At (step S302), in particular without moving the egg 1 and camera 14 with respect to each other.
  • a total number N of thermograms is obtained during the non-equilibrium period t transient .
  • the N thermograms are obtained from the same side of the egg 1, in particular a side that faces respective detecting means 14.
  • the N thermograms are transformed into N phase images (step S303).
  • the transformation step S303 proceeds as follows.
  • the N thermograms form a data cube having dimensions x c y c N, where x is the number of pixels along a first spatial dimension of a thermogram (for example, the horizontal dimension of the thermogram) and y is the number of pixels along a second, perpendicular spatial dimension of the thermogram (for example, the vertical dimension of the thermogram).
  • a pixel at x, y, in the /rth thermogram may be uniquely specified by the coordinates (x,> y s -, k).
  • a pixel / located at position x,> y has a temperature value for each of the N thermograms, that is, a temperature value at each of the N times f, at which a thermogram is obtained by detecting means 14.
  • each position x and y is associated with a sequence, or one-dimensional array, of N temperature values T .
  • a discrete Fourier transform is applied individually to each array T of temperature values, producing an array Z of complex numbers associated with each position x,> y,.
  • the nth element of an array Zy is given by
  • phase data cube comprises N phase images, each having dimensions x c y-
  • phase image of an egg 1 is shown.
  • the phase image is at a Fourier frequency of 0.25Flz with a Is pulse duration and a cooling down phase of 3s duration.
  • the egg 1 is supported by an egg holder 20 and is placed against a background 21 of a wooden plate.
  • the air cell 7 can be seen as an area of high phase contrast in comparison with the rest of the egg 1 and the background 21.
  • the size of the air cell 7 may be determined by, for example, counting the number of pixels in the high-contrast air cell region. Other methods of determining the size of the air cell 7 are possible. For example, a contour fitting method may be performed.
  • a segmentation step may also be performed in order to segment the egg from the background.
  • the segmentation step uses a first, 'cold' thermogram, taken before thermal excitation of the egg by an excitation pulse, and a second thermogram taken immediately after thermal excitation of the egg, that is, immediately after the excitation pulse finishes.
  • the second thermogram may also be referred to as an early recorded thermogram, or ERT.
  • thermogram a cold thermogram is shown.
  • Figure 6b an early recorded thermogram is shown.
  • the temperature in Celcius is indicated by temperature scale to the right of each image.
  • the egg surface was heated by 2 to 3 degrees Celcius after a heating pulse of duration 5 seconds.
  • the cold thermogram is subtracted from the early recorded thermogram to obtain a corrected thermal image.
  • the position of the air cell 7 is visible in this corrected thermal image.
  • the egg holder 20 can also be seen. Flowever, the contrast between the air cell and the egg, and between the air cell and the background, is low. Therefore, phase images are used to enhance the contrast between the air cell and the egg.
  • an egg mask is obtained using the corrected thermal image ( Figure 6c).
  • Otsu's segmentation method is applied to the corrected thermal image ( Figure 6c) and the result is shown in Figure 6d.
  • any image segmentation method is contemplated.
  • the image is separated into three regions, namely the background 21, the egg 1, and the egg holder 20. From this image, a mask of the egg is extracted using a blob extraction method. The second largest blob of connected pixels is assigned to the egg and is used as the mask of the egg.
  • the background is the largest region while the egg holder is the smallest region, and the middle region is thus assigned to the egg.
  • an ellipse detection method can also be used.
  • the largest blob of interconnected pixels will be assigned as the egg.
  • Image closing techniques may additionally be used to correct for noise.
  • a further segmentation step may be performed to separate the air cell. Due to the low contrast between the egg and the air cell in the corrected temperature image ( Figure 6c), phase images are preferably used to extract a mask of the air cell.
  • FIG. 7a an unprocessed phase image is shown.
  • the egg holder 20 can be seen at the base of the egg 1.
  • the egg mask is applied to the unprocessed image of Figure 7a and the result is shown in Figure 7b.
  • the histogram has been rescaled to take values between 0 and 1 and the background is set to value zero.
  • Otsu's segmentation method is applied to the masked image of 7b.
  • the segmentation either separates the egg 1 from both the background 21 and the air cell 7, or separates the air cell 7 from both the egg 1 and the background 21.
  • the mask of the egg, obtained as described in the preceding, is then applied to the result of the segmentation, and the resulting image is shown in Figure 7c.
  • the air cell 7 is extracted from the image of Figure 7c by a blob extraction method and the resulting image is shown in Figure 7d.
  • the largest region of interconnected pixels is assigned to the egg and the second-largest region of interconnected pixels is assigned to the air cell.
  • An air cell mask is thus provided as the region of interconnected pixels assigned to the air cell. Any noise present may be assigned to the egg region, which can be a preferred method of cancelling out noise as compared with, for example, using a median filter.
  • a restriction on distance of the air cell from the blunt end of the egg may be implemented in order to cancel out false positive detections of air cells.
  • the centroid of the air cell needs to be within a specified range of distances from the blunt end of the egg in order to be recognized as an air cell.
  • a restriction on the minimum number of pixels required to be assigned to the air cell may also be imposed. This can help to avoid noise being selected as being the air cell.
  • a restriction on the ratio of the size of the blob to the size of the egg can also be imposed. For example, if the number of pixels assigned to the air cell is greater than 50% of the number of pixels assigned to the egg, the segmentation may be determined to be unsuccessful. An unsuccessful result does not necessarily indicate that no air cell is present in the egg; it may be the case that an air cell is present and is not visible from the side of the egg which is facing the detection means.
  • the size of the air cell may be determined by counting the number of pixels assigned to the air cell in the blob extraction method described in relation to Figure 7d, that is, by counting the number of pixels comprising the air cell mask. Based on a spatial calibration, this number of pixels can be translated into a measure of air cell size. For example, a spatial calibration may result in a conversion factor between a single thermogram or phase image pixel and a distance.
  • phase images tend to be less influenced by geometrical features and less sensitive to non-uniform heating than thermal images and amplitude images.
  • This can be of use when analyzing thermograms of objects with complex geometries, such as eggs: the directional emissivity of a patch of surface on the object will decrease with increased viewing angle (that is, the angle between the thermal camera and the normal of the surface patch). It has been shown that phase images are much less affected by the unwanted effects of shape variation than amplitude images and thermal images.
  • step S303 of the above-described method is related to transformation of thermograms into phase images
  • the present invention is not limited to transformation into phase images.
  • Step S303 of the present invention comprises calculating at least one parameter based on the N thermograms.
  • the parameter is the phase, but in other embodiments of the present invention, the parameter may be different to the phase, for example in thermal signal reconstruction (TSR) methods.
  • TSR thermal signal reconstruction
  • the parameter is preferably calculated such that undesired influences such as uneven excitation or curvature effects do not significantly disturb the result. It will be clear to the person skilled in the art that other mathematical operations can deliver the same desired properties as the phase images.
  • TSR the time-temperature profile of each individual pixel is modelled. This is typically performed using a polynomial relating log(T) to log(t).
  • Step S304 then comprises identifying air cell size based on one or more parameters calculated in dependence upon the acquired thermograms.
  • phase image in the phase data cube may be selected for performing the segmentation steps to provide the air cell mask and thus estimate the size of the air cell.
  • the contrast between the air cell and the rest of the egg may vary between phase images in the sequence of N phase images. Choice of the phase image to be used can thus affect the accuracy of the determination of the air cell size.
  • phase profile 30 of an air cell pixel and an example phase profile 31 of an egg pixel are shown. Both phase profiles 30 and 31 take a negative value at lower frequencies and a positive value at higher frequencies.
  • the Fourier frequency at which each phase profile 30, 31 changes from negative to positive is different and may depend on parameter settings such as the excitation pulse duration t pu/se , the sampling frequency at which the N thermograms are obtained, and the cooling down time.
  • the cooling down time is the time period during which thermal images are acquired after the end of the excitation pulse.
  • the phase profile 30 of the air cell pixel changes sign at a Fourier frequency f 1 that is smaller than the Fourier frequency f 2 at which the phase profile 31 of the egg pixel changes sign.
  • the frequency range between fa and f 2 is the frequency range in which the largest contrast is seen between the air cell 7 and the egg 1. Successful segmentation may be more likely when performed on a phase image selected from within this optimal frequency range.
  • the method of identifying the size of an air cell contained in an egg may include a step of selecting a phase image within a Fourier frequency range specified by a first frequency at which the phase profile of an air cell pixel changes sign, and a second frequency at which the phase profile of an egg pixel changes sign.
  • parameters may affect the contrast between the egg and the air cell.
  • parameters are the excitation pulse duration, the sampling frequency, and the cooling down time.
  • a preferred or optimal parameter combination may be determined in order to provide phase images with an increased contrast between the egg and the air cell.
  • thermograms five different pulse durations (PD) and seven different cooling down times (CDT) are considered and twelve eggs are used as test objects for obtaining thermograms.
  • the pulse durations are 0.5s, Is, 2s, 3s, 5s.
  • the cool down times are 0, 0.5, 1, 2, 3, 7, 10.
  • different numbers and/or combinations and/or values of parameters may be compared in a similar manner.
  • the number of test objects is not limited to twelve.
  • a range of sampling frequency values may be assessed.
  • the sampling frequency is set to 60Hz.
  • a high sampling frequency for example 60Hz or higher, is preferable if the available computing power is adequate for processing the resulting larger datacube (for higher sampling frequencies, the number N of thermograms will be higher). This is because datasets taken at higher sampling frequencies tend to exhibit less noise.
  • the air cell segmentation process described hereinbefore is performed for each Fourier frequency within each of the 12 data cubes for each of the 35 parameter combinations.
  • the segmentation method described hereinbefore results in two groups of pixels within the mask of the egg: pixels that belong to the air cell and pixels that do not belong to the air cell. Each pixel within each group has an associated phase value. It is important for successful segmentation that there is a clear distinction between the two groups.
  • the groups may be compared with each other by using a one-way ANOVA F test.
  • the phase value of each pixel within a group is considered a sample of that group.
  • the F test compares the variance between each group with the variance within each group for one phase image which has been segmented.
  • the F test results in an F value for each Fourier frequency within each of the 12 data cubes for each of the 35 parameter combinations.
  • the F value quantifies the contrast between the air cell and egg in phase images. Thus, a high F value is preferred.
  • each group of F values (a group of F values being the set of 12 F values at a particular Fourier frequency, PD value, and CDT value) are removed, resulting in a modified group of F values.
  • a group of F values being the set of 12 F values at a particular Fourier frequency, PD value, and CDT value
  • a boxplot is constructed for each modified group of F values. This can help to enable comparison between modified groups of F values in order to determine the Fourier frequency, for each combination of PD and CDT, at which the segmentation is most successful.
  • the comparison can include comparing the median and the interquartile range of each modified group of F values. A high median and a small interquartile range is preferred.
  • Figure 10 shows F value boxplots for each of 10 Fourier frequencies for a pulse duration of 1 second and a cool down time of 7 seconds.
  • the boxplot at 0.25Flz has both a high median F value and a small interquartile range in comparison with the boxplots at other Fourier frequencies.
  • the 0.25Flz frequency may thus be determined to be the optimal Fourier frequency for this parameter combination.
  • one boxplot at a specific Fourier frequency is chosen for use in further analysis.
  • the Fourier frequency of the chosen boxplot for one combination of PD and CDT will not be the same as the Fourier frequency of the chosen boxplot for a different combination of PD and CDT.
  • step S901 a specific Fourier frequency is determined for each combination of PD and CDT.
  • the CDT value with the optimal segmentation is determined (step S902).
  • the seven boxplots having specific Fourier frequency values are compared, each boxplot being associated with a different CDT value.
  • the median and interquartile range are compared within each PD value and the CDT value showing the highest median and smallest interquartile range is selected for further analysis.
  • Figure 11 shows boxplots of the ANOVA F test for pulse duration Is, for six CDT values, each CDT value being associated with a Fourier frequency as selected in step S901.
  • the tick labels on the x axis are in the format: cool down time (Fourier frequency).
  • the CDT is given in seconds and the Fourier frequencies are given in Hz.
  • Figure 11 does not show a boxplot for a CDT of 0 seconds. This is because it was not possible to detect the air cell at this value of cool down time.
  • the cool down time of 3 seconds and corresponding Fourier frequency of 0.25 Hz were determined to be optimal for the pulse duration of 1 second.
  • step S902 an optimal CDT value is determined for each PD value.
  • the final step is to compare the PD value boxplots (step S903).
  • the median and interquartile range are compared between the five PD values and the PD value showing the highest median and smallest interquartile range is selected.
  • This combination of PD value, CDT value, and Fourier frequency is then determined to be the optimal combination for successful segmentation of the air cell.
  • Figure 12 shows boxplots of the ANOVA F test for five PD values, where each PD value is associated with a Fourier frequency and a CDT value as selected in steps S901 and S902 respectively.
  • the tick labels on the x axis are in the format: pulse duration (cool down time, Fourier frequency).
  • the PD and CDT are given in seconds and the Fourier frequencies are given in Hz.
  • the boxplots in Figure 12 it can be seen that it is advantageous to increase the pulse duration from 0.5s to Is.
  • the total recording time of the pulse groups in Figure 12 are from left to right 3.5s, 4s, 12s, 10s and 12s. The total recording time thus increases substantially for the three highest tested pulse durations. Therefore, it may be preferred to select the pulse duration of Is.
  • the boxplots of the five parameter combinations shown in Figure 12 may be analyzed in more detail as a function of the Fourier frequencies. Parameter combinations in which the two intersection points f f 2 are very close to each other should be avoided.
  • the physical attributes of an egg depend on various parameters which cannot be controlled, and the optimal parameters in this example were determined for a group of 12 eggs. It cannot be guaranteed that the optimal frequency will be the same for a second, different group of eggs. A slight shift of the optimal frequency is possible if a different group of eggs is being tested. Therefore, a broad peak or a platform of optimal frequencies is preferred, so that when moving between groups of eggs the pulse duration and cool down time do not need to be adjusted.
  • the Is pulse duration time is chosen as the ideal setting for pulse duration since it gives a slightly higher overall F-value, and since the total recording time is only increased by 0.5s in comparison to the pulse duration of 0.5s.
  • the ideal Fourier frequency is taken at 0.225 Hz instead of 0.25Flz as the 0.225 Hz is surrounded by two frequencies for which the F-value is not significantly different.
  • the optimal parameter settings were determined to be Is pulse duration, 3s cooling time, 60 Hz sampling rate, and the phase images at 0.225 Hz were used for the extraction of the air cell.
  • Determining an optimal parameter combination can be useful in applications wherein continual monitoring of air cell size is required. Continued reliable segmentation of air cells of eggs without requiring adjustment of parameters can provide an efficient and effective method of monitoring quality parameters, for example the air cell size, of a plurality of eggs.
  • Each egg was manually turned and monitored on four sides of the egg at angles of 0°, 90°, 180° and 270° relative to the initial orientation of the egg. Each egg was rotated around its axis of rotational symmetry. At the same time, the eggs were weighed on a O.lmg precision digital scale (FR-300 MK2, A&D Company, Ltd., Japan). The increase in volume of the air cell corresponds with a moisture loss and thus with a weight loss. The weight loss was used as the 'gold standard' in this experiment. On the basis of the water density, the weight loss can be converted to a volume increase. On the last day of the experiment the dimensions (i.e.
  • the air cell size was monitored during storage in two different ways.
  • the number of pixels of the air cell was determined for each of the four views of the egg, and these values were summed to obtain an overall air cell size.
  • the cumulative increase in a number of pixels can then be linked to the cumulative weight loss of each egg.
  • the volume and height of the air cell are estimated using geometry considerations. These two variables can be estimated if we assume that the blunt end of an egg approaches a spherical shape and that the air cell can be approached as a spherical cap of this sphere. Approximating the air cell as a spherical cap implicitly assumes that the inner shell membrane by which the air cell is formed has no curvature. Referring to Figure 13, the dimensions of a spherical cap 40 can be described by the parameters a, r and h for which a is the base of the spherical cap, r is the radius of the sphere, and h is the height of the spherical cap.
  • V —(3a 2 + h 2 )
  • the a parameter is extracted from the mask of the air cell.
  • the longest length of this mask was set to equal 2a.
  • the length of the mask was measured by the function "regionprops.m" within Matlab. To estimate the volume of the air cell as accurately as possible, it is important to estimate carefully the radius of the sphere that surrounds the air cell. A simple estimation of the curvature of the blunt end of the egg would not be sufficiently accurate as the air cell is not always centered at the blunt end.
  • One method of estimating r is to first extract the contour of the egg mask, and calculate the centroid coordinates of the air cell.
  • the contour of a mask is extracted by Matlab Function "bwboundaries.m” while the centroid of a mask can also be calculated by using "regionprops.m”.
  • the nearest point (xl,yl) of the egg contour to the centroid of the air cell was calculated with the Euclidian distance. This point was extracted from the contour, together with a segment of the contour around this point.
  • the two ends of this contour segment supplied two further points (x2,y2) and (x3,y3).
  • An example of this is illustrated in Figure 14 where the 3 points 41, 42, 43 are indicated with three red crosses. These three points belong to the circle 44 of which the radius is to be estimated.
  • a system of three equations may be formed and solved to estimate the radius r of the circle and to calculate the center (x0,y0) of the circle:
  • the volume V of the air cell can now be estimated by using parameters r and a as described in the above equations.
  • the final step is the conversion of dimensions.
  • Dimensions of the egg for example the width and/or the length of the egg as measured with a caliper, may be used to convert the parameters V, r, a, h from units of pixels to SI units.
  • Estimated volume increase of an air cell may be correlated with weight loss of the egg.
  • the excitation means may be any means configured to provide the egg in a non-equilibrium thermal condition for a period of time.
  • an optical excitation source may be used such as a photographic flash, IR lamp, or laser excitation by point or line scanning.
  • a mechanical excitation source may be used such as sound or ultrasound excitation, preferably causing internal stresses and thus heat (vibrothermography).
  • An inductive excitation source may be used such as one which causes eddy currents (eddy current thermography). Any appropriate means is contemplated including hot air or microwave heating.
  • An excitation means can be a heating or cooling source which can create a temperature difference or provide the egg in a non-equilibrium thermal condition for a period of time.
  • thermograms are described as being acquired at equally, or regularly, spaced points in time. However, the thermograms may be acquired at points in time which are not equally spaced.
  • the time spacing may be relatively shorter in the initial, fast part of cooling following excitation and the time space may be relatively longer as the cooling rate decreases.
  • the time spacing may increase or decrease monotonically as time elapsed since the excitation increases.
  • an egg may undergo a decrease in temperature as a result of excitation.
  • the excitation means is a source of cool air at a first temperature T 1 and the egg is initially in thermal equilibrium with surroundings at a second temperature T 2 , where T 2 is greater than T lt then the excitation may result in the temperature of the egg being decreased.
  • the at least one egg may comprise a plurality of eggs, in which case the active thermography method comprises obtaining a thermogram at a plurality of time intervals (N) of each egg of the plurality of eggs, resulting in N respective thermograms for each egg of the plurality of eggs.
  • a single thermogram may relate to, e.g. show, a single egg or a plurality of eggs, as long as at least one respective thermogram is obtained at a plurality of time intervals for each egg.
  • reference is made to "internal or subsurface” reference is made to existing or occurring within a body under its surface.
  • the biological sample is an egg
  • pulsed thermography reference is made to pulsed or flash thermography, whereby there is no restriction on the length of the pulse or flash; the duration of the pulse can be very short (flash) or longer.
  • the duration of the pulse can be at least 1 second, at least 5 seconds, at least 10 seconds, at least 20 seconds, or at least 30 seconds.
  • the duration of the pulse may be no more than 2 second, no more than 3 seconds, or no more than 5 seconds.
  • active thermography reference is made to using a means for excitation (e.g. an energy source) to produce a thermal contrast between different features of interest/areas of interest within the sample (e.g. egg).
  • excitation e.g. an energy source
  • the inspected eggs are usually in thermal equilibrium with the surroundings before excitation; the excitation disturbs the thermal equilibrium with the surroundings.
  • hot air reference is made to using air which has been heated and which is configured to heat up at least the internal or subsurface of the egg.
  • the heating of at least the internal of the egg is preferably kept to a minimum to avoid quality loss. Therefore, in preferred embodiments, a maximum heating of 2°C, preferably of 1°C is provided for the internal or subsurface of the egg.
  • thermoogram reference is made to an image (2D) of the temperature calculated or determined based on the radiation emitted by an object.

Landscapes

  • Engineering & Computer Science (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • Health & Medical Sciences (AREA)
  • Physics & Mathematics (AREA)
  • General Physics & Mathematics (AREA)
  • Chemical & Material Sciences (AREA)
  • Food Science & Technology (AREA)
  • Pathology (AREA)
  • Environmental Sciences (AREA)
  • Analytical Chemistry (AREA)
  • Biochemistry (AREA)
  • General Health & Medical Sciences (AREA)
  • Immunology (AREA)
  • Theoretical Computer Science (AREA)
  • Computer Vision & Pattern Recognition (AREA)
  • Medicinal Chemistry (AREA)
  • Animal Husbandry (AREA)
  • Biodiversity & Conservation Biology (AREA)
  • Quality & Reliability (AREA)
  • Investigating Or Analyzing Materials Using Thermal Means (AREA)
  • Investigating Or Analysing Materials By Optical Means (AREA)
  • Radiation Pyrometers (AREA)

Abstract

L'invention concerne un procédé de thermographie active permettant d'identifier au moins un paramètre interne ou de sous-surface d'au moins un œuf. Le procédé comprend l'excitation thermique d'au moins une surface dudit œuf à l'aide d'un moyen d'excitation de telle sorte que l'œuf soit soumis à une condition thermique de non-équilibre pendant une période de non-équilibre présentant une durée ttransient ; l'obtention d'un thermogramme à une pluralité d'intervalles temporels (N) dudit œuf pendant la période de non-équilibre, permettant l'obtention de N thermogrammes ; la détermination d'au moins un paramètre en fonction des N thermogrammes ; et l'identification desdits paramètres de qualité interne ou de sous-surface dudit œuf en fonction desdits paramètres.
PCT/EP2019/084077 2018-12-07 2019-12-06 Procédé d'identification d'un paramètre interne d'un œuf WO2020115316A1 (fr)

Applications Claiming Priority (2)

Application Number Priority Date Filing Date Title
GBGB1819994.3A GB201819994D0 (en) 2018-12-07 2018-12-07 Method for identifying internal parameter of an egg
GB1819994.3 2018-12-07

Publications (1)

Publication Number Publication Date
WO2020115316A1 true WO2020115316A1 (fr) 2020-06-11

Family

ID=65030183

Family Applications (1)

Application Number Title Priority Date Filing Date
PCT/EP2019/084077 WO2020115316A1 (fr) 2018-12-07 2019-12-06 Procédé d'identification d'un paramètre interne d'un œuf

Country Status (2)

Country Link
GB (1) GB201819994D0 (fr)
WO (1) WO2020115316A1 (fr)

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN117256514A (zh) * 2023-11-17 2023-12-22 河北玖兴农牧发展有限公司 一种用于分选系统的剔除不合格种蛋的方法

Citations (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
WO2002061858A2 (fr) * 2000-11-17 2002-08-08 Thermogenic Imaging, Inc. Dispositif et procedes de mesure calorimetrique infrarouge
US20080149033A1 (en) 2006-12-21 2008-06-26 Embrex, Inc. Methods and apparatus for candling avian eggs via thermal cameras
US20100074515A1 (en) * 2008-02-05 2010-03-25 Kla-Tencor Corporation Defect Detection and Response
US20130027547A1 (en) * 2010-04-13 2013-01-31 Christian Homma Apparatus and method for projecting information onto an object in thermographic investigations
JP5967766B2 (ja) 2012-09-21 2016-08-10 株式会社ナベル 卵検査装置
WO2017137837A1 (fr) 2016-02-08 2017-08-17 Egg-Chick Automated Technologies Appareil et procédé de détection d'œufs disposés à l'envers
WO2018212087A1 (fr) * 2017-05-15 2018-11-22 三菱電機株式会社 Appareil d'inspection de défaut et procédé d'inspection de défaut

Patent Citations (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
WO2002061858A2 (fr) * 2000-11-17 2002-08-08 Thermogenic Imaging, Inc. Dispositif et procedes de mesure calorimetrique infrarouge
US20080149033A1 (en) 2006-12-21 2008-06-26 Embrex, Inc. Methods and apparatus for candling avian eggs via thermal cameras
US20100074515A1 (en) * 2008-02-05 2010-03-25 Kla-Tencor Corporation Defect Detection and Response
US20130027547A1 (en) * 2010-04-13 2013-01-31 Christian Homma Apparatus and method for projecting information onto an object in thermographic investigations
JP5967766B2 (ja) 2012-09-21 2016-08-10 株式会社ナベル 卵検査装置
WO2017137837A1 (fr) 2016-02-08 2017-08-17 Egg-Chick Automated Technologies Appareil et procédé de détection d'œufs disposés à l'envers
WO2018212087A1 (fr) * 2017-05-15 2018-11-22 三菱電機株式会社 Appareil d'inspection de défaut et procédé d'inspection de défaut

Non-Patent Citations (1)

* Cited by examiner, † Cited by third party
Title
F. FRENI ET AL: "Assessment of eggs freshness by means of pulsed infrared thermography", PROCEEDINGS OF THE 14TH INTERNATIONAL CONFERENCE ON QUANTITATIVE INFRARED THERMOGRAPHY, 1 January 2018 (2018-01-01), pages 1 - 8, XP055675649, DOI: 10.21611/qirt.2018.068 *

Cited By (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN117256514A (zh) * 2023-11-17 2023-12-22 河北玖兴农牧发展有限公司 一种用于分选系统的剔除不合格种蛋的方法
CN117256514B (zh) * 2023-11-17 2024-04-23 河北玖兴农牧发展有限公司 一种用于分选系统的剔除不合格种蛋的方法

Also Published As

Publication number Publication date
GB201819994D0 (en) 2019-01-23

Similar Documents

Publication Publication Date Title
US11353439B2 (en) Method for bone scan in meat
US9905014B2 (en) Method and system for the examination of a sample by means of thermography
KR101686766B1 (ko) 레이저 스페클을 이용한 세균 및 미생물 탐지 장치 및 방법
BRPI1103714A2 (pt) Aparelho para determinar a espessura variável de um revestimento disposto sobre um substrato e método para determinar a espessura de um revestimento em vários pontos ao longo da superfície de um sustrato
CN103196564B (zh) 一种通过图像分割修正表面发射率的红外热成像测温方法
US7822268B2 (en) Advanced processing of active thermography signals
US10687713B2 (en) Diagnostic apparatus
US8055054B2 (en) Method and apparatus for thermographic nondestructive evaluation of an object
US11226301B2 (en) Thermographic examination means and method for non-destructive examination of a near-surface structure at a test object
US20120256777A1 (en) Method for Identifying Materials Using Dielectric Properties through Active Millimeter Wave Illumination
CZ20013674A3 (cs) Měření tlouą»ky kovové vrstvy pomocí tranzientní termografie
Usamentiaga et al. A quantitative comparison of stimulation and post-processing thermographic inspection methods applied to aeronautical carbon fibre reinforced polymer
KR101434720B1 (ko) 3ｄ 스캐너
US11193756B2 (en) Object shape measurement apparatus and method
WO2020115316A1 (fr) Procédé d'identification d'un paramètre interne d'un œuf
CN110456010A (zh) 一种无损快速检测鸡蛋新鲜度的方法
RU2753620C1 (ru) Способ определения кинетических теплофизических свойств анизотропных композитных материалов
Vavilov et al. Detecting corrosion in thick metals by applying active IR thermography
Sárosi et al. Evaluation of reflectivity of metal parts by a thermo-camera
JP4517044B2 (ja) 欠陥検査方法およびその装置
WO2022091702A1 (fr) Dispositif et procédé d'analyse d'images thermiques infrarouges, et programme
Ivanescu et al. Vision system for human body infrared thermography
Sun et al. Online Measurement of Content in Block Frozen Shrimp by Fusion of X-ray Imaging and Binocular Vision
Tsarenko et al. Valuation techniques of freshness of poultry eggs
WO2023157004A1 (fr) Profilage, modélisation et surveillance de température et de flux de chaleur dans de la viande ou des produits aliments pendant un processus de cuisson

Legal Events

Date Code Title Description
121 Ep: the epo has been informed by wipo that ep was designated in this application

Ref document number: 19816693

Country of ref document: EP

Kind code of ref document: A1

NENP Non-entry into the national phase

Ref country code: DE

122 Ep: pct application non-entry in european phase

Ref document number: 19816693

Country of ref document: EP

Kind code of ref document: A1