Nghiên Cứu Phương Pháp Khử Nhiễu Tín Hiệu Điện Tâm Đồ Bằng Wavelet Packet

Luận văn tốt nghiệp giáo dục nghiên cứu tốt nghiệp sư phạm vật lý cải tiến xử lý nhiễu của tín hiệu điện tim bằng phép biến đổi wavelet, điều tra thực trạng, phân tích số liệu, đề

• Tác giả: Nguyễn Trung Hiếu
• Người hướng dẫn: ThS. Trần Đặng Bảo Ân
• Trường học: Trường Đại học Sư phạm Thành phố Hồ Chí Minh
• Chuyên ngành: Sư phạm Vật lý
• Đề tài: Nghiên Cứu Phương Pháp Khử Nhiễu Tín Hiệu Điện Tâm Đồ Bằng Wavelet Packet
• Loại tài liệu: luận văn tốt nghiệp
• Năm xuất bản: 2016
• Địa điểm: Tp. Hồ Chí Minh

Mục lục chi tiết

• LỜI CAM ĐOAN
• LỜI CẢM ƠN
• MỤC LỤC
• DANH MỤC MỘT SỐ CHỮ VIẾT TẮT
• DANH MỤC HÌNH VẼ
• DANH MỤC BẢNG BIỂU
• 1. CHƯƠNG 1: TỔNG QUAN
• 1.1. Cơ sở lý thuyết tín hiệu điện tim
• 1.2. Cấu trúc giải phẫu và chức năng của tim
• 1.3. Quá trình hình thành tín hiệu điện tim
• 1.4. Nhiễu Gauss trắng trong tín hiệu điện tim
• 2. CHƯƠNG 2: CÁC PHƯƠNG PHÁP XỬ LÝ TÍN HIỆU
• 2.1. Phép biến đổi Fourier (FT - Fourier Transform)
• 2.2. Phép biến đổi Fourier liên tục (CFT — Continous Fourier Transform)
• 2.3. Phép biến đổi Fourier thời gian ngắn (STFT - Short Time Fourier Transform)
• 2.4. Phép biến đổi Wavelet (WT — Wavelet Transform)
• 2.5. Phép biến đổi Wavelet liên tục (CWT - Continous Wavelet Transform)
• 2.6. Phép biến đổi Wavelet rời rạc (DWT — Discrete Wavelet Transform)
• 2.7. Phép biến đổi Wavelet Packet (WPT)
• 3. CHƯƠNG 3: PHƯƠNG PHÁP THỰC HIỆN
• 3.1. Xử lý nhiễu tín hiệu bằng phép biến đổi Wavelet Packet
• 3.2. Mô hình xử lý nhiễu cơ bản
• 3.3. Nguyên tắc khử nhiễu
• 3.4. Phương pháp đặt ngưỡng
• 3.5. Lý thuyết ngưỡng và quy tắc chọn ngưỡng
• 3.6. Đánh giá chất lượng tín hiệu
• 3.7. Khảo sát quy tắc entropy tối ưu
• 4. CHƯƠNG 4: KẾT QUẢ
• 4.1. Đánh giá entropy tối ưu trong phép biến đổi Wavelet Packet
• 4.2. Đánh giá phép biến đổi Wavelet Packet và phép biến đổi Wavelet rời rạc
• 4.3. Về hiệu quả khử nhiễu
• 4.4. Về độ phân giải tín hiệu
• 4.5. Xây dựng giao diện xử lý nhiễu tín hiệu điện tim bằng phép biến đổi Wavelet Packet
• 4.5.1. Thiết kế giao diện
• 4.5.2. Thực thi hàm trong GUIDE
• 4.5.3. Sử dụng giao diện xử lý nhiễu tín hiệu điện tim bằng phép biến đổi Wavelet Packet
• 4.5.4. Đánh giá giao diện xử lý nhiễu tín hiệu điện tim bằng phép biến đổi Wavelet Packet
• 4.5.5. Xây dựng giao diện ngân hàng tín hiệu điện tim
• 4.5.6. Tiền xử lý tín hiệu
• 4.5.7. Nhập dữ liệu (Thực thi hàm trong GUIDE)
• KẾT LUẬN VÀ HƯỚNG PHÁT TRIỂN
• DANH MỤC CÔNG TRÌNH ĐÃ CÔNG BỐ
• TÀI LIỆU THAM KHẢO

Tóm tắt

I. Tổng quan về Nghiên Cứu Phương Pháp Khử Nhiễu Tín Hiệu Điện Tâm Đồ Bằng Wavelet Packet

Nghiên cứu về khử nhiễu tín hiệu điện tâm đồ (ECG) là một lĩnh vực quan trọng trong y học hiện đại. Tín hiệu ECG chứa đựng thông tin quý giá về hoạt động của tim, nhưng thường bị ảnh hưởng bởi nhiễu từ nhiều nguồn khác nhau. Phương pháp Wavelet Packet đã được chứng minh là một công cụ hiệu quả trong việc xử lý và khử nhiễu tín hiệu này. Bài viết sẽ trình bày tổng quan về phương pháp này và những ứng dụng của nó trong thực tiễn.

1.1. Khái niệm về tín hiệu điện tâm đồ và tầm quan trọng của nó

Tín hiệu điện tâm đồ (ECG) là các xung điện phản ánh hoạt động của tim. Việc phân tích tín hiệu này giúp chẩn đoán các bệnh lý tim mạch. Tuy nhiên, nhiễu tín hiệu có thể làm giảm độ chính xác trong chẩn đoán.

1.2. Giới thiệu về phương pháp Wavelet Packet

Phương pháp Wavelet Packet cho phép phân tích tín hiệu ở nhiều mức độ khác nhau, giúp tách biệt nhiễu và tín hiệu hữu ích. Điều này làm tăng khả năng khử nhiễu và giữ lại thông tin quan trọng trong tín hiệu ECG.

