Analysis of co-occurrence texture statistics as a function of gray-level quantization for classifying breast ultrasound.

In this paper, we investigated the behavior of 22 co-occurrence statistics combined to six gray-scale quantization levels to classify breast lesions on ultrasound (BUS) images. The database of 436 BUS images used in this investigation was formed by 217 carcinoma and 219 benign lesions images. The region delimited by a minimum bounding rectangle around the lesion was employed to calculate the gray-level co-occurrence matrix (GLCM). Next, 22 co-occurrence statistics were computed regarding six quantization levels (8, 16, 32, 64, 128, and 256), four orientations (0° , 45° , 90° , and 135°), and ten distances (1, 2,...,10 pixels). Also, to reduce feature space dimensionality, texture descriptors of the same distance were averaged over all orientations, which is a common practice in the literature. Thereafter, the feature space was ranked using mutual information technique with minimal-redundancy-maximal-relevance (mRMR) criterion. Fisher linear discriminant analysis (FLDA) was applied to assess the discrimination power of texture features, by adding the first m-ranked features to the classification procedure iteratively until all of them were considered. The area under ROC curve (AUC) was used as figure of merit to measure the performance of the classifier. It was observed that averaging texture descriptors of a same distance impacts negatively the classification performance, since the best AUC of 0.81 was achieved with 32 gray levels and 109 features. On the other hand, regarding the single texture features (i.e., without averaging procedure), the quantization level does not impact the discrimination power, since AUC = 0.87 was obtained for the six quantization levels. Moreover, the number of features was reduced (between 17 and 24 features). The texture descriptors that contributed notably to distinguish breast lesions were contrast and correlation computed from GLCMs with orientation of 90° and distance more than five pixels.

Computerized lesion segmentation of breast ultrasound based on marker-controlled watershed transformation.

PURPOSE: This paper presents a computerized segmentation method for breast lesions on ultrasound (US) images. METHODS: It consists of first applying a contrast-enhanced approach, i.e., a contrast-limited adaptive histogram equalization. Then, aiming at removing speckle and enhancing the lesion boundary, an anisotropic diffusion filter, guided by texture descriptors derived from a set of Gabor filters, is applied. To eliminate the distant pixels that do not belong to the tumor, the resulting filtered image is multiplied by a constraint Gaussian function. By doing so, both the segmentation and the marker functions are generated and could be used in the marker-controlled watershed transformation algorithm to create potential lesion boundaries. Finally, to determine the lesion contour, the average radial derivative function is evaluated. The proposed method was tested with 50 breast US images and 60 simulated "ultrasound-like" images. Accuracy and precision of the segmentation method were then assessed. For the accuracy, three parameters were used: Overlap ratio (OR), normalized residual value (nrv), and proportional distance (PD) between contours. RESULTS: The average results for US images were OR = 0.86 +/- 0.05, nrv = 0.16 +/- 0.06, and PD = 6.58 +/- 2.52%. For simulated ultrasound-like images, a better performance (OR = 0.92 +/- 0.01, nrv = 0.08 +/- 0.01, and PD = 3.20 +/- 0.53%) was achieved. CONCLUSIONS: The segmentation method proposed was capable of delineating the lesion contours with high accuracy in comparison to both the radiologists' delineations and the true delineations of simulated images. Moreover, this method was also found to be robust to human-dependent parameters variations.

Topographic distribution of the tibial somatosensory evoked potential using coherence

