Tổng quan nghiên cứu

Trong bối cảnh khoa học thần kinh và kỹ thuật y sinh hiện đại, công nghệ giao diện não - máy tính (Brain-Computer Interface - BCI) đang đạt tốc độ tăng trưởng ấn tượng với khoảng 15% đến 20% mỗi năm trên toàn cầu. Việc phát triển các phương pháp giải mã ý định vận động không xâm nhập đóng vai trò then chốt trong việc hỗ trợ phục hồi chức năng cho hàng triệu bệnh nhân bị tổn thương thần kinh hoặc đột quỵ. Tuy nhiên, các tín hiệu não thu thập từ thực nghiệm thường xuyên bị ảnh hưởng bởi nhiễu chuyển động và nhiễu sinh lý, khiến việc phân loại chính xác các trạng thái vận động tinh vi gặp nhiều thách thức.

Đề tài tập trung giải quyết vấn đề nhận dạng chính xác hoạt động gõ tay trái hoặc tay phải của con người thông qua phân tích biến thiên huyết động học vỏ não. Mục tiêu cụ thể là xây dựng một chuỗi giải thuật xử lý tín hiệu toàn diện, bao gồm khử nhiễu tín hiệu bằng bộ lọc thích nghi, trích xuất đặc trưng hình thái nồng độ oxy-hemoglobin và tự động phân loại bằng các mô hình học máy hiện đại. Nghiên cứu được triển khai thực nghiệm trên hệ thống đo quang phổ cận hồng ngoại chức năng fNIRS tại Bộ môn Kỹ thuật Y sinh, Trường Đại học Quốc tế thuộc Đại học Quốc gia Thành phố Hồ Chí Minh. Phạm vi nghiên cứu khảo sát trên các tình nguyện viên nam trong độ tuổi từ 22 đến 25 tuổi, tập trung thu thập dữ liệu tại vùng vỏ não điều khiển vận động với độ sâu thấu quang khoảng 3 cm xuyên qua hộp sọ. Kết quả nghiên cứu góp phần nâng cao độ chính xác nhận dạng cử động vận động đạt trên 90%, giảm thiểu số chiều dữ liệu ngõ vào cho các thuật toán phân loại và tạo nền tảng vững chắc cho việc ứng dụng công nghệ fNIRS trong chẩn đoán y khoa cũng như thiết bị trợ thị vận động 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

Nghiên cứu được xây dựng trên nền tảng quang phổ cận hồng ngoại chức năng (fNIRS) kết hợp với định luật Beer-Lambert mở rộng. Kỹ thuật fNIRS sử dụng chùm ánh sáng trong dải bước sóng cận hồng ngoại từ 650 nm đến 950 nm để xuyên qua mô xương sọ đến bề mặt vỏ não. Tại đây, mật độ quang hấp thụ thay đổi tương ứng với nồng độ của hai chất hấp thụ chính là oxy-hemoglobin (Oxy-Hb) và deoxy-hemoglobin (Deoxy-Hb). Thiết bị sử dụng 3 bước sóng tiêu chuẩn gồm 780 nm, 805 nm và 830 nm để tính toán định lượng sự biến thiên huyết động dựa trên phương trình hấp thụ vi phân.

Về mặt sinh lý học thần kinh, nghiên cứu dựa trên lý thuyết ghép nối thần kinh - mạch máu tại vùng vỏ não điều khiển chuyển động thuộc thùy đỉnh. Khi một khu vực não bộ được kích hoạt bởi hành vi vận động gõ tay, lưu lượng máu cục bộ sẽ tăng mạnh, dẫn đến sự gia tăng nồng độ Oxy-Hb và giảm nồng độ Deoxy-Hb tương ứng. Các khái niệm toán học và xử lý tín hiệu cốt lõi bao gồm: bộ lọc làm phẳng đa thức Savitzky-Golay bảo toàn biên dạng tín hiệu, mô hình hồi quy đa thức bình phương tối thiểu để biểu diễn động học nồng độ oxy, không gian siêu phẳng phân cách tối ưu của máy học vector hỗ trợ (SVM) và mạng nơ-ron truyền thẳng đa lớp lan truyền ngược (Backpropagation ANN).

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

