Articulo de referencia

Filtro bilateral

Izquierda: imagen original. Derecha: imagen procesada con filtro bilateral. Un filtro bilateral es un filtro de suavizado no lineal , que preserva los bordes y reduce el ruido e...

Izquierda: imagen original. Derecha: imagen procesada con filtro bilateral.

Un filtro bilateral es un filtro de suavizado no lineal , que preserva los bordes y reduce el ruido en las imágenes . Reemplaza la intensidad de cada píxel con un promedio ponderado de los valores de intensidad de los píxeles cercanos. Este peso puede basarse en una distribución gaussiana. Es fundamental que los pesos dependan no solo de la distancia euclidiana entre los píxeles, sino también de las diferencias radiométricas (por ejemplo, variaciones en la intensidad del color o la profundidad). Esta doble dependencia preserva los bordes nítidos a la vez que suprime el ruido.

Definición

El filtro bilateral se define como [ 1 ] [ 2 ]

Ifiltrado(incógnita)=1WpagincógnitaiΩI(incógnitai)Fr(I(incógnitai)I(incógnita))gramos(incógnitaiincógnita),{\displaystyle I^{\text{filtrado}}(x)={\frac {1}{W_{p}}}\sum _{x_{i}\in \Omega }I(x_{i})f_{r}(\|I(x_{i})-I(x)\|)g_{s}(\|x_{i}-x\|),}

y término de normalización,Wpag{\displaystyle {W_{p}}}, se define como

Wpag=incógnitaiΩFr(I(incógnitai)I(incógnita))gramos(incógnitaiincógnita){\displaystyle W_{p}=\sum _{x_{i}\in \Omega }{f_{r}(\|I(x_{i})-I(x)\|)g_{s}(\|x_{i}-x\|)}}

dónde

Ifiltrado{\displaystyle I^{\text{filtrado}}}es la imagen filtrada;
I{\displaystyle I}es la imagen de entrada original que se va a filtrar;
incógnita{\displaystyle x}son las coordenadas del píxel actual que se va a filtrar;
Ω{\displaystyle \Omega }¿Está la ventana centrada en?incógnita{\displaystyle x}, entoncesincógnitaiΩ{\displaystyle x_{i}\in \Omega }es otro píxel;
Fr{\displaystyle f_{r}}es el núcleo de rango para suavizar las diferencias de intensidades (esta función puede ser una función gaussiana );
gramos{\displaystyle g_{s}}es el núcleo espacial (o de dominio) para suavizar las diferencias en las coordenadas (esta función puede ser una función gaussiana).

El pesoWpag{\displaystyle W_{p}}se asigna utilizando la cercanía espacial (utilizando el núcleo espacial)gramos{\displaystyle g_{s}}) y la diferencia de intensidad (usando el kernel de rango)Fr{\displaystyle f_{r}}). [ 2 ] Consideremos un píxel ubicado en(i,j){\displaystyle (i,j)}que necesita ser desruido en la imagen usando sus píxeles vecinos y uno de sus píxeles vecinos está ubicado en(k,l){\displaystyle (k,l)}. Luego, suponiendo que los núcleos de rango y espaciales son núcleos gaussianos , el peso asignado para el píxel(k,l){\displaystyle (k,l)}para eliminar el ruido del píxel(i,j){\displaystyle (i,j)}es dado por

w(i,j,k,l)=exp((ik)2+(jl)22σd2I(i,j)I(k,l)22σr2),{\displaystyle w(i,j,k,l)=\exp \left(-{\frac {(ik)^{2}+(jl)^{2}}{2\sigma _{d}^{2}}}-{\frac {\|I(i,j)-I(k,l)\|^{2}}{2\sigma _{r}^{2}}}\right),}

dóndeσd{\displaystyle \sigma _{d}}yσr{\displaystyle \sigma _{r}}son parámetros de suavizado yI(i,j){\displaystyle I(i,j)}yI(k,l){\displaystyle I(k,l)}son la intensidad de los píxeles(i,j){\displaystyle (i,j)}y(k,l){\displaystyle (k,l)}respectivamente.