The objective of the present study was to determine the adequate cortical regions based on the signal-to-noise ratio (SNR) for somatosensory evoked potential (SEP) recording. This investigation was carried out using magnitude-squared coherence (MSC), a frequency domain objective response detection technique. Electroencephalographic signals were collected (International 10-20 System) from 38 volunteers, without history of neurological pathology, during somatosensory stimulation. Stimuli were applied to the right posterior tibial nerve at the rate of 5 Hz and intensity slightly above the motor threshold. Response detection was based on rejecting the null hypothesis of response absence (significance level α= 0.05 and M = 500 epochs). The best detection rates (maximum percentage of volunteers for whom the response was detected for the frequencies between 4.8 and 72 Hz) were obtained for the parietal and central leads mid-sagittal and ipsilateral to the stimulated leg: C4 (87 percent), P4 (82 percent), Cz (89 percent), and Pz (89 percent). The P37-N45 time-components of the SEP can also be observed in these leads. The other leads, including the central and parietal contralateral and the frontal and fronto-polar leads, presented low detection capacity. If only contralateral leads were considered, the centro-parietal region (C3 and P3) was among the best regions for response detection, presenting a correspondent well-defined N37; however, this was not observed in some volunteers. The results of the present study showed that the central and parietal regions, especially sagittal and ipsilateral to the stimuli, presented the best SNR in the gamma range. Furthermore, these findings suggest that the MSC can be a useful tool for monitoring purposes.

Topographic distribution of the tibial somatosensory evoked potential using coherence.

The objective of the present study was to determine the adequate cortical regions based on the signal-to-noise ratio (SNR) for somatosensory evoked potential (SEP) recording. This investigation was carried out using magnitude-squared coherence (MSC), a frequency domain objective response detection technique. Electroencephalographic signals were collected (International 10-20 System) from 38 volunteers, without history of neurological pathology, during somatosensory stimulation. Stimuli were applied to the right posterior tibial nerve at the rate of 5 Hz and intensity slightly above the motor threshold. Response detection was based on rejecting the null hypothesis of response absence (significance level alpha= 0.05 and M = 500 epochs). The best detection rates (maximum percentage of volunteers for whom the response was detected for the frequencies between 4.8 and 72 Hz) were obtained for the parietal and central leads mid-sagittal and ipsilateral to the stimulated leg: C4 (87%), P4 (82%), Cz (89%), and Pz (89%). The P37-N45 time-components of the SEP can also be observed in these leads. The other leads, including the central and parietal contralateral and the frontal and fronto-polar leads, presented low detection capacity. If only contralateral leads were considered, the centro-parietal region (C3 and P3) was among the best regions for response detection, presenting a correspondent well-defined N37; however, this was not observed in some volunteers. The results of the present study showed that the central and parietal regions, especially sagittal and ipsilateral to the stimuli, presented the best SNR in the gamma range. Furthermore, these findings suggest that the MSC can be a useful tool for monitoring purposes.

Use of magnitude-squared coherence to identify the maximum driving response band of the somatosensory evoked potential.

The present study proposes to apply magnitude-squared coherence (MSC) to the somatosensory evoked potential for identifying the maximum driving response band. EEG signals, leads [Fpz'-Cz'] and [C3'-C4'], were collected from two groups of normal volunteers, stimulated at the rate of 4.91 (G1: 26 volunteers) and 5.13 Hz (G2: 18 volunteers). About 1400 stimuli were applied to the right tibial nerve at the motor threshold level. After applying the anti-aliasing filter, the signals were digitized and then further low-pass filtered (200 Hz, 6th order Butterworth and zero-phase). Based on the rejection of the null hypothesis of response absence (MSC(f) > 0.0060 with 500 epochs and the level of significance set at a = 0.05), the beta and gamma bands, 15-66 Hz, were identified as the maximum driving response band. Taking both leads together ("logical-OR detector", with a false-alarm rate of a = 0.05, and hence a = 0.0253 for each derivation), the detection exceeded 70% for all multiples of the stimulation frequency within this range. Similar performance was achieved for MSC of both leads but at 15, 25, 35, and 40 Hz. Moreover, the response was detected in [C3'-C4'] at 35.9 Hz and in [Fpz'-Cz'] at 46.2 Hz for all members of G2. Using the "logical-OR detector" procedure, the response was detected at the 7th multiple of the stimulation frequency for the series as a whole (considering both groups). Based on these findings, the MSC technique may be used for monitoring purposes.

Statistical analysis of event-related potential elicited by verb-complement merge in Brazilian Portuguese.

