Image filtering (gaussian_filter) in SciPy

Introduction

We use image filtering to smooth images and reduce noise. The gaussian_filter helps blur an image gently by averaging nearby pixels with a bell-shaped curve.

To remove small specks or noise from a photo taken in low light.

To prepare an image before detecting edges or shapes.

To create a soft blur effect for artistic purposes.

To smooth data in scientific images like medical scans.

To reduce detail before compressing an image.

Syntax

SciPy

scipy.ndimage.gaussian_filter(input, sigma, order=0, output=None, mode='reflect', cval=0.0, truncate=4.0)

input: The image or array to filter.

sigma: Controls how much to blur. Bigger sigma means more blur.

Examples

Apply a mild blur with sigma 1 to the image.

SciPy

from scipy.ndimage import gaussian_filter

blurred = gaussian_filter(image, sigma=1)

Apply a stronger blur with sigma 3 for more smoothing.

SciPy

blurred = gaussian_filter(image, sigma=3)

Apply different blur amounts for each axis in a 2D image.

SciPy

blurred = gaussian_filter(image, sigma=[1, 2])

Sample Program

This code creates a small noisy image with a bright square in the center. Then it smooths the image using gaussian_filter with sigma=1. Finally, it prints both the original and blurred images as arrays rounded to 2 decimals.

SciPy

import numpy as np
from scipy.ndimage import gaussian_filter
import matplotlib.pyplot as plt

# Create a simple 2D image with noise
np.random.seed(0)
image = np.zeros((10, 10))
image[4:6, 4:6] = 10  # bright square in the middle
image += np.random.normal(0, 1, image.shape)  # add noise

# Apply gaussian filter with sigma=1
blurred_image = gaussian_filter(image, sigma=1)

# Print original and blurred images
print('Original image:\n', np.round(image, 2))
print('\nBlurred image:\n', np.round(blurred_image, 2))

Important Notes

The sigma value controls the blur strength. Try small values first.

Gaussian filter works well to reduce noise but keeps edges smoother than simple averaging.

Use mode parameter to control how edges are handled (default is 'reflect').

Summary

Gaussian filter smooths images by averaging pixels with a bell curve.

It helps reduce noise and prepare images for further analysis.

Adjust sigma to control how much blur you want.

Practice

1. What is the main purpose of using gaussian_filter in image processing?

easy

A. To detect edges in the image

B. To smooth the image by reducing noise

C. To convert the image to grayscale

D. To increase the image resolution

Solution

1. Step 1: Understand the function's role

   gaussian_filter applies a blur effect that smooths the image by averaging nearby pixels weighted by a Gaussian curve.

2. Step 2: Identify the effect on image quality

   This smoothing reduces noise and small details, making the image less sharp but cleaner.

3. Final Answer:

   To smooth the image by reducing noise -> Option B

4. Quick Check:

   Gaussian blur = noise reduction [OK]

Hint: Gaussian filter smooths images by blurring noise away [OK]

Common Mistakes:

• Thinking it increases resolution
• Confusing it with edge detection
• Assuming it changes color format

2. Which of the following is the correct way to import and apply a Gaussian filter with sigma=2 to a 2D numpy array named image?

easy

A. from scipy.ndimage import gaussian_filter
filtered = gaussian_filter(image, sigma=2)

B. import scipy
filtered = scipy.gaussian_filter(image, 2)

C. from scipy import gaussian_filter
filtered = gaussian_filter(image, 2)

D. import gaussian_filter from scipy.ndimage
filtered = gaussian_filter(image, sigma=2)

Solution

1. Step 1: Check the correct import statement

   The Gaussian filter is in scipy.ndimage, so import it with from scipy.ndimage import gaussian_filter.

2. Step 2: Verify function usage

   Apply it by calling gaussian_filter(image, sigma=2) to blur with sigma 2.

3. Final Answer:

   from scipy.ndimage import gaussian_filter filtered = gaussian_filter(image, sigma=2) -> Option A