Nghiên cứu sử dụng phương pháp thực nghiệm chọn mẫu có chủ đích trên 3 đối tượng nam giới khỏe mạnh, có độ tuổi lần lượt là 22, 24 và 25 tuổi (độ tuổi trung bình khoảng 23,7 tuổi) và đều có đặc điểm thuận tay phải nhằm kiểm soát tính đồng nhất của phản xạ bán cầu não. Nguồn dữ liệu thu thập từ máy fNIRS model FOIRE-3000 của hãng Shimadzu Nhật Bản với chu kỳ lấy mẫu 0,1 giây. Hệ thống bố trí ma trận đầu đo 4x2 tạo thành 20 kênh thu tín hiệu trên hai bán cầu não, trong đó nghiên cứu trích lọc 8 kênh trọng điểm tại khu vực điều khiển vận động, gồm các kênh 2, 5, 6, 9 ở bán cầu não trái và các kênh 12, 15, 16, 19 ở bán cầu não phải với khoảng cách giữa đầu phát và đầu thu cố định là 3 cm.

Quy trình phân tích dữ liệu trải qua ba giai đoạn chặt chẽ:

  • Tiền xử lý dữ liệu: Sử dụng bộ lọc Savitzky-Golay với chiều dài cửa sổ 11 mẫu và đa thức bậc 3. Phương pháp này được lựa chọn vì khả năng làm phẳng tín hiệu vượt trội, triệt tiêu nhiễu chuyển động đột biến có tần số cao mà không làm suy hao đỉnh biên độ sinh lý của Oxy-Hb, đạt tần số cắt tối ưu tại 0,1977 Hz.
  • Trích xuất đặc trưng: Áp dụng hồi quy đa thức bậc 3 theo phương pháp bình phương tối thiểu để biểu diễn chuỗi thời gian Oxy-Hb thành một vector các hệ số hồi quy cô đọng, giúp giảm số lượng biến ngõ vào nhưng vẫn giữ trọn vẹn đặc tính động học.
  • Huấn luyện và nhận dạng: Triển khai thuật toán PR-SVM kết hợp điều kiện tối ưu Karush-Kuhn-Tucker cùng thuật toán tối ưu tuần tự SMO, song song với mô hình PR-ANN có cấu trúc 3 lớp gồm 48 nút ngõ vào, 100 nút ẩn và 2 nút ngõ ra để so sánh hiệu năng.

Giao thức thời gian thực nghiệm được thiết kế theo chu kỳ chuẩn hóa 60 giây cho mỗi lượt đo, bao gồm 20 giây nghỉ ban đầu, 20 giây thực hiện tác vụ gõ tay (khoảng 10 lần gõ nhịp nhàng) và 20 giây nghỉ phục hồi, lặp lại nhiều phiên đo độc lập trên từng chủ thể.

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

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

Quá trình phân tích thực nghiệm trên tín hiệu fNIRS đã mang lại nhiều phát hiện khoa học quan trọng về mặt định lượng:

Thứ nhất, bộ lọc Savitzky-Golay cấu hình cửa sổ 11 và bậc đa thức 3 đã loại bỏ hoàn toàn các gai nhiễu đột biến do dịch chuyển đầu đo gây ra, giúp làm mượt đường cong tín hiệu Oxy-Hb trong khung thời gian vận động 7 giây. So với tín hiệu thô ban đầu, độ biến thiên sai số ngẫu nhiên giảm hơn 40%, phản ánh chính xác xu hướng tăng dần nồng độ oxy khi bắt đầu gõ tay và giảm dần khi bước vào pha nghỉ.

Thứ hai, các hệ số hồi quy đa thức thể hiện rõ tính bất đối xứng chức năng giữa hai bán cầu não đối với người thuận tay phải. Khi đối tượng thực hiện gõ tay phải, nồng độ Oxy-Hb tại các kênh 2, 5, 6, 9 trên bán cầu não trái tăng cao hơn từ 30% đến 35% so với bán cầu não phải. Ngược lại, khi gõ tay trái, bán cầu não phải có sự gia tăng huyết động nhưng bán cầu não trái vẫn duy trì mức kích hoạt nền đáng kể, chứng minh vai trò bán cầu não ưu thế ở người thuận tay phải.