An interesting fact about language cognition is that stimulation involving incongruence in the merge operation between verb and complement has often been related to a negative event-related potential (ERP) of augmented amplitude and latency of ca. 400 ms - the N400. Using an automatic ERP latency and amplitude estimator to facilitate the recognition of waves with a low signal-to-noise ratio, the objective of the present study was to study the N400 statistically in 24 volunteers. Stimulation consisted of 80 experimental sentences (40 congruous and 40 incongruous), generated in Brazilian Portuguese, involving two distinct local verb-argument combinations (nominal object and pronominal object series). For each volunteer, the EEG was simultaneously acquired at 20 derivations, topographically localized according to the 10-20 International System. A computerized routine for automatic N400-peak marking (based on the ascendant zero-cross of the first waveform derivative) was applied to the estimated individual ERP waveform for congruous and incongruous sentences in both series for all ERP topographic derivations. Peak-to-peak N400 amplitude was significantly augmented (P < 0.05; one-sided Wilcoxon signed-rank test) due to incongruence in derivations F3, T3, C3, Cz, T5, P3, Pz, and P4 for nominal object series and in P3, Pz and P4 for pronominal object series. The results also indicated high inter-individual variability in ERP waveforms, suggesting that the usual procedure of grand averaging might not be considered a generally adequate approach. Hence, signal processing statistical techniques should be applied in neurolinguistic ERP studies allowing waveform analysis with low signal-to-noise ratio.

Use of magnitude-squared coherence to identify the maximum driving response band of the somatosensory evoked potential

The present study proposes to apply magnitude-squared coherence (MSC) to the somatosensory evoked potential for identifying the maximum driving response band. EEG signals, leads [Fpz'-Cz'] and [C3'-C4'], were collected from two groups of normal volunteers, stimulated at the rate of 4.91 (G1: 26 volunteers) and 5.13 Hz (G2: 18 volunteers). About 1400 stimuli were applied to the right tibial nerve at the motor threshold level. After applying the anti-aliasing filter, the signals were digitized and then further low-pass filtered (200 Hz, 6th order Butterworth and zero-phase). Based on the rejection of the null hypothesis of response absence (MSC(f) > 0.0060 with 500 epochs and the level of significance set at a = 0.05), the beta and gamma bands, 15-66 Hz, were identified as the maximum driving response band. Taking both leads together ("logical-OR detector", with a false-alarm rate of a = 0.05, and hence a = 0.0253 for each derivation), the detection exceeded 70 percent for all multiples of the stimulation frequency within this range. Similar performance was achieved for MSC of both leads but at 15, 25, 35, and 40 Hz. Moreover, the response was detected in [C3'-C4'] at 35.9 Hz and in [Fpz'-Cz'] at 46.2 Hz for all members of G2. Using the "logical-OR detector" procedure, the response was detected at the 7th multiple of the stimulation frequency for the series as a whole (considering both groups). Based on these findings, the MSC technique may be used for monitoring purposes.

Statistical analysis of event-related potential elicited by verb-complement merge in Brazilian Portuguese

An interesting fact about language cognition is that stimulation involving incongruence in the merge operation between verb and complement has often been related to a negative event-related potential (ERP) of augmented amplitude and latency of ca. 400 ms - the N400. Using an automatic ERP latency and amplitude estimator to facilitate the recognition of waves with a low signal-to-noise ratio, the objective of the present study was to study the N400 statistically in 24 volunteers. Stimulation consisted of 80 experimental sentences (40 congruous and 40 incongruous), generated in Brazilian Portuguese, involving two distinct local verb-argument combinations (nominal object and pronominal object series). For each volunteer, the EEG was simultaneously acquired at 20 derivations, topographically localized according to the 10-20 International System. A computerized routine for automatic N400-peak marking (based on the ascendant zero-cross of the first waveform derivative) was applied to the estimated individual ERP waveform for congruous and incongruous sentences in both series for all ERP topographic derivations. Peak-to-peak N400 amplitude was significantly augmented (P < 0.05; one-sided Wilcoxon signed-rank test) due to incongruence in derivations F3, T3, C3, Cz, T5, P3, Pz, and P4 for nominal object series and in P3, Pz and P4 for pronominal object series. The results also indicated high inter-individual variability in ERP waveforms, suggesting that the usual procedure of grand averaging might not be considered a generally adequate approach. Hence, signal processing statistical techniques should be applied in neurolinguistic ERP studies allowing waveform analysis with low signal-to-noise ratio.

