- Image Processing
- What happens if we increase the number of bits in an image?
- What are the different colour spaces available?
- Where is Grey-colour Space useful?
- The Fourier Transform
- What are the advantages derived in the frequency domain?
- How do we use the Fourier Transform in digital image processing?
- High-pass filter
- What information do high-frequency signals convey?
- What is a kernel?
- What is a mask?
- Why is the noise identified as an edge?
- How do we prevent misclassification of noise as an edge?
- Contour detection
- What is the return type of findContours function?
The interest in computer vision started with a quest to make machines comprehend images to enhance the possibilities of image processing.
In this article, we cover basic image processing techniques, with ready-to-use code, and examples to understand concepts better. Various low-level vision algorithms dealing with image features like edges, contours, etc are considered.
Python and OpenCV will be our go-to resources. You may get stuck during installation, coding, compiling, or error-handling. You can refer to StackOverflow, Python docs, etc, to clarify to get solutions to the problems at hand. We find instructions for installation of OpenCV and its dependencies on the website mentioned in the references.
Let us understand the basic building block of an image – the pixels. The pixels are the most fundamental blocks in an image. A high definition (HD) image has higher volume of these pixels and a lower resolution image would have fewer pixels. The world is wrapped in a continuous field of reality, and the thought of pixelating the world is mind-boggling – which can you an understanding of how difficult it would be to apply concepts of computer vision for real world applications
Consider the two images: grayscale and colour. In the former, each pixel ranges between 0 and 255, where 0 is black and 255 is white. In the latter, three-colour channels together present the various colours present in the vast spectrum.
Let us look at the code to read and write an image.
1. import cv2 2. image =cv2.imread("messi.jpg") 3. cv2.imshow("Image", image)
The first line imports the library OpenCV, and we are accessing a file named “messi.jpg” present in the same directory as the code. Download an image and name it as messi.jpg, and store it in the same directory as the one with the code in it. Only then will the above code work. Finally, we print the image onto the output screen with the dialog box named “Image”.
Images are arrays of integers ranging between two numbers, depending on the number of bits used to represent each pixel. For example, in a 1-bit image, the number of colour choices for each pixel is white(1) or black(0).
What happens if we increase the number of bits in an image?
An increase in the number of bits results a larger variation in the representation of colours, and thus the contrast, sharpness and brightness standards increase. These images, as mentioned earlier, are 2-D arrays with grayscale and 3-D arrays, with colour images. Thus, processing of images means manipulating these images strategically to get the desired results. We’ve mentioned a few image transformations below for your reference.
There are common techniques applied to images, including translation, rotation, re-sizing, flipping, and cropping. Adobe Photoshop implements low-level and mid-level feature transformations. A few tools that Photoshop provides are easy to code. Let us explore a few of these.
The translation is a technique for shifting of an image along the x-axis and y-axis. Used to shift an image up, down, left, or right, along with any combination of the above translation is a good place to start.
1.image = cv2.imread("messi.jpg") 2.cv2.imshow("Original", image) 3.translation_matrix = np.float32([[1, 0, 25], [0, 1, 50]]) 4.shifted = cv2.warpAffine(image, translation_matirx, (image.shape, image.shape)) 5.cv2.imshow("Shifted Down and Right", shifted) 6.M = np.float32([[1, 0, -50], [0, 1, -90]]) 7.shifted = cv2.warpAffine(image, translation_matrix, (image.shape, image.shape)) 8.cv2.imshow("Shifted Up and Left", shifted)
The actual translation takes place on lines 3-5, where we initially define our translation matrix translation_matirx. This matrix handles the number of pixels that shift in any of the directions. The translation matrix “translation_matrix“ is a floating-point array. OpenCV expects the matrix to be of floating-point type. The first row of the matrix is [1, 0, shift_num_x], where shift_num_x is the number of pixels of shift for the image. Negative values of shift_num_x will shift the image to the left and positive values will shift the image to the right. The second row of the matrix is [0, 1, shift_num_y], where shift_num_y is the number of pixels the image shifts, up or down. The negative value of shift_num_y will shift the image up and positive values will shift the image down.
Now that we define the translation matrix, the actual translation takes place on Line 4 using the cv2.warpAffine function.
What are the parameters for the cv2.warpAffine function?
The first argument is the image of interest and the second argument is the translation matrix translation_matrix.
The dimensions (width and height) of our image are provided in place of the third argument. Line 5 shows the results of the translation.
Re-sizing is the next transformation. The inbuilt function provided by OpenCV accomplishes this task. An example always makes it easier to understand.
1.image = cv2.imread("messi.jpg") 2.cv2.imshow("Original", image) 3.r = 150.0 / image.shape 4.dim = (150, int(image.shape * r)) 5.resized = cv2.resize(image, dim, interpolation = cv2.INTER_AREA) 6.cv2.imshow("Resized (Width)", resized)
In the above code, re-sizing happens in lines 3-5. In line 3, we calculate the aspect ratio of the image.
Rotation is the transformation for rotating the image by an angle theta. It is useful to rotate images, as implemented in the photo editors on our phones and computers. Consider implementing this on your own, with the knowledge from the earlier two examples. Given this problem of rotation, how do you approach it? The solution always lies in the questions we ask ourselves, right? We mention the code below, in case you get stuck.
1 import numpy as np 2 import argparse 3 import imutils 4 import cv2 5 image = cv2.imread("messi.jpg") 6 cv2.imshow("Original", image) 7 (h, w) = image.shape[:2] 8 center = (w // 2, h // 2) 9 matrix = cv2.getRotationMatrix2D(center, 60, 1.0) 10 rotated = cv2.warpAffine(image, matrix, (w, h)) 11 cv2.imshow("Rotated by 60 Degrees", rotated) 12 matrix = cv2.getRotationMatrix2D(center, -90, 1.0) 13 rotated = cv2.warpAffine(image, matrix, (w, h)) 14 cv2.imshow("Rotated by -90 Degrees", rotated)
The point around which anything rotates decides everything about the rotation and hence is of utmost importance. The same applies for images, and the point of rotation is the point of concern. We usually rotate the image around the center of an image; However, OpenCV allows us to specify any arbitrary point of rotation. We shall go ahead with rotation around the center of the image. Lines 8 and 9 calculates the center of the image by dividing the length and width of the image by 2. Integer division represented by “//“ ensures integer-type results. A matrix like the one used to translate an image will be required here. Instead of manually constructing the matrix using NumPy, a call to the cv2.getRotationMatrix2D method on Line 9 does the job.
The cv2.getRotationMatrix2D function works with three arguments:
l the point of rotation
l theta, the angle of rotation(60 in this case)
l the scale of the image.
Similarly, there are various other transformations that you can experiment with.
Colour Spaces: Colour spaces refer to the coordinate system of the geometry of colours in the image. For example, BGR and grayscale are two examples of the same. There are multiple colour spaces and many methods available in OpenCV that do the job of conversion between colour spaces.
What are the different colour spaces available?
Three-colour spaces are in major use in the modern-day computer vision society:
l Hue, Saturation, Value (HSV).
Gray works with one feature channel and thus eliminates the role of colour, translating to shades of grey.
Where is Grey-colour Space useful?
It is used heavily during the intermediate processing phase, such as face detection.
BGR is the blue-green-red colour space. Web developers work with this colour space with a slight difference in the order of colours: RGB. The ordering of BGR is unique to OpenCV. Finally, with HSV, the hue is the colour tone, saturation is the intensity of a colour, and the value is its darkness. This is a colour space that is counter-intuitive and working with this colour space shall deepen our understanding of the subject.
The Fourier Transform
Image processing is a digital task because of the involvement of processors. The processors need to complete the tasks fast, and for high efficiency and throughput, one must have highly efficient algorithms. Processing images in their raw form of multidimensional arrays is not an efficient choice as attaining top processing speeds is a challenge. Thus, the conversion of these digital signals into the frequency domain is of utmost importance.
What are the advantages derived in the frequency domain?
In the frequency domain, we alter the characteristics of the image in a way that enables parallel processing and a lesser number of computations to arrive at the same results. We represent a signal as a sum of intermediate signals, thus enabling parallel processing. The Fourier Transform converts images in the spatial domain to the frequency domain.
Who discovered the Fourier Transform?
Joseph Fourier was an 18th-century French mathematician. He discovered that we can express waveforms as the sum of simple waveforms of various frequencies. The waveforms we see are the sum of waveforms of different frequencies equal to or smaller than the actual waveform’s frequency.
How do we use the Fourier Transform in digital image processing?
This concept is useful in manipulating images because it allows us to find regions in images where an image pixel changes a lot and regions where the change is less. Using these characteristics, we can then term these as the region of interest(ROI), or noise, or the foreground of an image or background. These divisions can help in faster processing and extrapolation of interesting patterns and results.
A high-pass filter (HPF) allows only higher frequencies to pass through and it blocks the lower frequencies. The threshold set for the filter is called the cut-off point. An HPF is a filter that examines a region of an image and increases the intensity of pixels.
How does an HPF work?
It increases the intensity of pixels in a region based on the difference in the values of the intensity in the ROI.
Consider the following kernel:
[[0, -1, 0], [-1, 4, -1], [0, -1, 0]]
This kernel, when convoluted with an image, gives a mask that has information about all the high-frequency signals in the image.
What information do high-frequency signals convey?
High-frequency signals add details to the image, for example, edges. Thus, the addition of this mask to the original image results in highlighting the edges present and enhancement of the image.
What is a kernel?
A kernel is a set of weights when applied to ROI results in a single pixel in the destination image. For example, a kernel size of 7 implies that 49 (7 x 7) source pixels are considered in generating each destination pixel. The intensity of the central pixel is boosted (or not), after finding large variations in the sum of differences of the intensities. If a pixel stands out from the surrounding pixels, it will get boosted.
Both high pass and low-pass filters use a property called the radius. Radius is the reach of the kernel in the image in one stride. Below is an example for an HPF:
1.import cv2 2.import numpy as np 3.from scipy import ndimage 4.kernel_3x3 = np.array([[-1, -1, -1], [-1, 8, -1], [-1, -1, -1]])
The kernel mentioned here is convoluted with the original image to get a mask.
What is a mask?
The mask has the information that the filter has derived from the image via convolution.
While HPF works with the intensity difference in images, a low-pass filter (LPF) smoothens the pixels by averaging the ROI. The output is a single-pixel value which is the average of the pixels in the kernel. Applications for LPF are blurring and denoising. An example of the same is the Gaussian filter used to blur the images.
LPF kernel example
1/8*([0 , 1 ,,0], [1,4,1], [0,1,0])
Edges are the lowest level of features, and they are also the first piece of information that we as humans infer about through our visual system. Usually, most MOOCs start with the human visual system. It creates a sense of intrigue in anyone who reads about it. Software, too, can reason about edges, poses, and archetypes.
How do we implement filters on OpenCV?
OpenCV lists many edge detecting filters.
These filters are faulty, not in terms of the implementation, but by their algorithms itself. The expected outcome is to turn non-edge regions to black while turning regions having edges to saturated colours. The faulty nature arises when the algorithms work with noisy images, wherein it identifies noise as edges.
Why is the noise identified as an edge?
Noise can be of many types, like uniform noise, impulsive noise, and Gaussian noise. With impulsive noise, the values of the pixels transform and thus the difference in the intensity of pixels is high, thus getting misclassified as an edge.
How do we prevent misclassification of noise as an edge?
We can eradicate this flaw. Logically thinking, we need to remove the impulsive increase in the intensity values occurring in the image. This can be done by passing the image through an LPF, which ends up blurring the image before trying to find its edges. OpenCV provides a range of blurring filters, including blur(), medianBlur(), and GaussianBlur(). The arguments are implementation-dependent, and we suggest you read the documentation to get a better understanding of the arguments in each of these filters.
For blurring, let us use medianBlur(), which is effective in removing digital video noise, especially in colour images. For edge-finding, let us use Laplacian(), which produces bold edge lines, especially in grayscale images. After applying medianBlur(), but before applying Laplacian(), we should convert the image from BGR to grayscale.
Once we have the result of Laplacian(), we can invert it to get black edges on a white background. Then, we can normalize it (so that its values range from 0 to 1) and multiply it with the source image to darken the edges. We carry out the function that gets the task done:
1.def strokeEdges(src, dst, blurKsize = 7, edgeKsize = 5): 2. if blurKsize >= 3: 3. blurredSrc = cv2.medianBlur(src, blurKsize) 4. graySrc = cv2.cvtColor(blurredSrc, cv2.COLOR_BGR2GRAY) 5. else: 6. graySrc = cv2.cvtColor(src, cv2.COLOR_BGR2GRAY) 7. cv2.Laplacian(graySrc, cv2.CV_8U, graySrc, ksize = edgeKsize) 8. normalizedInverseAlpha = (1.0 / 255) * (255 - graySrc) 9. channels = cv2.split(src) 10. for channel in channels: 11. channel[:] = channel * normalizedInverseAlpha 12. cv2.merge(channels, dst)
Note that we allow kernel sizes to be specified as arguments for strokeEdges(). We use the blurKsize argument as kernel size for medianBlur(), while edgeKsize is used as kernel size for Laplacian().
If you meet performance problems when running strokeEdges(), try decreasing the blurKsize value. To turn off blur, set it to a value less than 3.
Another important computer vision technique is called contour detection that has gained popularity among scientists and researchers. It deals with the aspect of detecting contours in ROI in an image. It gains importance because of its derivative operations connected with identifying contours.
What are the operations related to contour detection?
1. generating bounding boxes
2. approximating shapes
3. Calculating ROIs
ROIs simplify the complexity involved with the interaction with images. We use these concepts in the fields of object detection and face detection.
We shall familiarize ourselves with an example:
1.import cv2 2.import numpy as np 3.img = np.zeros((200, 200), dtype=np.uint8) 4.img[50:150, 50:150] = 255 5.ret, thresh = cv2.threshold(img, 127, 255, 0) 6.image, contours, hierarchy = cv2.findContours(thresh,cv2.RETR_TREE, cv2.CHAIN_APPROX_SIMPLE) 8.color = cv2.cvtColor(img, cv2.COLOR_GRAY2BGR) 9.img = cv2.drawContours(color, contours, -1, (0,255,0), 2) 10.cv2.imshow("contours", color) 11.cv2.waitKey() 12.cv2.destroyAllWindows()
First, create an empty black image that is 200×200 pixels in size. Then, a black square is placed in the center of the square. Thresholding is done, and we call the findContours() function.
What are the parameters of findContours() function in OpenCV?
1. the selected image to be fed in as input
2. type of hierarchy
3. the contour approximation method.
Several aspects are of particular interest in this function. The hierarchy tree returned by the function is important: cv2.RETR_TREE will find all the contours present in the image and thus enable us to find relationships between contours. If we only want to retrieve the most external contours, we make use of cv2.RETR_EXTERNAL. This is useful when we want to drop contours that are a subset of other contours.
What is the return type of findContours function?
It returns three elements:
1. the changed image
We use the contours to draw on the colour version of the image (so we can draw contours in green) and display it. The result is a white square with its contour drawn in green.
We covered various algorithms related to the basics of image processing. We hope you enjoy implementing the same and have fun with the many algorithms learned in the article. These basics learned well shall take us a long way in our journey in image processing and computer vision. This is part 1 in a 4-article series on Computer Vision. Be sure to read the next article that deals with the basics of deep learning applied in the domain of computer vision.
Check out Great Learning’s PGP-Artificial Intelligence and Machine Learning which covers all AIML and deep learning concepts including computer vision.