Después de calcular los pesos, normalícelos:

ID(i,j)=k,lI(k,l)w(i,j,k,l)k,lw(i,j,k,l),{\displaystyle I_{D}(i,j)={\frac {\sum _{k,l}I(k,l)w(i,j,k,l)}{\sum _{k,l}w(i,j,k,l)}},}

dóndeID{\displaystyle I_{D}}es la intensidad sin ruido del píxel(i,j){\displaystyle (i,j)}.

Parámetros

  • A medida que aumenta el parámetro de rango σ r , el filtro bilateral se vuelve gradualmente equivalente al desenfoque gaussiano ; es decir, σ d se vuelve dominante, porque el núcleo gaussiano es entonces mayormente plano.
  • A medida que aumenta el parámetro espacial σ d , las características más grandes se suavizan.

Limitaciones

El filtro bilateral en su forma directa puede introducir varios tipos de artefactos en la imagen:

  • Efecto de escalera: mesetas de intensidad que hacen que las imágenes parezcan dibujos animados [ 3 ].
  • Inversión de gradiente: introducción de bordes falsos en la imagen. [ 4 ]

Existen varias extensiones del filtro que abordan estos artefactos, como el filtro bilateral escalado que utiliza una imagen reducida para calcular los pesos. [ 5 ] También se han propuesto filtros alternativos, como el filtro guiado , [ 6 ] como una alternativa eficiente sin estas limitaciones.

Implementaciones

Adobe Photoshop implementa un filtro bilateral en su herramienta de desenfoque de superficie .

GIMP implementa un filtro bilateral en sus herramientas Filtros → Desenfoque ; y se llama Desenfoque gaussiano selectivo . El complemento gratuito G'MIC Reparación → Suavizado [bilateral] para GIMP añade más control. [ 7 ] Un truco sencillo para implementar eficientemente un filtro bilateral es aprovechar el submuestreo de disco de Poisson . [ 1 ]

OpenCV implementa la función: .bilateralFilter( source, destination, Ω{\displaystyle \Omega }, fr{\displaystyle f_{r}}, gs{\displaystyle g_{s}} )

Se ha demostrado que el filtro bilateral es una aplicación del núcleo de tiempo corto del flujo de Beltrami [ 8 ] [ 9 ] [ 10 ] que se introdujo como un mecanismo de suavizado selectivo que preserva los bordes antes del filtro bilateral.

Otros filtros de suavizado que preservan los bordes incluyen: difusión anisotrópica , [ 11 ] mínimos cuadrados ponderados, [ 12 ] ondículas que evitan los bordes, [ 13 ] edición geodésica, [ 14 ] filtrado guiado, [ 15 ] y transformaciones de dominio. [ 16 ]

Véase también

