US 8,731,295 B2Grant
Visual segmentation of lawn grass
Issue Date:2014-05-20
•20 Claims
•10 Drawing Sheets
Abstract
Disclosed is a method for identifying lawn grass which includes capturing an image of the terrain in front of a mower, segmenting the image into neighborhoods, calculating at least two image statistics for each of the neighborhoods, generating a binary representation of each image statistic. The binary representation of each image statistic is generated by comparing the calculated image statistic values to predetermined image statistic values for grass. The method further includes weighting each of the binary representations of each image statistic, and summing corresponding neighborhoods for all image statistics. A binary threshold is applied to each of the summed neighborhoods to generate a binary map representing grass containing areas and non-grass containing areas.
Metadata
Assignee
- MTD Products Inc.
Inventors
- Alexander Schepelmann
- Kathryn A. Daltorio
- Amaury D. Rolin
- Jonathan Beno
- Bradley E. Hughes
- James M. Green
- Michael S. Branicky
- Roger D. Quinn
- Henry H. Snow
- Francis L. Merat
- Richard E. Hudson
Application Information
Application Number:US 13/380,413
Filing Date:2010-07-01
Priority Date:2009-07-01
Art Unit:2666
Classifications
IPC:
G06K9/34
Patent Drawings (10 sheets)
Description
Cross Reference to Related Application
[0001] This application is a national phase entry of International Application No. PCT/US2010/001881 filed Jul. 1, 2010, which claims the priority benefit of U.S. Provisional Patent Application Ser. No. 61/269,962 filed Jul. 1, 2009; U.S. Provisional Patent Application Ser. No. 61/339,868 filed Mar. 10, 2010; and U.S. Provisional Patent Application Ser. No. 61/355,935 filed Jun. 17, 2010.
Field of the Invention
[0002] This invention relates generally to a method of identifying lawn grass. In particular, it relates to a method of identifying lawn grass for use with an autonomous lawnmower.
Background of the Invention
[0003] Autonomous robots are increasingly becoming part of our daily lives. Autonomous lawnmowers have recently gained popularity, yet, currently available consumer versions do not possess the ability to sense obstacles from a distance. It is known that LIDAR can sense obstacles from a distance; however, it is cost prohibitive to place a LIDAR unit on a consumer lawnmower. Accordingly, there is a need to develop a cost-effective obstacle detection method for an autonomous lawnmower.
Summary of the Invention
[0004] A method for identifying lawn grass comprising capturing an image of the terrain in front of a mower, segmenting the image into neighborhoods, calculating at least two image statistics for each of the neighborhoods, and generating a binary representation of each image statistic. The binary representation of each image statistic is generated by comparing the calculated image statistics to predetermined image statistics for grass. The method further comprises weighting each of the binary representations of each image statistic and summing corresponding neighborhoods for all image statistics. A binary threshold is applied to each of the summed neighborhoods to generate a binary map representing grass containing areas and non-grass containing areas.
Brief Description of the Drawings
[0005] These and other aspects of the invention will be understood from the description and claims herein, taken together with the drawings showing details of construction and illustrative embodiments, wherein:
[0006] FIG. 1 illustrates a side view of an autonomous robotic mower in accordance with an embodiment of the present invention.
[0007] FIG. 2 illustrates a front view of an autonomous robotic mower in accordance with an embodiment of the present invention.
[0008] FIGS. 3 a-b contain a flow chart in accordance with an embodiment of the present invention.
[0009] FIG. 4 illustrates the creation of a binary representation of an image statistic by comparing neighborhood values of the image statistic to a predetermined image statistic for grass in accordance with an embodiment of the present invention.
[0010] FIGS. 5 a-b illustrate the multiplication of resulting binary images for the neighborhood statistics created in block 375 by their normalized coefficients in accordance with an embodiment of the present invention.
[0011] FIG. 6 contains a flow chart in accordance with an embodiment of the present invention.
[0012] FIG. 7 contains a block diagram of a vision system in accordance with an embodiment of the present invention.
Detailed Description of the Invention
[0013] Approximating language, as used herein throughout the specification and claims, may be applied to modify any quantitative representation that could permissibly vary without resulting in a change in the basic function to which it is related. Accordingly, a value modified by a term or terms, such as “about”, is not limited to the precise value specified. In at least some instances, the approximating language may correspond to the precision of an instrument for measuring the value. Range limitations may be combined and/or interchanged, and such ranges are identified and include all the sub-ranges stated herein unless context or language indicates otherwise. Other than in the operating examples or where otherwise indicated, all numbers or expressions referring to quantities of ingredients, reaction conditions and the like, used in the specification and the claims, are to be understood as modified in all instances by the term “about”.
[0014] “Optional” or “optionally” means that the subsequently described event or circumstance may or may not occur, or that the subsequently identified material may or may not be present, and that the description includes instances where the event or circumstance occurs or where the material is present, and instances where the event or circumstance does not occur or the material is not present.
[0015] As used herein, the terms “comprises”, “comprising”, “includes”, “including”, “has”, “having”, or any other variation thereof, are intended to cover a non-exclusive inclusion. For example, a process, method, article or apparatus that comprises a list of elements is not necessarily limited to only those elements, but may include other elements not expressly listed or inherent to such process, method, article, or apparatus.
[0016] The singular forms “a”, “an”, and “the” include plural referents unless the context clearly dictates otherwise.
[0017] One embodiment of an obstacle detection method that can be used for vision processing in an autonomous lawnmower is to identify driveable terrain around the lawnmower through the identification of grass in images segmented into regions, or “neighborhoods”. Observed image regions that do not match the measured grass statistics are marked as “obstacle containing regions” and observed image regions that do match the measured grass statistics are marked as “grass containing regions.” This creates a binary freespace map around the robot for every camera frame, which is abstracted into range images.
[0018] Turning now to an embodiment of the invention as shown in FIGS. 3 a-b, in block 300, image data is captured of the terrain in front of the lawnmower. In the preferred embodiment, the data is acquired using a camera mounted on the top of the lawnmower elevated 0.75 meters at an angle of −45° to the horizontal. In test runs of the preferred embodiment, logged images included the following objects and conditions: 640×480 resolution images of fully illuminated grass, as well as 50×50 image regions taken at 640×480 resolution of each of the following obstacle surfaces: shaded grass, flowerbed edging, flowers, blue jeans, and a soccer ball. Statistical measurements were collected for a mixture of both illuminated and shaded obstacles to simulate common conditions an autonomous lawnmower may encounter. In one embodiment, the camera is an ImagingSource DFK21AF04 FireWire camera with a Kowa LM4PBR 4.0 mm F1.2 CS-Mount lens. The camera captures 640×480 resolution images at 30 frames/sec and is connected via firewire to a 1.83 GHz Mac Mini running Windows XP with 4 GB of memory. It is contemplated that a person having ordinary skill in the art can choose to use a different camera and/or hardware to process the images captured by the camera. FIGS. 1-2 show different views of an autonomous lawnmower 100 having camera 105 and an optional second camera 110. It is contemplated that some embodiments can have more than one camera.
[0019] In block 305, the hue component of the HSL color vector is calculated and extracted from the RGB image captured in block 300. Unlike RGB color values, hue values at corresponding saturation values have been found to be insensitive to shadows and changing lighting conditions. Since the raw camera image captured in block 300 are represented in the RGB color space, it is first necessary to calculate the hue component of the HSL color vector according to:
[0020]
where h is the normalized hue component and r, g, and b are the normalized components of the original RGB vector. The normalized RGB components are given by:
where h is the normalized hue component and r, g, and b are the normalized components of the original RGB vector. The normalized RGB components are given by:
[0021]
[0022] Since HSL color space can be represented by a cone where all possible hue values lie on the cone's plane base, hue values can be represented by an integer in the range [0, 360]. These values are represented as integers due to the way colors are represented in digital images. Since h is the normalized hue component, it is necessary to scale the normalized hue component h into the [0, 360] range through the equation:
[0023]
where (int) represents a typecasting operation of the scaled value into an integer and H is the scaled hue component of the HSL color vector. It is important to note that singularities occur when r, g, and b lie along the RGB color cube's grayline and therefore have the same numerical value. This results from color representation in the HSL color model where white, black, and grays have an undefined hue.
where (int) represents a typecasting operation of the scaled value into an integer and H is the scaled hue component of the HSL color vector. It is important to note that singularities occur when r, g, and b lie along the RGB color cube's grayline and therefore have the same numerical value. This results from color representation in the HSL color model where white, black, and grays have an undefined hue.
[0024] In block 310, the intensity plane is extracted from the image captured in block 310 by converting each pixel of the RGB image captured by the camera in block 300 into an unblurred grayscale image by extracting the intensity component I, of the HSL color model where I is given by the mean value of the RGB components for individual pixels, or:
[0025]
[0026] In block 315, the vertical edge response is extracted from the unblurred grayscale image created in block 310 by convolving the image within a vertical (Gy) Prewitt convolution kernel, respectively given by:
[0027]
[0028] In block 325, the horizontal edge response is extracted from the unblurred grayscale image created in block 310 by convolving the image with a horizontal (Gx) Prewitt convolution kernel, respectively given by:
[0029]
[0030] In both blocks 315 and 325, the convolution operation computes the edge strength based on grayscale intensity values of a pixel and its eight adjacent neighbors and plots the magnitude of the edge strength in the center pixel. The 2D spatial convolution operation is given by:
[0031]
where E(x,y) is the computed edge strength of a pixel located in the center of a 3×3 neighborhood located at (x,y) in the image, j→J−1 and k→K−1 are the indices of the convolution kernel G, and I(a,b) is the intensity value of a pixel located at the input coordinate pair (a,b). Convolving the image with Gx and Gy generates two new images that indicate edge response for the given filter, wherein the pixel intensity given by an integer value indicates the edge strength at that pixel location. We refer to the resulting images as “horizontal-” and “vertical texture images,” respectively. In one embodiment, the integer value is between [0, 255], but it is contemplated that a person having ordinary skill in the art can choose to use other ranges. In this application, we define “visual texture” as a collection of edges within an image region or “neighborhood.”
where E(x,y) is the computed edge strength of a pixel located in the center of a 3×3 neighborhood located at (x,y) in the image, j→J−1 and k→K−1 are the indices of the convolution kernel G, and I(a,b) is the intensity value of a pixel located at the input coordinate pair (a,b). Convolving the image with Gx and Gy generates two new images that indicate edge response for the given filter, wherein the pixel intensity given by an integer value indicates the edge strength at that pixel location. We refer to the resulting images as “horizontal-” and “vertical texture images,” respectively. In one embodiment, the integer value is between [0, 255], but it is contemplated that a person having ordinary skill in the art can choose to use other ranges. In this application, we define “visual texture” as a collection of edges within an image region or “neighborhood.”
[0032] In block 320, using the grayscale image produced in block 315, a vertical binary texture image is created by discarding edge strengths below an empirical threshold and setting corresponding pixel locations of edge strength above the threshold in a binary image equal to 1. Pixel groupings with an area of 1 in the binary image are removed. In one embodiment, the empirically determined threshold is 21, but it is contemplated that a person having ordinary skill in the art can choose to use a different threshold value. Unconnected pixels in the vertical binary texture image are also removed.
[0033] In block 325, using the grayscale image produced in block 330, a horizontal binary texture image is created by discarding edge strengths below an empirical threshold and setting corresponding pixel locations of edge strength above the threshold in a binary image equal to 1. Pixel groupings with an area of 1 in the binary image are removed. In one embodiment, the empirically determined threshold is 21, but it is contemplated that a person having ordinary skill in the art can choose to use a different threshold value. Unconnected pixels in the horizontal binary texture image are also removed.
[0034] In block 335, the corresponding horizontal and binary vertical texture images created in blocks 320 and 330 are combined to yield a directionally ambiguous binary texture image corresponding to a grayscale version of the image captured in block 300.
[0035] In block 340, the images produced in blocks 305, 315, 320, 325, 330, and 335 are each divided up into neighborhoods. In one embodiment, these neighborhoods range from 11×11 to 41×41 pixels in increments of 10 pixels, however it is contemplated that a person having ordinary skill in the art can choose to use a larger or smaller neighborhood. The operations that take place in block 340 are detailed below in blocks 345-375, in which the image statistics are calculated for each neighborhood of the image captured in step 300. The image statistics include the mean hue (block 345) and the following edge based texture statistics: the mean grayscale horizontal intensity value (block 350), mean grayscale vertical intensity value (block 350), vertical binary variance (block 355), mean vertical binary edge response area (block 355), horizontal and vertical neighborhood centroid location in directionally invariant binary texture image (block 360), horizontal binary variance (block 365), and mean horizontal binary edge response area (block 365).
[0036] Although this embodiment is discussed with reference to the neighborhood statistics listed in the paragraph above, it is understood that a person having ordinary skill in the art can choose to use more, less, or other image statistics.
[0037] In blocks 345, the mean hue of each neighborhood is calculated using the image produced in block 305. The mean hue value is calculated individually for each neighborhood by summing up all the mean hue pixel values for each neighborhood in the image generated in block 305, and dividing by the total number of pixels within the neighborhood. This is done for each neighborhood in the image. Accordingly, a mean hue image statistic value is calculated for each neighborhood in the image.
[0038] In block 350, the mean grayscale horizontal and vertical intensity value of each neighborhood is calculated. The mean grayscale horizontal intensity value is calculated individually for each neighborhood by summing up all pixel values for each neighborhood in the image generated in block 325 and dividing by the total number of pixels within the neighborhood. This is done for each neighborhood in the image. Accordingly, a mean grayscale horizontal intensity image statistic value is calculated for each neighborhood in the image.
[0039] Further, in block 350, the mean grayscale vertical intensity value is calculated individually for each neighborhood by summing up all the pixel values for each neighborhood in the image generated in block 315 and dividing by the total number of pixels within the neighborhood. This is done for each neighborhood in the image. Accordingly, a mean grayscale vertical intensity image statistic value is calculated for each neighborhood in the image.
[0040] In block 355 and 365, the variance of each neighborhood, Var(R), in the vertical binary texture image produced in block 320 and the horizontal binary texture image produced in block 330 are computed by:
Var(R)=(ΣX−X )/(n−1)
where X is the binary value of the current pixel,X is the mean value of the pixels within the square neighborhood, and n is the area of the neighborhood, given by the total number of pixels within the neighborhood. Since the variance of a binary neighborhood is being computed, this reduces to a polynomial in one dimension
Var(R)=(w−(w2/n))/(n−1)
where w is the number of pixels in the neighborhood above the binary threshold. Accordingly, a vertical binary variance image statistic value and horizontal binary variance image statistic value is calculated for each neighborhood in the image.
Var(R)=(ΣX−
where X is the binary value of the current pixel,
Var(R)=(w−(w2/n))/(n−1)
where w is the number of pixels in the neighborhood above the binary threshold. Accordingly, a vertical binary variance image statistic value and horizontal binary variance image statistic value is calculated for each neighborhood in the image.
[0041] Further, in blocks 355 and 365, the mean vertical binary edge response area and the mean horizontal binary edge response area is calculated for each neighborhood in the vertical binary texture image produced in block 320 and the horizontal binary texture image produced in block 330. These image statistics values are computed by:
[0042]
where the calculated moment Mij is given by the summation over the x and y dimensions of the pixel neighborhood, and i and j correspond to the index of the calculated moment. The binary area is given by M00. Accordingly, a mean vertical binary edge response area statistic value and horizontal binary edge response area statistic value is calculated for each neighborhood in the image.
where the calculated moment Mij is given by the summation over the x and y dimensions of the pixel neighborhood, and i and j correspond to the index of the calculated moment. The binary area is given by M00. Accordingly, a mean vertical binary edge response area statistic value and horizontal binary edge response area statistic value is calculated for each neighborhood in the image.
[0043] In block 360, the horizontal and vertical centroid location image statistic values for each neighborhood within the neighborhood, x and y , respectively, of the directionally ambiguous texture image generated in block 335 is calculated by:
x =M10/M00
and
y =M01/M00.
and
[0044] Accordingly, a horizontal and vertical centroid location statistic value is calculated for each neighborhood in the image.
[0045] Next, in blocks 370 and 375, the image statistics values for each neighborhood generated in blocks 345-365 are compared to predetermined neighborhood statistic values for grass. During this process, the value of each statistic in each individual neighborhood is compared to corresponding predetermined neighborhood statistics for grass. For each statistic, a binary (grass/not grass) representation of neighborhoods having mowable terrain is created. Accordingly, for each statistic, a neighborhood having a statistic value within a predetermined number of standard deviations of the corresponding predetermined image statistic value for grass is marked as grass (or assigned a “1”) and neighborhoods that fall out of this range are marked as obstacle containing (or assigned a “0”). In some embodiments, the predetermined number of standard deviations is 3, and in other embodiments, the predetermined number of standard deviations is 2, but it is contemplated that a person having ordinary skill in the art can choose another predetermined number of standard deviations. Accordingly, this creates a binary (grass/not grass) representation of driveable terrain of the neighborhoods.
[0046] This process is demonstrated in FIG. 4 , in which neighborhood values for statistic “a” are displayed for an image segmented into nine neighborhoods. The predetermined neighborhood statistic value for statistic “a” is 10, the standard deviation for “a” is 1, and the number of standard deviations being used is 3. Accordingly, a binary image is produced for statistic “a” in which all of the neighborhoods having values between 7 and 13 are assigned a “1” and all of the neighborhoods having values falling outside of the 7-13 range are assigned a “0”. For the purposes of this example, statistic “a” can be any statistic, such as mean grayscale vertical intensity value, vertical binary variance, or mean vertical binary edge response area.
[0047] The predetermined neighborhood statistic values used in one embodiment of the invention can be found in Table 1. However, it is contemplated a person having ordinary skill in the art can choose to use other statistical values.
[0048] The predetermined neighborhood statistic value for hue is not listed in the table, but in one embodiment of the invention that uses a 31×31 pixel neighborhood, the mean hue value is h =180.04, with a standard deviation of σ( h )=30.08.
[0049]
| TABLE I | ||||||
|---|---|---|---|---|---|---|
| NEIGHBORHOOD STATISTICS FOR NEIGHBORHOOD SIZE OF 31 | ||||||
| Measurement | Illuminated Grass | Shaded Grass | Flowerbed Edging | Flowers | Jeans | Soccer Ball |
| Ī(hTx,, vTx) | (−0.922, −10.169) | (−.027, −4.685) | (1.408, −5.644) | (−4.334, −13.333) | (0.251, −7.427) | (−4.739, −17.794) |
| σ(Ī)(hTx, vTx) | (5.053, 3.288) | (1.637, 1.225) | (5.749, 6.764) | (17.224, 12.480) | (8.535, 4.655) | (15.448, 12.688) |
| Var(R)(hTx, vTx) | (0.231, 0.225) | (0.164, 0.144) | (0.088, 0.124) | (0.172, 0.191) | (0.090, 0.065) | (0.068, 0.098) |
| σ (Var(R))(hTx, vTx) | (0.010, 0.014) | (0.039, 0.045) | (0.075, 0.066) | (0.0519, 0.0419) | (0.053, 0.067) | (0.056, 0.062) |
| Ā(hTx, vTx) | (350.590, 332.810) | (203.90, 172.80) | (103.80, 150.70) | (221.30, 261.20) | (101.40, 76.10) | (74.90, 112.40) |
| σ(Ā)(hTx, vTx) | (31.976, 39.559) | (59.056, 57.965) | (94.373, 92.161) | (76.288, 83.751) | (64.547, 87.778) | (65.829, 73.552) |
| ( x, y) | (15.368, 15.015) | (15.415, 14.570) | (7.465, 13.008) | (16.691, 15.120) | (15.527, 13.126) | (18.407, 14.952) |
| σ( x, y) | (0.606, 0.599) | (0.951, 0.722) | (2.564, 1.067) | (2.469, 2.000) | (2.473, 6.364) | (6.152, 5.322) |
[0050] The predetermined neighborhood statistic values were calculated from randomly selected neighborhoods that were randomly selected from images in two data sets taken between 9:00 AM and 11:59 AM on Aug. 31, 2009. All the neighborhoods contained only one type of surface. For example: neighborhoods containing illuminated grass contained only illuminated grass, shaded grass contained only shaded grass, and plastic fence neighborhoods contained only plastic fencing. Images from both sets were recorded on a plot of land which contained a mixture of rye, blue, and creeping red fescue grasses. These data sets were chosen for analysis for several reasons.
[0051] First, this grass mixture is the most common type of lawn grass found in the Midwest, and therefore results represent an accurate performance indicator of hue- and texture-based grass identification in a real-world environment. Second, data collection times are representative of typical times people tend to mow their lawns. Both data sets maintained constant lighting conditions, where the sun was visible at all times and never obscured by cloud cover when images were recorded.
[0052] By randomly selecting neighborhoods from both data sets, calculated statistics are valid for multiple times of day and are therefore not tuned to function only under specific, time dependent conditions.
[0053] Grass identification via some of these statistics was tested on 40 randomly selected samples of 31×31 neighborhoods containing: illuminated grass, shaded grass, artificial obstacles, and flowers. These results are tabulated in Table 2, which lists the number of correctly identified texture regions. Hue correctly identified 21/40 samples (52% accuracy) of shaded grass, 32/40 (80%) of artificial obstacles, and 16/40 (40% of flowers correctly. Illuminated grass was correctly identified with at least (38/40) 95% accuracy for all texture measurements except mean vertical grayscale intensity. Shaded grass was identified with at least 90% accuracy for all texture measurements except the mean vertical grayscale intensity. Binary horizontal variance and binary horizontal and vertical area identified obstacles correctly with at least 87.5% accuracy.
[0054]
| TABLE 2 | ||||||||
|---|---|---|---|---|---|---|---|---|
| NUMBER OF CORRECTLY IDENTIFIED TEXTURE REGIONS | ||||||||
| Measurement | ||||||||
| Hue | Var(R)hTx | Var(R)vTx | ĀhTx | ĀvTx | ( x, y) | ĪhTx | ĪvTx | |
| Illuminated Grass | 40 | 40 | 40 | 38 | 38 | 40 | 38 | 14 |
| Shaded Grass | 21 | 36 | 36 | 36 | 36 | 39 | 40 | 8 |
| Artificial Obstacles | 32 | 37 | 31 | 35 | 35 | 28 | 10 | 33 |
| Flowers | 16 | 11 | 5 | 9 | 9 | 21 | 14 | 27 |
[0055] As is demonstrated in Table 2 above, some image statistics exhibit better performance than others for recognizing grass. Accordingly, a weighting coefficient is applied to the neighborhoods of image statistics produced in block 375 prior to the weighted statistics being added together in block 385.
[0056] In one embodiment of the invention, the performance of each neighborhood statistic is estimated by its ability to correctly classify 40 sample neighborhoods each of illuminated grass, shaded grass, and obstacles. These estimates were then used to generate a weighting coefficient for individual neighborhood statistic:
[0057]
where α is the normalized weighting coefficient and Cillium, Cshaded, and Cobs are the number of correctly identified neighborhoods of illuminated grass, shaded grass, and obstacles respectively. Each weighting coefficient is normalized by dividing the coefficient by the total number of image measurements (40×3=120 in this embodiment) to generate a normalized weighting coefficient between [0,1].
where α is the normalized weighting coefficient and Cillium, Cshaded, and Cobs are the number of correctly identified neighborhoods of illuminated grass, shaded grass, and obstacles respectively. Each weighting coefficient is normalized by dividing the coefficient by the total number of image measurements (40×3=120 in this embodiment) to generate a normalized weighting coefficient between [0,1].
[0058] Table 3 contains pre-weight coefficients for each image statistic based on Table 2, and Table 4 contains the normalized pre-weight coefficients for six statistics:
[0059]
| TABLE 3 | ||||||||
|---|---|---|---|---|---|---|---|---|
| Pre-weight coefficients for each statistic. | ||||||||
| Weight Coeff. | Hue | Var(R)hTx | Var(R)vTx | ĀhTx | ĀvTx | ( x, y) | ĪhTx | ĪvTx |
| a | 0.7750 | 0.9416 | 0.8916 | 0.9083 | 0.9083 | 0.8916 | 0.7333 | 0.45833 |
[0060]
| TABLE 4 | ||||||
|---|---|---|---|---|---|---|
| Normalized pre-weight coefficients for six identifiers. | ||||||
| Normal- | ||||||
| ized | ||||||
| Coeff. | Hue | Var(R)hTx | Var(R)vTx | ĀhTx | ĀvTx | ( x, y) |
| aNorm | 0.1458 | 0.1771 | 0.1677 | 0.1703 | 0.1703 | 0.1677 |
[0061] In one embodiment, the resulting neighborhood statistics binary images created in block 375 are multiplied by their normalized coefficients and added together. More specifically, the resulting binary images for the neighborhood statistics are multiplied by their normalized coefficients. Then, corresponding neighborhoods are added with their applied weights to generate a grayscale representation of mowable terrain in the image, known as a probability map. The likelihood of a neighborhood containing grass is determined by the sum of the weighted statistics for the neighborhood.
[0062] In this embodiment, if a binary neighborhood is identified as grass containing for a neighborhood statistic, and the statistic's normalized weighting coefficient is 0.1, the likelihood of that neighborhood containing grass based on the neighborhood statistic increases by 0.1. This is repeated for all statistics for that neighborhood, such that the normalized weighting coefficients for the corresponding neighborhood statistics indicating “grass” are added together and multiplied by 1 and the corresponding normalized weighting coefficients for the corresponding neighborhood statistics indicating “not grass” are added together and multiplied by 0. The weighted “grass” and “not grass” results are added together and indicate the likelihood that the neighborhood contains grass.
[0063] This concept can be illustrated by looking at FIGS. 5 a-b in which five neighborhood statistics are used to identify grass in an image that is segmented into nine neighborhoods. Four of the neighborhood statistics are equally weighted and one has a weight twice that of the others. The neighborhood of interest is located at the upper left hand corner of the image and designated by the bold number.
[0064] For simplicity purposes, the statistics are identified as a through e. Looking at FIG. 5 a, statistics a through d have been assigned a weight of ⅙, and e has been assigned a weight of ⅓. Now, assume that statistics a, c, and e identify the neighborhood as grass, indicated by a binary value of 1. Neighborhood statistics b and d state that the neighborhood of interest is not grass, indicated by a 0.
[0065] FIG. 5 b illustrates the neighborhood statistics in the segmented binary texture images after the neighborhood statistics are multiplied by their normalized coefficients. Once the statistics for the neighborhood of interest are added up, the likelihood L, or probability, of that neighborhood being grass is:
[0066]
[0067] In another embodiment, also illustrated by FIG. 5 a, the likelihood L of the neighborhood of interest containing grass is calculated by adding up the normalized coefficients of the neighborhood statistics which identify the neighborhood of interest as grass:
[0068]
[0069] In block 390, using the neighborhood probability information produced in block 385, a binary grass image map is created by assigning a “0” to all of the pixels in corresponding neighborhood locations having a probability of containing grass below an empirical threshold as not containing grass (assigned a “0”), and assigning a “1” to all of the pixels in corresponding neighborhood locations having a probability of containing grass above the threshold. Pixel groupings with an area of 1 in the binary image are removed. Accordingly, areas assigned a “1” represent grass containing areas, and areas assigned a “0” represent obstacle containing areas where no grass is present.
[0070] In one embodiment, a scale of [0,1] and an empirically determined probability threshold of 0.50 were used. In another embodiment, a scale of [0,1] and an empirically determined probability threshold of 0.75 were used. However, it is contemplated that a person having ordinary skill in the art can choose to use a different threshold value and/or scale.
[0071] Grass identification using combined statistics was tested on 40 randomly selected samples of 31×31 neighborhoods containing the same objects that were used for individual texture measurements. The combined statistical measurements relied on the six neighborhood statistics that are contained in Table 4 below. The weighting coefficients were normalized for the six neighborhood statistics.
[0072]
| TABLE 4 | ||||||
|---|---|---|---|---|---|---|
| Normalized pre-weight coefficients for six identifiers. | ||||||
| Normal- | ||||||
| ized | ||||||
| Coeff. | Hue | Var(R)hTx | Var(R)vTx | ĀhTx | ĀvTx | ( x, y) |
| aNorm | 0.1458 | 0.1771 | 0.1677 | 0.1708 | 0.1708 | 0.1677 |
[0073] In this test, correct neighborhood identification was deemed to occur if the neighborhood was correctly identified as grass or obstacle containing with a voting agreement greater than 0.5 (or 50%). Since each of the utilized neighborhood statistics had approximately the same individual accuracy, this corresponds to the majority of the neighborhoods identifying the neighborhood as either grass containing or obstacle containing. Neighborhood identification rates were also calculated for a voting agreement greater than 0.75 (or 75%). The number of correctly identified neighborhoods based on this criterion of the combined statistics for the neighborhoods are tabulated below is Table 5.
[0074]
| TABLE 5 | |||
|---|---|---|---|
| Number of correctly identified 31 × 31 neighborhoods for combined | |||
| measurement (n = 40 for each measurement). | |||
| Combined (50% | Combined (75% | ||
| Measurement | Voting Threshold) | Voting Threshold) | |
| Grass | 40 | 40 | |
| Shadow | 38 | 33 | |
| Artificial | 32 | 38 | |
| Obstacles | |||
| Flowers | 6 | 15 | |
[0075] For 0.5 (50%) voting agreement or greater, the ability to correctly identify shaded grass increases compared to nearly all statistics to 38/40 (95%) when combining multiple neighborhood statistics. Illuminated grass is still recognized with 40/40 (100%) accuracy. Artificial obstacle recognition performance also increases compared to stand-alone hue and binary vertical variance measurements. For 0.75 (75%) voting agreement or greater, artificial obstacle recognition performance improves significantly to 38/40 (95%) and illuminated grass is still recognized with 40/40 (100%) accuracy. In one embodiment, MATLAB, produced by The MathWorks, Inc., was used to carry out the processes in blocks 300-390. However, it is contemplated that a person having ordinary skill in the art can use other software or hardware to carry out the operations.
[0076] In block 395, the binary grass image map generated in block 390 is abstracted to a range image. In one embodiment, the range image conversion occurs in the following way. First, the binary grass map, which is an array, is converted from rectangular real-world coordinates to a polar representation using the “Rectangular-to-Polar Array” function in LabVIEW, produced by National Instruments. This creates a binary range image of the mowable terrain and obstacle locations in the current camera frame, where the (0,0) location of the camera is located in the middle element of the bottom row of the array. Corresponding (r,θ) pairs are then sorted using LabVIEW's “Sort 1D Array” method, and the shortest range for each integer value of theta in the camera's field of view is inserted into a 1-D Array shifted so that the ranges are relative to the (0,0) location of the robot body. Ranges are inserted into the array sequentially, which creates a data structure identical to what is output by the LIDAR, where the first element corresponds to the range to the nearest obstacle at θ=1°, the second element corresponds to the range to the nearest obstacle at θ=2° degrees, etc. If an obstacle is not observed at a θ value, a range of 4 meters is entered into the corresponding 1-D array element, a distance equivalent to the maximum range of the SICK LMS291, indicating that no obstacle is present at that angle. This 1-D range array is referred to as a “pseudo-LIDAR” scan.
[0077] Optionally, to visualize the pseudo-LIDAR scan, a polar function is fit to the 1-D array and plotted in LabVIEW's front-panel. This allows the user to monitor what obstacles the mower was or was not able to see.
[0078] This range image processing occurs at a minimum of 10 Hz during mower operation. The resulting 1-D array is passed to the freespace observer. The freespace observer accomplishes two things. First, it combines the new pseudo-LIDAR observations with previously observed freespace to create an accurate estimate of traversable terrain around the robot. Second, it shifts observed ranges in front of the mower to behind the mower in accordance with how the mower moves through the environment. This shifted freespace is simply the previous freespace estimation, shifted according to the change in position and orientation of the mower's center. This creates a 360° representation of the currently traversable space around the mower, which is passed to the navigation system.
[0079] In one embodiment, the navigation systems is an onboard National Instruments compactRIO programmable automation controller; however, it is contemplated that a person having ordinary skill in the art can choose to use another navigation system.
[0080] Turning now to another embodiment of the invention as shown in FIG. 6 , in block 600, an image is captured of the terrain in front of the lawnmower. In block 605, the image is segmented into neighborhoods. In block 610, image statistic values are calculated for each neighborhood. In one embodiment, at least two image statistics are calculated for each neighborhood. In block 615, the image statistic values calculated in block 610 are compared to predetermined neighborhood grass statistics. In block 620, a binary representation of each image statistic is generated, wherein each neighborhood statistic value that compares favorably in block 615 is marked as grass (assigned a “1”), and each neighborhood statistic value that compares unfavorably in block 615 is marked as not containing grass (assigned a “0”). This is illustrated in FIG. 4 .
[0081] In block 625, each binary representation is weighted using predetermined normalized weighting coefficients. In block 630, the corresponding neighborhoods for all of the image statistics are summed. In block 635, a binary threshold is applied to the image statistics summed in block 630, thereby creating a grass/no grass binary neighborhood map. In block 640, the binary map generated in block 635 is abstracted into a range image.
[0082] Turning to FIG. 7 , which discloses a vision processing system 700 for an autonomous lawnmower having a mounted camera 705 and a processing arrangement 710. The processing arrangement 710, such the Mac Mini discussed above, connected to camera 705 join together as a vision system for the lawnmower. The processing arrangement provides structure for segmenting each image captured by cameras 705 into neighborhoods. The processing arrangement further provides structure for calculating at least two image statistics for each neighborhood. Additionally, the processing arrangement provides structure for generating a binary representation of each image statistic, wherein a neighborhood for an image statistic is marked as grass containing when the neighborhood image statistic compares favorably to a predetermined image statistics for grass. The processing arrangement further provides structure of applying a weight to each of the binary representations of each image statistic. The processing arrangement further providing structure of summing corresponding neighborhoods for all image statistics. Further, the processing arrangement further providing structure for applying a binary threshold to each of said summed neighborhoods, wherein a binary map is produced representing grass containing areas and non-grass containing areas. Further, the binary map representing grass containing areas and non-grass containing areas can be abstracted to a range image.
[0083] While this invention has been described in conjunction with the specific embodiments described above, it is evident that many alternatives, combinations, modifications and variations are apparent to those skilled in the art. Accordingly, the preferred embodiments of this invention, as set forth above, are intended to be illustrative only, and not in a limiting sense. Various changes can be made without departing from the spirit and scope of this invention. Therefore, the technical scope of the present invention encompasses not only those embodiments described above, but also all that fall within the scope of the appended claims.
[0084] This written description uses examples to disclose the invention, including the best mode, and also to enable any person skilled in the art to practice the invention, including making and using any devices or systems and performing any incorporated processes. The patentable scope of the invention is defined by the claims, and may include other examples that occur to those skilled in the art. These other examples are intended to be within the scope of the claims if they have structural elements that do not differ from the literal language of the claims, or if they include equivalent structural elements with insubstantial differences from the literal language of the claims.
[0085] Of course, it is to be understood that not necessarily all such objects or advantages described above may be achieved in accordance with any particular embodiment. Thus, for example, those skilled in the art will recognize that the systems and techniques described herein may be embodied or carried out in a manner that achieves or optimizes one advantage or group of advantages as taught herein without necessarily achieving other objects or advantages as may be taught or suggested herein.
[0086] Furthermore, the skilled artisan will recognize the interchangeability of various features from different embodiments. The various features described, as well as other known equivalents for each feature, can be mixed and matched by one of ordinary skill in this art to construct additional systems and techniques in accordance with principles of this disclosure.
Claims
What is claimed is:
1. A method for identifying lawn grass comprising:
capturing an image of terrain in front of a mower;
segmenting said image into neighborhoods;
calculating values for at least two image statistics for each of said neighborhoods;
comparing said image statistic values of said neighborhoods to predetermined image statistic values for grass;
generating a binary representation of each image statistic value, wherein a neighborhood for an image statistic value is marked as grass containing when said image statistic value of said neighborhood compares favorably to a predetermined image statistic value for grass;
weighting each of said binary representations of each image statistic;
summing corresponding neighborhoods for all image statistics; and
applying a binary threshold to each of said summed neighborhoods, wherein said binary threshold produces a binary map representing grass containing areas and non-grass containing areas.
2. The method of claim 1, wherein one of said image statistics is mean hue.
3. The method of claim 1, wherein said image statistic value of said neighborhood compares favorably to a predetermined image statistic for grass when said image statistic value of said neighborhood is within a predetermined number of standard deviations of said predetermined image statistic value for grass.
4. The method of claim 1, wherein at least one of said image statistics is calculated using an intensity plane extracted from said captured image.
5. The method of claim 4, wherein said image statistics include at least one of: mean grayscale horizontal intensity, mean grayscale vertical intensity, vertical binary variance, mean vertical binary edge response area, horizontal and vertical neighborhood centroid location, horizontal binary variance, or mean horizontal binary edge response area.
6. The method of claim 1, wherein said binary representations of each image statistic are weighted by multiplying said neighborhoods in said binary representations of each image statistic by predetermined normalizing weighting coefficients.
7. The method of claim 1, wherein said binary map representing grass containing areas and non-grass containing areas is abstracted to a range image.
8. A vision processing system for an autonomous lawnmower comprising:
an image capturing camera mounted on said lawnmower for capturing an image of terrain in front of said lawnmower;
a processing arrangement connected to said camera for segmenting said image into neighborhoods;
said processing arrangement further providing structure for calculating values for at least two image statistics for each neighborhood;
said processing arrangement further providing structure for generating a binary representation of each image statistic value, wherein a neighborhood for an image statistic value is marked as grass containing when said image statistic value of said neighborhood compares favorably to a predetermined image statistic for grass;
said processing arrangement further providing structure for applying a weight to each of said binary representations of each image statistic;
said processing arrangement further providing structure for summing corresponding neighborhoods for all image statistics; and
said processing arrangement further providing structure for applying a binary threshold to each of said summed neighborhoods, wherein a binary map is produced representing grass containing areas and non-grass containing areas.
9. The vision system of claim 8, wherein said image statistic of said neighborhood compares favorably to a predetermined image statistic for grass when said image statistic of said neighborhood is within a predetermined number of standard deviations of the predetermined image statistic for grass.
10. The vision system of claim 8, wherein one of said image statistics is mean hue.
11. The vision system of claim 8, wherein at least one of said image statistics is calculated using an intensity plane extracted from said captured image.
12. The vision system of claim 11, wherein said image statistics include at least one of: mean grayscale horizontal intensity, mean grayscale vertical intensity, vertical binary variance, mean vertical binary edge response area, horizontal and vertical neighborhood centroid location, horizontal binary variance, or mean horizontal binary edge response area.
13. The vision system of claim 8, wherein said binary representations of each image statistic are weighted by multiplying said neighborhoods in said binary representations of each image statistic by predetermined normalizing weighting coefficients.
14. The vision system of claim 8, wherein said binary map representing grass containing areas and non-grass containing areas is abstracted to a range image.
15. A method for identifying lawn grass by an autonomous lawnmower comprising:
capturing an image of terrain in front of said mower;
segmenting said image into neighborhoods;
calculating values for at least two image statistics for each of said neighborhoods;
generating a binary representation of each image statistic, wherein said binary representation is generated by comparing said calculated image statistic values to predetermined image statistic values for grass;
weighting each of said binary representations of each image statistic;
summing corresponding neighborhoods for all image statistics; and
applying a binary threshold to each of said summed neighborhoods to generate a binary map representing grass containing areas and non-grass containing areas.
16. The method of claim 15, wherein one of said image statistics is mean hue.
17. The method of claim 15, wherein at least one of said image statistics is calculated using an intensity plane extracted from said captured image.
18. The method of claim 17, wherein said image statistics include at least one of: mean grayscale horizontal intensity, mean grayscale vertical intensity, vertical binary variance, mean vertical binary edge response area, horizontal and vertical neighborhood centroid location, horizontal binary variance, or mean horizontal binary edge response area.
19. The method of claim 15, wherein said binary representations of each image statistic are weighted before summing said corresponding neighborhoods of all image statistics.
20. The method of claim 15, wherein said image statistic of said neighborhood compares favorably to a predetermined image statistic for grass when said image statistic of said neighborhood is within a predetermined number of standard deviations of the predetermined image statistic for grass.
Patent Citations (9)
| Patent | Date | Inventor | Cited By |
|---|---|---|---|
| US5134661(A) | 1992-07-01 | Reinsch | Applicant |
| US5911669(A) | 1999-06-01 | Stentz et al. | Examiner |
| US7400957(B2) | 2008-07-01 | Hofer et al. | Applicant |
| US2008/0219508(A1) | 2008-09-01 | Ganguli et al. | Applicant |
| US2009/0319170(A1) | 2009-12-01 | Madsen et al. | Applicant |
| US2010/0053593(A1) | 2010-03-01 | Bedros et al. | Applicant |
| US2011/0157221(A1) | 2011-06-01 | Ptucha et al. | Examiner |
| US2011/0157226(A1) | 2011-06-01 | Ptucha et al. | Examiner |
| US2012/0158236(A1) | 2012-06-01 | Chung et al. | Examiner |
Non-Patent Literature (92)
- Bacha et al., “Odin: Team VictorTango's Entry in the DARPA Urban Challenge” “Special Issue on the 2007 DARPA Urban Challenge, Part 1” Journal of Field Robotics, 25:467-492 (Jan. 2008).Applicant
- Miller et al., “Team Cornell's Skynet: Robust Perception and Planning in an Urban Environment” “Special Issue on the 2007 DARPA Urban Challenge, Part 1” Journal of Field Robotics, 25:493-527 (Jan. 2008).Applicant
- Patz et al., “A Practical Approach to Robotic Design for the DARPA Urban Challenge” “Special Issue on the 2007 DARPA Urban Challenge, Part 1” Journal of Field Robotics, 25:528-566 (Jan. 2008).Applicant
- Forsyth, D.A. and Ponce, J., “Computer Vision: A Modern Approach” Prentice Hall (Jan. 2003).Applicant
- Gonzalez, R.C. and Woods, R.E., “Digital Image Processing” Third Edition, Prentice Hall (Jan. 2008).Applicant
- International Search Report and Written Opinion mailed Sep. 24, 2010 for PCT/US2010/001881 filed Jul. 1, 2010.Applicant
- International Preliminary Report on Patentability mailed Jan. 4, 2012 for PCT/US2010/001881 filed Jul. 1, 2010.Applicant
- Smith, “Color Gamut Transform Pairs”, Computer Graphics Lab, New York Institute of Technology, Old Westbury, NY, pp. 12-19, 1978.Applicant
- Marr et al., “Theory of Edge Detection”, Proceedings of the Royal Society of London, Series B, vol. 207, No. 1167, pp. 187-217, Feb. 29, 1980.Applicant
- Chateau et al., “An Original Correlation and Data Fusion Based Approach to Detect a Reap Limit into a Gray Level Image”, IEEE, pp. 1258-1263, 1997.Applicant
- Bertozzi et al., “GOLD: A Parallel Real-Time Stereo Vision System for Generic Obstacle and Lane Detection”, IEEE Transactions on Image Processing, vol. 7, No. 1, pp. 62-80, Jan. 1998.Applicant
- Hicks, II, et al., “A Survey of Robot Lawn Mowers”, 8 pages, 2000.Applicant
- Ulrich et al., “Appearance-Based Obstacle Detection with Monocular Color Vision”, Proceedings of the AAAI National Conference on Artificial Intelligence, Austin, Texas, Jul./Aug. 2000, 6 pages.Applicant
- Bellutta et al., “Terrain Perception for DEMO III”, Proceedings of the 2000 Intelligent Vehicles Conference, Dearborn, Michigan, Oct. 4-5, 2000, 6 pages.Applicant
- Hemming et al., “Computer-Vision-Based Weed Identification Under Field Conditions Using Controlled Lighting”, J. Agric. Engng. Res. vol. 78, No. 3, pp. 233-243, 2001.Applicant
- Batavia et al., “Obstacle Detection Using Adaptive Color Segmentation and Color Stereo Homography”, 6 pages, 2001.Applicant
- Desouza et al., “Vision for Mobile Robot Navigation: A Survey”, IEEE Transactions on Pattern Analysis and Machine Intelligence, vol. 24, No. 2, pp. 237-267, Feb. 2002.Applicant
- Philipp et al., “Improving Plant Discrimination in Image Processing by Use of Different Colour Space Transformations”, Computers and Electronics in Agriculture, vol. 35, pp. 1-15, 2002.Applicant
- Salvador et al., “Spatio-Temporal Shadow Segmentation and Tracking”, Signal Processing Institute (ITS), Swiss Federal Institute of Technology (EPFL) Lausanne, Switzerland, 12 pages, 2003.Applicant
- Karcher et al., “Quantifying Turfgrass Color Using Digital Image Analysis”, Turfgrass Science, Crop Sci., vol. 43, pp. 943-951, 2003.Applicant
- Karklin et al., “Learning Higher-Order Structures in Natural Images”, Institute of Physics Publishing, Network: Computation in Neural Systems, vol. 14, pp. 483-499, 2003.Applicant
- Dima et al., “Classifier Fusion for Outdoor Obstacle Detection”, Carnegie Mellon University, The Robotics Institute, Apr. 26, 2004, 7 pages.Applicant
- Taylor et al., “Molecular Vision as a Range Sensor”, CIMCA 2004 Proceedings pp. 566-575, Jul. 12-14, 2004, Gold Coast, Australia.Applicant
- Montemerlo et al., “Winning the DARPA Grand Challenge with an AI Robot”, pp. 982-987, 2006.Applicant
- Sheh et al., “A Low-Cost, Compact, Lightweight 3D Range Sensor”, 8 pages, 2006.Applicant
- Watchareeruetai et al., “Computer Vision Based Methods for Detecting Weeds in Lawns”, Machine Vision and Applications, vol. 17, pp. 287-296, 2006.Applicant
- Aziz et al., “Evaluation of Visual Attention Models for Robots”, Proceedings of the Fourth IEEE International Conference on Computer Vision Systems (ICVS 2006), 6 pages, 2006.Applicant
- Kim et al., “Traversability Classification for UGV Navigation: A Comparison of Patch and Superpixel Representations”, Oct. 29, 2007.Applicant
- Ali, “Tree Detection Using Color, and Texture Cues for Autonomous Navigation in Forest Environment”, Master's Thesis in Computing Science, Umea University, Department of Computing Science Sweden, Jun. 2006.Applicant
- Zafarifar et al., “Grass Field Detection for TV Picture Quality Enhancement”, 2 pages, Aug. 5, 2008.Applicant
- McMichael, “Lane Detection for DEXTER, and Autonomous Robot, in the Urban Challenge”, Master's Thesis, Submitted for the Degree of Master of Science, Case Western Reserve University, May 2008.Applicant
- Beno et al., “2008 Cutter Technical Report”, Case Western Reserve University, 30 pages, May 7, 2008.Applicant
- Beno et al., CWRU Cutter 2.0 Technical Report, Case Western Reserve University, 31 pages, May 14, 2009.Applicant
- “The Sixth! Annual Robotic Lawnmower Competition Rulebook”, Revision 1.2.2009, May 28-31, 2009, Dayton, Ohio, 21 pages.Applicant
- Bohren et al., “Little Ben: The Ben Franklin Racing Team's Entry in the 2007 DARPA Urban Challenge”, M. Bueher et al., (Eds.): The DARPA Urban Challenge, STAR 56, pp. 231-255, 2009.Applicant
- Schepelmann et al., “Vision-Based Obstacle Detection and Avoidance for the CWRU Cutter Autonomous Lawnmower”, 6 pages, Nov. 9, 2009.Applicant
- Wang, “Autonomous Machine Vision for Off-Road Vehicles in Unstructured Fields,” Urbana, Illinois, Dissertation, Submitted for the Degree of Doctor of Philosophy in Agricultural and Biological Engineering, University of Illinois at Urbana-Champaign, 2009.Applicant
- Hughes et al., “C3-CWRU Cutter C, Technical Report”, Case Western Reserve University, 7th Autonomous Lawnmower Competition Group, 31 pages, May 14, 2010.Applicant
- Wybo et al., “Obstacles Avoidance by Monocular Multi-Cue Image Analysis”, 12 pages, Nov. 15, 2008, 15th World Congress on ITS.Applicant
- Daltorio et al., “An Obstacle-Edging Reflex for an Autonomous Lawnmower”, Case Western Reserve University, 14 pages, May 2010.Applicant
- Borenstein et al., “The Vector Field Histogram-Fast Obstacle Avoidance for Mobile Robots” IEEE Journal of Robots and Automation, 7:278-288 (Jun. 1991).Applicant
- Broggi et al., “The Single Frame Stereo Vision System for Reliable Obstacle Detection Used During the 2005 DARPA Grand Challenge on TerraMax” Proceedings of the IEEE ITSC 2006 IEEE Intelligent Transportation Systems Conference, 745-752 (Sep. 2006).Applicant
- Camhi et al., “High-Frequency Steering Maneuvers Mediated by Tactile Cues: Antennal Wall-Following in the Cockroach” J. Exp. Bio., 202:631-643 (Feb. 1999).Applicant
- Chen et al., “A Crash Avoidance System Based Upon the Cockroach Escape Response Circuit” Proceedings of the 1997 IEEE International Conference on Robotics and Automation, 2007-2012 (Apr. 1997).Applicant
- Frazzoli, “Real-Time Motion Planning for Agile Autonomous Vehicles” J. Guidance, Control, and Dynamics, 25:116-129 (Jan. 2002).Applicant
- Kalman, “A New Approach to Linear Filtering and Prediction Problems” Transactions of the ASME—Journal of Basic Engineering 82(Series D):35-45 (Jan. 1960).Applicant
- Kelly, “An Intelligent Predictive Control Approach to the High-Speed Cross-Country Autonomous Navigation Problem” PHD Thesis, Carnegie Mellon (Sep. 1995).Applicant
- Lumelski et al., “Dynamic Path Planning for a Automaton with Limited Information on the Environment” IEEE Transactions on Automatic Control, AC-31:1058-1063 (Nov. 1986).Applicant
- Maxwell, “On the Theory of Compound Colours, and the Relations of the Colours of the Spectrum” Phil. Trans. R. Soc., 150:57-84 (Mar. 1860).Applicant
- Minguez et al., “Abstracting Vehicle Shape and Kinematic Constraints from Obstacle Avoidance Methods” Autonomous Robots, 20:43-59 (Feb. 2006).Applicant
- Newman, “Team CASE and the 2007 DARPA Urban Challenge” DARPA Grand Challenge Tech Report, http/www.darpa.mil/Jun. 2007.Applicant
- Newman, “Roberto” IGVC Design Report retrieved from http:/igvc.org/design/reports/dr148.pdf (Jan. 2006).Applicant
- Nourani-Vatani et al., “Practical Path Planning and Obstacle Avoidance for Autonomous Mowing” Australasian Conference on Robotics and Automation (Jan. 2006).Applicant
- Schepelmann et al., “Visual Segmentation of Lawn Grass for a Mobile Robotic Lawnmower” IEEE/RSJ International Conference on Intelligent Robots and Systems, 18-22 (Submitted) (Mar. 2010).Applicant
- Smith et al., “Design and Implementation of a Control Algorithm for an Autonomous Lawnmower” Midwest Symposium on Circuits and Systems (Jan. 2005).Applicant
- Thrun et al., “Stanley: The Robot that Won the DARPA Grand Challenge” Journal of Field Robotics, 23:661-692 (Jun. 2006).Applicant
- Wikman et al., “Reflexive Collision Avoidance: A Generalized Approach” IEEE ICRA, pp. 31-36 (Jan. 1993).Applicant
- Yata et al., “Wall Following Using Angle Information Measured by a Single Ultrasonic Transducer” Proceedings of the IEEE International Conference on Robotics & Automation, pp. 1590-1596 (May 1998).Applicant
- Ritzmann, R.E. “The Neural Organization of Cockroach Escape and Its Role in Context Dependent Orientation”, Chapter VI in “Biological Neural Networks in Invertebrate Neuroethology and Robotics” Ed: Beer et al., Academic Press, pp. 113-137 (Jan. 1993).Applicant
- Chen, G., “Texture Based Road Surface Detection” M.S. Thesis, Case Western Reserve University (Aug. 2008).Applicant
- Bemporad et al., “Sonar-Based Wall-Following Control of Mobile Robots” submitted to Journal of Dynamic, Measurement and Control (Jan. 1997).Applicant
- Garcia-Alegre et al., “Autonomous Robot in Agricultural Tasks” European Conference on Precision Agriculture, France (Jan. 2001).Applicant
- Kanayama et al., “A New Line Tracking Method for Nonholonomic Vehicles” Proceedings of the IEEE International Conference on Robots and Automation, New Mexico (Apr. 1997).Applicant
- Oh et al., “Complete Coverage Navigation of Cleaning Robots Using Triangular-Cell-Based Map” IEEE Transactions on Industrial Electronics, 51:718-726 (Jun. 2004).Applicant
- Prassler et al., “A Short History of Cleaning Robots” Autonomous Robots, 9:211-226 (Jan. 2000).Applicant
- Schafer et al., “3D Obstacle Detection and Avoidance in Vegetated Off-Road Terrain” IEEE International Conference on Robotics and Automation, California, pp. 923-928 (May 2008).Applicant
- Tsugawa, S., “A History of Automated High Systems in Japan and Future Issues” Proceedings of the 2008 IEEE International Conference on Vehicular Electronics and Safety, Ohio (Sep. 2008).Applicant
- Ulrich et al., “Autonomous Vacuum Cleaner” Robotics and Autonomous Systems, 19:233-245 (Jan. 1997).Applicant
- Yokoi et al., “Development of a Whisker Sensor System and Simulation of Active Whisking for Agent Navigation” Proceedings of 2004 IEEE/RSJ International Conference on Intelligent Robots and Systems, Japan (Sep. 2004).Applicant
- Kyodo America, “LB3510 Manual”, http:/www.lawnbott.com/pdf/LB3510-manual.pdf (Jan. 2009).Applicant
- Robomower Review, http://bamabots.com/rl1000review.htm (Nov. 2005).Applicant
- Kyodo America, “How LawnBott Works”, http://lawnbott.com/buyers/how-it-works/ (Jan. 2012).Applicant
- Mahalanobis, P.C., “On the Generalized Distance in Statistics” Proceedings of the National Institute of Sciences in India (Apr. 1936).Applicant
- Grewal et al., “Global Positioning Systems, Inertial Navigation, and Integration” Hoboken, NJ: Wiley, p. 394 (Jan. 2008).Applicant
- Liu et al., “Hue Image RGB Colour Composition. a Simple technique to suppress shadow and enhance special signature” International Journal of Remote Sensing, 11:1521-1530 (Jan. 1990).Applicant
- Carelli et al., “Corridor navigation and wall-following stable control for sonar-based mobile robots” Robotics and Autonomous Systems, 45:235-247 (Jan. 2003).Applicant
- Lamperski et al., “Dynamical Wall Following for a Wheeled Robot Using a Passive Tactile Sensor” Proceedings of the IEEE International Conference on Robotics and Automation, Barcelona, Spain, pp. 3838-3843 (Apr. 2005).Applicant
- Brooks, R., “A robust layered control system for a mobile robot” IEEE J. Robotics Autom., vol. RA-2:14-23 (Mar. 1986).Applicant
- Noborio et al., “Sensor-based path-planning algorithms for a nonholonomic mobile robot” International Conference on Intelligent Robots and Systems, pp. 917-924 (Jan. 2000).Applicant
- Kanayama et al., “A locomotion control method for autonomous vehicles” Proceedings of the IEEE International Conference on Robotics and Automation, Philadelphia, PA, pp. 1315-1317 (Jan. 1988).Applicant
- Lavalle et al., “Rapidly-exploring random trees: Progress and Prospects” Algorithmic and Computational Robotics, Eds. Donald, Lynch, and Russ; AK Peters Wellesley, MA, pp. 1-17 (Jan. 2001).Applicant
- Kavraki et al., “Probablistic roadmaps for path planning in high-dimensional configuration spaces” IEEE Trans. Robot. & Auton., 12:566-580 (Aug. 1996).Applicant
- Simmons, R., “The Curvature-Velocity Method for Local Obstacle Avoidance” International Conference on Robotics and Automation, pp. 3375-3382 (Apr. 1996).Applicant
- Jung et al., “Whisker-based mobile robot navigation” Proceedings of the IEEE/RSJ International Conference on Intelligence Robots and Systems (IROS), 2:497-504 (Jan. 1996).Applicant
- Sekhavat et al., “Multilevel path planning for nonholonomic robots using semiholonomic subsystems” Int. J. Robot. Res., 17:840-857 (Jan. 1998).Applicant
- Khatib et al., “Real-Time Obstacle Avoidance for Manipulators and Mobile Robots” In Proceedings of the IEEE International Conference on Robotics and Automation, St. Louis, MO, pp. 500-505 (Jan. 1985).Applicant
- Stachniss et al., “An Integrated Approach to Goal-Directed Obstacle Avoidance under Dynamic Constraints for Dynamic Environments” IROS (Oct. 2002).Applicant
- Chatterjee et al., “Use of single side reflex for autonomous navigation of mobile robots in unknown environments” Robotics and Autonomous Systems, 35:77-96 (Jan. 2001).Applicant
- Buehler et al., “Editorial” “Special Issue on the 2007 DARPA Urban Challenge, Part 1” Journal of Field Robotics, 25:423-424 (Jan. 2008).Applicant
- Urmson et al., “Autonomous Driving in Urban Environments: Boss and the Urban Challenge” “Special Issue on the 2007 DARPA Urban Challenge, Part 1” Journal of Field Robotics, 25:425-466 (Jan. 2008).Applicant
- Dima et al., “Sensor and Classifier Fusion for Outdoor Obstacle Detection: an Application of Data Fusion to Autonomous Off-Road Navigation” Proceedings of the 32nd Applied Imagery Pattern Recognition Workshop (2003).Applicant
- Communication from European Patent Office, dated Nov. 13, 2012, for corresponding European Patent Application No. 10736866.4.Applicant