2.4.1 실험
대학병원 소아과 진료실에서 1~4회 내방한 40여명의 소아를 대상으로 심음 데이터를 수집하였고, 이중에서 질환 치료를 통해 정상이 되었거나, 잡음이
심한 경우 등의 데이터를 제외하고, 본 연구의 연령별, 질환별 관찰을 위한 15인의 소아를 대상으로 분석 실험을 진행하였고, 이는 표 2에 나타내었다.
심음 데이터 수집 시스템은 필터 및 심음 증폭기가 내장된 전자청진기를 제작하여 사용하였고, 노트북 컴퓨터로 구성되어 있다(6). 수집된 소아 심음은 증폭 및 필터링되고 사운드 카드의 A/D 변환기를 통해서 8[kHz]의 샘플링 주파수와 16비트 모노의 해상도로 웨이브 파일로
디지털 변환되어 저장된다. 청진음 데이터에서 음성, 기계음, 피부마찰음 등 아티팩트(artifact)는 잘라내고, 약 2 주기의 심박이 포함되는 1초를
분석 윈도우(window) 크기로 설정하였다.
표 2. 피검자 정보
Table 2 Subject information
피검자
|
나이
|
성별
|
진단
|
JDK
|
0.41
|
M
|
small VSD
|
BJW
|
0.83
|
M
|
functional murmur
|
SKJ
|
0.91
|
M
|
VSD
|
KMJ
|
1.00
|
F
|
VSD
|
KTS
|
1.50
|
M
|
ASD
|
KSM
|
1.50
|
M
|
small ASD
|
JES
|
2.00
|
F
|
정상(functional murmur)
|
LDM
|
2.00
|
M
|
VSD
|
PDJ
|
2.33
|
M
|
정상
|
CHH
|
2.83
|
M
|
정상
|
YHS
|
3.00
|
M
|
정상
|
LHJ
|
3.16
|
M
|
VSD
|
PHS
|
3.75
|
F
|
정상
|
SSH
|
4.00
|
M
|
정상
|
LSJ
|
4.00
|
F
|
functional murmur
|
심음의 선형 예측 분석을 통해 구한 특징 파라미터의 성능을 분석하기 위하여 심음의 AR 스펙트럼을 구하고, 특징 파라미터로는 심음 발생 모델의 전달함수
극점에 해당하는 형성음 주파수를 추출하였다. 이를 위한 실험 절차는 다음 그림과 같다.
그림 3 형성음 검출 실험 절차
Fig. 3 Procedure for formant detection
디지털화된 심음은 격자 필터 구조의 10차 Burg 알고리즘 으로 선형 예측 분석이 이루어지고, AR 모델의 스펙트럼 추정법에 의해 스펙트럼을 추출한다.
스펙트럼 파형의 1차 도함수를 이용해 기울기가 양에서 음으로 변하는 지점을 검출하여 1단계의 형성음 주파수를 구한다. 추출된 스팩트럼을 이동 평균과
2차 도함수를 이용하여 0교차를 검출하여 2단계로 형성음을 구하여 결과를 취합하여 형성음 주파수를 결정한다.
2.4.2 결과 고찰
그림 4~7은 정상 소아의 스펙트럼 파형의 예를 1세부터 4세까지 나타낸 것이고, 그림 8~10은 세 질환 소아의 파형을 나타낸 것이다. 그림에서 파형과 수평선 사이에 저주파수 쪽에서부터 검출된 형성음 주파수를 수직선으로 나타내었다.
그림 4 1세 피검자(VSD)의 스펙트럼의 예 : 피검자 KMJ
Fig. 4 Example of the spectra for 1-year old child(VSD) : subject KMJ
그림 5 2세 피검자(정상)의 스펙트럼의 예 : 피검자 JES
Fig. 5 Example of the spectra for 2-year old child(normal) : subject JES
그림 6 3세 피검자(정상)의 스펙트럼의 예 : 피검자 YHS
Fig. 6 Example of the spectra for 3-year old child(normal) : subject YHS
그림 7 4세 피검자(정상)의 스펙트럼의 예 : 피검자 SSH
Fig. 7 Example of the spectra for 4-year old child(normal) : subject SSH
형성음 주파수는 극점에 해당하는 위치이지만, 영점의 영향 등으로 검출이 되지 않는 경우도 빈번히 있으며, 제1, 2 형성음 주파수(F1, F2)는
저주파 부근에서 비교적 일정하게 검출되나, 제3~5 형성음 주파수(F3~F5)는 선형 예측 분석의 특성상 위치가 변동이 심하거나, 경우에 따라 검출되지
않기도 한다.
그림 8 VSD 질환자의 스펙트럼의 예 : 피검자 LDM
Fig. 8 Example of the spectra for VSD patient : subject LDM
그림 9 ASD 질환자의 스펙트럼의 예 : 피검자 KTS
Fig. 9 Example of the spectra for ASD patient : subject KTS
그림에서 질환별 스펙트럼마다 유사해 보이지만 약간씩 차이가 있다는 것을 알 수 있다. VSD 질환 스펙트럼의 경우 F1 주변 경사가 급하고, 이후에도
굴곡(F2, F3등)이 있다. 그리고 정상보다는 F1-F2 간격이 좁은 것을 알 수 있다. ASD 질병 스펙트럼의 경우에는 F1 주변이 정상과 비슷하게
주변에서 완만하고, 전체적으로도 파형이 완만함을 볼 수 있으며, F1-F2 간격은 VSD 질병처럼 간격이 좁은 것을 알 수가 있다. 기능성 심잡음
질환 어린이의 경우에는 스펙트럼의 형태가 정상과 거의 비슷한데 특정 프레임에서 F3 이후의 고주파 형성음이 크고 뚜렷하게 나타나는 경우가 나타난다.
그림 10 기능성심잡음 질환자의 스펙트럼의 예 : 피검자 BJW
Fig. 10 Example of the spectra for functional murmur patient : subject BJW
표 3은 소아의 5개 프레임에서 추출된 형성음주파수 평균값을 연령별로 나타낸 것이고, 표 4는 같은 방법으로 질병별로 결과를 나타낸 것이며, 단위는 Hz이다. F1, F2,...는 각각 제1 형성음주파수, 제2 형성음주파수, ...를 의미한다.
표 3. 소아 연령별 형성음 주파수
Table 3 Formant frequencies for child age
(단위: Hz)
subj.
|
age
|
F1
|
F2
|
F3
|
F4
|
F5
|
JDK
|
0.41
|
274.8
|
720.6
|
1610.2
|
2753.3
|
3487.5
|
BJW
|
0.83
|
237.4
|
824.2
|
1732.3
|
2668.5
|
3293.5
|
SKJ
|
0.91
|
342.4
|
746.0
|
1479.8
|
2308.0
|
3293.0
|
KMJ
|
1.00
|
249.4
|
849.0
|
1360.8
|
2660.0
|
3513.6
|
average
|
276.0
|
785.0
|
1545.8
|
2597.5
|
3396.9
|
KTS
|
1.50
|
229.4
|
604.0
|
1567.7
|
2525.8
|
3423.3
|
KSM
|
1.50
|
233.0
|
734.3
|
1637.6
|
2714.8
|
3699.5
|
JES
|
2.00
|
243.2
|
795.2
|
1671.3
|
2353.4
|
3521.8
|
LDM
|
2.00
|
307.4
|
658.8
|
×
|
2517.2
|
3670.0
|
average
|
253.3
|
698.1
|
1625.5
|
2527.8
|
3578.7
|
PDJ
|
2.33
|
207.6
|
577.0
|
1341.8
|
2337.5
|
3284.8
|
CHH
|
2.83
|
204.2
|
602.7
|
1669.8
|
2569.5
|
3404.6
|
YHS
|
3.00
|
239.8
|
748.0
|
1683.0
|
2470.4
|
3609.5
|
average
|
217.2
|
642.6
|
1564.9
|
2459.1
|
3433.0
|
LHJ
|
3.16
|
246.0
|
650.5
|
1745.0
|
2650.0
|
3409.3
|
PHS
|
3.75
|
188.0
|
757.4
|
1787.6
|
2677.4
|
3529.3
|
SSH
|
4.00
|
331.0
|
944.8
|
1768.0
|
2617.3
|
3468.6
|
LSJ
|
4.00
|
199.0
|
798.5
|
1787.3
|
2622.0
|
3425.7
|
average
|
241.0
|
787.8
|
1772.0
|
2641.7
|
3458.2
|
AR 모델링의 스펙트럼 특성상 형성음 주파수 위치의 바이어스가 나타나며, 고주파로 갈수록 위치의 편차가 크다. 따라서 F1, F2가 F3 이후의 것들보다
중요한 파라미터로 볼 수 있을 것이다. 표 3에서 보면 연령이 높아짐에 따라 F1, F2 등의 평균 주파수 값이 낮아지는 것을 볼 수 있다. 이는 성장하면서 체격이 커짐에 따라 심박도 느려지고,
따라서 형성음 주파수도 값이 낮아지는 것으로 추정된다. 그러나 피검자 SKJ, LDM, SSH 등은 F1이 심하게 높아서 상대적으로 성장이 느려 체격이
작았을 것으로 예측할 수 있다. 표 4의 질환별 형성음 주파수 평균치를 살펴 보면, F3 이후의 주파수 값은 특별한 차이를 나타내지 않지만, 위에서 언급한 특별히 편차가 있는 어린이를
제외하면 F1이 VSD 질환 어린이의 주파수 평균값이 높은 경향을 나타내는 것을 알 수 있다. ASD와 심잡음 어린이의 경우는 정상 소아와 큰 차이가
없는 것으로 나타남을 알 수 있다. 심잡음의 경우 정상 소아의 경우에도 빈번히 나타나고, 성장하면서 사라지기도 한다.
또한 정상 어린이의 형성음 주파수 평균값의 경우 F1과 F2의 간격이 482Hz이며, VSD와 ASD보다 약 42Hz 넓고, 기능성심잡음보다 약 80Hz
작게 나타난 것을 알 수가 있다. VSD 질환의 평균값의 경우 평균적으로 다른 질병보다 F1이 284Hz로 높은 주파수를 갖으나, F3은 상대적으로
낮은 주파수를 갖음을 알 수가 있다. ASD 질환의 평균값의 경우 F1은 정상과 거의 비슷한 주파수 분포를 갖으며, F1-F2 간격은 VSD 질환과
유사하게 나타났다.
표 4 소아 심장 질환별 형성음 주파수
Table 4 Formant frequencies for child heart disease
(단위: Hz)
class
|
subj.
|
age
|
F1
|
F2
|
F3
|
F4
|
F5
|
Nor-mal
|
JES
|
2.00
|
243.2
|
795.2
|
1,671.3
|
2,353.4
|
3,521.8
|
PDJ
|
2.33
|
207.6
|
577.0
|
1,341.8
|
2,337.5
|
3,284.8
|
CHH
|
2.83
|
204.2
|
602.7
|
1,669.8
|
2,569.5
|
3,404.6
|
YHS
|
3.00
|
239.8
|
748.0
|
1,683.0
|
2,470.4
|
3,609.5
|
PHS
|
3.75
|
188.0
|
757.4
|
1,787.6
|
2,677.4
|
3,529.3
|
SSH
|
4.00
|
331.0
|
944.8
|
1,768.0
|
2,617.3
|
3,468.6
|
average
|
235.6
|
737.5
|
1653.6
|
2504.3
|
3469.8
|
VSD
|
JDK
|
0.41
|
274.8
|
720.6
|
1,610.2
|
2,753.3
|
3,487.5
|
SKJ
|
0.91
|
342.4
|
746.0
|
1,479.8
|
2,308.0
|
3,293.0
|
KMJ
|
1.00
|
249.4
|
849.0
|
1,360.8
|
2,660.0
|
3,513.6
|
LDM
|
2.00
|
307.4
|
658.8
|
×
|
2,517.2
|
3,670.0
|
LHJ
|
3.16
|
246.0
|
650.5
|
1,745
|
2,650.0
|
3,409.3
|
average
|
284.0
|
725.0
|
1549.0
|
2577.7
|
3474.7
|
ASD
|
KTS
|
1.50
|
229.4
|
604.0
|
1,567.7
|
2,525.8
|
3,423.3
|
KSM
|
1.50
|
233.0
|
734.3
|
1,637.6
|
2,714.8
|
3,699.5
|
average
|
231.2
|
669.2
|
1602.7
|
2620.3
|
3561.4
|
Func. mur.
|
BJW
|
0.83
|
237.4
|
824.2
|
1,732.3
|
2,668.5
|
3,293.5
|
LSJ
|
4.00
|
199.0
|
798.5
|
1,787.3
|
2,622.0
|
3,425.7
|
average
|
218.2
|
811.4
|
1759.8
|
2645.3
|
3359.6
|