II. Vấn đề và Thách thức trong Khử Nhiễu Tín Hiệu Điện Tâm Đồ

Khử nhiễu tín hiệu điện tâm đồ gặp nhiều thách thức do tính chất phức tạp của tín hiệu và sự đa dạng của nhiễu. Các yếu tố như sự thay đổi tần số nguồn điện, tiếp xúc không tốt giữa điện cực và da, và các yếu tố sinh lý khác có thể gây ra nhiễu tín hiệu. Việc hiểu rõ các vấn đề này là cần thiết để phát triển các phương pháp khử nhiễu hiệu quả.

2.1. Các nguồn gốc gây nhiễu trong tín hiệu ECG

Nhiễu trong tín hiệu ECG có thể đến từ nhiều nguồn như điện từ trường, chuyển động cơ thể, và chất lượng tiếp xúc của điện cực. Những yếu tố này làm cho việc phân tích tín hiệu trở nên khó khăn hơn.

2.2. Hạn chế của các phương pháp khử nhiễu hiện tại

Các phương pháp khử nhiễu truyền thống như biến đổi Fourier có những hạn chế trong việc xử lý tín hiệu không ổn định như ECG. Điều này dẫn đến việc cần thiết phải phát triển các phương pháp mới như Wavelet Packet.

III. Phương Pháp Khử Nhiễu Tín Hiệu Điện Tâm Đồ Bằng Wavelet Packet

Phương pháp Wavelet Packet được áp dụng để khử nhiễu tín hiệu ECG bằng cách phân tích tín hiệu thành các thành phần khác nhau. Phương pháp này cho phép tách biệt nhiễu và tín hiệu hữu ích một cách hiệu quả hơn so với các phương pháp truyền thống.

3.1. Nguyên tắc hoạt động của Wavelet Packet

Phương pháp Wavelet Packet sử dụng các hàm wavelet để phân tích tín hiệu ở nhiều mức độ khác nhau, từ đó giúp tách biệt nhiễu và tín hiệu chính xác hơn.

3.2. Quy trình khử nhiễu bằng Wavelet Packet

Quy trình khử nhiễu bao gồm các bước như phân tích tín hiệu, xác định ngưỡng khử nhiễu, và tái tạo tín hiệu. Mỗi bước đều quan trọng để đảm bảo chất lượng tín hiệu sau khi khử nhiễu.

IV. Ứng Dụng Thực Tiễn và Kết Quả Nghiên Cứu

Kết quả từ việc áp dụng phương pháp Wavelet Packet cho thấy khả năng khử nhiễu tốt hơn so với các phương pháp khác. Nghiên cứu đã chỉ ra rằng phương pháp này không chỉ giúp cải thiện chất lượng tín hiệu mà còn giữ lại các thông tin lâm sàng quan trọng.

4.1. Đánh giá hiệu quả khử nhiễu

Các chỉ số như tỷ lệ tín hiệu trên nhiễu (SNR) và độ lệch bình quân (MSE) cho thấy phương pháp Wavelet Packet có hiệu quả cao trong việc khử nhiễu tín hiệu ECG.

4.2. Ứng dụng trong chẩn đoán y tế

Kết quả từ nghiên cứu cho thấy tín hiệu sau khử nhiễu có thể được sử dụng trong chẩn đoán và điều trị các bệnh lý tim mạch, góp phần nâng cao chất lượng chăm sóc sức khỏe.

V. Kết Luận và Tương Lai của Nghiên Cứu Khử Nhiễu Tín Hiệu Điện Tâm Đồ

Nghiên cứu về phương pháp khử nhiễu tín hiệu điện tâm đồ bằng Wavelet Packet đã mở ra nhiều hướng đi mới trong việc cải thiện chất lượng tín hiệu. Tương lai của nghiên cứu này có thể bao gồm việc phát triển các thuật toán khử nhiễu tiên tiến hơn và ứng dụng trong các lĩnh vực khác.

5.1. Tóm tắt kết quả nghiên cứu

Kết quả nghiên cứu cho thấy phương pháp Wavelet Packet có khả năng khử nhiễu hiệu quả, giữ lại thông tin lâm sàng quan trọng trong tín hiệu ECG.

5.2. Hướng phát triển trong tương lai

Nghiên cứu có thể mở rộng sang các lĩnh vực khác như xử lý tín hiệu y sinh, từ đó nâng cao khả năng chẩn đoán và điều trị bệnh lý tim mạch.

Nội dung chính

Tổng quan nghiên cứu

Bệnh lý tim mạch hiện là nguyên nhân gây tử vong hàng đầu trên toàn cầu. Thống kê từ Tổ chức Y tế Thế giới chỉ ra rằng mỗi năm có khoảng 17,5 triệu người tử vong do các biến chứng tim mạch, và con số này được dự báo vượt mức 25 triệu người. Trong thực hành y khoa, điện tâm đồ là phương pháp chẩn đoán không xâm lấn hàng đầu nhờ chi phí thấp và cung cấp dữ liệu lâm sàng quan trọng. Tín hiệu điện tim phản ánh hoạt động điện học của các buồng tim với biên độ rất nhỏ, dao động trong khoảng từ 0,1 mV đến 5 mV và phân bố chủ yếu ở dải tần số từ 0 Hz đến 100 Hz.