Topographic aspects of photic driving in the electroencephalogram of children and adolescents.

The electroencephalogram amplitude spectra at 11 fixed frequencies of intermittent photic stimulation of 3 to 24 Hz were combined into driving "profiles" for 14 scalp points in 8 male and 7 female normal subjects aged 9 to 17 years. The driving response varied over frequency and was detected in 70 to 100% of cases in the occipital areas (maximum) and in 27 to 77% of cases in the frontal areas (minimum) using as a criterion peak amplitude 20% higher than those of the neighbors. Each subject responded, on average, to 9.7 +/- 1.15 intermittent photic stimulation frequencies in the right occipital area and to 6.8 +/- 1.97 frequencies in the right frontal area. Most of the driving responses (in relation to the previous background) were significant according to the spectral F-test (alpha = 0.05), which also detected changes in some cases of low amplitude responses not revealed by the peak criterion. The profiles had two maxima in the alpha and theta bands in all leads. The latter was not present in the background spectra in the posterior areas and was less pronounced in the anterior ones. The weight of the profile theta maximum increased towards the frontal areas where the two maxima were similar, while the profile amplitudes decreased. The profiles repeated the shape of the background spectra, except for the theta band. The interhemispheric correlation between profiles was high. The theta driving detected in all areas recorded suggests a generalized influence of the theta generators in prepubertal and pubertal subjects.

Topographic aspects of photic driving in the electroencephalogram of children and adolescents

The electroencephalogram amplitude spectra at 11 fixed frequencies of intermittent photic stimulation of 3 to 24 Hz were combined into driving "profiles" for 14 scalp points in 8 male and 7 female normal subjects aged 9 to 17 years. The driving response varied over frequency and was detected in 70 to 100 percent of cases in the occipital areas (maximum) and in 27 to 77 percent of cases in the frontal areas (minimum) using as a criterion peak amplitude 20 percent higher than those of the neighbors. Each subject responded, on average, to 9.7 ñ 1.15 intermittent photic stimulation frequencies in the right occipital area and to 6.8 ñ 1.97 frequencies in the right frontal area. Most of the driving responses (in relation to the previous background) were significant according to the spectral F-test (a = 0.05), which also detected changes in some cases of low amplitude responses not revealed by the peak criterion. The profiles had two maxima in the alpha and theta bands in all leads. The latter was not present in the background spectra in the posterior areas and was less pronounced in the anterior ones. The weight of the profile theta maximum increased towards the frontal areas where the two maxima were similar, while the profile amplitudes decreased. The profiles repeated the shape of the background spectra, except for the theta band. The interhemispheric correlation between profiles was high. The theta driving detected in all areas recorded suggests a generalized influence of the theta generators in prepubertal and pubertal subjects.

Modeling extracellular space electrodiffusion during Leão's spreading depression.

Computational modeling of spreading depression (SD) has been used increasingly to study the different mechanisms that are involved in this phenomenon. One of them that is still under discussion involves the mechanisms that originate the extracellular electrical field responsible for the dc potential shift. The main goal of this paper is to present a mathematical derivation for the extracellular electric field that is incorporated in a SD model that has the basic structure of Tuckwell and Miura's model, but with the ionic variations calculated electrochemically. Electrodiffusion equations were used to describe the ionic movement of the four ions Na+, K+, Cl-, and Ca2+. These are mutually coupled by the electric field within the extracellular space (ECS). The results from the simulations show that the model is able to calculate the effect of the ionic changes along the ECS on the electric field, and to reproduce the SD in respect to the most important features that characterize the phenomenon experimentally in the retina or hippocampus. It is suggested that the extracellular negative field-potential shift during SD is due to an electrical field generated by a Goldman-Hodgkin-Katz equation acting within the ECS.

Functional imaging of the retinal layers by laser scattering: an approach for the study of Leão's spreading depression in intact tissue.

