Tổng quan nghiên cứu
Rối loạn động kinh là một bệnh lý thần kinh phổ biến với khoảng 50 triệu người mắc trên toàn cầu, trong đó phần lớn nằm ở các nước nghèo hoặc đang phát triển. Việc chuẩn đoán và theo dõi bệnh động kinh dựa trên dữ liệu điện não đồ (EEG) đóng vai trò quan trọng trong việc xác định chính xác các cơn co giật và dự đoán cơn động kinh. Tuy nhiên, việc phân tích EEG truyền thống phụ thuộc nhiều vào kinh nghiệm của bác sĩ, dẫn đến hiệu quả thấp và tốn nhiều thời gian. Nghiên cứu này nhằm mục tiêu xây dựng một hệ thống tự động phát hiện xung động kinh trong dữ liệu EEG bằng phương pháp phân tích tensor, cụ thể là thuật toán phân tích HOSVD (Higher-Order Singular Value Decomposition).
Phạm vi nghiên cứu tập trung vào dữ liệu EEG được thu thập từ hai bệnh nhân với tần số lấy mẫu 5000Hz, trong đó các đoạn dữ liệu kéo dài 6 giây, tổng cộng khoảng 3 triệu mẫu dữ liệu trên mỗi kênh. Hệ thống được phát triển nhằm tăng độ chính xác nhận diện, giảm thiểu tỷ lệ sai sót trong phát hiện cơn động kinh, giúp nâng cao hiệu quả chẩn đoán và điều trị. Kết quả nghiên cứu có ý nghĩa lớn trong việc áp dụng trí tuệ nhân tạo và các kỹ thuật xử lý tín hiệu hiện đại vào y học, đặc biệt trong lĩnh vực thần kinh học và chẩn đoán tự động.
Cơ sở lý thuyết và phương pháp nghiên cứu
Khung lý thuyết áp dụng
Nghiên cứu này dựa trên hai lý thuyết chính:
-
Phân tích Tensor và HOSVD: Tensor là cấu trúc toán học mở rộng từ ma trận cho dữ liệu đa chiều, được sử dụng để biểu diễn các tín hiệu EEG theo các chiều thời gian, kênh đo và tần số. HOSVD là một dạng mở rộng của phân tích giá trị riêng (SVD) cho tensor, giúp phân tách và trích xuất đặc trưng một cách hiệu quả. HOSVD phù hợp với dữ liệu EEG đa chiều, cho phép nhận diện các yếu tố đặc trưng riêng biệt của xung động kinh.
-
Biến đổi Sóng liên tục (CWT): CWT được sử dụng để phân tích tín hiệu EEG ở cả miền thời gian và tần số, nhờ đó có thể trích xuất thông tin đặc trưng chính xác hơn. Hàm sóng mẹ Mexican Hat được chọn nhờ ưu điểm làm nổi bật tín hiệu bất thường trên sóng EEG.
Các khái niệm chính gồm: tensor bậc 3, lát cắt tensor, chế độ tensor, chuẩn Frobenius, ma trận hóa tensor, và phân tích HOSVD với ba loại ma trận thành phần đại diện cho thời gian, không gian và kênh tín hiệu EEG.
Phương pháp nghiên cứu
Nguồn dữ liệu sử dụng từ bộ dữ liệu Kaggle, chứa tín hiệu EEG ghi lại bằng 15 kênh với tần số lấy mẫu 5000 Hz, tổng thời gian ghi khoảng một giờ cho hai bệnh nhân. Dữ liệu được chia thành các đoạn 6 giây, với khoảng 900 đoạn chứa xung động kinh (Preictal) và 2500 đoạn không có xung động kinh (Interictal). Thời gian lấy mẫu mới sau khi xử lý được giảm xuống còn 200 Hz để phù hợp với tiêu chuẩn Nyquist và tăng tốc quá trình tính toán.
Quá trình thu thập và xử lý dữ liệu gồm các bước:
-
Biến đổi Wavelet liên tục (CWT): Dữ liệu EEG được biến đổi sang dạng tensor bậc 3 qua phân tích sóng mẹ Mexican Hat, biểu diễn chi tiết thông tin thời gian-tần số-kênh.
-
Xây dựng tensor bậc 3 cho xung động kinh và hoạt động bình thường: Các đoạn dữ liệu EEG với phân loại rõ ràng được biểu diễn tensor để sắp xếp thành tensor bậc 4 phục vụ cho phân tích HOSVD.
-
Phân tích HOSVD: Thực hiện phân rã tensor để trích xuất các ma trận thành phần đặc trưng (A_time, B_scale, C_channel). Việc xác định hạng đa tuyến được tiến hành dựa trên biến thiên eigen gap và độ tương quan giữa các thành phần liên tiếp.
-
Chia dữ liệu huấn luyện và kiểm thử: Sử dụng phương pháp cross-validation với tỷ lệ 4:1, huấn luyện bộ phân loại SVM và KNN dựa trên các đặc trưng trích xuất từ tensor.
-
Đánh giá mô hình: Sử dụng các chỉ số chính gồm độ nhạy (Sensitivity), độ chính xác (Accuracy), độ đặc hiệu (Specificity) và diện tích dưới đường cong ROC (AUC) để đánh giá hiệu năng.
Phân tích được thực hiện trên phần mềm Matlab với cấu hình máy tính Intel Core i5, RAM 16GB, GPU GeForce GTX 1060.
Kết quả nghiên cứu và thảo luận
Những phát hiện chính
-
Hiệu quả phân tích HOSVD: Kết quả phân tích dữ liệu EEG cho thấy việc sử dụng hạng đa tuyến R = [200, 5, 15] là tối ưu, với 200 thành phần thời gian, 5 thành phần theo miền tần số và 15 thành phần kênh tín hiệu. Độ tương quan tổng thể giữa tensor nguyên gốc và tensor xấp xỉ đạt gần 100%.
-
Độ nhạy và độ chính xác của hệ thống: Với bộ phân loại SVM, hệ thống đạt độ nhạy lên tới 98%, độ chính xác trung bình 78%, độ đặc hiệu 72%, trong khi KNN có độ nhạy và độ chính xác thấp hơn nhưng chênh lệch không đáng kể, cho thấy SVM phù hợp hơn trong bài toán này.
-
So sánh với phương pháp CP/PARAFAC: Phương pháp HOSVD vượt trội hơn phương pháp CP/PARAFAC với độ nhạy đạt 99%, so với 75% của CP, chứng tỏ HOSVD hiệu quả hơn trong xử lý dữ liệu tensor cao chiều từ EEG.
-
Biểu đồ ROC: Diện tích dưới đường cong ROC của phương pháp HOSVD kết hợp SVM đạt trên 86%, khẳng định khả năng phân biệt tốt giữa các thời điểm có và không có xung động kinh.
Thảo luận kết quả
Sự vượt trội của HOSVD trong phân tích xung động kinh EEG có thể được giải thích bởi khả năng trích xuất đa chiều đặc trưng chính xác từ dữ liệu tensor đa chiều phức tạp, trong khi CP/PARAFAC dễ mất mát thông tin khi dữ liệu có nhiều yếu tố không tuyến tính. Sự lựa chọn hàm sóng mẹ Mexican Hat trong CWT càng giúp làm nổi bật các đặc trưng bất thường trên miền tần số-thời gian.
Kết quả cũng phù hợp với các báo cáo ngành và một số nghiên cứu sử dụng phân tích tensor trong xử lý EEG, đồng thời nhấn mạnh vai trò của SVM với cơ chế superplane giúp phân tách tốt hơn so với KNN trong dữ liệu có nhiều tính chất phức tạp và nhiễu như EEG.
Việc giảm tần số lấy mẫu xuống còn 200Hz làm tăng tốc độ xử lý mà không làm mất đi lượng thông tin có giá trị nhờ vào việc sử dụng stratified sampling và biến đổi Wavelet hiệu quả.
Các biểu đồ ROC và đồ thị tương quan eigen gap cho thấy dữ liệu rất thích hợp cho phân tích tensor với hạng đa tuyến phù hợp. Ngoài ra, kết quả cho thấy cách tiếp cận này có thể triển khai trên các hệ thống chẩn đoán tự động thực tế, gián tiếp cải thiện chất lượng sống bệnh nhân.
Đề xuất và khuyến nghị
-
Ứng dụng mô hình phát hiện xung động kinh tự động: Triển khai mô hình phân tích HOSVD kết hợp SVM vào các thiết bị chẩn đoán EEG tại các cơ sở y tế từ năm 2024 để nâng cao hiệu quả chẩn đoán, giảm thiểu sai sót và tốn kém thời gian của bác sĩ.
-
Mở rộng tập dữ liệu và đa dạng hóa kênh EEG: Thu thập dữ liệu nhiều bệnh nhân với các đặc trưng khác nhau và bổ sung thêm kênh để hoàn thiện mô hình, cải thiện khả năng tổng quát của hệ thống trong vòng 2-3 năm tiếp theo.
-
Phát triển thêm thuật toán phân tích tensor nâng cao: Thử nghiệm các kỹ thuật phân tích tensor không âm (NTD) và phân tích thành phần chính tuyến tính phi tuyến để tăng cường độ chính xác và khả năng phân loại các dạng bất thường đa dạng hơn.
-
Tích hợp mô hình vào hệ thống hỗ trợ quyết định lâm sàng: Phát triển giao diện thân thiện, tích hợp vào hệ thống y tế điện tử, hỗ trợ bác sĩ trong việc đánh giá dữ liệu EEG theo thời gian thực nhằm phát hiện sớm cơn động kinh, triển khai trong vòng 1-2 năm.
Đối tượng nên tham khảo luận văn
-
Chuyên gia và kỹ sư trong lĩnh vực xử lý tín hiệu y sinh: Luận văn cung cấp phương pháp và mô hình phân tích tensor nâng cao cùng các thuật toán học máy, hữu ích trong nghiên cứu và phát triển ứng dụng xử lý điện não đồ.
-
Bác sĩ và nhà nghiên cứu thần kinh học: Các mô hình phát hiện tự động giúp cải thiện chẩn đoán và theo dõi bệnh động kinh, góp phần nâng cao chất lượng điều trị và chăm sóc bệnh nhân.
-
Sinh viên và học viên cao học ngành công nghệ thông tin, kỹ thuật y sinh, trí tuệ nhân tạo: Luận văn trình bày chi tiết khung lý thuyết, phương pháp và thực nghiệm, là tài liệu tham khảo quan trọng trong đào tạo và nghiên cứu.
-
Doanh nghiệp phát triển thiết bị y tế và phần mềm chẩn đoán: Các công nghệ mới như phân tích tensor và học máy có thể ứng dụng vào sản phẩm tự động phát hiện cơn động kinh, nâng cao giá trị sản phẩm trên thị trường.
Câu hỏi thường gặp
-
Phân tích tensor là gì và tại sao lại chọn HOSVD cho dữ liệu EEG?
Phân tích tensor mở rộng từ ma trận cho dữ liệu đa chiều. HOSVD là thuật toán phân tích tensor tổng quát, cho khả năng trích xuất đặc trưng tốt hơn trong dữ liệu EEG đa chiều với các chiều thời gian, kênh và tần số, giúp phát hiện xung động kinh hiệu quả hơn. -
Biến đổi sóng liên tục (CWT) có vai trò gì trong xử lý tín hiệu EEG?
CWT cho phép phân tích tín hiệu theo thời gian-tần số đồng thời, giúp trích xuất các đặc trưng tín hiệu phức tạp của EEG, như xung động và các biến đổi trong tần số, làm tăng độ nhạy của hệ thống phát hiện. -
SVM và KNN khác nhau như thế nào trong việc phân loại dữ liệu EEG?
SVM tạo siêu phẳng phân tách tối ưu giữa hai lớp dữ liệu, phù hợp với dữ liệu có nhiều đặc trưng phức tạp và có nhiễu, còn KNN dựa trên khoảng cách đến các điểm lân cận nên dễ bị ảnh hưởng bởi nhiễu và độ phân bố dữ liệu kém đồng nhất. -
Thế nào là hạng đa tuyến trong phân tích HOSVD và cách lựa chọn nó?
Hạng đa tuyến là số thành phần trích xuất theo từng chiều tensor. Lựa chọn dựa trên phân tích các khoảng cách giữa giá trị riêng liên tiếp (eigen gap) và độ tương quan các thành phần tạo ra mô hình xấp xỉ tốt nhất với tensor gốc. -
Hệ thống này có thể áp dụng thực tế ra sao trong y tế?
Hệ thống có thể tích hợp vào phần mềm chẩn đoán EEG tự động, giúp bác sĩ phát hiện chính xác thời điểm cơn động kinh xuất hiện, giảm phụ thuộc vào phân tích thủ công, tiết kiệm thời gian và nâng cao chất lượng điều trị.
Kết luận
- Nghiên cứu thành công ứng dụng phân tích HOSVD trên dữ liệu tensor EEG để phát hiện xung động kinh tự động với độ nhạy đến 99%.
- Hệ thống kết hợp HOSVD và SVM cho hiệu quả phân loại cao hơn so với các phương pháp truyền thống và CP/PARAFAC.
- Kết quả thực nghiệm với bộ dữ liệu đa kênh và tần số lấy mẫu cao chứng minh tính khả thi và áp dụng của phương pháp.
- Các chỉ số đánh giá như độ nhạy, độ chính xác và AUC đều đạt mức cao, phù hợp cho ứng dụng chẩn đoán y khoa.
- Đề xuất phát triển mô hình đa dạng hơn và tích hợp vào hệ thống hỗ trợ chẩn đoán thực tế nhằm cải thiện chất lượng cuộc sống bệnh nhân.
Luận văn mở ra hướng nghiên cứu hiệu quả kết hợp phân tích tensor và học máy trong chăm sóc sức khỏe thần kinh, khuyến khích nghiên cứu sâu rộng hơn vào các mô hình nâng cao và ứng dụng công nghệ AI trong y sinh. Hành động tiếp theo là triển khai thử nghiệm tiền lâm sàng và mở rộng dữ liệu đa dạng hơn để hoàn thiện hệ thống.