Thứ ba, cả hai mô hình nhận dạng đề xuất đều đạt độ chính xác phân loại vượt trội. Thuật toán PR-SVM đạt độ chính xác trung bình ấn tượng từ 88% đến 93% trên cả 3 chủ thể thử nghiệm. Mô hình PR-ANN cũng cho kết quả tương đồng với độ chính xác dao động từ 86% đến 91%, khẳng định tính hiệu quả và độ ổn định cao của bộ đặc trưng hệ số hồi quy đa thức.

Thảo luận kết quả

Hiện tượng tăng mạnh nồng độ Oxy-Hb ở bán cầu não đối bên hoàn toàn phù hợp với các quy luật sinh lý học thần kinh kinh điển về đường dẫn truyền vận động chéo qua tủy sống. Kết quả so sánh giữa các kênh cho thấy kênh số 6 ở não trái và kênh số 16 ở não phải là hai vị trí ghi nhận mức độ kích hoạt tín hiệu mạnh nhất, tương ứng với tâm điểm phân vùng điều khiển cử động bàn tay và ngón tay trên thùy đỉnh.

Khi so sánh với các nghiên cứu trước đây trong cùng lĩnh vực, việc đưa trực tiếp toàn bộ chuỗi tín hiệu thô Oxy-Hb và Deoxy-Hb vào bộ phân loại SVM thường dẫn đến hiện tượng quá khớp do số chiều dữ liệu quá lớn, hoặc phương pháp hồi quy tuyến tính đơn biến chỉ đạt độ chính xác hạn chế do không mô tả được độ cong phi tuyến của đáp ứng huyết động. Giải thuật PR-SVM trong nghiên cứu này đã tối ưu hóa không gian đặc trưng bằng cách thu gọn dữ liệu chuỗi thời gian thành các hệ số đa thức đại diện, giúp tăng tốc độ huấn luyện mô hình lên khoảng 3 lần và cải thiện độ chính xác phân loại thêm 12% so với mô hình hồi quy tuyến tính cổ điển.

Dữ liệu kết quả nghiên cứu được biểu diễn trực quan thông qua các biểu đồ phổ đáp ứng tần số biên độ với tần số cắt 0,1977 Hz, đồ thị phân tích tín hiệu huyết động theo thang thời gian thực nghiệm và các bảng ma trận kiểm tra chéo độ chính xác giữa từng chủ thể, mang lại độ tin cậy khoa học cao cho toàn bộ kết luận của đề tài.

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

Dựa trên những kết quả thực nghiệm và tiềm năng phát triển của đề tài, các giải pháp và khuyến nghị kỹ thuật được đề xuất như sau:

  • Phát triển thuật toán xử lý luồng dữ liệu thời gian thực: Nhóm nghiên cứu kỹ thuật y sinh cần xây dựng mô-đun cửa sổ trượt trực tuyến có khả năng tự động lọc nhiễu Savitzky-Golay và tính toán hệ số hồi quy tức thời, mục tiêu giảm độ trễ phản hồi hệ thống xuống dưới 100 mili-giây trong vòng 6 đến 12 tháng tới.
  • Mở rộng quy mô mẫu thử nghiệm đa dạng: Các đơn vị nghiên cứu cần mở rộng cỡ mẫu lên tối thiểu 50 tình nguyện viên trong vòng 12 đến 18 tháng, bổ sung thêm nhóm người thuận tay trái, các nhóm tuổi khác nhau và bệnh nhân đang điều trị di chứng tai biến mạch máu não nhằm nâng cao tính đại diện và khả năng khái quát hóa của mô hình học máy.
  • Tối ưu hóa kiến trúc học sâu kết hợp: Các kỹ sư công nghệ thông tin và trí tuệ nhân tạo nên nghiên cứu tích hợp mạng nơ-ron tích chập kết hợp mạng hồi quy dài hạn (CNN-LSTM) để tự động hóa hoàn toàn khâu khai phá đặc trưng chuỗi thời gian, hướng tới mục tiêu nâng tỷ lệ phân loại chính xác các tác vụ vận động phức tạp lên trên 95% trong lộ trình 12 tháng.
  • Nghiên cứu chế tạo thiết bị fNIRS không dây dạng đeo: Các viện nghiên cứu cần chủ động liên kết với các doanh nghiệp sản xuất trang thiết bị y tế trong thời gian 24 tháng để thiết kế các mũ đo fNIRS thu nhỏ 8 kênh với kết nối không dây, tạo điều kiện thuận lợi cho việc ứng dụng giao diện não - máy tính tại các trung tâm phục hồi chức năng và bệnh viện chuyên khoa thần kinh.

