Articulo de referencia

Ranura de caja

En los campos matemáticos del análisis numérico y la teoría de la aproximación , las splines de caja son funciones polinómicas por partes de varias variables. [ 1 ] Las splines ...

En los campos matemáticos del análisis numérico y la teoría de la aproximación , las splines de caja son funciones polinómicas por partes de varias variables. [ 1 ] Las splines de caja se consideran una generalización multivariada de las splines de base (B-splines) y se utilizan generalmente para la aproximación/interpolación multivariada. Geométricamente, una spline de caja es la sombra (rayos X) de un hipercubo proyectado a un espacio de menor dimensión. [ 2 ] Las splines de caja y las splines simplex son casos especiales bien estudiados de splines poliédricas que se definen como sombras de politopos generales .

Definición

Una spline de caja es una función multivariada.(RdR){\displaystyle (\mathbb {R} ^{d}\to \mathbb {R} )}definido para un conjunto de vectores,ξRd,{\displaystyle \xi \in \mathbb {R} ^{d},}generalmente reunidos en una matrizΞ:=[ξ1ξnorte].{\displaystyle \mathbf {\Xi } :=\left[\xi _{1}\dots \xi _{N}\right].}

Cuando el número de vectores es igual a la dimensión del dominio (es decir,norte=d{\displaystyle N=d}) entonces la spline de caja es simplemente la función indicadora (normalizada) del paralelepípedo formado por los vectores enΞ{\displaystyle \mathbf {\Xi } }:

