High Quality Integrated Data Reconstruction for Medical Applications

In this paper, the implementation of a high quality integrated data reconstruction model and algorithm has been proposed, especially for medical applications. Patients' Information acquired at the sending end and reconstructed at the receiving end by using a technique that would be high quality for the signal reconstruction process. A method is proposed in which the reconstruction of data like ECG, audio and other patients' vital parameters that are acquired in the time-domain and operated in the frequency-domain. Further the data will be reconstructed in the time-domain from the frequency domain where high quality data is required. In this particular case, high quality ensures the distortion less and noiseless recovered baseband signal. This would usually require the application of Fast Fourier Transform (FFT) and Inverse Fast Fourier Transform (IFFT) return the data to the spatial domain. The simulation is performed using Matlab. The Composite baseband signal has been generated by developing a program as well as by acquiring to the workspace. The feature of the method is that it can achieve high-quality integrated data reconstruction and can be associated easily with spatial domain.


I. INTRODUCTION
High quality integrated data reconstruction model and algorithm are used, within an electronic system, to extract the desired time-domain signal from the frequency-domain signal acquired from the human body (in offline), especially Electrocardiogram (ECG), audio and other vital parameters [1].Each stage automatically generates a template of a source from the candidate events in the initialization period, and thereafter performs classification of the remaining candidate events based on a template matching technique.Matlab simulation results on offline demonstrate the effectiveness of the proposed method.
In recent literature [2,3,4], the perception of conveying vital information reconstruction used by medical practitioners has had some concentration.Li, Mueller and Ernst [2] emphasized on the methods for efficient, high quality volume resampling in the frequency domain.It was described in the use of frequency-domain filters for the accurate resampling of images and volumes at arbitrary magnification factor.It was also investigated the Freq approach in relation to higher-quality filters.
Chazan, Hoory, Cohen, and Zibulski [3] illustrated speech reconstruction from Mel frequency cepstral coefficients and pitch frequency.They presented a novel low complexity, frequency domain algorithm for reconstruction of speech from the Mel frequency cepstral coefficient.The construction technique was based on the sinusoidal speech representation.
Kong [4] focused on GPS modeling in Frequency Domain.The author presented a frequency domain modeling approach to model GPS errors and increase GPS positioning shaping filter.The application of approach was mainly used for vehicle navigation system.
In this paper, most of the information contained in the baseband signal is found below 100 Hz.High frequencies random noise may corrupt the reconstructed time domain baseband signal.To remedy the situation, Filter Design Tool (FDA) has been used to eliminate the high frequency component.There is no distortion appearing in that particular spatial domain instead of attenuation amplified by the gain of the signal.
Bioelectrical signals are typically very small in amplitude (mV) and an amplifier is required to accurately depending on the hardware and software used, the biological amplifier serves to amplify the signal.It is also known that the frequency of heart signals is very low, approximately 5 to 10 H Z [5,6].

II. MATERIALS AND METHODS
The reconstruction process would be in order of data collection, Discrete Fourier transform (DFT), Fast Fourier Transform (FFT), Inverse Fourier Transform (IFFT), and finally, noise cancellation.These topics are discussed below: The discrete Fourier transform, or DFT, is the primary tool of digital signal processing.The foundation of the Signal Processing Toolbox is the fast Fourier transform http://ijacsa.thesai.org/(FFT), a method for computing the DFT with reduced execution time [7,8].Many of the toolbox functions (including z-domain frequency response, spectrum and cepstrum analysis, and some filter design and implementation functions) incorporate the FFT.MATLAB provides the functions FFT and IFFT to compute the discrete Fourier transform and it"s inverse, respectively.For the input sequence x and its transformed version X (the discrete-time Fourier transform at equally spaced frequencies around the unit circle), the two functions implement the relationships.
In these equations, the series subscripts begin with 1 instead of 0 because of the MATLAB vector indexing scheme, and Working with the Fourier transform on a computer usually involves a form of the transform known as the discrete Fourier transform (DFT).There are two principal reasons for using this form: The input and output of the DFT are both discrete, which makes it convenient for computer manipulations.There is a fast algorithm for computing the DFT known as the fast Fourier transform (FFT).
The DFT is usually defined for a discrete function f (m, n) that is nonzero only over the finite region 0<= m <= M and 0<= n <= N. The two-dimensional M-by-N DFT and inverse M-by-N DFT relationships are given by The values F (p , q) are the DFT coefficients of f (m , n).
The zero-frequency coefficient, F (0, 0), is often called the "DC component."DC is an electrical engineering term that stands for direct current.(Note that matrix indices in MATLAB always start at 1 rather than 0; therefore, the matrix elements f(1,1) and F(1,1) correspond to the mathematical quantities f (0, 0) and F (0 , 0) , respectively.) The DFT coefficients F (p , q) are samples of the Fourier transform F(ω 1 , ω 2 ) .

