Image restoration, Noise, Segmentation and Contours
So far, you have done some very cool things with your image processing skills! In this chapter, you will apply image restoration to remove objects, logos, text, or damaged areas in pictures! You will also learn how to apply noise, use segmentation to speed up processing, and find elements in images by their contours. This is the Summary of lecture "Image Processing in Python", via datacamp.
import numpy as np
import matplotlib.pyplot as plt
import pandas as pd
plt.rcParams['figure.figsize'] = (10, 8)
Let's restore a damaged image
In this exercise, we'll restore an image that has missing parts in it, using the inpaint_biharmonic() function.
We'll work on an image with damaged. Some of the pixels have been replaced by 1s using a binary mask, on purpose, to simulate a damaged image. Replacing pixels with 1s turns them totally black.
The mask is a black and white image with patches that have the position of the image bits that have been corrupted. We can apply the restoration function on these areas.
Remember that inpainting is the process of reconstructing lost or deteriorated parts of images and videos.
def show_image(image, title='Image', cmap_type='gray'):
plt.imshow(image, cmap=cmap_type)
plt.title(title)
plt.axis('off')
def plot_comparison(img_original, img_filtered, img_title_filtered):
fig, (ax1, ax2) = plt.subplots(ncols=2, figsize=(10, 8), sharex=True, sharey=True)
ax1.imshow(img_original, cmap=plt.cm.gray)
ax1.set_title('Original')
ax1.axis('off')
ax2.imshow(img_filtered, cmap=plt.cm.gray)
ax2.set_title(img_title_filtered)
ax2.axis('off')
from skimage.restoration import inpaint
from skimage.transform import resize
from skimage import color
defect_image = plt.imread('./dataset/damaged_astronaut.png')
defect_image = resize(defect_image, (512, 512))
defect_image = color.rgba2rgb(defect_image)
mask = pd.read_csv('./dataset/astronaut_mask.csv').to_numpy()
# Apply the restoration function to the image using the mask
restored_image = inpaint.inpaint_biharmonic(defect_image, mask, multichannel=True)
# Show ther defective image
plot_comparison(defect_image, restored_image, 'Restored image')
Removing logos
As we saw in the video, another use of image restoration is removing objects from an scene. In this exercise, we'll remove the Datacamp logo from an image.
You will create and set the mask to be able to erase the logo by inpainting this area.
Remember that when you want to remove an object from an image you can either manually delineate that object or run some image analysis algorithm to find it.
image_with_logo = plt.imread('./dataset/4.2.06_w_logo_2_2.png')
# Initialize the mask
mask = np.zeros(image_with_logo.shape[:-1])
# Set the pixels where the logo is to 1
mask[210:272, 360:425] = 1
# Apply inpainting to remove the logo
image_logo_removed = inpaint.inpaint_biharmonic(image_with_logo,
mask,
multichannel=True)
# Show the original and logo removed images
plot_comparison(image_with_logo, image_logo_removed, 'Image with logo removed')
from skimage.util import random_noise
fruit_image = plt.imread('./dataset/fruits_square.jpg')
# Add noise to the image
noisy_image = random_noise(fruit_image)
# Show th original and resulting image
plot_comparison(fruit_image, noisy_image, 'Noisy image')
from skimage.restoration import denoise_tv_chambolle
noisy_image = plt.imread('./dataset/miny.jpeg')
# Apply total variation filter denoising
denoised_image = denoise_tv_chambolle(noisy_image, multichannel=True)
# Show the noisy and denoised image
plot_comparison(noisy_image, denoised_image, 'Denoised Image')
from skimage.restoration import denoise_bilateral
landscape_image = plt.imread('./dataset/noise-noisy-nature.jpg')
# Apply bilateral filter denoising
denoised_image = denoise_bilateral(landscape_image, multichannel=True)
# Show original and resulting images
plot_comparison(landscape_image, denoised_image, 'Denoised Image')
Superpixels & segmentation
- Superpixel
- A group of connected pixels with similar colors or gray levels
- Benefits
- meaningful regions
- Confutational efficiency
- Segmentation
- supervised
- unsupervised
- Simple Linear Iterative Clustering (CLIC): takes all the pixel values of the image and tries to separate them into a predifined number of sub-regions
from skimage.segmentation import slic
from skimage.color import label2rgb
face_image = plt.imread('./dataset/chinese.jpg')
# Obtain the segmentation with 400 regions
segments = slic(face_image, n_segments=400)
# Put segments on top of original image to compare
segmented_image = label2rgb(segments, face_image, kind='avg')
# Show the segmented image
plot_comparison(face_image, segmented_image, 'Segmented image, 400 superpixels')
Contouring shapes
In this exercise we'll find the contour of a horse.
For that we will make use of a binarized image provided by scikit-image in its data module. Binarized images are easier to process when finding contours with this algorithm. Remember that contour finding only supports 2D image arrays.
Once the contour is detected, we will display it together with the original image. That way we can check if our analysis was correct!
def show_image_contour(image, contours):
plt.figure()
for n, contour in enumerate(contours):
plt.plot(contour[:, 1], contour[:, 0], linewidth=3)
plt.imshow(image, interpolation='nearest', cmap='gray_r')
plt.title('Contours')
plt.axis('off')
from skimage import measure, data
# Obtain the horse image
horse_image = data.horse()
# Find the contours with a constant level value of 0.8
contours = measure.find_contours(horse_image, level=0.8)
# Shows the image with contours found
show_image_contour(horse_image, contours)
Find contours of an image that is not binary
Let's work a bit more on how to prepare an image to be able to find its contours and extract information from it.
We'll process an image of two purple dices loaded as image_dices and determine what number was rolled for each dice.
In this case, the image is not grayscale or binary yet. This means we need to perform some image pre-processing steps before looking for the contours. First, we'll transform the image to a 2D array grayscale image and next apply thresholding. Finally, the contours are displayed together with the original image.
from skimage.io import imread
from skimage.filters import threshold_otsu
image_dices = imread('./dataset/dices.png')
# Make the image grayscale
image_dices = color.rgb2gray(image_dices)
# Obtain the optimal thresh value
thresh = threshold_otsu(image_dices)
# Apply thresholding
binary = image_dices > thresh
# Find contours at a constant value of 0.8
contours = measure.find_contours(binary, level=0.8)
# Show the image
show_image_contour(image_dices, contours)
Count the dots in a dice's image
Now we have found the contours, we can extract information from it.
In the previous exercise, we prepared a purple dices image to find its contours:
This time we'll determine what number was rolled for the dice, by counting the dots in the image.
Create a list with all contour's shapes as shape_contours. You can see all the contours shapes by calling shape_contours in the console, once you have created it.
Check that most of the contours aren't bigger in size than 50. If you count them, they are the exact number of dots in the image.
shape_contours = [cnt.shape[0] for cnt in contours]
# Set 50 as the maximum size of the dots shape
max_dots_shape = 50
# Count dots in contours excluding bigger than dots size
dots_contours = [cnt for cnt in contours if np.shape(cnt)[0] < max_dots_shape]
# Shows all contours found
show_image_contour(binary, contours)
# Print the dice's number
print('Dice`s dots number: {}.'.format(len(dots_contours)))