4. Quick Check:

   Correct import and sigma usage = from scipy.ndimage import gaussian_filter filtered = gaussian_filter(image, sigma=2) [OK]

Hint: Import from scipy.ndimage, use sigma as keyword [OK]

Common Mistakes:

• Wrong import path for gaussian_filter
• Passing sigma as positional without keyword
• Using incorrect import syntax

3. Given the code below, what will be the output array after applying the Gaussian filter?

import numpy as np
from scipy.ndimage import gaussian_filter

image = np.array([[0, 0, 0],
                  [0, 10, 0],
                  [0, 0, 0]])
filtered = gaussian_filter(image, sigma=1)
print(np.round(filtered, 2))

medium

A. [[0.46 1.23 0.46]
[1.23 2.99 1.23]
[0.46 1.23 0.46]]

B. [[0 0 0]
[0 10 0]
[0 0 0]]

C. [[2.5 2.5 2.5]
[2.5 2.5 2.5]
[2.5 2.5 2.5]]

D. [[0.1 0.1 0.1]
[0.1 10 0.1]
[0.1 0.1 0.1]]

Solution

1. Step 1: Understand Gaussian filter effect on sparse image

   The single bright pixel (10) will blur into neighbors, spreading intensity smoothly.

2. Step 2: Calculate approximate blurred values

   Using sigma=1, the center pixel reduces from 10 to about 3, neighbors get values around 1.2 to 0.46.

3. Final Answer:

   [[0.46 1.23 0.46] [1.23 2.99 1.23] [0.46 1.23 0.46]] -> Option A

4. Quick Check:

   Blur spreads intensity smoothly = [[0.46 1.23 0.46] [1.23 2.99 1.23] [0.46 1.23 0.46]] [OK]

Hint: Gaussian blur spreads bright pixels softly to neighbors [OK]

Common Mistakes:

• Expecting no change in array
• Assuming uniform average instead of weighted blur
• Misreading sigma effect as sharpening

4. Identify the error in the following code snippet that applies a Gaussian filter:

import numpy as np
from scipy.ndimage import gaussian_filter

image = np.ones((5,5))
filtered = gaussian_filter(image, sigma='2')
print(filtered)

medium

A. gaussian_filter cannot be applied to numpy arrays

B. The array shape is invalid for gaussian_filter

C. Missing import for numpy

D. sigma should be a number, not a string

Solution

1. Step 1: Check parameter types

   The sigma parameter must be a numeric value (int or float), not a string.

2. Step 2: Identify the cause of error

   Passing '2' as a string causes a type error when the filter tries to compute the blur.

3. Final Answer:

   sigma should be a number, not a string -> Option D

4. Quick Check:

   Numeric sigma required, string causes error [OK]

Hint: Use numeric sigma, not string, to avoid errors [OK]

Common Mistakes:

• Passing sigma as string
• Assuming gaussian_filter needs special array types
• Ignoring import errors

5. You have a noisy grayscale image stored as a 2D numpy array. You want to smooth the image but keep edges relatively sharp. Which approach using gaussian_filter and its parameters is best?

hard

A. Apply gaussian_filter multiple times with sigma=3

B. Use a large sigma value like 5 to blur the image heavily

C. Use a small sigma value like 0.5 to reduce noise but preserve edges

D. Do not use gaussian_filter; use median filter instead

Solution

1. Step 1: Understand sigma effect on smoothing

   Small sigma values cause light blur, preserving edges better than large sigma which blurs heavily.

2. Step 2: Choose best sigma for noise reduction and edge preservation

   Using sigma=0.5 smooths noise but keeps edges sharper than larger sigma values or repeated blurring.

3. Final Answer:

   Use a small sigma value like 0.5 to reduce noise but preserve edges -> Option C

4. Quick Check:

   Small sigma = smooth noise + keep edges [OK]

Hint: Small sigma blurs less, preserving edges better [OK]

Common Mistakes:

• Using large sigma blurs edges too much
• Applying filter multiple times causes over-blur
• Confusing gaussian_filter with median filter