Tuy nhiên, quá trình ghi nhận điện tâm đồ thường xuyên bị ảnh hưởng bởi nhiều nguồn tạp âm phức tạp như nhiễu nguồn điện 50 Hz hoặc 60 Hz, nhiễu điện cơ do rung cơ, và hiện tượng trôi đường đẳng điện do tiếp xúc điện cực không ổn định. Nhiễu Gauss trắng với tỷ lệ tín hiệu trên nhiễu nằm trong khoảng từ 0 dB đến 18 dB thường làm lu mờ các đỉnh sóng quan trọng, dẫn đến nguy cơ chẩn đoán sai lệch các bệnh lý nguy hiểm như nhồi máu cơ tim hay rối loạn nhịp tim. Các phương pháp xử lý truyền thống như biến đổi Fourier hoặc biến đổi Wavelet rời rạc bộc lộ nhược điểm lớn khi xử lý tín hiệu không dừng, làm suy giảm biên độ đỉnh sóng và xóa mất thông tin tần số cao.

Đề tài tập trung nghiên cứu giải pháp cải tiến xử lý nhiễu điện tâm đồ bằng phép biến đổi Wavelet Packet. Mục tiêu cụ thể là khảo sát quy tắc entropy tối ưu, xác định hàm Wavelet mẹ và mức phân tách phù hợp nhất, đồng thời xây dựng giao diện xử lý trực quan và ngân hàng dữ liệu tín hiệu điện tim. Nghiên cứu được triển khai thực nghiệm trên các bản ghi tín hiệu chuẩn từ cơ sở dữ liệu y sinh quốc tế, kết hợp tham vấn chuyên môn từ bác sĩ tim mạch tại Bệnh viện huyện Bình Chánh. Kết quả nghiên cứu mang ý nghĩa thực tiễn to lớn khi giúp nâng cao tỷ lệ tín hiệu trên nhiễu, giảm thiểu độ lệch bình phương trung bình xuống dưới 24 x 10^-4 mV^2 và bảo toàn nguyên vẹn cấu trúc hình học của sóng điện tim, tạo tiền đề kỹ thuật vững chắc cho việc phát triển thiết bị y tế số tại Việt Nam.

Cơ sở lý thuyết và phương pháp nghiên cứu

Khung lý thuyết áp dụng

Cơ sở sinh lý học của điện tâm đồ dựa trên hệ thống dẫn truyền tự kích của tim. Nút xoang nhĩ đóng vai trò chủ nhịp phát xung với tần số khoảng 120 lần/phút, lan truyền làm khử cực tâm nhĩ tạo nên sóng P với biên độ khoảng 0,25 mV. Xung động tiếp tục qua nút nhĩ thất với tốc độ làm chậm dẫn truyền ở tần số 50 đến 60 lần/phút, đi qua bó His với tần số phát xung 30 đến 40 lần/phút và phân nhánh vào mạng Purkinje để khử cực tâm thất, tạo nên phức bộ sóng QRS nhọn và cao, kết thúc bằng quá trình tái cực tạo nên đoạn ST và sóng T.

Về mặt toán học xử lý tín hiệu, biến đổi Fourier liên tục và biến đổi Fourier thời gian ngắn bị giới hạn bởi nguyên lý bất định Heisenberg, không thể tối ưu hóa đồng thời độ phân giải thời gian và tần số cho các tín hiệu không ổn định. Biến đổi Wavelet liên tục và biến đổi Wavelet rời rạc khắc phục hạn chế này bằng cách co giãn hàm Wavelet mẹ, nhưng biến đổi Wavelet rời rạc chỉ phân tách dải tần số thấp (thành phần xấp xỉ) mà bỏ qua dải tần số cao (thành phần chi tiết), dẫn đến việc loại bỏ nhầm các thành phần hữu ích mang tính cục bộ.

Phép biến đổi Wavelet Packet khắc phục hoàn toàn nhược điểm trên bằng cách phân tách đệ quy đồng thời cả thành phần xấp xỉ lẫn thành phần chi tiết, tạo nên cấu trúc cây nhị phân hoàn chỉnh. Để lựa chọn cấu trúc cây phân tách tốt nhất, phương pháp áp dụng 5 tiêu chuẩn entropy toán học bao gồm: entropy Shannon, Logarit energy, Sure, Threshold và Norm. Tiêu chuẩn entropy đo lường mức độ mất mát thông tin của hệ thống, chỉ cho phép phân tách nút nhánh khi tổng entropy của hai nút con nhỏ hơn nút cha. Kỹ thuật lọc nhiễu sử dụng phương pháp đặt ngưỡng mềm theo tiêu chuẩn Birgé - Massart với hệ số alpha bằng 2, kết hợp ước lượng độ lệch chuẩn của nhiễu Gauss trắng dựa trên trung vị của các hệ số chi tiết chia cho hằng số 0,6745.

Phương pháp nghiên cứu

Nghiên cứu sử dụng nguồn dữ liệu thứ cấp chuẩn hóa từ cơ sở dữ liệu tín hiệu y sinh quốc tế PhysioNet MIT-BIH. Mẫu thử nghiệm chính là bản ghi tín hiệu điện tim có tên rec1m với quy mô cỡ mẫu gồm 5000 điểm mẫu được thu thập liên tục trong 10 giây ở tần số lấy mẫu 500 mẫu/giây. Phương pháp chọn mẫu tập trung vào các đoạn tín hiệu có đầy đủ các phức bộ hình thái sóng chuẩn, sau đó tiến hành cộng thêm nhiễu Gauss trắng nhân tạo với tỷ lệ tín hiệu trên nhiễu đầu vào biến thiên từ 0 dB đến 18 dB với bước nhảy 0,5 dB để mô phỏng chính xác các cấp độ nhiễu lâm sàng trong thực tế.