Đố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ị học thuật và ứng dụng thực tiễn cao cho 4 nhóm đối tượng trọng tâm:

  • Các nhà nghiên cứu và giảng viên chuyên ngành Kỹ thuật Y sinh, Khoa học Thần kinh: Tài liệu cung cấp cơ sở lý thuyết chuẩn mực về quang phổ cận hồng ngoại fNIRS, mô hình biến thiên huyết động học não bộ và phương pháp thiết kế giao thức thực nghiệm BCI đạt tiêu chuẩn quốc tế.
  • Kỹ sư phát triển giải pháp Trí tuệ nhân tạo và Xử lý tín hiệu số: Cung cấp tài liệu tham khảo chi tiết về quy trình lập trình bộ lọc thích nghi Savitzky-Golay, kỹ thuật trích xuất đặc trưng hồi quy đa thức OLS và thuật toán tối ưu hóa siêu phẳng SVM trên dữ liệu y sinh nhiều chiều.
  • Bác sĩ chuyên khoa Phục hồi chức năng và Thần kinh học: Nắm bắt cơ chế kích hoạt huyết động vỏ não trong các tác vụ vận động, từ đó ứng dụng để đánh giá định lượng mức độ phục hồi chức năng của người bệnh sau đột quỵ hoặc chấn thương sọ não.
  • Học viên cao học và sinh viên các ngành Kỹ thuật Điện tử, Tự động hóa và Y sinh: Là nguồn tài liệu tham khảo mẫu mực về phương pháp nghiên cứu khoa học, quy trình vận hành thiết bị đo đạc chuyên dụng FOIRE-3000 và kỹ năng trình bày báo cáo học thuật chuyên sâu.

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

Kỹ thuật fNIRS có ưu thế gì vượt trội so với EEG và MRI trong nghiên cứu giao diện não - máy tính?

fNIRS sử dụng ánh sáng cận hồng ngoại với độ sâu thấu quang 3 cm, có khả năng định vị không gian vượt trội so với điện não đồ EEG và ít bị ảnh hưởng bởi nhiễu điện cơ. So với máy cộng hưởng từ MRI cồng kềnh và đắt đỏ, fNIRS sở hữu thiết bị nhỏ gọn, chi phí vận hành thấp, môi trường đo thân thiện và cho phép người tham gia cử động tự nhiên trong suốt quá trình thử nghiệm.

Tại sao bộ lọc Savitzky-Golay lại hiệu quả hơn các bộ lọc số truyền thống trong việc khử nhiễu fNIRS?

Bộ lọc Savitzky-Golay với cửa sổ 11 điểm và đa thức bậc 3 sử dụng nguyên lý làm phẳng bình phương tối thiểu cục bộ. Nhờ đó, bộ lọc triệt tiêu hoàn toàn các gai nhiễu chuyển động tần số cao nhưng vẫn giữ nguyên vẹn độ cao đỉnh, độ rộng và diện tích tích phân của đáp ứng huyết động Oxy-Hb, tránh hiện tượng làm méo dạng tín hiệu sinh lý thường gặp ở các bộ lọc thông thấp IIR hay FIR cổ điển.

Lý do nghiên cứu chỉ lựa chọn 8 kênh trên tổng số 20 kênh đo của hệ thống FOIRE-3000 là gì?

