Digamos que tengo una función compleja (por ejemplo, una imagen de resonancia magnética) que tiene una magnitud constante casi por piezas, pero una fase no constante.
Si tengo un problema de optimización para encontrar y configurar una función objetivo con un término de variación total (por ejemplo, para detección de ruido o compresión), generalmente tiene la siguiente forma:
Sin embargo, dado que supongo que tiene una magnitud constante por partes, creo que podría ser mejor usar:
Sin embargo, para un solucionador basado en gradiente, uno debería conocer el gradiente de obj2. El gradiente para es: . ¿Cuál es el gradiente de ?
Actualizar:
Intuitivamente, supondría algo como lo siguiente (dado que la fase no tiene influencia en , deje la fase "intacta"):
Sin embargo, mi conocimiento en análisis complejo es muy limitado y no estoy seguro de si esto tiene sentido.
Respuestas:
El problema conEl | FEl | es que, dado que no es analítico, no se aplica la definición estándar de derivada compleja. Una solución es usar derivados de Wirtinger:
http://en.wikipedia.org/wiki/Wirtinger_derivatives
Una cuenta detallada del cálculo de Wirtinger para problemas de procesamiento de señal es
http://arxiv.org/abs/0906.4835
Otra opción (probablemente más simple) es tratar la imagen compleja como una imagen de dos canales (real, imag) y usar la definición de derivada para campos vectoriales. Este documento tiene una explicación muy clara sobre cómo hacer esto:
Lee, H.-C .; Cok, DR; "Detección de límites en un campo vectorial" (IEEE Transactions on Signal Processing, vol.39, no.5, pp.1181-1194, mayo de 1991)
fuente