medfilt3

Filtrado de mediana de 3D

Sintaxis

B = medfilt3(A)

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

B = medfilt3(___,padopt)

Descripción

B = medfilt3(A) filtra la imagen 3D A con un filtro de 3 por 3 por 3. Por defecto, medfilt3 rellena la imagen replicando los valores de forma reflejada en los bordes.

B = medfilt3(A,[m n p]) realiza un filtrado de mediana de la imagen A 3D en tres dimensiones. Cada vóxel de salida en B contiene el valor de la mediana del entorno de m por n por p alrededor del vóxel correspondiente en A.

B = medfilt3(___,padopt) controla cómo medfilt3 rellena los límites del arreglo.

Ejemplos

Utilizar el filtrado de mediana para eliminar valores atípicos en datos 3D

Cree una superficie 3D con ruido.

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

Aplique el filtrado de mediana.

filteredV = medfilt3(noisyV);

Muestre las superficies con ruido y filtradas juntas.

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

Argumentos de entrada

A — Imagen de entrada

Arreglo numérico 3D | Arreglo lógico 3D

Imagen de entrada, especificada como arreglo numérico o lógico 3D.

Tipos de datos: single | double | int8 | int16 | int32 | int64 | uint8 | uint16 | uint32 | uint64 | logical

[m n p] — Tamaño del entorno

[3 3 3] (predeterminado) | Vector de 3 elementos

Tamaño del entorno, especificado como vector de 3 elementos de enteros positivos impares.

Tipos de datos: single | double | int8 | int16 | int32 | uint8 | uint16 | uint32

padopt — Opción de relleno

"symmetric" (predeterminado) | "zeros" | "replicate"

Opción de relleno, especificada como uno de los siguientes valores:

Valor	Descripción	Ejemplo
"symmetric"	Rellenar el arreglo con reflejos de espejo de sí mismo	$$\left[\begin{array}{ccc}3& 1& 4\\ 1& 5& 9\\ 2& 6& 5\end{array}\right]\to \left[\begin{array}{ccccccc}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{array}\right]$$
"replicate"	Rellenar el arreglo repitiendo los elementos del borde	$$\left[\begin{array}{ccc}3& 1& 4\\ 1& 5& 9\\ 2& 6& 5\end{array}\right]\to \left[\begin{array}{ccccccc}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{array}\right]$$
"zeros"	Rellenar el arreglo con el valor 0	$$\left[\begin{array}{ccc}3& 1& 4\\ 1& 5& 9\\ 2& 6& 5\end{array}\right]\to \left[\begin{array}{ccccccc}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{array}\right]$$

Tipos de datos: char | string

Argumentos de salida

B — Imagen de salida

Arreglo numérico 3D

Imagen de salida, devuelta como arreglo numérico 3D del mismo tamaño y tipo de datos que la imagen de entrada A.

Capacidades ampliadas

Entorno basado en subprocesos

Ejecute código en segundo plano con MATLAB® backgroundPool o acelere código con Parallel Computing Toolbox™ ThreadPool.

Esta función es totalmente compatible con los entornos basados en hilos. Para obtener más información, consulte Ejecutar funciones de MATLAB en entornos basados en subprocesos.

Historial de versiones

Introducido en R2016b

R2022b: Compatibilidad con entornos basados en hilos

Ahora medfilt3 es compatible con los entornos basados en hilos.

R2021b: Rendimiento mejorado para tamaños de entorno pequeños

La función medfilt3 muestra un rendimiento mejorado para los tamaños de entorno desde [3, 3, 3] hasta [31, 31, 31].

Por ejemplo, este código es aproximadamente 3 veces más rápido que el de la versión anterior.

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

Los tiempos de ejecución aproximados son:

R2021a: 0,24 s

R2021b: 0,08 s

El código se ha cronometrado en un sistema de prueba Intel^® Xeon^® Gold 5220 CPU @ 2.2 GHz con Windows^® 10, llamando a la función timingTestMedfilt3.

Consulte también

medfilt2

Temas

• Opciones de relleno de límites de imfilter