Phương pháp phân tích thực nghiệm được thiết kế thông qua quy trình khảo sát toàn diện 54 hàm Wavelet mẹ thuộc các họ toán học Daubechies, Symlets, Coiflets và Biorthogonal trong thư viện Wavelet. Nhằm triệt tiêu hoàn toàn tính ngẫu nhiên của mô hình nhiễu Gauss trắng, thuật toán được lập trình chạy lặp lại 1000 chu kỳ cho mỗi cấu hình thử nghiệm và lấy giá trị trung bình toán học của các hệ số phẩm chất. Ba chỉ số định lượng then chốt được áp dụng để đánh giá chất lượng tín hiệu sau tái tạo gồm: Tỷ lệ tín hiệu trên nhiễu SNR (đo lường bằng đơn vị dB), Độ lệch bình phương trung bình MSE, và Phần trăm căn bậc hai trung bình bình phương sai số PRD. Lý do lựa chọn phương pháp phân tích này là nhằm cung cấp bằng chứng thực nghiệm khách quan, định lượng hóa chính xác mối tương quan giữa dạng sóng Wavelet mẹ và cấu trúc sinh lý của sóng điện tim người.

Kết quả nghiên cứu và thảo luận

Những phát hiện chính

Quá trình khảo sát thực nghiệm trên 54 hàm Wavelet mẹ với tiêu chuẩn entropy Shannon ở mức phân tách bậc 3 đã mang lại các phát hiện định lượng quan trọng. Thứ nhất, nghiên cứu phát hiện hàm Wavelet rbio5.5 (Reverse Biorthogonal 5.5) đạt hiệu quả khử nhiễu tối ưu nhất toàn diện. Cụ thể, hàm rbio5.5 đạt tỷ lệ SNR đầu ra cao nhất, đồng thời kéo giảm chỉ số sai số MSE xuống mức thấp nhất khoảng 23,2 x 10^-4 mV^2 và chỉ số sai số tương đối PRD đạt mức tối thiểu 0,256%. Các hàm thuộc nhóm hiệu năng cao tiếp theo lần lượt là Daubechies db10, Symlet sym8 và Coiflet coif4. Ngược lại, hàm rbio3.1 cho kết quả kém nhất với SNR thấp và độ méo dạng tín hiệu cao.

Thứ hai, khi tiến hành biến thiên số mức phân tách từ mức 1 đến mức 10, nghiên cứu xác định mức phân tách tối ưu cho hàm rbio5.5 là mức N = 4. Trong khi các hàm db10, sym8 và coif4 đạt điểm cực trị chất lượng tại mức phân tách N = 3, hàm rbio5.5 tại mức phân tách N = 4 cho phép triệt tiêu tối đa nhiễu nền mà không làm biến dạng các pha khử cực nhanh.

Thứ ba, khảo sát so sánh giữa 5 quy tắc entropy cho thấy entropy Threshold và Shannon vượt trội nhất về khả năng khử nhiễu ở các dải nhiễu nặng có SNR đầu vào từ 0 dB đến 6 dB. Tuy nhiên, về phương diện duy trì độ phân giải và bảo toàn thông tin lâm sàng khi tăng số mức phân tách từ N = 4 lên N = 10, quy tắc entropy Norm chứng minh tính ổn định xuất sắc nhất khi giá trị SNR hầu như không bị suy giảm theo chiều sâu phân tách cây nhị phân.

Thứ tư, kết quả so sánh đối đầu trực tiếp giữa phép biến đổi Wavelet rời rạc (DWT) và biến đổi Wavelet Packet (WPT) ở cùng mức nhiễu nặng SNR = 5 dB khẳng định WPT vượt trội hơn rõ rệt. WPT giúp tăng tỷ lệ SNR đầu ra cao hơn DWT từ 1,5 dB đến 2,8 dB trên toàn bộ các mẫu thử nghiệm, loại bỏ hoàn toàn các gợn sóng nhiễu tần số cao nhưng vẫn bảo toàn trọn vẹn 100% hình thái sóng P và biên độ đỉnh nhọn sóng R.

Thảo luận kết quả

Hiệu năng vượt trội của phép biến đổi Wavelet Packet bắt nguồn từ cơ chế phân tích đa độ phân giải toàn diện trên toàn bộ dải tần số. Trong khi biến đổi Wavelet rời rạc truyền thống bỏ qua việc lọc ở các nút chi tiết dẫn đến hiện tượng thất thoát năng lượng của các thành phần có ích mang tần số cao, Wavelet Packet phân nhánh chi tiết hóa toàn bộ không gian tín hiệu thành cây nhị phân cân bằng. Khi kết hợp với thuật toán chọn cây tối ưu bằng entropy, các hệ số chứa thông tin giải phẫu quan trọng được giữ lại tối đa trước khi áp dụng ngưỡng mềm Birgé - Massart.

Dữ liệu nghiên cứu có thể được trực quan hóa thông qua hệ thống đồ thị và bảng biểu chi tiết. Đồ thị phân bố 2 chiều thể hiện rõ mối tương quan giữa 54 hàm Wavelet với trục tung biểu diễn SNR và trục hoành biểu diễn MSE, làm nổi bật vị trí vượt trội của nhánh hàm Reverse Biorthogonal. Đồ thị đường đa trục biểu diễn sự phụ thuộc của PRD vào số mức phân tách từ N = 1 đến N = 10 minh họa rõ nét điểm uốn tối ưu tại N = 4. Bên cạnh đó, bảng số liệu đối sánh giữa DWT và WPT cung cấp các tham số định lượng minh chứng rằng khi mức phân tách tăng lên các bậc cao từ N = 6 đến N = 10, tín hiệu sau lọc của DWT bị biến dạng nghiêm trọng làm mất đỉnh R và phẳng hóa sóng P, trong khi tín hiệu xử lý bằng WPT vẫn duy trì sự tương đồng hình học gần như tuyệt đối so với tín hiệu gốc sạch nhiễu ban đầu.

