Leer Data File
Leer Data File
Leer Data File
ans =
'C:\Users\Luis Sanchez\Documents\UNI Ciclo 2020-II\Procesamiento de Señales y Datos\Semana 04'
filename1 = fullfile(pwd,'Electrocardiograma_basico.jpg')
filename1 =
'C:\Users\Luis Sanchez\Documents\UNI Ciclo 2020-II\Procesamiento de Señales y Datos\Semana 04\Electrocardiograma_bas
1
typicalECG = imread('https://upload.wikimedia.org/wikipedia/commons/b/bd/12leadECG.jpg');
imshow(typicalECG)
Esta imagen (cargada en Wikipedia por MoodyGroove el 24-01-2007) muestra un ECG típico. Tenga en cuenta
que podemos leer y mostrar imágenes de Internet directamente en Matlab, sin necesidad de descargarlas en
su computadora.
Cargando y mostrando un ECG desde un archivo en su computadora. El primer paso es leer el archivo que
contiene el electrocardiograma, esto se puede hacer con el comando 'load' de la siguiente forma:
load Electrocardiogram
Este archivo contiene 250 segundos de una grabación de ECG, ya que la frecuencia de muestreo es de 200
Hz, el archivo contiene 50,000 puntos. Ahora podemos mostrar el ECG en el tiempo con comando plot:
plot(ecg)
2
Esta señal de datos ha sido leída del sitio web del profesor Rangaraj M. Rangayyan (http://people.ucalgary.ca/
~ranga/). Si sigue la sección "Notas de clase, ejercicios de laboratorio y archivos de datos de señales",
encontrará una gran cantidad de señales de diferentes modalidades de adquisición (ECG, EEG, EMF, etc.).
Puede guardar los archivos de un navegador en su disco duro y luego abrirlos con el comando cargar como se
indicó anteriormente. Alternativamente, también puede leer los conjuntos de datos directamente desde Matlab.
Las siguientes líneas le mostrarán cómo hacer esto.
file = 'http://people.ucalgary.ca/~ranga/enel563/SIGNAL_DATA_FILES/ecgpvc.dat';
ecg_string = urlread(file);
Debe notar que el comando 'urlread' ha proporcionado una variable de un tipo de datos diferente del archivo
que leímos anteriormente. Podemos enumerar las variables y sus características con el comando 'whos'
whos
3
k 1x1 8 double
s 1x10 2305344 struct
typicalECG 407x1040x3 1269840 uint8
typicalECG01 479x638x3 916806 uint8
vidFrame 240x320x3 230400 uint8
vidObj 1x1 8 VideoReader
videoReader 1x1 8 VideoReader
vidframes 240x320x3x42 9676800 uint8
xylophone 240x320x3x141 32486400 uint8
Dado que el archivo fue leído desde un sitio web, se supuso que sería una cadena de caracteres, por lo tanto,
debemos convertir esos caracteres en un número. Esto es sencillo, con el comando 'str2num'.
ecg2 = str2num(ecg_string);
plot(ecg2)
Hemos leído el archivo directamente desde el sitio web. Intente leer y mostrar graficamente otros archivos.
Zooming in
Como es evidente en el gráfico, el número de muestras puede dificultar la observación de los detalles de
cada ciclo cardíaco. Podemos acercarnos manualmente con la lupa de la figura, y también desde la línea de
comando con el comando axis, de la siguiente forma:
4
Leer Video
videoReader = VideoReader('xylophone.mp4');
implay(xylophone);
5
The audioread function can support WAVE, OGG, FLAC, AU, MP3, and MPEG-4 AAC files.
[y,Fs] = audioread('FemaleSpeech-16-4-mono-20secs.wav');
Play el audio.
info = audioinfo('FemaleSpeech-16-4-mono-20secs.wav')
sound(y,Fs)
t = 0:seconds(1/Fs):seconds(info.Duration);
t = t(1:end-1);
plot(t,y)
xlabel('Time')
ylabel('Audio Signal')