III. DATA COLLECTION
Data collection process can be acquired by real time or offline.In real time, data, especially ECG are acquired by data acquisition tools.Other vital parameters like audio, video can be acquired by existing audio-video editing software [9].If the existing software does not support to acquire the above parameters, it can be considered for developing the software to perform the acquisition properly.The array of electrodes used for the capture of Body Surface Potentials (BSP) comprises a flexible plastic anterior and posterior electrode harness and a portable recording unit.The anterior harness contains 64 electrodes, including 3 proximal bipolar limb leads (Mason-Likar position), and a posterior harness with 16 electrodes.This lead configuration enables recording of 77 unipolar ECG signals with respect to the Wilson central terminal.Recordings are taken over 4198 ms, the sampling rate is 1 kHz and the potentials are recorded using a 12-bit resolution to represent a 5 mV range.The procedure for attaching the harnesses to a patient and taking measurements of electrical potentials lasts about 3-4 minutes [10].
In this regard, data have been taken off-line generated by the software Matlab.It is worth mentioning that different data carries the different information that would be recovered by the method given in the next part of this paper.

IV. PROPOSED METHOD
The algorithm is based on using a sinusoidal model, where audio signal and ECG are represented by workspace block converted the captured file, especially .wavfile into 2dimentional data file.signals are characterized by the frequencies, amplitudes and phases of the sine-wave components.The complete reconstruction scheme is presented in Figure .1.It is applied in every 10 ms.frame using the simulink block system"s rules and regulations.
The algorithm is comprised of the following stages: 1. Sum block -Add or subtract inputs.Specify one of the following: a) String containing + orfor each input port.b) Scalar >=1.A value >1 sums all inputs; 1 sums elements of a single input vector.http://ijacsa.thesai.org/ 2. FFT -Outputs the complex Fast Fourier Transform (FFT) of real or complex input by computing radix-2 decimation-in-time (DIT) or decimation-in-frequency (DIF), depending on block options.Uses half-length and doublesignal algorithm for real inputs where possible.Computes the FFT Figure 1.Block Diagram of proposed simulation method along the vector dimension for sample-based vector inputs, which must have a power-of-2 length.Computes the FFT along each column for all other inputs, where columns must be a power-of-2 length.
3. IFFT -Outputs the Inverse Fast Fourier Transform (FFT) of real or complex input by computing radix-2 decimation-in-time (DIT) or decimation-in-frequency (DIF), depending on block options.
Outputs are real if we select "Input is conjugate symmetric" option; otherwise, outputs are complex.Computes the IFFT along the vector dimension for samplebased vector inputs, which must have a power-of-2 length.Computes the IFFT along each column for all other inputs, where columns must be a power-of-2 length.
4. FDA Tool -The block applies the specified filter to each channel of a discrete-time input signal, and outputs the result.The outputs of the block numerically match the outputs of the Digital Filter block, the filter function in MATLAB, and the filter function in the Filter Design Toolbox.The sampling frequency, Fs, which specified in the FDA tool"s Graphical User Interface (GUI), should be identical to the sampling frequency of the digital filter design block's input block.When the sampling frequencies of these blocks do not match, the digital filter design block returns a warning message and inherits the sampling frequency of the input block.
The block accepts inputs that are sample-based or framebased vectors and matrices.The block filters each input channel independently over time, where  Each column of a frame-based vector or matrix is an independent channel.
 Each element of a sample-based vector or matrix is an independent channel. 5. Workspace -Simulink allows us to import input signal and initial state data from the MATLAB workspace and export output signal and state data to the MATLAB workspace during simulation.This capability allows us to use standard or custom MATLAB functions to generate a simulated system's input signals and to graph, analyze, or otherwise postprocess the system's outputs STFT reconstruction (Fig. 1) -The frequencies, phases and [8192x1] [8192x1] [8192x1] http://ijacsa.thesai.org/amplitudes are combined to form a sine-wave representation.The final reconstructed STFT is constructed from the sine waves by a convolution procedure.