So với các nghiên cứu sử dụng thuật toán lọc thích nghi LMS hoặc phép biến đổi Cosin rời rạc DCT trước đây, phương pháp cải tiến bằng Wavelet Packet trong luận văn này giải quyết triệt để bài toán triệt nhiễu trong điều kiện tỷ lệ tạp âm cực lớn mà không làm xuất hiện hiện tượng trễ pha hay méo dạng cục bộ. Nghiên cứu đã hiện thực hóa thuật toán thành công cụ phần mềm với giao diện đồ họa GUIDE thân thiện trên môi trường MATLAB, tích hợp đầy đủ chức năng nạp tín hiệu, tùy biến tham số lọc và quản lý ngân hàng dữ liệu điện tâm đồ phục vụ công tác nghiên cứu chuyên sâu.

Đề xuất và khuyến nghị

Nhằm chuyển hóa các kết quả nghiên cứu lý thuyết và thực nghiệm vào thực tiễn y tế và kỹ thuật y sinh, bốn khuyến nghị hành động cụ thể được đề xuất như sau:

Thứ nhất, tích hợp thuật toán Wavelet Packet cải tiến sử dụng hàm mẹ rbio5.5 và entropy Threshold vào vi xử lý của các thiết bị đo điện tim kỹ thuật số. Mục tiêu là nâng cao tỷ lệ SNR đầu ra của tín hiệu đo đạt trên mức 15 dB ngay cả trong môi trường nhiều tạp âm điện từ, với lộ trình triển khai trong vòng 6 tháng do các kỹ sư và chuyên gia nghiên cứu phát triển thiết bị y tế chủ trì thực hiện.

Thứ hai, chuẩn hóa quy trình lọc nhiễu tín hiệu điện tâm đồ tại các bệnh viện và phòng khám tim mạch bằng kỹ thuật đặt ngưỡng mềm Birgé - Massart đa mức. Giải pháp này hướng tới việc giảm tỷ lệ chẩn đoán sai lệch hình thái sóng xuống dưới 1%, hoàn thành trong thời gian 9 tháng dưới sự phối hợp của phòng trang thiết bị y tế và các bác sĩ chuyên khoa tim mạch.

Thứ ba, nâng cấp và phát triển giao diện phần mềm GUIDE hiện tại sang nền tảng ứng dụng di động và nền tảng điện toán đám mây đa phương thức. Mục tiêu kỹ thuật là đạt tốc độ xử lý lọc nhiễu thời gian thực với độ trễ phản hồi dưới 50 mili-giây trên mỗi bản ghi 5000 mẫu, thực hiện trong khung thời gian 12 tháng bởi nhóm kỹ sư phần mềm y tế và xử lý tín hiệu số.

Thứ tư, mở rộng xây dựng ngân hàng dữ liệu tín hiệu điện tim số hóa đồng bộ quy mô lớn, tích hợp các công cụ phân tích tự động hỗ trợ giảng dạy đại học và nghiên cứu khoa học chuyên sâu về vật lý kỹ thuật y sinh. Mục tiêu là thu thập và chuẩn hóa hơn 1000 mẫu bệnh lý tim mạch đa dạng trong lộ trình 18 tháng do các cơ sở đào tạo đại học và viện nghiên cứu chuyên ngành đảm trách.

Đối tượng nên tham khảo luận văn

Nội dung và kết quả nghiên cứu của luận văn mang lại giá trị thiết thực cho bốn nhóm đối tượng chính:

Nhóm thứ nhất là các bác sĩ chuyên khoa tim mạch, bác sĩ hồi sức cấp cứu và kỹ thuật viên điện sinh lý tim. Luận văn cung cấp hiểu biết sâu sắc về bản chất toán học của các dạng nhiễu và cơ chế bảo toàn hình thái sóng điện tim, giúp y bác sĩ tự tin phân tích chính xác các chỉ dấu bệnh lý vi mô như đoạn ST chênh lệch, sóng T đảo ngược hay biên độ sóng P bị mờ trong các ca bệnh phức tạp.

Nhóm thứ hai là các kỹ sư nghiên cứu và phát triển thiết bị y tế tại các doanh nghiệp công nghệ sinh học. Luận văn cung cấp toàn bộ quy trình thuật toán từ lựa chọn hàm Wavelet mẹ rbio5.5, cấu trúc cây phân tách bậc 4 đến phương pháp tính toán ngưỡng Birgé - Massart, đóng vai trò như một tài liệu thiết kế kỹ thuật hoàn chỉnh để lập trình nhúng trực tiếp vào các máy đo điện tim cầm tay hoặc thiết bị theo dõi bệnh nhân liên tục.

Nhóm thứ ba là giảng viên, nghiên cứu sinh và học viên cao học thuộc các chuyên ngành Vật lý Y sinh, Kỹ thuật Điện tử - Viễn thông, Tự động hóa và Xử lý Tín hiệu Số. Luận văn là tài liệu tham khảo học thuật chuẩn mực về lý thuyết Wavelet Packet nâng cao, phương pháp thiết kế thực nghiệm kiểm chứng 1000 chu kỳ và hệ thống tiêu chí đánh giá định lượng SNR, MSE, PRD.

Nhóm thứ tư là các nhà khoa học dữ liệu y tế và lập trình viên phát triển hệ thống trí tuệ nhân tạo chẩn đoán hình ảnh y sinh. Cấu trúc ngân hàng dữ liệu và thuật toán tiền xử lý làm sạch tín hiệu trong luận văn đóng vai trò là tầng xử lý dữ liệu đầu vào chuẩn xác, giúp tối ưu hóa độ chính xác cho các mô hình học sâu nhận diện tự động rối loạn nhịp tim.

Câu hỏi thường gặp

Nguyên nhân nào khiến phép biến đổi Wavelet Packet khử nhiễu điện tâm đồ tốt hơn biến đổi Wavelet rời rạc? Biến đổi Wavelet Packet phân tách đệ quy đồng thời cả thành phần xấp xỉ tần số thấp và thành phần chi tiết tần số cao, tạo nên cấu trúc cây nhị phân hoàn chỉnh. Điều này giúp thuật toán giữ lại được các đặc tính địa phương sắc nét của tín hiệu như đỉnh sóng R và sóng P, trong khi biến đổi Wavelet rời rạc bỏ qua việc lọc chi tiết nên làm mất thông tin hữu ích khi gặp mức nhiễu nặng có SNR khoảng 5 dB.

