Prewitt operator · Lablypedia

In a nutshell

The Prewitt operator is 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 Prewitt operator is either the corresponding gradient vector or the norm of this vector. The Prewitt operator is based on convolving the image with a small, separable, and integer valued filter in horizontal and vertical directions and is therefore relatively inexpensive in terms of computations like Sobel and Kayyali operators. On the other hand, the gradient approximation which it produces is relatively crude, in particular for high frequency variations in the image. The Prewitt operator was developed by Judith M. S. Prewitt.

Related on Lably

Through Digital image processing, Recurrence relation and Gradient

Lab · Class 10–12Handwriting Recognizer Lab · Class 10–12Pixel Art Explorer Lab · Class 6–12AI & Python Builds (Medha AI)Lab · Class 7–11How the Monsoon Works

Simplified description

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 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 Prewitt 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.

3 more chapters are free with a Lably account.

Formulation

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 $A$ as the source image, and $G_{x}$ and $G_{y}$ are two images which at each point contain the horizontal and vertical derivative approximations, the latter are computed as:

G_{y} = [\begin{matrix} + 1 & + 1 & + 1 \\ 0 & 0 & 0 \\ - 1 & - 1 & - 1 \end{matrix}] * A and G_{x} = [\begin{matrix} + 1 & 0 & - 1 \\ + 1 & 0 & - 1 \\ + 1 & 0 & - 1 \end{matrix}] * A

where $*$ here denotes the 2-dimensional convolution operation.

Since the Prewitt kernels can be decomposed as the products of an averaging and a differentiation kernel, they compute the gradient with smoothing. Therefore, it is a separable filter. For example, $G_{x}$ can be written as

[\begin{matrix} + 1 & 0 & - 1 \\ + 1 & 0 & - 1 \\ + 1 & 0 & - 1 \end{matrix}] = [\begin{matrix} 1 \\ 1 \\ 1 \end{matrix}] [\begin{matrix} + 1 & 0 & - 1 \end{matrix}]

The x-coordinate is defined here as increasing in the "left"-direction, and the y-coordinate is defined as increasing in the "up"-direction. At each point in the image, the resulting gradient approximations can be combined to give the gradient magnitude, using:

G = \sqrt{{G_{x}}^{2} + {G_{y}}^{2}}

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

Θ = atan2 (G_{y}, G_{x})

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

Example

Code example

