Home Backend Development Python Tutorial How to implement basic morphological filtering in Python

How to implement basic morphological filtering in Python

May 21, 2023 pm 11:29 PM
python

There are four most basic morphological operations, namely corrosion, expansion, open calculation and closed calculation. `scipy.ndimage implements these four operations for binary arrays and grayscale arrays respectively

##binary_erosiongrey_erosionCorrosionbinary_dilationgrey_dilationInflationbinary_closinggrey_closingClose (expand first and then corrode)binary_openinggrey_openingOpen (corrode first and then expand)
Binary Grayscale
Binary Morphology

The so-called corrosion is expressed in mathematical symbols as

How to implement basic morphological filtering in Python

where B

ij represents the set of all values ​​​​that are 1 in B when the origin of B BB is at (i, j).

This formula means that when structure B is used to corrode A, when the origin of B is translated to the pixel (i,j) of image A, if B is completely surrounded by the overlapping area of ​​the two, then Assigned a value of 1, otherwise assigned a value of 0. When a certain element in B is 1, if the corresponding position in A is also 1, then the value at (i, j) is 1. This is a more intuitive example.

Inflation is the opposite, and can be expressed as

How to implement basic morphological filtering in Python

#In other words, as long as the overlapping area of ​​B and A is not an empty set, then (i,j) point Just set it to 1.

An example is as follows

import numpy as np
import matplotlib.pyplot as plt
import scipy.ndimage as sn

x = np.zeros([20,20])
x[5:15, 5:15] = 1
x_ero = sn.binary_erosion(x)
x_dil = sn.binary_dilation(x)

fig = plt.figure()
ax = fig.add_subplot(1,3,1)
ax.imshow(x)
plt.title("original")
ax = fig.add_subplot(1,3,2)
ax.imshow(x_ero)
plt.title("erosion")
ax = fig.add_subplot(1,3,3)
ax.imshow(x_dil)
plt.title("dilation")
plt.show()
Copy after login

The effect is as follows

How to implement basic morphological filtering in Python

The opening operation is to corrode first and then expand; the closing operation is to expand first and then corrode. The example is as follows

x = np.zeros([20,20])
x[5:15, 5:15] = 1
x[10:12,10:12] = 0
x[2:4, 2:4] = 1

x_open = sn.binary_opening(x)
x_close = sn.binary_closing(x)

fig = plt.figure()
ax = fig.add_subplot(1,3,1)
ax.imshow(x)
plt.title("original")
ax = fig.add_subplot(1,3,2)
ax.imshow(x_open)
plt.title("opening")
ax = fig.add_subplot(1,3,3)
ax.imshow(x_close)
plt.title("closing")
plt.show()
Copy after login

The effect is as follows. It can be seen that the opening operation will remove the isolated 1, and the closing operation will remove the isolated 0.

How to implement basic morphological filtering in Python

Grayscale morphology

The corrosion, expansion and opening and closing operations of grayscale images are an extension of its binary situation, using similar For the logic of convolution, the staircase image is directly retrieved from

scipy, and the erosion, expansion, and opening and closing operations are performed in sequence.

from scipy.misc import ascent
img = ascent()

funcs = {
    "original": lambda x, tmp:x,
    "erosion" : sn.grey_erosion,
    "dilation" : sn.grey_dilation,
    "opening" : sn.grey_opening,
    "closing" : sn.grey_closing
}

fig = plt.figure()
for i, key in enumerate(funcs):
    ax = fig.add_subplot(2,3,i+1)
    plt.imshow(funcs[key](img, (10,10)), cmap=plt.cm.gray)
    plt.title(key)

plt.show()
Copy after login

Parameter list

The parameters of the binary function and the grayscale function are not the same. The following takes the closing operation as an example. All parameters of the binary and grayscale functions , in addition to the input input, the parameters shared by both are

  • structure is an array type, indicating the construction element, which can be understood as a convolution template

  • output is an array with the same dimensions as the input, and can save the result

  • orgin filter setting, the default is 0

Binary form Other parameters for learning filtering are as follows

binary_closing(input, iterations=1, mask=None, border_value=0, brute_force=False)
Copy after login

  • iterations is the number of executions

  • mask mask array, which is bool type Array, the position corresponding to False will not change

  • border_value The value at the edge

  • brute_force If it is False, only the last iteration Only the changed values ​​in will be updated

  • grey_closing(input, size=None, footprint=None, mode='reflect', cval=0.0)
    Copy after login

    • size is the filter template

    • mode optional reflect,constant, nearest, mirror, wrap, edge filling method

    • cval edge filling value

    The above is the detailed content of How to implement basic morphological filtering in Python. For more information, please follow other related articles on the PHP Chinese website!

Statement of this Website
The content of this article is voluntarily contributed by netizens, and the copyright belongs to the original author. This site does not assume corresponding legal responsibility. If you find any content suspected of plagiarism or infringement, please contact admin@php.cn

Hot AI Tools

Undresser.AI Undress

Undresser.AI Undress

AI-powered app for creating realistic nude photos

AI Clothes Remover

AI Clothes Remover

Online AI tool for removing clothes from photos.

Undress AI Tool

Undress AI Tool

Undress images for free

Clothoff.io

Clothoff.io

AI clothes remover

AI Hentai Generator

AI Hentai Generator

Generate AI Hentai for free.

Hot Article

R.E.P.O. Energy Crystals Explained and What They Do (Yellow Crystal)
2 weeks ago By 尊渡假赌尊渡假赌尊渡假赌
Repo: How To Revive Teammates
4 weeks ago By 尊渡假赌尊渡假赌尊渡假赌
Hello Kitty Island Adventure: How To Get Giant Seeds
4 weeks ago By 尊渡假赌尊渡假赌尊渡假赌

Hot Tools

Notepad++7.3.1

Notepad++7.3.1

Easy-to-use and free code editor

SublimeText3 Chinese version

SublimeText3 Chinese version

Chinese version, very easy to use

Zend Studio 13.0.1

Zend Studio 13.0.1

Powerful PHP integrated development environment

Dreamweaver CS6

Dreamweaver CS6

Visual web development tools

SublimeText3 Mac version

SublimeText3 Mac version

God-level code editing software (SublimeText3)

Can the Python interpreter be deleted in Linux system? Can the Python interpreter be deleted in Linux system? Apr 02, 2025 am 07:00 AM

Regarding the problem of removing the Python interpreter that comes with Linux systems, many Linux distributions will preinstall the Python interpreter when installed, and it does not use the package manager...

How to solve the problem of Pylance type detection of custom decorators in Python? How to solve the problem of Pylance type detection of custom decorators in Python? Apr 02, 2025 am 06:42 AM

Pylance type detection problem solution when using custom decorator In Python programming, decorator is a powerful tool that can be used to add rows...

How to solve permission issues when using python --version command in Linux terminal? How to solve permission issues when using python --version command in Linux terminal? Apr 02, 2025 am 06:36 AM

Using python in Linux terminal...

Python 3.6 loading pickle file error ModuleNotFoundError: What should I do if I load pickle file '__builtin__'? Python 3.6 loading pickle file error ModuleNotFoundError: What should I do if I load pickle file '__builtin__'? Apr 02, 2025 am 06:27 AM

Loading pickle file in Python 3.6 environment error: ModuleNotFoundError:Nomodulenamed...

Do FastAPI and aiohttp share the same global event loop? Do FastAPI and aiohttp share the same global event loop? Apr 02, 2025 am 06:12 AM

Compatibility issues between Python asynchronous libraries In Python, asynchronous programming has become the process of high concurrency and I/O...

What should I do if the '__builtin__' module is not found when loading the Pickle file in Python 3.6? What should I do if the '__builtin__' module is not found when loading the Pickle file in Python 3.6? Apr 02, 2025 am 07:12 AM

Error loading Pickle file in Python 3.6 environment: ModuleNotFoundError:Nomodulenamed...

How to ensure that the child process also terminates after killing the parent process via signal in Python? How to ensure that the child process also terminates after killing the parent process via signal in Python? Apr 02, 2025 am 06:39 AM

The problem and solution of the child process continuing to run when using signals to kill the parent process. In Python programming, after killing the parent process through signals, the child process still...

See all articles