Clase 1 – 5-Septiembre-2014

[Link]

[Link]

James E. Holden. Lectures on Image Science, Medical Physics, University of Wisconsin.

Arfken, Weber. Essential Mathematical Methods for Physicists

Barret- Myers, “Foundations of Imaging Science”

Steven Smith, “The Scientist and Engineers Guide to Digital Signal Processing”

*Feeman T.G. ‘’The mathematics of medical imaging. A beginners guide’’

*Epstein, Charles. ‘’The Mathematics of Medical Imaging’’

Contenido:

1.- Introducción

2.- Math Background: El plano complejo, funciones pares e impares. La ‘’función’’ Delta de
Dirac.

1.- Introducción

Encontramos imágenes en todos lados en el mundo actual. Ya sea para mantenerse

informados, entretenimiento, comunicaciones, así como en medicina, ciencia y tecnología.

Varios ejemplos:

(a) imágenes médicas son vitales para el diagnóstico y tratamiento en enfermedades.

(b) imágenes satelitales nos dan información del clima

1
(c) en investigación biomédicas y ciencia de los materiales se usan imágenes microscópicas
de muchos tipos.

(d) reservas de petróleo usualmente son encontradas a través de imágenes sísmicas.

(e) en sistemas de búsquedas submarina son usadas imágenes sónicas.