<span id="mwUg"></span><span class="c" id="mwUw">% MATLAB Code | Prewitt Operator from Scratch </span> <span class="c" id="mwVA">% Read Input Image </span> <span class="n" id="mwVQ">input_image</span><span class="w" id="mwVg"> </span><span class="p" id="mwVw">=</span><span class="w" id="mwWA"> </span><span class="nb" id="mwWQ">imread</span><span class="p" id="mwWg">(</span><span class="s" id="mwWw">'[name of input image file].[file format]'</span><span class="p" id="mwXA">);</span><span class="w" id="mwXQ"> </span> <span class="c" id="mwXg">% Displaying Input Image </span> <span class="n" id="mwXw">input_image</span><span class="w" id="mwYA"> </span><span class="p" id="mwYQ">=</span><span class="w" id="mwYg"> </span><span class="nb" id="mwYw">uint8</span><span class="p" id="mwZA">(</span><span class="n" id="mwZQ">input_image</span><span class="p" id="mwZg">);</span><span class="w" id="mwZw"> </span> <span class="nb" id="mwaA">figure</span><span class="p" id="mwaQ">,</span><span class="w" id="mwag"> </span><span class="nb" id="mwaw">imshow</span><span class="p" id="mwbA">(</span><span class="n" id="mwbQ">input_image</span><span class="p" id="mwbg">);</span><span class="w" id="mwbw"> </span><span class="nb" id="mwcA">title</span><span class="p" id="mwcQ">(</span><span class="s" id="mwcg">'Input Image'</span><span class="p" id="mwcw">);</span><span class="w" id="mwdA"> </span> <span class="c" id="mwdQ">% Convert the truecolor RGB image to the grayscale image </span> <span class="n" id="mwdg">input_image</span><span class="w" id="mwdw"> </span><span class="p" id="mweA">=</span><span class="w" id="mweQ"> </span><span class="nb" id="mweg">rgb2gray</span><span class="p" id="mwew">(</span><span class="n" id="mwfA">input_image</span><span class="p" id="mwfQ">);</span><span class="w" id="mwfg"> </span> <span class="c" id="mwfw">% Convert the image to double </span> <span class="n" id="mwgA">input_image</span><span class="w" id="mwgQ"> </span><span class="p" id="mwgg">=</span><span class="w" id="mwgw"> </span><span class="nb" id="mwhA">double</span><span class="p" id="mwhQ">(</span><span class="n" id="mwhg">input_image</span><span class="p" id="mwhw">);</span><span class="w" id="mwiA"> </span> <span class="c" id="mwiQ">% Pre-allocate the filtered_image matrix with zeros </span> <span class="n" id="mwig">filtered_image</span><span class="w" id="mwiw"> </span><span class="p" id="mwjA">=</span><span class="w" id="mwjQ"> </span><span class="nb" id="mwjg">zeros</span><span class="p" id="mwjw">(</span><span class="nb" id="mwkA">size</span><span class="p" id="mwkQ">(</span><span class="n" id="mwkg">input_image</span><span class="p" id="mwkw">));</span><span class="w" id="mwlA"> </span> <span class="c" id="mwlQ">% Prewitt Operator Mask </span> <span class="n" id="mwlg">Mx</span><span class="w" id="mwlw"> </span><span class="p" id="mwmA">=</span><span class="w" id="mwmQ"> </span><span class="p" id="mwmg">[</span><span class="o" id="mwmw">-</span><span class="mi" id="mwnA">1</span><span class="w" id="mwnQ"> </span><span class="mi" id="mwng">0</span><span class="w" id="mwnw"> </span><span class="mi" id="mwoA">1</span><span class="p" id="mwoQ">;</span><span class="w" id="mwog"> </span><span class="o" id="mwow">-</span><span class="mi" id="mwpA">1</span><span class="w" id="mwpQ"> </span><span class="mi" id="mwpg">0</span><span class="w" id="mwpw"> </span><span class="mi" id="mwqA">1</span><span class="p" id="mwqQ">;</span><span class="w" id="mwqg"> </span><span class="o" id="mwqw">-</span><span class="mi" id="mwrA">1</span><span class="w" id="mwrQ"> </span><span class="mi" id="mwrg">0</span><span class="w" id="mwrw"> </span><span class="mi" id="mwsA">1</span><span class="p" id="mwsQ">];</span><span class="w" id="mwsg"> </span> <span class="n" id="mwsw">My</span><span class="w" id="mwtA"> </span><span class="p" id="mwtQ">=</span><span class="w" id="mwtg"> </span><span class="p" id="mwtw">[</span><span class="o" id="mwuA">-</span><span class="mi" id="mwuQ">1</span><span class="w" id="mwug"> </span><span class="o" id="mwuw">-</span><span class="mi" id="mwvA">1</span><span class="w" id="mwvQ"> </span><span class="o" id="mwvg">-</span><span class="mi" id="mwvw">1</span><span class="p" id="mwwA">;</span><span class="w" id="mwwQ"> </span><span class="n" id="mwwg">0</span><span class="w" id="mwww"> </span><span class="s" id="mwxA">0</span><span class="w" id="mwxQ"> </span><span class="s" id="mwxg">0</span><span class="p" id="mwxw">;</span><span class="w" id="mwyA"> </span><span class="n" id="mwyQ">1</span><span class="w" id="mwyg"> </span><span class="s" id="mwyw">1</span><span class="w" id="mwzA"> </span><span class="s" id="mwzQ">1]</span><span class="p" id="mwzg">;</span><span class="w" id="mwzw"> </span> <span class="c" id="mw0A">% Edge Detection Process </span> <span class="c" id="mw0Q">% When i = 1 and j = 1, then filtered_image pixel </span> <span class="c" id="mw0g">% position will be filtered_image(2, 2) </span> <span class="c" id="mw0w">% The mask is of 3x3, so we need to traverse </span> <span class="c" id="mw1A">% to filtered_image(size(input_image, 1) - 2 </span> <span class="c" id="mw1Q">%, size(input_image, 2) - 2) </span> <span class="c" id="mw1g">% Thus we are not considering the borders. </span> <span class="k" id="mw1w">for</span><span class="w" id="mw2A"> </span><span class="nb" id="mw2Q">i</span><span class="w" id="mw2g"> </span><span class="p" id="mw2w">=</span><span class="w" id="mw3A"> </span><span class="mi" id="mw3Q">1</span><span class="p" id="mw3g">:</span><span class="nb" id="mw3w">size</span><span class="p" id="mw4A">(</span><span class="n" id="mw4Q">input_image</span><span class="p" id="mw4g">,</span><span class="w" id="mw4w"> </span><span class="mi" id="mw5A">1</span><span class="p" id="mw5Q">)</span><span class="w" id="mw5g"> </span><span class="o" id="mw5w">-</span><span class="w" id="mw6A"> </span><span class="mi" id="mw6Q">2</span><span class="w" id="mw6g"> </span> <span class="w" id="mw6w"> </span><span class="k" id="mw7A">for</span><span class="w" id="mw7Q"> </span><span class="nb" id="mw7g">j</span><span class="w" id="mw7w"> </span><span class="p" id="mw8A">=</span><span class="w" id="mw8Q"> </span><span class="mi" id="mw8g">1</span><span class="p" id="mw8w">:</span><span class="nb" id="mw9A">size</span><span class="p" id="mw9Q">(</span><span class="n" id="mw9g">input_image</span><span class="p" id="mw9w">,</span><span class="w" id="mw-A"> </span><span class="mi" id="mw-Q">2</span><span class="p" id="mw-g">)</span><span class="w" id="mw-w"> </span><span class="o" id="mw_A">-</span><span class="w" id="mw_Q"> </span><span class="mi" id="mw_g">2</span><span class="w" id="mw_w"> </span> <span class="w" id="mwAQA"> </span><span class="c" id="mwAQE">% Gradient approximations </span> <span class="w" id="mwAQI"> </span><span class="n" id="mwAQM">Gx</span><span class="w" id="mwAQQ"> </span><span class="p" id="mwAQU">=</span><span class="w" id="mwAQY"> </span><span class="nb" id="mwAQc">sum</span><span class="p" id="mwAQg">(</span><span class="nb" id="mwAQk">sum</span><span class="p" id="mwAQo">(</span><span class="n" id="mwAQs">Mx</span><span class="o" id="mwAQw">.*</span><span class="n" id="mwAQ0">input_image</span><span class="p" id="mwAQ4">(</span><span class="nb" id="mwAQ8">i</span><span class="p" id="mwARA">:</span><span class="nb" id="mwARE">i</span><span class="o" id="mwARI">+</span><span class="mi" id="mwARM">2</span><span class="p" id="mwARQ">,</span><span class="w" id="mwARU"> </span><span class="nb" id="mwARY">j</span><span class="p" id="mwARc">:</span><span class="nb" id="mwARg">j</span><span class="o" id="mwARk">+</span><span class="mi" id="mwARo">2</span><span class="p" id="mwARs">)));</span><span class="w" id="mwARw"> </span> <span class="w" id="mwAR0"> </span><span class="n" id="mwAR4">Gy</span><span class="w" id="mwAR8"> </span><span class="p" id="mwASA">=</span><span class="w" id="mwASE"> </span><span class="nb" id="mwASI">sum</span><span class="p" id="mwASM">(</span><span class="nb" id="mwASQ">sum</span><span class="p" id="mwASU">(</span><span class="n" id="mwASY">My</span><span class="o" id="mwASc">.*</span><span class="n" id="mwASg">input_image</span><span class="p" id="mwASk">(</span><span class="nb" id="mwASo">i</span><span class="p" id="mwASs">:</span><span class="nb" id="mwASw">i</span><span class="o" id="mwAS0">+</span><span class="mi" id="mwAS4">2</span><span class="p" id="mwAS8">,</span><span class="w" id="mwATA"> </span><span class="nb" id="mwATE">j</span><span class="p" id="mwATI">:</span><span class="nb" id="mwATM">j</span><span class="o" id="mwATQ">+</span><span class="mi" id="mwATU">2</span><span class="p" id="mwATY">)));</span><span class="w" id="mwATc"> </span> <span class="w" id="mwATg"> </span> <span class="w" id="mwATk"> </span><span class="c" id="mwATo">% Calculate magnitude of vector </span> <span class="w" id="mwATs"> </span><span class="n" id="mwATw">filtered_image</span><span class="p" id="mwAT0">(</span><span class="nb" id="mwAT4">i</span><span class="o" id="mwAT8">+</span><span class="mi" id="mwAUA">1</span><span class="p" id="mwAUE">,</span><span class="w" id="mwAUI"> </span><span class="nb" id="mwAUM">j</span><span class="o" id="mwAUQ">+</span><span class="mi" id="mwAUU">1</span><span class="p" id="mwAUY">)</span><span class="w" id="mwAUc"> </span><span class="p" id="mwAUg">=</span><span class="w" id="mwAUk"> </span><span class="nb" id="mwAUo">sqrt</span><span class="p" id="mwAUs">(</span><span class="n" id="mwAUw">Gx</span><span class="o" id="mwAU0">.^</span><span class="mi" id="mwAU4">2</span><span class="w" id="mwAU8"> </span><span class="o" id="mwAVA">+</span><span class="w" id="mwAVE"> </span><span class="n" id="mwAVI">Gy</span><span class="o" id="mwAVM">.^</span><span class="mi" id="mwAVQ">2</span><span class="p" id="mwAVU">);</span><span class="w" id="mwAVY"> </span> <span class="w" id="mwAVc"> </span> <span class="w" id="mwAVg"> </span><span class="k" id="mwAVk">end</span> <span class="k" id="mwAVo">end</span> <span class="c" id="mwAVs">% Displaying Filtered Image </span> <span class="n" id="mwAVw">filtered_image</span><span class="w" id="mwAV0"> </span><span class="p" id="mwAV4">=</span><span class="w" id="mwAV8"> </span><span class="nb" id="mwAWA">uint8</span><span class="p" id="mwAWE">(</span><span class="n" id="mwAWI">filtered_image</span><span class="p" id="mwAWM">);</span><span class="w" id="mwAWQ"> </span> <span class="nb" id="mwAWU">figure</span><span class="p" id="mwAWY">,</span><span class="w" id="mwAWc"> </span><span class="nb" id="mwAWg">imshow</span><span class="p" id="mwAWk">(</span><span class="n" id="mwAWo">filtered_image</span><span class="p" id="mwAWs">);</span><span class="w" id="mwAWw"> </span><span class="nb" id="mwAW0">title</span><span class="p" id="mwAW4">(</span><span class="s" id="mwAW8">'Filtered Image'</span><span class="p" id="mwAXA">);</span><span class="w" id="mwAXE"> </span> <span class="c" id="mwAXI">% Define a threshold value </span> <span class="n" id="mwAXM">thresholdValue</span><span class="w" id="mwAXQ"> </span><span class="p" id="mwAXU">=</span><span class="w" id="mwAXY"> </span><span class="mi" id="mwAXc">100</span><span class="p" id="mwAXg">;</span><span class="w" id="mwAXk"> </span><span class="c" id="mwAXo">% varies between [0 255] </span> <span class="n" id="mwAXs">output_image</span><span class="w" id="mwAXw"> </span><span class="p" id="mwAX0">=</span><span class="w" id="mwAX4"> </span><span class="nb" id="mwAX8">max</span><span class="p" id="mwAYA">(</span><span class="n" id="mwAYE">filtered_image</span><span class="p" id="mwAYI">,</span><span class="w" id="mwAYM"> </span><span class="n" id="mwAYQ">thresholdValue</span><span class="p" id="mwAYU">);</span><span class="w" id="mwAYY"> </span> <span class="n" id="mwAYc">output_image</span><span class="p" id="mwAYg">(</span><span class="n" id="mwAYk">output_image</span><span class="w" id="mwAYo"> </span><span class="o" id="mwAYs">==</span><span class="w" id="mwAYw"> </span><span class="nb" id="mwAY0">round</span><span class="p" id="mwAY4">(</span><span class="n" id="mwAY8">thresholdValue</span><span class="p" id="mwAZA">))</span><span class="w" id="mwAZE"> </span><span class="p" id="mwAZI">=</span><span class="w" id="mwAZM"> </span><span class="mi" id="mwAZQ">0</span><span class="p" id="mwAZU">;</span><span class="w" id="mwAZY"> </span> <span class="c" id="mwAZc">% Displaying Output Image </span> <span class="n" id="mwAZg">output_image</span><span class="w" id="mwAZk"> </span><span class="p" id="mwAZo">=</span><span class="w" id="mwAZs"> </span><span class="n" id="mwAZw">im2bw</span><span class="p" id="mwAZ0">(</span><span class="n" id="mwAZ4">output_image</span><span class="p" id="mwAZ8">);</span><span class="w" id="mwAaA"> </span> <span class="nb" id="mwAaE">figure</span><span class="p" id="mwAaI">,</span><span class="w" id="mwAaM"> </span><span class="nb" id="mwAaQ">imshow</span><span class="p" id="mwAaU">(</span><span class="n" id="mwAaY">output_image</span><span class="p" id="mwAac">);</span><span class="w" id="mwAag"> </span><span class="nb" id="mwAak">title</span><span class="p" id="mwAao">(</span><span class="s" id="mwAas">'Edge Detected Image'</span><span class="p" id="mwAaw">);</span>

More in Robotics

AnthropomorphismAttribution of human traits to non-human entities Meta PlatformsAmerican multinational technology conglomerate Space probeUnmanned robotic spacecraft StereoscopyTechnique for creating or enhancing the impression of depth in an image Control systemSystem that manages the behavior of other systems Inertial navigation systemContinuously computed dead reckoning