This paper presents a novel optical approach for the study of spreading depression in isolated retina. The method makes it possible to register the laser light scattered from each layer of the tissue, yielding a functional image of the retina during spreading depression. The tissue is kept intact, since histological cuts are not necessary. Measurements of other variables, such as extracellular potential, are also allowed by the described method. This is done simultaneously with the functional image in a high spatial resolution, with the positioning of the microelectrode tip being easily monitored. The information about temporal and spatial evolution of light was compacted in a single image. The image-processing technique used here enables the visualization of the light scattered by the inner plexiform layer (IPL), which is the most prominent scatter layer during spreading depression. The wavefront velocity and its increase as two wavefronts approach each other can then be determined, and it is also possible to observe the thickness variation of the tissue during the wave travel. The relationship between two peaks of light-scattering sequence during the phenomenon was studied at two wavelengths (632.8 and 543.5 nm). This relationship is shown to be dependent on the wavelength.

Coherence between one random and one periodic signal for measuring the strength of responses in the electro-encephalogram during sensory stimulation.

Coherence between a pulse train representing periodic stimuli and the EEG has been used in the objective detection of steady-state evoked potentials. This work aimed to quantify the strength of the stimulus responses based on the statistics of coherence estimate between one random and one periodic signal, focusing on the confidence limits and power of significance tests in detecting responses. To detect the responses in 95% of cases, a signal-to-noise ratio of about -7.9 dB was required when using 48 windows (M) in the coherence estimation. The ratio, however, increased to -1.2 dB when M was 12. The results were tested in Monte Carlo simulations and applied to EEGs obtained from 14 subjects during visual stimulation. The method showed differences in the strength of responses at the stimulus frequency and its harmonics, as well as variations between individuals and over cortical regions. In contrast to those from the parietal and temporal regions, results for the occipital region gave confidence limits (with M = 12) that were above zero for all subjects, indicating statistically significant responses. The proposed technique extends the usefulness of coherence as a measure of stimulus responses and allows statistical analysis that could also be applied usefully in a range of other biological signals.

Detecção da resposta auditiva no EEG de crianças utilizando técnicas no domínio da freqüência / Detection of EEG auditory response in children, using frequency domain techniques

A triagem auditiva, especialmente de recém nascidos (RN) que apresentem fatores de risco, tem sido recomendada por diversas entidades internacionais. O exame de potencial evocado auditivo de tronco cerebral (BERA - Brainstem Evoked Responses Audiometry) é atualmente considerado como referência para detecção do limiar auditivo em RNs, embora seja um procedimento demorado e de interpretação subjetiva. Para detecção objetiva da resposta auditiva em crianças foram aplicadas duas técnicas que utilizam a Transformada de Fourier: o CSM (Component Synchrony Meadure) e o MSC (Magnitude Squared Coherence). Sinais EEG de 32 crianças foram adquiridos durante BERA com estimulação por clicks a 20 Hz. Colheram-se ao todo 73 sinais EEG, cada um abrangendo 1024 estímulos de 85 dBHL. Destes 73 sinais, a interpretação clínica dos potenciais evocados indicou haver resposta (audição intacta) em 55. Um critério heurístico de detecção de resposta baseado em um conjunto de harmônicos da freqüência de estimulação (20 a 900 Hz) resultou em sensibilidade de 93 por cento e especificidade de 78 por cento para o MSC, e 98 por cento e 72 por cento para o CSM, respectivamente. Como o desempenho na detecção da surdez (resultado negativo) é o desejado, MSC seria mais adequado à detecção de resposta auditiva no EEG de crianças e recém-nascidos.

Monitorização do plano anestésico por técnicas estatísticas no domínio da freqüência / Monitoring depth of anaesthesia by frequency-domain statistical techniques

