medfilt3

三维中位数滤波

全页折叠

语法

B = medfilt3(A)

B = medfilt3(A,[m n p])

B = medfilt3(___,padopt)

说明

B = medfilt3(A) 使用 3×3×3 滤波器对三维图像 A 进行滤波。默认情况下，medfilt3 通过在边界处以镜像方式复制值来填充图像。

示例

B = medfilt3(A,[m n p]) 对三维图像 A 执行三维中位数滤波。B 中的每个输出体素包含 A 中对应体素周围 m×n×p 邻域的中位数值。

B = medfilt3(___,padopt) 控制 medfilt3 如何填充数组边界。

示例

全部折叠

使用中位数滤波去除三维数据中的离群值

打开实时脚本

创建一个含噪三维曲面。

[x,y,z,V] = flow(50);
noisyV = V + 0.1*double(rand(size(V))>0.95) - 0.1*double(rand(size(V))<0.05);

应用中位数滤波。

filteredV = medfilt3(noisyV);

一起显示含噪曲面和滤波曲面。

subplot(1,2,1)
hpatch1 = patch(isosurface(x,y,z,noisyV,0));
isonormals(x,y,z,noisyV,hpatch1)
set(hpatch1,FaceColor="red",EdgeColor="none")
daspect([1,4,4])
view([-65,20]) 
axis tight off
camlight left
lighting gouraud

subplot(1,2,2)
hpatch2 = patch(isosurface(x,y,z,filteredV,0));
isonormals(x,y,z,filteredV,hpatch2)
set(hpatch2,FaceColor="red",EdgeColor="none")
daspect([1,4,4])
view([-65,20])
axis tight off
camlight left 
lighting gouraud

Figure contains 2 axes objects. Hidden axes object 1 contains an object of type patch. Hidden axes object 2 contains an object of type patch.

输入参数

全部折叠

`A` — 输入图像
三维数值数组 | 三维逻辑数组

输入图像，指定为三维数值或逻辑数组。

`[m n p]` — 邻域大小
`[3 3 3]` (默认) | 三元素向量

邻域大小，指定为由正奇整数组成的三元素向量。

`padopt` — 填充选项
`"symmetric"` (默认) | `"zeros"` | `"replicate"`

填充选项，指定为下列值之一：

值	描述	示例
`"symmetric"`	用自身的镜像翻转填充数组	$[\begin{matrix} 3 & 1 & 4 \\ 1 & 5 & 9 \\ 2 & 6 & 5 \end{matrix}] \to [\begin{matrix} 5 & 1 & 1 & 5 & 9 & 9 & 5 \\ 1 & 3 & 3 & 1 & 4 & 4 & 1 \\ 1 & 3 & 3 & 1 & 4 & 4 & 1 \\ 5 & 1 & 1 & 5 & 9 & 9 & 5 \\ 6 & 2 & 2 & 6 & 5 & 5 & 6 \\ 6 & 2 & 2 & 6 & 5 & 5 & 6 \\ 5 & 1 & 1 & 5 & 9 & 9 & 5 \end{matrix}]$
`"replicate"`	通过重复边框元素填充数组	$[\begin{matrix} 3 & 1 & 4 \\ 1 & 5 & 9 \\ 2 & 6 & 5 \end{matrix}] \to [\begin{matrix} 3 & 3 & 3 & 1 & 4 & 4 & 4 \\ 3 & 3 & 3 & 1 & 4 & 4 & 4 \\ 3 & 3 & 3 & 1 & 4 & 4 & 4 \\ 1 & 1 & 1 & 5 & 9 & 9 & 9 \\ 2 & 2 & 2 & 6 & 5 & 5 & 5 \\ 2 & 2 & 2 & 6 & 5 & 5 & 5 \\ 2 & 2 & 2 & 6 & 5 & 5 & 5 \end{matrix}]$
`"zeros"`	用值 `0` 填充数组	$[\begin{matrix} 3 & 1 & 4 \\ 1 & 5 & 9 \\ 2 & 6 & 5 \end{matrix}] \to [\begin{matrix} 0 & 0 & 0 & 0 & 0 & 0 & 0 \\ 0 & 0 & 0 & 0 & 0 & 0 & 0 \\ 0 & 0 & 3 & 1 & 4 & 0 & 0 \\ 0 & 0 & 1 & 5 & 9 & 0 & 0 \\ 0 & 0 & 2 & 6 & 5 & 0 & 0 \\ 0 & 0 & 0 & 0 & 0 & 0 & 0 \\ 0 & 0 & 0 & 0 & 0 & 0 & 0 \end{matrix}]$

数据类型: char | string

输出参量

全部折叠

`B` — 输出图像
三维数值数组

输出图像，以与输入图像 A 大小和数据类型相同的数值三维数组形式返回。

扩展功能

全部展开

基于线程的环境
使用 MATLAB® `backgroundPool` 在后台运行代码或使用 Parallel Computing Toolbox™ `ThreadPool` 加快代码运行速度。

此函数完全支持基于线程的环境。有关详细信息，请参阅在基于线程的环境中运行 MATLAB 函数。

版本历史记录

在 R2016b 中推出

全部展开

R2022b: 支持基于线程的环境

medfilt3 现在支持基于线程的环境。

R2021b: 改进了对小邻域的处理性能

medfilt3 函数在处理从 [3, 3, 3] 到 [31, 31, 31] 的邻域大小时的性能得到改进。

例如，以下代码比以前的版本快 3 倍。

function timingTestMedfilt3
      load mristack;
      noisyV = imnoise(mristack,'salt & pepper',0.2);
      tic
      filteredV = medfilt3(noisyV);
      toc
end

大致的执行时间是：

R2021a：0.24 秒

R2021b：0.08 秒

代码是在运行 Windows^® 10 的Intel^® Xeon^® Gold 5220 CPU @ 2.2 GHz 测试系统上通过调用函数 timingTestMedfilt3 进行计时的。

另请参阅

medfilt2

主题

imfilter 边界填充选项

medfilt3

语法

说明

示例

使用中位数滤波去除三维数据中的离群值

输入参数

A — 输入图像 三维数值数组 | 三维逻辑数组

[m n p] — 邻域大小 [3 3 3] (默认) | 三元素向量

padopt — 填充选项 "symmetric" (默认) | "zeros" | "replicate"

输出参量

B — 输出图像 三维数值数组

扩展功能

基于线程的环境 使用 MATLAB® backgroundPool 在后台运行代码或使用 Parallel Computing Toolbox™ ThreadPool 加快代码运行速度。

版本历史记录

R2022b: 支持基于线程的环境

R2021b: 改进了对小邻域的处理性能

另请参阅

主题

`A` — 输入图像
三维数值数组 | 三维逻辑数组

`[m n p]` — 邻域大小
`[3 3 3]` (默认) | 三元素向量

`padopt` — 填充选项
`"symmetric"` (默认) | `"zeros"` | `"replicate"`

`B` — 输出图像
三维数值数组

基于线程的环境
使用 MATLAB® `backgroundPool` 在后台运行代码或使用 Parallel Computing Toolbox™ `ThreadPool` 加快代码运行速度。