DSP >> Digital Signal Processing (Ref: Steven Smith, “The Scientist and Engineers Guide to Digital Signal
Processing”

2
Steven Smith, “The Scientist and Engineers Guide to Digital Signal Processing”

3
4
Estos sistemas de imágenes que hemos mencionado incluyen lo último en alta tecnología,
sin embargo no resulta tan obvio pensar en los fundamentos teóricos en los cuales
descansan estas aplicaciones y que a su vez relacionan las diferentes tecnologías
permitiendo un diseño sistemático y optimización de los diferentes sistemas de imágenes.

Existen subdisciplinas de lo que denominaremos “ciencia de la imagen” que incluyen:

-óptica cuántica

-óptica de rayos

-propagación de ondas

-procesamiento e interpretación de imágenes

Para que la ciencia de las imágenes, sea considerada ciencia, se deben cumplir tres
elementos:

- el campo debe tener un lenguaje común, estar de acuerdo sobre un conjunto de

definiciones.

- debe tener un conjunto de procedimientos experimentales aceptados.

- y finalmente, debe tener una teoría con valor para predecir.

Además, es necesario recordar que el objetivo de los profesionales es obtener (y hacer)

cada vez mejores imágenes

Según Barrett-Myers: “la calidad de la imagen puede ser definida objetivamente, medida
experimentalmente y predicha y optimizada teóricamente”

5
En este curso trabajaremos los fundamentos matemáticos (en primer lugar, en los aspectos
determinísticos, y luego los estadísticos ) detrás de las imágenes.

La imagen de entrada o de “input” es el ‘’objeto’’ de estudio o escaneado

La imagen de salida o de “output” usualmente es la imagen procesada para mejorar la

visualización o para cuantificar las características del objeto.

6
2.- Los Números Complejos

Ref: Complex Variables- Series Schaum // Complex Variables and Applications – Brown, Churchill

Como sabemos el conjunto de números complejos surgen de la necesidad de dar solución a

la ecuación polinómica:

x2 + 1 = 0
Se definen como:
z ≡ a + ib
Donde a, b ∈ ℝ (son números reales) y la denominada unidad imaginaria i que tiene la
propiedad

i 2 = −1
Entonces en un número complejo z = a + ib , al número a se le denomina parte real de z ,
Re{z } y al número b , parte imaginaria de z , Im{z } . En general, z puede designar a
cualquier número complejo por lo que puede llamársele variable compleja.

Propiedades de los Números Complejos

(a) Dos números complejos z 1 = a + ib y z 2 = c + id son iguales si y solo si a = c y b = d .

(b) Consideremos un número complejo, z = a + ib . Si b = 0 , entonces z = a + i 0 , es un

número real. Esto es, podemos considerar a los números reales como un subconjunto de los
números complejos con b = 0 .

(c) Cuando en un número complejo, z = a + ib , a = 0 , esto es, z = 0 + ib = +ib , se habla

de un número complejo imaginario puro o puramente imaginario.

7
(d) El complejo conjugado (o el conjugado) de un número complejo z = a + ib , se define

como z = a − ib , y se designa por z o z * .

(e) Suma de números complejos:

z 1 + z 2 = (a + ib) + (c + id ) = (a + c) + i(b + d )

(f) Resta o sustracción de números complejos

z 1 − z 2 = (a + ib) − (c + id ) = (a − c) + i(b − d )

(g) Producto o multiplicación de números complejos

z 1 ⋅ z 2 = z 1z 2 = (a + ib)(c + id ) = ac + iad + ibc + (ib)(id ) =

z1 ⋅ z 2 = ac + iad + ibc + i 2 (bd ) = (ac − bd ) + i(ad + bc)

(h) División de números complejos

z1 (a + ib )
= , la cual para que esté bien definida ⇒ z 2 ≠ 0 + 0i, ó c, d ≠ 0
z2 (c + id )

Entonces,

z1 (a + ib ) (c − id ) ac − iad + icb − i 2bd (ac + bd ) + i(−ad + cb )

= = =
z 2 (c + id ) (c − id ) c 2 − icd + icd − i 2d 2 c2 + d 2

z1 (ac + bd ) (bc − ad )
= + i
z2 c2 + d 2 c2 + d 2

(i) El módulo o valor absoluto de un número complejo z = a + ib se define como

z = a + ib = a 2 + b 2

(j) Si z 1, z 2 ,..., z n son números complejos, se verifican las siguientes propiedades:

8
(*) z 1z 2 = z 1 z 2 o z 1z 2 ...z n = z 1 z 2 ... z n

z1 z1
(**) = , z2 ≠ 0
z2 z2

(***) z 1 + z 2 ≤ z 1 + z 2

(*v) z 1 ± z 2 ≥ z 1 − z 2

(k) si z = a − ib entonces , z = z y z = z para todo z .

(l) Dado z = a − ib , entonces zz = z

Representación gráfica de los números complejos

Un número complejo x + iy se puede considerar como un par ordenado de números reales

que puede representarse como un punto en un plano xy al cual se le denomina Plano
Complejo o diagrama de Argand. A este plano se le denomina también z -plano.

A cada número complejo le corresponde uno y solo un punto en el plano, y al revés, a cada
punto en el plano le corresponde uno y solo un número complejo. Al eje X se le denomina
eje real, y al eje Y , eje imaginario.

9
La distancia entre dos puntos z 1 = x 1 + iy1 y z 2 = x 2 + iy2 en el plano complejo está dada
por:

z1 − z 2 = (x 1 − x 2 )2 + (y1 − y2 )2 , el cual es un número positivo. Es fácil ver que la última

ecuación se cumple cuando escribimos al número complejo como un vector:

Ya que geométricamente, z = x 2 + y 2 es la distancia entre el punto (x , y ) y el origen

del plano.

10
Para el complejo conjugado, en el plano complejo, representa una reflexión con respecto al
eje real.

Representación gráfica de los números complejos

La distancia entre dos puntos z 1 = x 1 + iy1 y z 2 = x 2 + iy2 en el plano complejo está dada
por:

z1 − z 2 = (x 1 − x 2 )2 + (y1 − y2 )2 , el cual es un número positivo.

Es fácil ver que la última ecuación se cumple cuando escribimos al número complejo como
un vector:

11
Ya que geométricamente, z = x 2 + y 2 es la distancia entre el punto (x , y ) y el origen
del plano.

Forma Polar – Forma exponencial de los números complejos.

(1) Sea un número complejo z = x + iy , representado en la figura. Las coordenadas

polares del punto (x , y ) son (r , ϕ ) , donde r = z . La relación entre ambas

representaciones:

x = r cos ϕ , y = rsenϕ .

12
⇒ z = x + iy → z = r cos ϕ + i rsenϕ = r (cos ϕ + i senϕ )

Con ϕ en radianes. El ángulo ϕ tiene un número infinito de valores que difieren en

y
múltiplos enteros de 2π . Dichos valores se calculan con la ecuación tgϕ = .
x

Cada valor de ϕ es llamado argumento de z , y el conjunto de todos los valores se denota

con arg z . El valor principal de arg z , denotado por Arg z , es el único valor Ψ tal que
0 ≤ Ψ < 2π . Se pueden escoger otros intervalos, por ejemplo: −π < Ψ ≤ π . Luego,
arg z = Arg z + 2nπ , n = 0, ±1, ±2,.... .

iϕ
(2) Se define el símbolo e o exp(iϕ ) a través de la fórmula de Euler

e iϕ = (cos ϕ + i senϕ )

(3) Con esto, podemos escribir:

z = reiϕ
Esta ecuación la podemos interpretar como sigue:

e iϕ dice en cuál dirección está el punto, mientras que r dice, cuán lejos está el punto del
origen del plano complejo en esa dirección. A ϕ se le denomina fase del número complejo
y r , magnitud.

(4) Se define la operación:

iϕ iϕ2
e 1e = (cos ϕ1 + i senϕ1 )(cos ϕ2 + i senϕ2 )

13
= (cos ϕ1 cos ϕ2 − senϕ1senϕ2 ) + i (senϕ1 cos ϕ2 + cos ϕ1senϕ2 )

i (ϕ1 +ϕ2 )
= cos(ϕ1 + ϕ2 ) + isen(ϕ1 + ϕ2 ) = e

(5) Con esta puede demostrarse:

i (ϕ1 +ϕ2 )
z 1z 2 = r1r2 e

z1 r1 i (ϕ1 −ϕ2 )
= e
z2 r2

1 1 −iϕ
z −1 = = e ,z ≠0
z r
z n = r n einϕ n = 0, ±1, ±2,....
(6) La fórmula de Moivre

(cos ϕ + i senϕ )n = cos nϕ + i sennϕ con n = 0, ±1, ±2,....

14
Funciones Pares

Sea una función f (x ) de valor real de una variable real x . Entonces esta función será par si
se cumple, para todo x en el dominio de f (x ) , que:

f (x ) = f (−x )
El gráfico de una función par es simétrico con respecto al eje y = f (x ) , es decir, su gráfico
es el mismo después de hacer una reflexión alrededor del eje y = f (x ) . Ejemplos de
funciones pares: C , x 2 , cos x , x , cosh x .

15
 ƒ(x) = x2

16
* Los conceptos de paridad solo están definidos para funciones cuyo dominio y rango tienen
el inverso aditivo (esto es: dado un número x , el inverso aditivo es el número que cuando
se le suma, resulta 0 ). Estos razonamientos pueden aplicarse a funciones de valor
complejo.

Funciones Impares

Sea una función f (x ) de valor real de una variable real x . Entonces es impar si se cumple
que:

f (−x ) = − f (x )

para todo x en el dominio de f (x ) .

f (−x ) + f (x ) = 0

El gráfico de una función impar tiene simetría con respecto al origen en el sentido que el
gráfico es invariante luego de una rotación de 180 alrededor del origen. Ejemplos: x , x 3 ,
senx , senhx , tan x , erf (x ) , erf −1 (x ) , integrales de Fresnel C (x ) y S (x ) .

17
 ƒ(x) = x3

Función error o Función error de Gauss

18
Integrales de Fresnel

ƒ(x) = x3 + 1 no es ni par ni impar.

19
Si una función par es diferenciable, entonces su derivada es una función impar.
Analógamente, si una función impar es diferenciable, su derivada es una función par.
(diferenciabilidad quiere decir, en general, que la derivada cuando se calcula en un punto,
existe)

La serie de Maclaurin de una función par contiene solo potencias pares, mientras que las de
una función impar contiene solo potencias impares.

Recordemos que la serie de Maclaurin es una expansión en serie de Taylor evaluada en 0 :

f ′(0) f ′′(0) 2 f ′′′(0) 3

f (x ) = f (0) + x+ x + x + ...
1! 2! 3!
O
∞
f (n ) (0) n
f (x ) = ∑
0 n!
x

20
Como se sabe, mientras más alto es el grado de la aproximación del polinomio de Taylor,
obtenemos una curva cada vez más cercana a la función correcta. Esta imagen muestra las
curvas correspondientes a los polinomios de grado 1, 3, 5, 7, 9, 11 y 13.

21
Propiedades de las funciones pares e impares

(1) La función f (x ) = 0 , f (x ) = 0, ∀x , es tanto par como impar.

(2) Si f (x ), g(x ) son pares, entonces F (x ) = f (x ) ± g(x ) , también es par.

(3) Si f (x ), g(x ) son impares, entonces F (x ) = f (x ) ± g(x ) , también es impar.

(4) Si f (x ) es par, entonces F (x ) = af (x ) ( a contante) es par. El mismo razonamiento

sigue si es impar.

(5) Si f (x ), g(x ) son pares, entonces f (x )g(x ) , también es par.

(6) Si f (x ), g(x ) son impares, entonces f (x )g(x ) , es par.

(7) Si f (x ) es una función par, y g(x ) es una función impar, entonces f (x )g(x ) es una
función impar.

f (x )
(8) Si f (x ), g(x ) son pares, entonces , también es par.
g (x )

f (x )
(9) Si f (x ), g(x ) son impares, entonces , es par.
g (x )

f (x )
(10) Si f (x ) es una función par, y g(x ) es una función impar, entonces , es impar.
g (x )

(11) Si f (x ) es una función par, f ′(x ) es impar.

(12) Si f (x ) es una función impar, f ′(x ) es par.

(13) Si f (x ), g(x ) son pares, entonces f (x ) g(x ) = f (g(x )) (la composición) es par.

(14) Si f (x ), g(x ) son impares, entonces f (x ) g(x ) = f (g(x )) es impar.

(15) Si f (x ) es una función par, y g(x ) es una función impar, entonces

f (x ) g(x ) = f (g(x )) es par.
22
(16) Si f (x ) es una función par o impar, y g(x ) es una función par, entonces
f (x ) g(x ) = f (g(x )) es par.

+a

(17) Si f (x ) es una función impar, la integral ∫

−a
f (x )dx = 0 (con a finito, y f (x ) no

tenga asíntotas verticales en el intervalo.)

+a +a

(18) Si f (x ) es una función par, la integral ∫

−a
f (x )dx = 2 ∫ f (x )dx
0
, (con a finito, y

f (x ) no tenga asíntotas verticales en el intervalo. Es válido también cuando a → ∞ pero

solo si la integral converge)

(19) La suma de una función par y una función impar no es ni par ni impar, a menos que una
de las funciones sea igual a cero sobre el dominio dado.

(20) Una función cualquiera f (x ) puede ser expresada como la suma de una función par y
una función impar.

f (x ) f (x ) f (−x ) f (−x )
f (x ) = + + −
2 2 2 2

f (x ) f (−x ) f (x ) f (−x )
f (x ) = + + − = g(x ) + h(x )
2 2 2 2
Con

f (x ) f (−x ) f (x ) f (−x )
g(x ) ≡ + , h(x ) ≡ −
2 2 2 2

Luego, calculemos g (−x ) y h(−x )

23
 f (−x ) f (−(−x )) f (−x ) + f (x )
g(−x ) ≡ + = = g(x ) Función par
2 2 2

f (−x ) f (−(−x )) f (−x ) − f (x )

h(−x ) ≡ − = = −h(x ) Función Impar
2 2 2

Esto es verdad también para funciones de valor complejo. En ese caso, todas las funciones
pueden descomponerse en:

f (x ) = g(x ) + i j (x )

con g(x ),j (x ) a valores reales. Luego, se puede expresar cada función g(x ), j (x ) como:

g(x ) = G par + Gimpar y j(x ) = J par + J impar .

24
La Delta de Dirac Ref: Arfken, Weber. Essential Mathematical Methods for Physicists

Un objeto puntual puede representarse a través de la denominada “función” delta de Dirac

que tiene las siguientes propiedades (nos concentramos en problemas con una variable):

δ (x ) = 0, x ≠0

+∞ +∞
f (0) = ∫
−∞
f (x )δ (x )dx , como caso especial ∫
−∞
δ (x )dx = 1 .

Con f (x ) una función bien comportada y la integración incluye al origen.

δ (x ) debe ser infinitamente ‘’alta’’ en x = 0 . Por esta razón, los matemáticos para ser más
exactos, denominan a la delta de Dirac, una distribución. Es decir, una secuencia de
+∞
funciones que el límite verifican la ecuación f (0) = ∫
−∞
f (x )δ (x )dx . Por ejemplo: la delta

de Dirac puede ser aproximada por la secuencia de funciones en n para n → ∞

 1
 0, x < −
 2n
 1 1
δ n (x ) = 0, − <x <
 2n 2n
 1
0, x >
 2n

25
 n
δ n (x ) = exp(−n 2x 2 )
π

n 1
δ n (x ) =
π 1 + n 2x 2

26