Tatiana Valentina Polanco Vélez - 20201573017
Jonathan David Aguirre Barrios - 20201573014
BITACORA
Ajustes de nivel de intensidad
Para controlar la saturación de los colores, es decir, que tanta intensidad de los colores hay en la
imagen, se usó la función imadjust la cual tiene como parámetros un rango de intensidad el cual
será reemplazado por los segundos valores.
En el código primero se leyó la imagen mano1.png con el comando imread, luego se le cambio el
nivel de intensidad y por último se mostró tanto la imagen original como la imagen modificada
por medio de imshow.
En caso de ser imágenes a color, son necesarios 3 rangos, uno para cada color del RGB, y por
ultimo los 3 rangos finales de saturación.
Ecualización e histograma
La ecualización de una imagen es un proceso de procesamiento de imágenes que ajusta la
distribución de la intensidad de los píxeles en una imagen para mejorar su contraste y claridad;
Este primer código, se creó una imagen desde cero, esta imagen es creada en escala de grises
La función imhist devuelve los recuentos del histograma de counts y las ubicaciones de los bins de
binLocations, el funcionamiento de esta función se ve reflejada en las gráficas de la derecha de cada
imagen, o sea en los subplot 222 y 224
En este segundo código se utiliza igualmente la función imhist que se ve en la imagen 322, pero, esto ya
a una imagen del archivo que nos compartió el profesor y no una creada.
En la imagen 323 se utilizó el comando histeq, esta función ecualiza la imagen directamente
Interpolación
La interpolación el cambio de la resolución de la imagen al aumentar o disminuir la cantidad de
pixeles, y la función para esta era interp2.
Rotación
La rotación es girar la imagen alrededor de un punto fijo, para esto se usa el comando
imrotate(imagen, x), requiriendo como parámetros la imagen a girar, y el valor de cuantos grados
sexagesimales se va a girar la imagen.
Traslación
La traslación se refiere a una transformación geométrica que se aplica a una imagen para
moverla de su posición original a una nueva posición, la función que se utiliza para hacer la
traslación de una imagen es imtranslate, donde se utilizan como parámetros la imagen en escala
de grises y el vector de traslación.
Binarización
La binarizacion de una imagen es reemplazar los valores de los pixeles que estén de 0 a 255 que
es el caso de escala de grises a valores 0 y 1, es decir, negro y blanco (bajar la resolución de 8
bits a 1 bit), para esto toma un umbral donde divide los valores de 0 hasta el umbral (0 ≤ x <
umbral) y los reemplaza por 0, y luego del valor umbral hasta 255 (umbral ≤ x ≤ 255) los
reemplaza por 1.
Primero se visualizó la imagen original en el subplot 221, luego se realizó una binarizacion por
medio de la función for para un único umbral y se imprimió en el subplot 222, luego se cambió
el valor del umbral y se presentó en el subplot 223, para por último se realizó el mismo ejercicio,
pero con dos umbrales, donde los valores que estén entre 0 y el primer umbral (0 ≤ x < umbral 1)
los reemplazó por 0, los valores que estén entre el primer y segundo umbral (umbral 1 ≤ x ≤
umbral 2) los reemplazó por 0.5 (gris) y por último el resto de valores (umbral2 < x ≤ 255) se
reemplazaron por 1, y este proceso se mostró en el subplot 224.
Aunque este proceso se puede resumir por medio de la función im2bw (imagen, umbral).
Mapeo RGB A YIQ
En el procesamiento de imágenes, el mapeo RGB a YIQ es una técnica utilizada para convertir
una imagen en formato RGB (rojo, verde, azul) en una imagen en formato YIQ (luminancia, in-
phase, quadrature-phase), en Matlab el comando que hace automáticamente esta conversión es
rgb2ntsc, ya al tener la imagen en YIQ se puede modificar la luminosidad, saturación o matiz
Filtrado básico
El primer filtro básico, fue el filtro pasa bajos el cual se encarga de quitar partículas blancas que
se encuentren en la imagen, lo hace por medio de un kernel o por medio de un filtro gaussiano
con la función fspecial donde se le indica el tamaño (k_size) y su desviación estándar (sigma),
para aplicar estos filtros a la imagen original se usa la función filter2 y para pasar esa nueva
matriz a una imagen en escala de grises se usa la función mat2gray.
El segundo filtro básico es el filtro mediano, el cual se encarga de enfocar o desenfocar una
imagen, para esto se usa la función medfilt2 en donde la mediana del entorno de los dos
parámetros solicitados es tomada por cada pixel alrededor de la imagen de entrada y lo coloca a
la salida.
Uno puede aplicar filtros para poder detectar el borde de una imagen, para esto se aplica al igual
que el filtro pasa bajo con kernel, pero con un kernel horizontal o vertical.
Operadores morfológicos
Los operadores morfológicos se pueden dividir en operadores lógicos y ventanas. Las ventanas
se refieren a como se van a tomar los pixeles.
- Fila: Row (2P+1)
- Columna: Col (2P+1)
- Cuadrados: Square (2P+1) ^ 2
- Cruces: Cross (2P+1)
Y entre los operadores lógicos se encuentran la dilatación que consiste en ensanchar las zonas
blancas de una imagen y adelgazar las zonas negras y su función es imdilate, está la erosión que
es su inverso (ensanchas las zonas negras y adelgazar las zonas blancas) su función es imerode,
está el majority que es un filtro mediano, el filtro open que es la dilatación de una erosión y su
función es imopen y por último el filtro close que es la erosión de una dilatación y su función es
imclose.
Así mismo se pueden mezclar los operadores lógicos con las ventanas.
Análisis de Fourier
El análisis de Fourier se utiliza en el procesamiento de imágenes para analizar las características
de frecuencia de una imagen y descomponerla en sus componentes espectrales
La transformada discreta de Fourier (DTF) se utiliza para analizar y descomponer una imagen en sus
componentes de frecuencia, la función fft2 devuelve la transformada bidimensional de Fourier de una
matriz utilizando un algoritmo de la transformada rápida de Fourier, la función bar es este código es
para modelar la transformada bidimensional de Fourier en una gráfica de barras y por último la función
stem convierte la transformada bidimensional de Fourier en una gráfica secuencia de datos P ( en este
caso) como puntos con líneas verticales que se extienden desde una línea de referencia a lo largo del eje
x.
Visualización DFT 2D
La transformada discreta de Fourier para imágenes, nos permite encontrar las frecuencias
espaciales escalizadas que están en ciclos/pixel, la función para crear esto es fft2, y para
reorganizar la transformada para que el componente con frecuencia cero quede en el centro es
por medio de fftshift.
Por medio de una transformada inversa de Fourier y una máscara pasa bajo, se puede crear una
convolución cíclica, la función para una transformada inversa es ifft2.
Filtrado con funciones
Se pueden generar filtrados de imágenes en el dominio de la frecuencia, entre estos está el filtro
pasa bajos ideal LPF el cual en la práctica es un círculo blanco en fondo negro, para este filtro se
usaron las funciones de la transformada de Fourier discreta, su transformada inversa y por último
la función distmatrix que nos permite hallas la distancia de la matriz de la imagen.
También está el filtro LPF gaussiano, el cual entre mayor se la n, este será más cercano al LPF
ideal.
Así como existe el filtro ideal pasa bajos, también está el filtro pasa altos ideal HPF, este en
cambio usa un círculo negro con fondo blanco.
Luego es esta el filtro pasa altos gaussiano.
El pasa altos Butterworth
Y por último el pasa altos de énfasis
DWT
DWT se refiere a la transformada wavelet discreta la cual descompone la señal en varios
conjuntos, en la práctica se generó la suma de dos señales senoidales, la primera donde la
amplitud es 1 y su frecuencia es 5, y la segunda es de amplitud 0.2 con frecuencia 500.
Después se aplicó su transformada por medio del comando DWT devolviendo dos datos, el
coeficiente de aproximación 1 y el coeficiente de detalle 1, para por último graficarlos.
Por último, se volvió a sacar la DWT del coeficiente de aproximación, para obtener los segundos
coeficientes de aproximación y detalle.
Después se realizó la DWT, pero con una imagen para esto el comando es dwt2, al realizar la
transformada se reciben cuatro valores que corresponden al primer coeficiente de aproximación,
coeficiente horizontal, coeficiente vertical y coeficiente diagonal, y se visualizaron por medio del
comando image y la conversión de la función wcodemat que convierte la escala del coeficiente a
un rango de 1 hasta 255.
Luego se procedió a reconstruir una imagen por medio de una wavelet boirtogonal, para esto se
usó el comando wfilters(‘bior3.5’), la cual devuelve los dos filtros pasa bajos y los dos filtros
pasa altos asociados a la wavelet bior3.5, luego se usó la función wavedec2 que usa los filtros de
descomposición pasa bajos y el filtro pasa altos para descomponer la imagen X.
Después se organizaron los elementos de la descomposición C1, para más adelante con el
comando waverec2 reconstruir la imagen, con los filtros pasa bajos y pasa altos.
Por último, por medio de la familia wavelet haar con el comando haart2 se realizó su
transformada hasta el nivel 3, se aplicó esto tanto para la imagen original, como la imagen de la
marca de agua, para después por medio del comando ihaart2 realizar la transformada inversa con
los coeficientes de aproximación y detalle horizontal, vertical y diagonal.
Fourier vs Wavelet
La Transformada de Fourier y la Transformada Wavelet ambas son técnicas se utilizan para
analizar las características de frecuencia de una imagen, hay diferencias fundamentales entre
ellas:
La Transformada de Fourier transforma una imagen del dominio espacial al dominio de la
frecuencia, descomponiendo la imagen en una serie de componentes sinusoidales de diferentes
frecuencias. Muestra la distribución de amplitudes y fases en el dominio de la frecuencia.
Por otro lado, la Transformada Wavelet también realiza una descomposición de la imagen, pero
utiliza funciones llamadas wavelets que son localizadas en el tiempo y la frecuencia. Esto
permite una representación más localizada de la información de la imagen, ya que las wavelets
pueden capturar tanto las características de baja frecuencia como las de alta frecuencia en
diferentes regiones de la imagen.
Para este código se usó la función syms, para crear 2 variables y y t, con la y se crea una función
sinusoidal; La función subs(y) se utiliza para evaluar la función simbólica y en el vector t, generando así
la señal f en el dominio del tiempo; En la línea 10 f(7*N/16:9*N/16) = zeros; se establecen a cero los
valores de la señal f en un rango específico del vector t; Se aplica la Transformada de Fourier a la señal f
utilizando fft(f). El resultado se almacena en la variable F.
Luego se crea un vector freq que representa los valores de frecuencia correspondientes a la
Transformada de Fourier, en el rango de 0 a 30 veces la inversa de T.
En la subplot 221, se traza la señal f en función de t; En el subplot 223 se grafica la magnitud de los
coeficientes de Fourier correspondientes a diferentes frecuencias; En el subplot 222 se realiza una
transformada wavelet continua de la señal f utilizando la función cwt. Se utiliza la función de onda de
Morlet y se grafica el resultado.
En el subplot 224 se muestra una representación de imagen de los coeficientes de la transformada
wavelet continúa utilizando imagesc. Los coeficientes se representan en función de t y las diferentes
escalas.
Talleres de investigación
1. Equipo de rayos x estacionarios:
* Equipo de rayos x telecomandado digital con sustracción digital y detector plano
Características:
- Mesa de pacientes: Basulable mediante motor a +90° y - 17° en Trendelemburg,
con parada automática en 0°, permite soportar hasta 150 kg, consola de telemando
para todas las funciones de la mesa
-Detector de panel plano: Detector de plano dinámico digital de alta resolución
con rejilla desmontable integrada, Material semiconductor: silicio amorfo (a-Si),
tamaño de 43 cm x 43 cm, tamaño de pixel 148 mµ (6.76 pixeles por mm),
Tamaño de matriz: 2880 x 2880 (8.3 millones de pixeles), Resolución de detalles:
3.4 Pl/mm, profundidad de adquisición de 16 bits.
-Generador de rayos X: Generador de rayos X de alta frecuencia, controlado por
microprocesadores, de 65 KW de potencia con 100 KV, cuenta con un rango
desde 40 hasta 150 KV, Un tubo de rayos X de alta capacidad calorífica, de ánodo
giratorio, capacidad de almacenamiento de calor en el ánodo sobre 780000 HU,
con 2 focos, fino y grueso en un rango de 0.5 y 1.5 mm, Funciona con 380 hasta 4
80 Voltios a 60 Hz.
-Sistema digital: Disco duro para almacenar 50.000 imágenes en matriz 1024 x
1024, fluoroscopia digital de 30 cuadros/seg, con las siguientes funciones de
procesamiento de imágenes, binarización de imagen, zoom, roaming, inversor de
imagen derecha/izquierda, medidas de ángulo, distancias, Directorio de pacientes
por nombre, número, órganos, médico, etc. Substracción digital DSA en tiempo
real.
El equipo cuenta con la posibilidad de conectarse de forma remota con la fábrica o
centro de servicio remoto para detección y reparación de fallas.
2. Tomógrafo computarizado
* Unidad para tomografía computarizada para estudios avanzados de 128 cortes
Características:
- Con 64 hileras o filas de detectores, diagnóstico remoto con el proveedor local y
con fabrica (referido al soporte técnico) el router y/0 dispositivo de conexión
deberá será proporcionado por el proveedor del equipo, la línea y/0 conexión a
internet estará a cargo de la entidad, tecnología de reducción de dosis de radiación
por reconstrucción iterativa y técnica de modulación de dosis, fantomas propios
del equipo para su calibración, verificación y/0 control de calidad; necesarios para
su puesta en funcionamiento
- Gantry: apertura de 700 mm, ángulo de inclinación: de - 30° a +30°, con
incrementos de 1° (un grado) o menos controlado desde el gantry o desde la
consola, luz laser para posicionamiento del paciente
- Generador: Por Inversor de alta frecuencia, potencia nominal de 70 kW, voltaje
máximo de 135 KV, rango de corriente desde 20mA hasta 600mA
-Tubo de rayos X: características técnicas de voltaje y corriente iguales a las del
generador, tubos de rayos x convencionales con capacidad de almacenamiento de
calor en el ánodo de 7.0 mhu (millones de unidades de calor).
-Detectores: resolución espacial de 19.5 ip/cm (en pares de líneas por centímetro),
cobertura máxima del detector en el eje z, por rotación: 38 mm
-Software especializados en la estación de trabajo: imágenes tridimensionales 3d,
software cardiológico especializado (visualización y cuantificación de calcificación en
coronarias - calcium scoring). software con sistema avanzado de evaluación vascular y
traumas múltiples con segmentación automática vascular. visualización cardiaca
(software de evaluación cardiaca no invasivo (análisis funcional del corazón) con software
para uso del gatillador). software cardiológico funcional. Software de evaluación de
coronarias. programa de perfusión cerebral con evaluación de isquemias con mapas
cerebrales a color y visualización de flujo sanguíneo dinámico en volumen 3d, mip o vrt
programa de colonoscopia virtual (puede presentarse en la consola de adquisición).
software automático de cuantificación y evaluación de nódulos pulmonares. software de
sustracción automática ósea (auto eliminación de huesos) y volumetría, software de
segmentación vascular hepática automática y/0 manual. Fusión de imágenes entre ct - mr
- pet – spect, software para estudios dentales software para estudios oncológicos software
para traumatología y ortopedia (puede presentarse en la consola de adquisición), software
para estudios de vías biliares.
2.1 X-ray detectors
* Detector de rayos x digital, 14HK701G-WP
Características:
- Contiene tecnología TFT, convirtiendo con precisión la energía de rayos x que atraviesa
el cuerpo humano en imágenes digitales en un monitor. El dispositivo ofrece
procesamiento de imágenes de 16 bits, lo que mejora aún más el contraste de la imagen y
admite una distancia entre pixeles de 140µm para imágenes de alta resolución, este
contiene una dimensión externa de 14 x 17 pulgadas, junto con una duración de carga de
3 horas.
- El detector puede conectarse de manera inalámbrica y tiene un ciclo de 11 segundos
obteniendo una imagen completa en tan solo 6 segundos, y en caso de estar conectado de
manera inalámbrica los tiempos se reducen a 8 segundos para el ciclo y de 3 segundos.
- Aparte de su calidad de imagen excepcional, también está diseñado para ser duradero y
portátil, el dispositivo está fabricado con materiales ligeros de fibra de carbono, lo que
facilita su transporte, con un peso de tan solo 2,95 kg. A pesar de su ligereza, el
dispositivo puede soportar un peso de carga completa de hasta 300 kg y un peso de carga
local de hasta 100 kg, lo que lo hace extremadamente robusto.
3. Sistema de resonancia magnética
* Maquina IRC Toshiba Vantage Titan 1.5T
- Contiene un conjunto integral de aplicaciones disponibles en los sistemas de 16 y 32
canales de transmisión y recepción. Contiene cuatro generaciones de técnicas sin
contraste aumentan la seguridad del paciente al capturar imágenes de angiografías por
RM de alta calidad, eliminando la necesidad de pruebas de laboratorio y reduciendo la
interrupciones de los horarios de los pacientes debido a complicaciones.
- Tiene un consumo de energía reducido, con plataforma de software M-Power y
tecnología de bobina integrada que se ajusta automáticamente a la configuración del
sistema y al examen que se esté realizando.
- Usa una potencia de imán de 1.5 T con un peso magnético de 4 toneladas, contiene una
amplitud máxima de los gradientes de 30 mT/m, unos requisitos mínimos para la
superficie de 55 m², campo de visión: 55 х 55 х 50 cm, velocidad de aumento del gradiente de
120 mT/m/m, ruido acústico de 30 dB, peso corporal máximo de 230 kg, soporte de infusión
montado en el gantry, reposabrazos integrado para fácil aplicación de aguja, requiere
trabajos de instalación menores, todo esto bajo el uso de técnicas sin contraste, adquisición
isotrópica 3D, supresión de grasa, angiograma tiempo-espacio.