Hàm Wavelet mẹ nào đạt hiệu suất tối ưu nhất cho tín hiệu điện tâm đồ trong nghiên cứu? Qua khảo sát thực nghiệm 54 hàm Wavelet mẹ, hàm Reverse Biorthogonal rbio5.5 đạt hiệu quả cao nhất. Hàm này đạt giá trị SNR đầu ra lớn nhất, đồng thời kéo giảm độ lệch bình phương trung bình MSE xuống khoảng 23,2 x 10^-4 mV^2 và phần trăm sai số PRD xuống 0,256%, chứng tỏ mức độ tương quan cấu trúc cao nhất với dạng sóng điện tim người.

Nên lựa chọn quy tắc entropy nào khi xử lý tín hiệu điện tim nhiễm nhiễu nặng? Khi tín hiệu bị nhiễm nhiễu Gauss trắng nặng với tỷ lệ SNR đầu vào từ 0 dB đến 6 dB, quy tắc entropy Threshold và Shannon là lựa chọn tối ưu nhất để triệt tiêu tạp âm. Tuy nhiên, nếu cần bảo toàn độ phân giải và đặc tính hình học của tín hiệu qua nhiều mức phân tách sâu từ bậc 4 đến bậc 10, quy tắc entropy Norm là giải pháp ổn định nhất.

Mức phân tách tối ưu của thuật toán Wavelet Packet đối với tín hiệu điện tim là bao nhiêu? Nghiên cứu thực nghiệm chứng minh mức phân tách tối ưu cho hàm Wavelet rbio5.5 là mức N = 4. Tại mức phân tách này, các chỉ số phẩm chất SNR đạt giá trị cực đại trong khi MSE và PRD đạt giá trị cực tiểu. Việc phân tách sâu hơn mức 4 không làm tăng thêm chất lượng lọc mà còn làm gia tăng khối lượng tính toán của hệ thống.

Các chỉ số nào được sử dụng để đánh giá định lượng chất lượng tín hiệu sau khi khử nhiễu? Chất lượng tín hiệu được định lượng chính xác qua 3 hệ số phẩm chất gồm: Tỷ lệ tín hiệu trên nhiễu SNR tính theo đơn vị dB (càng lớn chất lượng càng cao), Độ lệch bình phương trung bình MSE (đo lường độ tương quan hình học, giá trị càng nhỏ càng tốt), và Phần trăm căn bậc hai trung bình bình phương sai số PRD (đo lường tỷ lệ sai số tổng thể giữa tín hiệu tái tạo và tín hiệu gốc).

Kết luận

Luận văn đã giải quyết trọn vẹn và khoa học bài toán nâng cao chất lượng tín hiệu điện tâm đồ trong môi trường nhiễu phức tạp thông qua các kết quả then chốt:

• Khẳng định tính ưu việt vượt trội của phép biến đổi Wavelet Packet so với các phương pháp cổ điển như biến đổi Fourier và biến đổi Wavelet rời rạc trong việc xử lý tín hiệu y sinh không dừng.
• Xác định chính xác hàm Wavelet mẹ Reverse Biorthogonal rbio5.5 kết hợp mức phân tách N = 4 là cấu hình tối ưu nhất giúp tối đa hóa tỷ lệ SNR và cực tiểu hóa các sai số MSE, PRD.
• Làm rõ cơ chế hoạt động của 5 tiêu chuẩn entropy, chứng minh tính vượt trội của entropy Threshold trong lọc nhiễu nặng và tính ổn định độ phân giải của entropy Norm.
• Xây dựng thành công bộ công cụ phần mềm với giao diện đồ họa người dùng GUIDE trên MATLAB và thiết lập ngân hàng dữ liệu tín hiệu điện tim phục vụ nghiên cứu lâm sàng.
• Đóng góp giải pháp kỹ thuật có độ tin cậy cao, hỗ trợ đắc lực cho các bác sĩ trong việc chẩn đoán sớm và chính xác các bệnh lý tim mạch nguy hiểm.

Lộ trình phát triển tiếp theo trong 12 đến 24 tháng tập trung vào việc nhúng thuật toán lên phần cứng vi xử lý DSP chuyên dụng và mở rộng thử nghiệm trên các dạng bệnh lý loạn nhịp phức tạp. Các nhà nghiên cứu, kỹ sư y sinh và đơn vị y tế quan tâm có thể tiếp cận khai thác toàn bộ mô hình thuật toán và cơ sở dữ liệu để thúc đẩy ứng dụng thực tiễn công nghệ xử lý tín hiệu y sinh hiện đại.

Trích đoạn nội dung tài liệu

Chương 1: Tông quan: Phân tích, đánh giá công trình nghiên cứu của các tác giả trong và ngoài nước liên quan mật thiết đến luận văn, nêu ra những vấn đề còn tồn đọng, chi ra những van đề mà luận văn tập trung nghiên cứu, giải quyết và khái quất nội dung luận văn. Ngoài ra, luận văn còn trình bảy các cơ sở hình thành tín hiệu điện tim và nhiễu trong tín hiệu điện tim. Chương 2: Các phương pháp xử lý tín hiệu: Nội dung chính của chương này chính là phép biến đổi Wavelet Packet — phép biến đổi được phát triển dựa trên nền tang của phép biến đôi Wavelet rời rac. Chương 3: Phương pháp thực hiện: trình bày phương pháp khử nhiễu trên tín hiện điện tim bang phép biến đổi Wavelet Packet.