METROΞ(incógnita):=1detΞχΞ(incógnita)={1detΞsi incógnita=norte=1dtnorteξnorte para algunos 0tnorte<10de lo contrario.{\displaystyle M_{\mathbf {\Xi } }(\mathbf {x} ):={\frac {1}{\mid {\det {\Xi }}\mid }}\chi _{\mathbf {\Xi } }(\mathbf {x} )={\begin{cases}{\frac {1}{\mid {\det {\Xi }}\mid }}&{\text{if }}\mathbf {x} =\sum _{n=1}^{d}{t_{n}\xi _{n}}{\text{ for some }}0\leq t_{n}<1\\0&{\text{otherwise}}.\end{cases}}}

Añadiendo una nueva dirección,ξ,{\displaystyle \xi ,}aΞ,{\displaystyle \mathbf {\Xi } ,}o generalmente cuandonorte>d,{\displaystyle N>d,}La spline de caja se define recursivamente: [ 1 ]

METROΞξ(incógnita)=01METROΞ(incógnitatξ)dt.{\displaystyle M_{\mathbf {\Xi } \cup \xi }(\mathbf {x} )=\int _{0}^{1}M_{\mathbf {\Xi } }(\mathbf {x} -t\xi )\,{\rm {d}}t.}
Ejemplos de splines de caja bivariados correspondientes a 1, 2, 3 y 4 vectores en 2-D.

La ranura de la cajaMETROΞ{\displaystyle M_{\mathbf {\Xi } }}puede interpretarse como la sombra de la función indicadora del hipercubo unitario enRnorte{\displaystyle \mathbb {R} ^{N}}cuando se proyecta hacia abajo enRd.{\displaystyle \mathbb {R} ^{d}.}Desde esta perspectiva, los vectoresξΞ{\displaystyle \xi \in \mathbf {\Xi } }son la proyección geométrica de la base estándar enRnorte{\displaystyle \mathbb {R} ^{N}}(es decir, las aristas del hipercubo) aRd.{\displaystyle \mathbb {R} ^{d}.}

Considerando distribuciones templadas, una spline de caja asociada con un vector de dirección única es una función generalizada de tipo Dirac soportada entξ{\displaystyle t\xi }para0t<1{\displaystyle 0\leq t<1}. Entonces, la spline de caja general se define como la convolución de las distribuciones asociadas a las splines de caja de vector único:

METROΞ=METROξ1METROξ2METROξnorte.{\displaystyle M_{\mathbf {\Xi } }=M_{\xi _{1}}\ast M_{\xi _{2}}\ast \dots \ast M_{\xi _{N}}.}

Propiedades

  • Dejarκ{\displaystyle \kappa }ser el número mínimo de direcciones cuya eliminación deΞ{\displaystyle \Xi }hace que las instrucciones restantes no se extiendanRd{\displaystyle \mathbb {R} ^{d}}. Entonces la spline de caja tieneκ2{\displaystyle \kappa -2}grados de continuidad:METROΞdoκ2(Rd){\displaystyle M_{\mathbf {\Xi } }\in C^{\kappa -2}(\mathbb {R} ^{d})}. [ 1 ]
  • Cuandonorted{\displaystyle N\geq d}(y vectores enΞ{\displaystyle \Xi }durarRd{\displaystyle \mathbb {R} ^{d}}) la spline de caja es una función con soporte compacto cuyo soporte es un zonotopo enRd{\displaystyle \mathbb {R} ^{d}}formada por la suma de Minkowski de los vectores de direcciónξΞ{\displaystyle \xi \in \mathbf {\Xi } }.
  • Dado que los zonotopos son simétricos respecto al centro, el soporte de la spline de caja es simétrico con respecto a su centro:doΞ:=12norte=1norteξnorte.{\displaystyle \mathbf {c} _{\Xi }:={\frac {1}{2}}\sum _{n=1}^{N}\xi _{n}.}
  • Transformada de Fourier de la spline de caja, end{\displaystyle d}dimensiones, está dado por
METRO^Ξ(ω)=exp(jdoΞω)norte=1nortedesde(ξnorteω).{\displaystyle {\hat {M}}_{\Xi }(\omega )=\exp(-j\mathbf {c} _{\Xi }\cdot \omega )\prod _{n=1}^{N}\operatorname {sinc} (\xi _{n}\cdot \omega ).}

Aplicaciones

Para diversas aplicaciones, se utilizan combinaciones lineales de desplazamientos de una o más splines de caja sobre una malla. Estas splines son más eficientes que las combinaciones lineales de splines simplex, ya que son refinables y, por definición, invariantes a los desplazamientos. Por lo tanto, constituyen el punto de partida para muchas construcciones de superficies de subdivisión .

Las splines de caja han sido útiles en la caracterización de arreglos de hiperplanos . [ 3 ] Además, las splines de caja se pueden usar para calcular el volumen de politopos. [ 4 ]

En el contexto del procesamiento de señales multidimensionales , las splines de caja pueden proporcionar núcleos de interpolación multivariados (filtros de reconstrucción) adaptados a redes de muestreo no cartesianas , [ 5 ] y redes cristalográficas (redes raíz) que incluyen muchas redes de muestreo óptimas desde el punto de vista de la teoría de la información. [ 6 ] Generalmente, las redes de empaquetamiento de esferas óptimas y las redes de recubrimiento de esferas [ 7 ] son ​​útiles para muestrear funciones multivariadas en 2-D, 3-D y dimensiones superiores. [ 8 ] En el entorno 2-D, la spline de caja de tres direcciones [ 9 ] se utiliza para la interpolación de imágenes muestreadas hexagonalmente. En el entorno 3-D, las splines de caja de cuatro direcciones [ 10 ] y seis direcciones [ 11 ] se utilizan para la interpolación de datos muestreados en las redes cúbicas centradas en el cuerpo y cúbicas centradas en las caras (óptimas) respectivamente. [ 5 ] La spline de caja de siete direcciones [ 12 ] se ha utilizado para modelar superficies y puede utilizarse para interpolación de datos en la red cartesiana [ 13 ] así como en la red cúbica centrada en el cuerpo . [ 14 ] La generalización de las splines de caja de cuatro [ 10 ] y seis direcciones [ 11 ] a dimensiones superiores [ 15 ] puede utilizarse para construir splines en redes raíz . [ 16 ] Las splines de caja son ingredientes clave de las splines hexagonales [ 17 ] y las splines de Voronoi [ 18 ] que, sin embargo, no son refinables.

Las splines de caja han encontrado aplicaciones en el filtrado de alta dimensión, específicamente para el filtrado bilateral rápido y los algoritmos de medias no locales . [ 19 ] Además, las splines de caja se utilizan para diseñar filtros eficientes que varían espacialmente (es decir, no convolucionales). [ 20 ]

Las splines de caja son funciones base útiles para la representación de imágenes en el contexto de problemas de reconstrucción tomográfica , ya que los espacios de splines generados por los espacios de splines de caja son cerrados bajo transformadas de rayos X y de Radon . [ 21 ] [ 22 ] En esta aplicación, mientras que la señal se representa en espacios invariantes a traslaciones, las proyecciones se obtienen, en forma cerrada, mediante traslaciones no uniformes de splines de caja. [ 21 ]

En el contexto del procesamiento de imágenes, se ha demostrado que los marcos de spline de caja son eficaces en la detección de bordes . [ 23 ]

Referencias

  1. 1 2 3 Aburrido, C.; Hollig, K.; Riemenschneider, S. (1993). Splines de caja . Ciencias Matemáticas Aplicadas. vol.  98.doi : 10.1007 /978-1-4757-2244-4 . ISBN 978-1-4419-2834-4.
  2. Prautzsch, H.; Boehm, W.; Paluszny, M. (2002). "Box splines". Técnicas de Bézier y B-Spline . Matemáticas y visualización. págs. 239–258 . doi : 10.1007/978-3-662-04919-8_17 . ISBN  978-3-642-07842-2.
  3. De Concini, C.; Procesi, C. (2010). Temas en arreglos hiperplanos, politopos y splines de caja . doi : 10.1007/978-0-387-78963-7 . ISBN 978-0-387-78962-0.
  4. Xu, Z. (2011). "Splines multivariados y politopos". Journal of Approximation Theory . 163 (3): 377– 387. arXiv : 0806.1127 . doi : 10.1016/j.jat.2010.10.005 . S2CID 10063913 . 
  5. 1 2 Entezari, Alireza. Retículos de muestreo óptimos y splines de caja trivariados. [Vancouver, BC.]: Universidad Simon Fraser, 2007. < http://summit.sfu.ca/item/8178 >.
  6. Kunsch, HR; Agrell, E.; Hamprecht, FA (2005). "Optimal Lattices for Sampling" . IEEE Transactions on Information Theory . 51 (2): 634. doi : 10.1109/TIT.2004.840864 . S2CID 16942177 . 
  7. JH Conway, NJA Sloane. Empaquetamientos de esferas, retículos y grupos. Springer, 1999.
  8. Petersen, DP; Middleton, D. (1962). "Muestreo y reconstrucción de funciones limitadas por número de onda en espacios euclidianos N-dimensionales" . Information and Control . 5 (4): 279. doi : 10.1016/S0019-9958(62)90633-2 .
  9. Condat, L.; Van De Ville, D. (2006). "Splines de caja tridimensionales: Caracterización y evaluación eficiente" (PDF) . IEEE Signal Processing Letters . 13 (7): 417. Bibcode : 2006ISPL...13..417C . doi : 10.1109/LSP.2006.871852 . S2CID 9023102 . 
  10. 1 2 Entezari, A.; Van De Ville, D.; Moller, T. (2008). "Splines de caja prácticos para reconstrucción en la red cúbica centrada en el cuerpo" ( PDF) . IEEE Transactions on Visualization and Computer Graphics . 14 (2): 313– 328. doi : 10.1109/TVCG.2007.70429 . PMID 18192712. S2CID 6395127 .  
  11. 1 2 Minho Kim, M.; Entezari, A.; Peters, Jorg (2008). "Reconstrucción de Box Spline en la red cúbica centrada en las caras". IEEE Transactions on Visualization and Computer Graphics . 14 (6): 1523– 1530. CiteSeerX 10.1.1.216.408 . doi : 10.1109/TVCG.2008.115 . PMID 18989005 . S2CID 194024 .   
  12. Peters, Jorg; Wittman, M. (1997). "Mezclas CSG basadas en splines de caja" . Actas del cuarto simposio ACM sobre modelado sólido y aplicaciones - SMA '97 . pp. 195. doi : 10.1145/267734.267783 . ISBN  0897919467. S2CID 10064302 . 
  13. Entezari, A.; Moller, T. (2006). "Extensiones de la spline de caja de Zwart-Powell para la reconstrucción de datos volumétricos en la red cartesiana". IEEE Transactions on Visualization and Computer Graphics . 12 (5): 1337– 1344. doi : 10.1109/TVCG.2006.141 . PMID 17080870 . S2CID 232110 .  
  14. Minho Kim (2013). "Reconstrucción de Box-Spline cuártica en la red BCC". IEEE Transactions on Visualization and Computer Graphics . 19 (2): 319– 330. doi : 10.1109/TVCG.2012.130 . PMID 22614329 . S2CID 7338997 .  
  15. Kim, Minho. Box-Splines simétricos en redes de raíces. [Gainesville, Florida]: Universidad de Florida, 2008. < http://uf.catalog.fcla.edu/permalink.jsp?20UF021643670 >.
  16. Kim, M.; Peters, Jorg (2011). "Splines de caja simétricas en retículos de raíz" . Journal of Computational and Applied Mathematics . 235 (14): 3972. doi : 10.1016/j.cam.2010.11.027 .
  17. Van De Ville, D.; Blu, T.; Unser, M.; Philips, W.; Lemahieu, I.; Van De Walle, R. (2004). "Hex-Splines: Una nueva familia de splines para redes hexagonales" (PDF) . IEEE Transactions on Image Processing . 13 (6): 758– 772. Bibcode : 2004ITIP...13..758V . doi : 10.1109/TIP.2004.827231 . PMID 15648867. S2CID 9832708 .  
  18. Mirzargar, M.; Entezari, A. (2010). "Voronoi Splines". IEEE Transactions on Signal Processing . 58 (9): 4572. Bibcode : 2010ITSP...58.4572M . doi : 10.1109/TSP.2010.2051808 . S2CID 9712416 . 
  19. Baek, J.; Adams, A.; Dolson, J. (2012). "Filtrado gaussiano de alta dimensión basado en retículos y el retículo permutoédrico". Journal of Mathematical Imaging and Vision . 46 (2): 211. doi : 10.1007/s10851-012-0379-2 . hdl : 1721.1/105344 . S2CID 16576761 . 
  20. Chaudhury, KN; MuñOz-Barrutia, A. ; Unser, M. (2010). "Filtrado elíptico rápido con variación espacial mediante splines de caja". IEEE Transactions on Image Processing . 19 (9): 2290– 2306. arXiv : 1003.2022 . Bibcode : 2010ITIP...19.2290C . doi : 10.1109/TIP.2010.2046953 . PMID 20350851 . S2CID 16383503 .  
  21. 1 2 Entezari, A.; Nilchian, M.; Unser, M. (2012). "Un cálculo de splines de caja para la discretización de problemas de reconstrucción de tomografía computarizada" ( PDF) . IEEE Transactions on Medical Imaging . 31 (8): 1532– 1541. doi : 10.1109/TMI.2012.2191417 . PMID 22453611. S2CID 3787118 .  
  22. Entezari, A.; Unser, M. (2010). "Un cálculo de splines de caja para tomografía computarizada". Simposio Internacional IEEE de 2010 sobre Imágenes Biomédicas: De lo Nano a lo Macro . pág. 600. doi : 10.1109/ISBI.2010.5490105 . ISBN  978-1-4244-4125-9. S2CID 17368057 . 
  23. Guo, W.; Lai, MJ (2013). "Box Spline Wavelet Frames for Image Edge Analysis" . SIAM Journal on Imaging Sciences . 6 (3): 1553. doi : 10.1137/120881348 . S2CID 2708993 .