O potencial evocado auditivo de média latência (MLAEP) tem sido empregado na avaliação do nível anestésico, sendo usual a análise morfológica (AEPidx - Indice de Potencial Evocado Auditivo, e DT - Distância Taxonômica). Neste trabalho, técnicas de detecção objetiva de respostas no domínio da freqüência (MSC - Magnitude Quadrada da Coerência, e CSM - Medida de Sincronismo de Componentes) foram avaliadas a partir de EEG sob estimulação auditiva. Dez voluntários foram estimulados com cliques repetitivos (7,645 Hz, nível de pressão sonora de 75 dB SPL), sendo EEG adquirido na derivação Vértex - Lobo Auricular durante 10 min sem administração de anestésico e dois períodos consecutivos de infusão contínua de propofol (3 mg/Kg/h, sedação, e 9 mg/Kg/h, anestesia) de 45 min cada. O instante de perda de consciência foi estimado com base no pressionamento de uma maniqueta em resposta a estímulos luminosos a cada 15 s. MSC e CSM permitem a inferência estatística de ausência de resposta auditiva (valor crítico, para e.g. alfa = 1 por cento) ao longo do processo anestésico. Além disso, MSC(Ÿ) e CSM(Ÿ) se reduzem durante a infusão (mais expressivamente na anestesia), tendo maior concordância (diagramas de Bland-Altman) com os resultados da maniqueta que AEPidx e DT, e maior equivalência (Wilcoxon Signed-Rank Test, alfa = 5 por cento) que DT.

Radiologia diagnóstica e programa de garantia de qualidade: avaliaçäo crítica / Diagnostic Radiology and quality assurance program: critical evaluation

O Programa de Garantia de Qualidade, implementado no Serviço de Radiologia do Hospital da Polícia Militar - RJ, propiciou a reduçäo de 70 por cento na taxa de rejeiçäo de filmes e de 75 por cento nos custos. Entretanto, a descontinuidade do Programa tem conduzido ao retorno a situaçäo anterior a implementaçäo.

A Quality Assurance Program, implemented at the Radiology Service of the Hospital of the Military Police in Rio de Janeiro, brought about a reduction of 70% in the number of rejected films and 75% in the costs. The discontinuation of the program has however lead to a return to the former conditions.

Avaliando serviços de radiologia diagnóstica / Evaluate service of the diagnostic radiology

Radiação ionizante tem sido, cada vez mais, utilizada no diagnóstico e tratamento de doenças, sem a adequada adoção de métodos que garantem a qualidade do exame como um todo. Um Sistema de Garantia de Qualidade que vise reverter este cenário, depende das características de cada unidade, na qual Programas de Qualidade serão implementados. A avaliação preliminar de Serviços de Radiologia Diagnóstica, no que concerne a complexidade, oferta e qualidade dos serviços e disponibilidade de pssoal, constataram a utilização de técnicas de baixa voltagem, a não realização de testes de aceitação de aprelhos e acessórios e a não existência de profissional responsável pela proteção radiológica, entre outras limitações e deficiências.

The use of ionizing radiation is steadily increasing in diagnosis and treatment of diseases, in many cases however without the adoption of suitable methods that would garantee the quality of the exam as a whole. Improvements in this respect must be based on the present characteristics of the units in which Quality Control is to be introduced. Preliminary evaluation of Diagnostic Radiology Services showed, among other problems, the preferential use of low-voltage X-rays, the absence of acceptance triais of equipment and acessories and the lack of professionals responsible for radiological protection.

Análise espectral e de coerência de padröes do siclo sono-vígilia do EEG neonatal / Spectral analysis of pastterns coherence of sllep-awake cicle of neonatal ECG

O espectro de potência e a coerência inter-hemisférica foram estimados para segmentos dos padrões do ciclo sono-vigélia Low Voltage Irregular, Misto e High Voltage Slow. Diferenças significativas entre estes padrões foram encontradas, através da ANOVA, em uma larga faixa de frequência, principalmente para o espectro de potência até 15 Hz.

Modelagem AR de sinais EEG multicanal usando o TMS320C30 / AR modelling of multicanal EEG signals using the TMS320C30

Uma placa processadora digital de sinais (DSP), com o TMS320C30, em conjunto com um PC foi empregada na modelagem autoregressiva (AR) de sinais EEG. O sistema PC-DSP mostrou-se capaz de realizar a análise espectral destes sinais com tempo de processamento compatível com monitoração em tempo real do espectro de 16 canais EEG.