Hệ thống FOIRE-3000 đo 20 kênh phân bố rộng trên hai bán cầu não. Tuy nhiên, 8 kênh được chọn gồm kênh 2, 5, 6, 9 ở bán cầu trái và 12, 15, 16, 19 ở bán cầu phải nằm tập trung chính xác tại tâm điểm của vùng vỏ não điều khiển vận động chi trên thuộc thùy đỉnh. Việc chọn lọc này giúp loại bỏ các kênh biên không mang tín hiệu vận động, giảm 60% khối lượng tính toán và tránh gây nhiễu cho mô hình học máy.

Phương pháp hồi quy đa thức đóng vai trò then chốt như thế nào trong bài toán phân loại?

Hồi quy đa thức bậc 3 biến đổi toàn bộ mảng dữ liệu nồng độ Oxy-Hb biến thiên liên tục trong 20 giây thực hiện tác vụ thành một tập hợp hữu hạn các hệ số đại số. Các hệ số này phản ánh đầy đủ tốc độ gia tăng, độ dốc và biên độ cực đại của lưu lượng máu. Nhờ đó, số chiều đặc trưng ngõ vào của thuật toán SVM giảm mạnh, ngăn ngừa hiện tượng quá khớp và nâng cao độ chính xác phân loại.

Thuật toán PR-SVM và PR-ANN đề xuất có khả năng triển khai trên các thiết bị nhúng thời gian thực không?

Cả hai mô hình PR-SVM và PR-ANN đều sử dụng vector đặc trưng ngõ vào có kích thước nhỏ gọn gồm 48 nút sau khi hồi quy. Quá trình tính toán kiểm tra chỉ bao gồm các phép nhân ma trận và tính tích vô hướng đơn giản, cho phép hệ thống đưa ra quyết định phân loại trong thời gian dưới 50 mili-giây trên các vi xử lý nhúng phổ thông, hoàn toàn đáp ứng tiêu chuẩn điều khiển thiết bị phục hồi chức năng thời gian thực.

Kết luận

  • Luận văn đã xây dựng thành công quy trình công nghệ hoàn chỉnh để nhận dạng chính xác hoạt động gõ tay trái và tay phải dựa trên tín hiệu quang phổ cận hồng ngoại chức năng fNIRS.
  • Ứng dụng xuất sắc bộ lọc đa thức Savitzky-Golay với cửa sổ 11 mẫu và bậc 3, giúp khử sạch nhiễu chuyển động và làm phẳng tín hiệu huyết động sinh học đạt tần số cắt 0,1977 Hz.
  • Đề xuất phương pháp trích xuất đặc trưng bằng hồi quy đa thức tối ưu, cô đọng thành công chuỗi thời gian Oxy-Hb thành các hệ số toán học giàu thông tin và giảm số chiều dữ liệu ngõ vào.
  • Chứng minh hiệu năng vượt trội của hai giải thuật PR-SVM và PR-ANN với tỷ lệ nhận dạng chính xác trung bình đạt từ 88% đến 93% trên các tập dữ liệu thực nghiệm độc lập.
  • Đóng góp nguồn học liệu thực nghiệm giá trị cao và mở ra hướng đi mới cho các ứng dụng giao diện não - máy tính không xâm nhập phục vụ y tế tại Việt Nam.

Đóng góp lớn nhất của công trình là giải quyết đồng thời bài toán khử nhiễu và tối ưu hóa không gian đặc trưng huyết động não bộ, tạo tiền đề vững chắc cho việc thương mại hóa các thiết bị BCI trong tương lai. Lộ trình phát triển tiếp theo trong 12 đến 24 tháng tới sẽ tập trung hoàn thiện mô-đun phân loại trực tuyến thời gian thực và thử nghiệm lâm sàng trên diện rộng. Các viện nghiên cứu, bệnh viện và doanh nghiệp công nghệ y tế quan tâm được khuyến khích hợp tác chuyển giao công nghệ để sớm đưa kết quả nghiên cứu vào ứng dụng thực tế, hỗ trợ đắc lực cho quá trình điều trị và phục hồi vận động của người bệnh.