Referencias

  1. 1 2 Banterle, F.; Corsini, M.; Cignoni, P.; Scopigno, R. (2011). "Un filtro bilateral rápido, sencillo y de baja memoria mediante submuestreo en el dominio espacial" . Computer Graphics Forum . 31 (1): 19– 32. doi : 10.1111/j.1467-8659.2011.02078.x . hdl : 11380/1170184 . S2CID 18288647 . 
  2. 1 2 Tomasi, C; Manduchi, R (1998). Filtrado bilateral para imágenes en escala de grises y en color (PDF) . Sexta Conferencia Internacional sobre Visión por Computadora. Bombay. pp. 839–846 . doi : 10.1109/ICCV.1998.710815 . Archivado del original (PDF) el 29 de agosto de 2008. Recuperado el 18 de noviembre de 2009 . 
  3. Kornprobst, Pierre (2007). "¿Limitaciones? - Una introducción sencilla al filtrado bilateral y sus aplicaciones" (PDF) . Recuperado el 7 de mayo de 2017 .
  4. He, Kaiming; Sun, Jian; Tang, Xiaoou. "Filtrado de imágenes guiado" (PDF) . Archivado del original (PDF) el 9 de diciembre de 2017. Recuperado el 7 de mayo de 2017 .
  5. Aswatha, Shashaank M.; Mukhopadhyay, Jayanta; Bhowmick, Partha (diciembre de 2011). "Reducción de ruido en imágenes mediante filtrado bilateral escalado". Tercera Conferencia Nacional de 2011 sobre Visión por Computadora, Reconocimiento de Patrones, Procesamiento de Imágenes y Gráficos . págs. 122–125 . doi : 10.1109/NCVPRIPG.2011.33 . ISBN  978-1-4577-2102-1. S2CID 25738863 . 
  6. He, Kaiming. "Filtrado de imágenes guiado" . Archivado del original el 19 de diciembre de 2017. Recuperado el 7 de mayo de 2017 .
  7. "G'MIC - La magia de GREYC para el procesamiento de imágenes: un marco de trabajo abierto y completo para el procesamiento de imágenes" . gmic.eu. Archivado del original el 27 de diciembre de 2014.
  8. R. Kimmel, R. Malladi y N. Sochen. Imágenes como mapas de incrustación y superficies mínimas: Películas, color e imágenes médicas volumétricas. IEEE CVPR'97, págs. 350-355, Puerto Rico, 17-19 de junio de 1997. https://www.cs.technion.ac.il/~ron/PAPERS/cvpr97.pdf
  9. R. Kimmel, R. Malladi y N. Sochen. Imágenes como mapas incrustados y superficies mínimas: películas, color, textura e imágenes médicas volumétricas. International Journal of Computer Vision, 39(2):111-129, septiembre de 2000. https://www.cs.technion.ac.il/~ron/PAPERS/KimMalSoc_IJCV2000.pdf
  10. N. Sochen, R. Kimmel y A. M. Bruckstein. Difusiones y confusiones en el procesamiento de señales e imágenes, Journal of Mathematical Imaging and Vision, 14(3):195-209, 2001. https://www.cs.technion.ac.il/~ron/PAPERS/SocKimBru_JMIV2001.pdf
  11. Durand, Frédo y Dorsey, Julie. «Filtrado bilateral rápido para la visualización de imágenes de alto rango dinámico». En ACM Transactions on Graphics, volumen 21, número 3, págs. 257-266. https://people.csail.mit.edu/fredo/PUBLI/Siggraph2002/DurandBilateral.pdf
  12. Farbman, Zeev, Raanan Fattal, Dani Lischinski y Richard Szeliski. «Descomposiciones que preservan los bordes para la manipulación de tonos y detalles a múltiples escalas». En ACM Transactions on Graphics, vol. 27, n.º 3 (2008): 67. http://www.cs.huji.ac.il/~danix/epd/
  13. Fattal, Raanan. "Ondículas que evitan bordes y sus aplicaciones." En ACM Transactions on Graphics vol. 28, n.º 3 (2009): 22. http://www.cs.huji.ac.il/~raananf/projects/eaw/
  14. Criminisi, Antonio, Toby Sharp, Carsten Rother y Patrick Pérez. «Edición geodésica de imágenes y vídeo». En ACM Transactions on Graphics (TOG), vol. 29, n.º 5 (2010): 134. http://research.microsoft.com/apps/pubs/default.aspx?id=81528
  15. He, Kaiming, Jian Sun y Xiaoou Tang. "Filtrado de imágenes guiado". En Computer Vision–ECCV 2010, págs. 1-14. Springer Berlin Heidelberg, 2010. http://kaiminghe.com/eccv10/index.html Archivado el 19/12/2017 en Wayback Machine
  16. Gastal, Eduardo SL y Manuel M. Oliveira. "Transformación de dominio para el procesamiento de imágenes y vídeo con detección de bordes". En ACM Transactions on Graphics, vol. 30, n.º 4 (2011): 69. http://inf.ufrgs.br/~eslgastal/DomainTransform/