## Introduction

In this post, we will see various ways of doing image subtraction in Python by using NumPy, OpenCV, and Pillow libraries. The subtraction of images may sound a bit strange to beginners but images consist of numeric pixels hence we can perform pixel-wise subtraction between two images.

**Subtraction of Images**

When one image is subtracted from another it does pixel-wise subtraction where the pixel of the 2nd image is subtracted from the 1st image resulting in a new 3rd image.

Mathematically, when image B is subtracted from image A, the pixel-wise subtraction operation to get the output image C can be represented as follows â€“

**C(x,y) = A(x,y) â€“ B(x,y)**

Where,

- A(x,y) represents the pixel value at position (x, y) in the image A.
- B(x, y) represents the pixel value at position (x, y) in the image B.
- C(x, y) represents the resulting pixel value at position (x, y) in the new image C.

A point to be noted here is that both images should be of the same size. Otherwise, we crop, resize, or pad to make both images the same size.

**Use of Image Subtraction**

Image subtraction has many practical use cases in computer vision, image processing, and computer graphics. Some of the regular use cases are â€“

- Detecting Changes in Images
- Background Subtraction
- Image Enhancement
- Image Registration

**Input Image for Our Examples**

We will use the images below for all our examples of image subtraction with Python.

**Image 1**

**Image 2**

**Image Subtraction with Python NumPy**

In this example, we first read the two images and convert them into NumPy arrays. Next, the 2nd array is subtracted from the 1st array by element-wise subtraction. We then take the absolute difference to ensure that the result is non-negative and then convert the result back to an unsigned 8-bit integer array. Finally the subtracted image is displayed as output.

Â

import cv2 import numpy as np # Load the two images image1 = cv2.imread('image1.jpg') image2 = cv2.imread('image2.jpg') # Convert images to NumPy arrays np_image1 = np.array(image1, dtype=np.int16) np_image2 = np.array(image2, dtype=np.int16) # Perform image subtraction subtracted_image = np.abs(np_image1 - np_image2).astype(np.uint8) # Display the result cv2.imshow('Subtracted Image (NumPy)', subtracted_image) cv2.waitKey(0) cv2.destroyAllWindows()

**Output**

Â

**Image Subtraction with Python OpenCV**

Here we we first read the two images with and then use cv2.subtract() function of the OpenCV library to subtract the 2nd image from the 1st image. Then the subtracted image is displayed as output.

Â

import cv2 # Load the two images image1 = cv2.imread('image1.jpg') image2 = cv2.imread('image2.jpg') # Perform image subtraction subtracted_image = cv2.subtract(image2, image1) # Display the result cv2.imshow('Subtracted Image (OpenCV)', subtracted_image) cv2.waitKey(0) cv2.destroyAllWindows()

**Output**

Â

**Image Subtraction with Python Pillow**

In this example, the two images are opened and then use ImageChops.subtract() function of the PIL library to subtract the 2nd image from the 1st image and then show the image as output.

Â

from PIL import Image, ImageChops # Load the two images image1 = Image.open('image1.jpg') image2 = Image.open('image2.jpg') # Perform image subtraction subtracted_image = ImageChops.subtract(image2, image1) # Display the result subtracted_image.show()

**Output**

Â

Â