To use this format, it selects input in the load from workspace pane and selects the Array option from the Format list on the data import/export pane.Selecting this option causes Simulink to evaluate the expression next to the Input check box and use the result as the input to the model.The expression must evaluate to a real (noncomplex) matrix of data type double.The first column of the matrix must be a vector of times in ascending order.The remaining columns specify input values.In particular, each column represents the input for a different Import block signal (in sequential order) and each row is the input value for the corresponding time point.Simulink linearly interpolates or extrapolates input values as necessary if the Interpolate data option is selected for the corresponding import.The total number of columns of the input matrix must equal n + 1, where n is the total number of signals entering the model's imports.
6. Vector Scope -The vector scope block is a comprehensive display tool similar to a digital oscilloscope.The block can display time-domain, frequency-domain, or user-defined signals.We can use the Vector Scope block to plot consecutive time samples from a frame-based vector, or to plot vectors containing data such as filter coefficients or spectral magnitudes.To compute and plot the periodogram of a signal with a single block, use the Spectrum Scope block.The input to the Vector Scope block can be any realvalued M-by-N matrix, column or row vector, or 1-D (unoriented) vector, where 1-D vectors are treated as column vectors.Regardless of the input frame status, the block treats each column of an M-by-N input as an independent channel of data with M consecutive samples.The block plots each sample of each input channel sequentially across the horizontal axis of the plot conversion to time domain (Fig. 1) -This is the final stage in which the time domain signal frame is computed.
7. Sine wave -The sine wave block may generate a multichannel real or complex sinusoidal signal, with independent amplitude, frequency, and phase in each output channel.A real sinusoidal signal is generated when the output complexity parameter is set to Real, and is defined by an expression of the type y = A sin (2πft + φ) where A in the amplitude parameter, f in hertz in the frequency parameter, and in radians in the phase offset parameter.A complex exponential signal is generated when the output complexity parameter is set to complex, and is defined by an expression of the type 8. Random source -Output a random signal with uniform or Gaussian (normal) distribution.Set output repeatability to nonrepeatable (block randomly selects initial seed every time simulation starts), repeatable (block randomly selects initial seed once and uses it every time simulation starts), or specify seed (block uses specified initial seed every time simulation, starts producing repeatable output).

V. ANALYSIS AND DISCUSSION
The first stage of the implementation is to pass the simulated integrated baseband signal corrupt by high frequency random noise.The method tries to produce the filtered simulated distortion less and noiseless integrated baseband signal in a relatively short period of time.Increasing the filter order above one slows down the convergence rate but makes the results more precise.The recovered signal closely resembles the original simulated signal minus the noise and hence indicates that the implementation of the algorithm functions correctly and efficiently.The second stage is to recover the original composite baseband attenuated signal using FDA tool.Finally, the original composite banseband signal has been recovered where noise, distortion and attenuation have been eliminated.
High frequency noise is used, as the signal corrupt by noise as well as the reference noise.With this setup, it is observed that the output signal has been attenuated a little bit but not distorted due to the affect of distortion.But the attenuation would be disappeared by amplifying the gain of the signal.All involved steps of the simulation are discussed below.
Step 1: Generation of different sinusoidal signals having amplitude of 1 mV, ECG pattern having amplitude of 3.5mV and audio signal having amplitude of 0-1 mV. Figure 2, 3, 4, 5, and 6 show the generation of input signal.Here ECG signal has been generated in offline and the time domain data imported from the workspace.The data of the audio signal may take from stored file or real time.In this particular case, it has been taken from the stored file and the time domain data also imported from the workspace.
Step 2: The composite baseband signal for the different input signals (Amplitude of 4.5 mV) has been showed in Fig. 7.
Step 3: Fig. 8 shows the original composite baseband signal in frequency domain.
Step 4: Generation of high frequency random noise signal amplitude of 3.5 mV is shown in Fig. 9.
Step 5: The noisy IFFT signal (recovered baseband signal + noise) is shown in Fig. 10.
Step 6: Fig. 11 shows the reconstructed filtered attenuated composite baseband signal using FDA tool.

VI. CONCLUSION
In this work, an attempt has been made to employ a method that can be implemented for reconstructing a high quality integrated data, especially for medical applications.Distortions of reconstructed filtered baseband signal may occur due to the affect of frequency offset in integrated data that have been eliminated.Generally, the reconstruction of data like ECG, audio and other patients" vital parameters that are acquired in the time-domain and operated in the frequency-domain.Data acquisition and the simulation are performed using Matlab.The Composite baseband signal, especially ECG, has been generated by developing a program.Filter Design Tool (FDA) has been used to eliminate the high frequency component and the original filtered baseband signal has been obtained.Distortions have been eliminated in that particular spatial domain and attenuation has been eliminated by amplifying the gain of the signal.This approach may be considered for reconstruction in the advanced technology which needs

Figure 8 .
Figure 8. Original Composite Baseband Signal in Frequency Domain.