Trong phần này, mô hình và nguyên tắc khử nhiễu được đưa ra. Ngoài ra, tôi còn trình bày các hệ số đánh giá chất lượng tín hiệu điện tim, từ đó đưa ra thuật toán khử nhiễu của mình. Chương 4: Kết quả: trình bày các kết quả thu được trong quá trình ứng dụng thuật toán khử nhiễu dé khảo sát quy tắc entropy tôi ưu và đưa ra đánh giá hiệu qua khử nhiễu giữa phép biến đôi Wavelet rời rac và phép biến đổi Wavelet Packet. Mặc khác, luận văn còn đưa ra giao điện xử lý nhiễu tín hiệu điện tim băng phép biến đôi Wavelet Packet và đánh giá so với giao diện GUI trong Wavelet toolbox.

Đông thời, giao diện ngân hàng tín hiệu điện tim được xây dựng bằng lập trình giao điện GUI. Giới thiệu tổng quan Trên thé giới, đã có rất nhiều công trình nghiên cứu về tín hiệu điện tim từ rất sớm. Galvani đã nhận thấy sự co cơ đùi của con ếch khi bị phóng điện, từ đó, ông đưa ra giả thuyết về dong điện trong động vật và đặt nền tảng ban đầu cho nghiên cứu về tín hiệu điện tim. Năm 1887 - 1888, Augustus Desiree Waller ghi nhận điện thế tương ứng với nhịp đập của tim trên bề mặt cơ the người.

Năm 1934, Frank Wilson đưa ra các đạo trình lay ý tưởng từ tam giác “Eithoven”. Sau đó, các nhà khoa học tiếp tục nghiên cứu nhằm hoàn thiện và phát triển lý thuyết vẻ tín hiệu điện tìm. Ngoài ra, phân tích và xử lý tín hiệu y sinh trong đó có tín hiệu điện tim là một trong những chủ dé nóng nhận được sự quan tâm của đông đảo các nhà khoa học trong những năm gan đây. Trong [11], Himanshu Gothwal đã sử dụng FFT (Fast Fourier Transform) và mang Neutron nhân tạo dé dò tìm các dấu hiệu bệnh lý của chứng réi loạn nhịp tim.

Trong [12], Mahajan, R đã sử dụng DCT (Discrete Cosin Transform) và FFT dé nén tín hiệu điện tim. Mặc khác, phép biến đổi Wavelet Packet đã được ứng dụng trong phân tích và xử lý các đối tượng tín hiệu khác nhau va bước đầu cho thay những hiệu quả đáng ké so với phép biến đổi Wavelet rời rac. Trong công trình [13], Wang Jungchen và cộng sự đã sử dụng biến đổi Wavelet Packet dé khử nhiễu tín hiệu âm thanh dưới nước và kết quả cho thấy khử nhiều bằng biến đôi Wavelet Packet thì tốt hon so với biến đổi Wavelet rời rac. Trong [6], Đỗ Huy Khôi và cộng sự đã dùng biến đôi Wavelet Packet kết hợp với thông kê bậc cao (HOS) dé khử nhiễu xung tín hiệu RE và kết quả vẫn tốt hơn so với biến đôi Wavelet rời rạc.

Có nhiều phương pháp nâng cao chất lượng tín hiệu điện tim, trong đó, phương pháp được sử dụng rộng rãi là thuật toán thích nghỉ LMS nhưng thuật toán này không thích hợp với những tín hiệu không ồn định như tín hiệu điện tim, nó gây ra sự dư thừa các thông số trung bình trong lọc thông thấp [3]. Mặc khác, biến đổi Fourier cũng được ứng dụng trong phân tích và xử lý tín hiệu, tuy nhiên, cũng giống như thuật toán thích nghi LMS, biến đôi Fourier chỉ thích hợp với những tín hiệu tuần hoàn và khó phát hiện các đột biến cũng như những tín hiệu không ồn định như tín hiệu điện tim. Hạn chế này chi có thê được khắc phục bởi phép biến đôi Wavelet vì biến đổi Wavelet thích hợp với những tín hiệu không ôn định như là tín hiệu điện tim. Tuy nhiên, từ kết quả của công trình [5] cho thấy phép biến đôi Wavelet rời rac có hiệu quả khử nhiễu chưa cao, vẫn còn tồn đọng một lượng lớn nhiễu trong tín hiệu sau tái tạo, đặc biệt là đối với các mức nhiễu nặng.

Trong luận văn nảy, tôi tập trung nghiên cứu nhằm cải tiền xử lý nhiều của tín hiệu điện tim bằng Wavelet Packet vì phép biến đổi này có thé hạn chế ảnh hưởng của nhiễu và quá trình khử nhiễu lên tín hiệu nhờ phân tích cả thành phân xap xi lẫn thành phan chỉ tiết trên tín hiệu. Mặc khác, phép biến đôi Wavelet Packet còn sử dụng tiêu chuân Entropy đề tính toán phân tích tối ưu do đó có thé hạn chế được sự mất mát thông tin và giữ lại được nhiều đặc tính địa phương trên tín hiệu điện tim. Cơ sở lý thuyết tín hiệu điện tim 1. Cấu trúc giải phẫu và chức năng của tim Động mach chú Tỉnh mạch chủ trên Tere Van động mach cha Tam nhi trai Tĩnh mạch phối phải eS „~~ Tinh mach phôi trái Tâm nhỉ phải Van 2 lá Van 3 lá (Van nhĩ - thất trái) KHI eee Tâm thất trải Tỉnh mạch chủ đưới Tâm thất phải Vách Bên thất Hình 1.

Cấu tao tim người [7]. Tim là một tô chức cơ rong gom 4 buông. Bên ngoài được bao bọc bởi một túi sợi gọi là bao tim, bên trong được cấu tao bằng cơ tim có vách ngăn chia tim thành hai nửa riêng biệt gọi là tim trái và um phải. Tim trai bơm máu ra ngoại vi, còn um phải bom máu lên phối.

