Sobel

From Wikipedia, the free encyclopedia

The Sobel operator is an operator used in image processing, particularly within edge detection algorithms. Technically, it is a discrete differentiation operator, computing an approximation of the gradient of the image intensity function. At each point in the image, the result of the Sobel operator is either the corresponding gradient vector or the norm of this vector.

Contents

In simple terms, the operator calculates the gradient of the image intensity at each point, giving the direction of the largest possible increase from light to dark and the rate of change in that direction. The result therefore shows how "abruptly" or "smoothly" the image changes at that point, and therefore how likely it is that that part of the image represents an edge, as well as how that edge is likely to be oriented. In practice, the magnitude (likelihood of an edge) calculation is more reliable and easier to interpret than the direction calculation.

Mathematically, the gradient of a two-variable function (here the image intensity function) is at each image point a 2D vector with the components given by the derivatives in the horizontal and vertical directions. At each image point, the gradient vector points in the direction of largest possible intensity increase, and the length of the gradient vector corresponds to the rate of change in that direction. This implies that the result of the Sobel operator at an image point which is in a region of constant image intensity is a zero vector and at a point on an edge is a vector which points across the edge, from darker to brighter values.

Mathematically, the operator uses two 3×3 kernels which are convolved with the original image to calculate approximations of the derivatives - one for horizontal changes, and one for vertical. If we define \mathbf{A} as the source image, and \mathbf{G_x} and \mathbf{G_y} are two images which at each point contain the horizontal and vertical derivative approximations, the latter are computed as:

\mathbf{G_x} = \begin{bmatrix}  -1 & 0 & +1 \\ -2 & 0 & +2 \\ -1 & 0 & +1  \end{bmatrix} * \mathbf{A} \quad \mbox{and} \quad  \mathbf{G_y} = \begin{bmatrix}  +1 & +2 & +1 \\ 0 & 0 & 0 \\ -1 & -2 & -1  \end{bmatrix} * \mathbf{A}

At each point in the image, the resulting gradient approximations can be combined to give the gradient magnitude, using:

\mathbf{G} = \sqrt{ \mathbf{G_x}^2 + \mathbf{G_y}^2 }

Using this information, we can also calculate the gradient's direction:

\mathbf{\Theta} = \operatorname{arctan}\left({ \mathbf{G_y} \over \mathbf{G_x} }\right)

where, for example, Θ is 0 for a vertical edge which is darker on the left side.

Since the intensity function of a digital image is only known at discrete points, derivatives of this function cannot be defined unless we assume that there is an underlying continuous intensity function which has been sampled at the image points. With some additional assumptions, the derivative of the continuous intensity function can be computed as a function on the sampled intensity function, i.e. the digital image. It turns out that the derivatives at any particular point are functions of the intensity values at virtually all image points. However, approximations of these derivative functions can be defined at lesser or larger degrees of accuracy.

The Sobel operator represents a rather inaccurate approximation of the image gradient, but is still of sufficient quality to be of practical use in many applications. More precisely, it uses intensity values only in a 3×3 region around each image point to approximate the corresponding image gradient, and it uses only integer values for the coefficients which weight the image intensities to produce the gradient approximation.

As a consequence of its definition, the Sobel operator can be implemented by simple means in both hardware and software: only eight image points around a point are needed to compute the corresponding result and only integer arithmetic is needed to compute the gradient vector approximation. Furthermore, the two discrete filters described above are both separable:

\begin{bmatrix}  -1 & 0 & +1 \\ -2 & 0 & +2 \\ -1 & 0 & +1  \end{bmatrix} = \begin{bmatrix}  1 \\ 2 \\ 1   \end{bmatrix} * \begin{bmatrix}  -1 & 0 & +1 \end{bmatrix} \quad \quad \begin{bmatrix}  +1 & +2 & +1 \\ 0 & 0 & 0 \\ -1 & -2 & -1  \end{bmatrix} = \begin{bmatrix}  +1 \\ 0 \\ -1   \end{bmatrix} * \begin{bmatrix}  1 & 2 & 1 \end{bmatrix}

and the two derivatives \mathbf{G_x} and \mathbf{G_y} can therefore be computed as

\mathbf{G_x} = \begin{bmatrix}  1 \\ 2 \\ 1 \end{bmatrix} * \begin{bmatrix}  -1 & 0 & +1   \end{bmatrix} * \mathbf{A} \quad \mbox{and} \quad  \mathbf{G_y} = \begin{bmatrix}  +1 \\ 0 \\ -1   \end{bmatrix} * \begin{bmatrix}  1 & 2 & 1 \end{bmatrix} * \mathbf{A}

In certain implementations, this separable computation may be advantageous since it implies fewer arithmetic computations for each image point.

Because the result of the Sobel operator is a 2-dimensional map of the gradient at each point, it can be processed and viewed as though it is itself an image, with the areas of high gradient (the likely edges) visible as white lines. The following images illustrate this, by showing the computation of the Sobel operator on a simple image..

Grayscale image of plastic figure of a penguin
Grayscale image of plastic figure of a penguin
Normalised sobel gradient image of penguin figure
Normalised sobel gradient image of penguin figure
Normalised sobel x-gradient image
Normalised sobel x-gradient image
Normalised sobel y-gradient image
Normalised sobel y-gradient image

Unpublished:

   Sobel, I., Feldman,G., "A 3x3 Isotropic Gradient Operator for Image
   Processing", presented at a talk at the Stanford Artificial Project in
   1968, unpublished but often cited, orig. in Pattern Classification and
   Scene Analysis, Duda,R. and Hart,P., John Wiley and Sons,'73, pp271-2

Advanced Search
Included Web Search Engines


Safe Search

close

Top Matching Results

Occasionally Search.com will highlight specialized results that are based on the context of your query. Examples of specialized results include specific links to news, images, or video.

Top Matching Results may highlight information from other Search.com pages, content from the CNET Network of sites, or third party content. The listings are based purely on relevance. Search.com does not receive payment for listings in this section but our partners that provide this data may get paid for listing these products.

Sponsored Links

This section contains paid listings which have been purchased by companies that want to have their sites appear for specific search terms and related content. These listings are administered, sorted and maintained by a third party and are not endorsed by Search.com.

Search Results

Search.com sends your search query to several search engines at one time and integrates the results into one list which has been sorted by relevance using Search.com's proprietary algorithm. You can customize the list of search engines included in your metasearch from the preferences.

The search engines that are used in your metasearch may allow companies to pay to have their Web sites included within the results. To view the Paid Inclusion policy for a specific search engine, please visit their Web site. Search.com does not accept payment or share revenue with any search engine partner for listings in this section.