Mỗi nửa tim lại được chia ra thành hai buồng: buông trên là tâm nhĩ có thành mong làm nhiệm vụ chứa máu; buồng dưới là tâm thất có thành dày, khối cơ lớn giúp cung cấp lực đây máu đi đến các bộ phận. Giữa tâm nhĩ và tâm thất có van nhĩ thất, giữa tâm thất trái và động mạch chú, tâm thất phải và động mạch phôi có van bán nguyệt. Các van này đảm bảo cho máu chỉ di chuyên theo một chiều từ tâm nhĩ xuống tâm thắt, từ tâm thất xuống động mạch chứ không cho đi ngược lại, nhờ vậy đảm bảo được sự tuần hoàn máu [7]. Ngoài ra, tim còn có một cau trúc đặc biệt thực hiện chức năng phát và dẫn truyền xung được gọi là hệ dẫn truyền.

Hệ thống dẫn truyền gồm: Nút xoang nhĩ (SAN): là nút tạo nhịp cho toàn bộ trái tim, nằm ở cơ tâm nhĩ phải, phát xung với tan số khoảng 120 làn/phút. Các đường liên nút: nằm ở giữa nút xoang nhĩ và nút nhĩ thất, thực hiện chức năng dan truyền các xung động giữa nút xoang nhĩ và nút nhĩ that, Nút nhĩ that (AVN): nằm ở bên phải vách liên nhĩ, giữ nhiệm vụ làm chậm dan truyền trước khi các xung động được truyền xuống thất với tần số khoảng 50 — 60 lần/phút. Bó Hís: bắt đầu từ nút nhĩ thất đến vách liên thất thì chía thành hai nhánh trái và phải chạy dưới nội tâm mạc hai that dé dẫn truyén xung động đến hai thất, tại đây, chúng phân nhánh thành mạng Purkinje chạy giữa các sợi cơ tim giúp dẫn truyền xung động xuyên qua các thành của thất. Bó His phát xung khoảng 30 — 40 lần/phút.

Hệ thống dan truyền tim [7]. Quá trình hình thành tín hiệu điện tim a. Nhĩ đồ Tim hoạt động được nhờ vào một xung động truyền qua một hệ thong than kinh tự kích của tim. Đầu tiên, nút xoang nhĩ sẽ phát xung tự động, xung động tỏa ra làm cơ nhĩ khử cực trước.

Sóng khử cực có hướng chung là từ trên xuống dưới, từ phải sang trái và hợp với phương ngang một góc 49° [7]. Dợt sóng này được máy ghi điện tìm ghi lại với dang một sóng dương, don, thấp, nhỏ và có biên độ khoảng 0,25mV gọi là sóng P. ° ~” ~~ oo na lees Leeda Hình 1. Khử cực ở tâm nhĩ va sự hình thành sóng P [7].

Thất đồ Ngay khí nhĩ còn đang khứ cực thì xung động đã bắt đầu truyền vào nút nhĩ thất xuống thất và hai nhánh bó His xuống khử cực thất. Sóng khử cực hướng từ giữa mặt trái đi xuyên qua mặt phải của vách liên thất, Máy sẽ ghi nhận được một sóng âm nho, gọn gọi là sóng Q. Sự khử cực ở vách liên thất và hình thành sóng Q [7]. Xung tiếp tục truyền xuống và tiên hành khử cực đồng thời cả hai tâm thất theo hướng xuyên qua bé mặt đảy cơ tim, từ dưới nội tâm mạc ra dưới thượng tâm mạc.

Vectơ khử cực hướng từ phải sang trái và máy ghi nhận được một làn sóng dương, cao và nhọn gọi là sóng R. Sau cùng, xung động truyền xuống và khử cực vùng đáy thất. Vectơ khử cực hướng từ trái sang phải, máy sẽ ghi nhận được một sóng âm, nhỏ và nhọn gọi là sóng S. Sự khứ cực ở tâm thất và hình thành sóng R, S [7].

Sau khi thất khử cực xong sẽ qua thời kỳ tái cực chậm. Giai đoạn nảy được thê hiện trên điện tâm đỗ bằng một đường đăng điện gọi là đoạn S — T. Sau đó là thời kỳ tái cực nhanh tạo nên sóng T. Tái cực có hướng xuyên qua cơ tim, từ lớp đưới thượng tâm mạc vào lớp dưới nội tâm mạc.

Vectơ tái cực có hướng từ trên xuống dưới và từ phải sang trái tạo ra một sóng dương, thấp, không đối xứng mà có sườn lên thoai thoải hơn và sườn xuống dốc đứng hơn gọi là sóng T. Sau khi kết thúc 10 sóng T còn có thé thấy được một sóng chậm nhỏ gọi là sóng U đặc trưng cho giai đoạn tát cực muộn. Sự tdi cực ở tâm thất và hình thành sóng T [7]. Nhiễu Gauss trắng trong tín hiệu điện tim Tín hiệu điện tim chứa đựng nhiều thông tin lâm sàng quý giá trong khoảng tan số 0 — 100Hz.

Tuy nhiên trong khoảng tần số nay, ECG dé bị ảnh hưởng bởi bởi các loại nhiều khác nhau như nhiễu nguồn điện, nhiễu tín hiệu điện cơ EMG, nhiễu do tiếp xúc điện cực và da không tốt gây khó khăn cho bác sĩ trong quá trình chân đoán các bệnh lý tim mạch. Trong quá trình khảo sát ảnh hưởng của nhiễu đến chất lượng tín hiệu điện tim, trong [3] đã chỉ rõ rằng lọc nhiễu từ mạng cung cấp điện là cần thiết nhất vì tính chất phố biến và khó kiểm soát của nó.