Bai tap lon báo cáo xử lý ảnh
|
(1) Show
ĐẠI HỌC BÁCH KHOA HÀ NỘIVIỆN CÔNG NGHỆ THÔNG TIN VÀ TRUYỀN THÔNG ------BÀI TẬP LỚN MÔN XỬ LÝ ẢNHĐỀ TÀI:Tìm hiểu các đặc trưng sinh trắc ảnh khuôn mặt, nghiên cứu ứngdụng của phép biến đổi KL và phân tích thành các thành phầnchính PCA trong trích chọn đặc trưng khuôn mặtGiảng viên hướng dẫn : PGS.TS Nguyễn Thị Hoàng Lan Sinh viên thực hiện : Nguyễn Văn Thành SHSV : 20072604 Lớp: Hệ thống thông tin và truyền thông – KSCLC-K52 (2) MỤC LỤCLỜI NÓI ĐẦU ... 3 I. ĐẶC TRƯNG SINH TRẮC ẢNH KHUÔN MẶT ... 4 II. NGHIÊN CỨU ỨNG DỤNG CỦA PHÉP BIẾN ĐỔI KL ... 7 1. Phép biến đổi KL ... 7 2. Ứng dụng của phép biến đổi KL ... 9 III. PHÂN TÍCH THÀNH CÁC THÀNH PHẦN CHÍNH PCA TRONG TRÍCH CHỌN ĐẶC TRƯNG KHUÔN MẶT ... 14 1. Phép biến đổi PCA... 14 2. Phân tích thành chính PCA trong trích chọn đặc trưng khuôn mặt ... 14 2.1. Tính toán các vector riêng ... 14 2.2. Biểu diễn khuôn mặt theo cơ sở tìm được ... 17 IV. KẾT LUẬN ... 18 (3) LỜI NÓI ĐẦU “Xử lý ảnh” là môn học quan trọng đối với sinh viên ngành công nghệ thôngtin. Đây là môn học khó đối với hầu hết nhiều sinh viên do yêu cầu kiến thức vềtoán và xác suất.Với mục tiêu có thêm kiến thức cơ bản trong lĩnh vực về xử lý ảnh, em chọnđề tài “Tìm hiểu các đặc trưng sinh trắc ảnh khuôn mặt, nghiên cứu ứngdụng của phép biến đổi KL và phân tích thành các thành phần chính PCAtrong trích chọn đặc trưng khuôn mặt”. Qua bài tập lớn, em đã có thêm đượccái nhìn chung về hệ thống nhận dạng khuôn mặt, các đặc trưng sinh trắc khuônmặt giúp phát hiện và trích rút được các đặc trưng cho việc nhận dạng. Đồngthời em được củng cố thêm về kiến thức toán và xác suất thông kê, có được hiểubiết tốt hơn về ứng dụng của kiến thức cơ bản trong các bài toán thực tế.Dù rất cố gắng trong việc tìm hiểu tài liệu, nhưng do thiếu sót về kiến thức cơbản, hạn chế tầm hiểu biết mà bài báo cáo còn rất nhiều thiếu khuyết. Kínhmong nhận được những ý kiến góp ý để em hoàn thiện hơn.Nhân dịp này, em xin gửi lời cảm ơn chân thành tới PGS.TS Nguyễn ThịHoàng Lan đã nhiệt tình hướng dẫn góp ý giúp em hoàn thành báo cáo môn họcnày.Em xin chân thành cảm ơn!(4) I. ĐẶC TRƯNG SINH TRẮC ẢNH KHUÔN MẶT Các đặc trưng khuôn mặt bao gồm: - Đặc trưng hình học: cấu trúc, hình dạng và các thành phần trên khuôn mặt: miệng, mắt, mũi, lông mày. Khoảng cách giữa mắt, mũi, miêng và hàm; đường bao các hốc mắt; các cạnh của miệng; vị trí của mũi, hai mắt và các vùng xung quanh.Các thành phần khuôn mặt được trích rút để hình thành vector đặc trưng biểu diễn hình học khuôn mặt. (5) Hình 1.1. Minh họa đặc trưng hình học của khuôn mặt - Đặc trưng về diện mạo biểu diễn sự thay đổi vẻ bề ngoài: kết cấu da như các nếp nhăn trên khuôn mặt; biểu đồ nhiệt của khuôn mặt: các mẫu nhiệt khuôn mặt là duy nhất với mỗi người và đặc trưng về nụ cười. Các đặc trưng về diện mạo có thể được trích rút trên cả khuôn mặt hoặc phần nào đó trên khuôn mặt. (6) (7) II. NGHIÊN CỨU ỨNG DỤNG CỦA PHÉP BIẾN ĐỔI KL 1. Phép biến đổi KL Xét không gian mẫu S = {x} gồm n vector dữ liệu mẫu. Trong đó x là vector các biến ngẫu nhiên N chiều: x = [ x1 x2 … xN ]T Phép biến đổi KL đối với các vector x được định nghĩa như sau: y = WTx (1.1) Trong đó x = [ x1 x2 … xN ]T, y = [ y1 y2 … yM ]T và ma trận W là ma trận phép biến đổi với kích thước NxM (M ≤ N) có dạng: MN M N T i T wwwwwW...............1 1 11 Ta định nghĩa ma trận hiệp tương quan của các vector x: [C]x = E [xxT ](1.2) Ma trận tương quan được ước lượng từ n mẫu trong không gian dữ liệu quan sát được xác định bởi biểu thức: n i T i i x x x n C 1 1 Phép biến đổi KL là bài toán tìm một ma trận biến đổi W thỏa mãn (1.1) Mỗi cột iN i i www...1 của W là vector cơ sở trực giao của không gian mới hay: 1 . Tk j w (8) Vì vậy mỗi phần tử yi của y sẽ được tính: N iN i i T i i wxwxwxwxy1 1 2 2 ...(1.4) Do wi là các vector cơ sở trực giao theo (1.3) nên W là ma trận trực chuẩn thỏa mãn: WTW = I = WWT (I là ma trận đơn vị) (1.5) Từ đó: WT = W-1 và ta có dạng biến đổi ngược của (1.1): x = Wy (1.6) Đặt Cy là ma trận đường chéo mong muốn của vector biến ngẫu nhiên y: Cy = N ... ... 0 ... 1 Trong đó các phần tử đường chéo là các phương sai của dữ liệu được biến đổi. Ma trận đường chéo này có thể tính toàn từ ma trận hiệp tương quan gốc như sau: Cy = T yy E = T T T x W x W E = E WT xxT W = WTCxW Hay WCy= CxW (1.7) Cy là ma trận đường chéo nên biểu thức (1.7) được đưa về sạng Cxwi = iwi (1.8) Từ (1.8) ta thấy λi và wi là cặp giá trị riêng và vector riêng của ma trận tương quan Cx trong biểu thức (1.2). Hay các cột wi của ma trận W là các vector riêng của ma trận Cx. Để đảm bảo W là trực chuẩn cần áp dụng trực giao và chuẩn hóa Gram-Schmidt với các vector riêng tìm được. Như vậy, phép biến đổi KL chính là tìm các vector riêng wi ứng với M giá trị riêng lớn nhất của ma trận hiệp tương quan của các vector ngẫu nhiên quan sát được. Phép biển đổi KL là cách tối ưu nhằm giảm thứ nguyên từ không gian dữ liệu có số chiều lớn (9) thành không gian mới có số chiều bé hơn rất nhiều với sai số là bé nhất nằm ở các vector riêng ứng với các giá trị riêng nhỏ nhất. 2. Ứng dụng của phép biến đổi KL Hình 2.1 biểu diễn một ảnh đa mức xám kích thước 512 x 512 với giá trị mức xám của mỗi điểm ảnh được biểu diễn bằng 8 bit (giá trị trong [0 – 255]). Giá trị mức xám của các điểm ảnh kề nhau có xu hướng tương tự nhau. Hình 2.2 minh họa giá trị mức xám của các cặp điểm ảnh liền kề. Hình 2.1. Ảnh đa mức xám kích thước 512 x 512 Trong hình 2.2, mỗi chấm biểu diễn một điểm ảnh trong bức ảnh hình 2.1 với hoành độ x là giá trị của điểm ảnh đó và tung độ y là giá trị của điểm ảnh liền kề bên phải nó.Từ đồ thị, ta thấy quan hệ x = y thể hiện độ tương quan mạnh giữa các điểm ảnh liền kề. Chia bức ảnh thành các khối 1x2 rời nhau như hình 2.3, chúng ta biểu diễn (10) một bức ảnh thành tập các vector 2 chiều xi. Phân bố giá trị mức xám của mỗi thành phần được vẽ như hình 2.4. Chúng ta thấy rằng phân bố mức xám của mỗi thành phần tương đối rộng và phủ hầu hết dải 0 – 255. Hơn nữa, hai phân bố này khá giống với phân bố chung của mỗi điểm ảnh trong bức ảnh hình 2.1. (11) (12) Hình 2.4. Phân bố giá trị mức xám của mỗi thành phần trong các khối Bây giờ, chúng ta quay phân bố trong hình 2.2 đi một góc 45o. Kết quả trong hình 2.5 cho thấy, hai thành phần mới không tương quan nghĩa là biết giá trị của thành phần thứ nhất sẽ không giúp ta xác định được giá trị của thành phần thứ hai. Phân bố của hai thành phần mới được vẽ như hình 2.6. Thành phần thứ nhất vẫn khá giống với phân bố trước, tức là phân bố rộng, trải hầu hết khoảng giá trị. Tuy nhiên thành phần thứ hai thì khác, nó hẹp hơn rất nhiều và giá trị đỉnh tại 0. Do nó có khoảng biến đổi nhỏ hơn nhiều nên chúng ta cần ít số bit để mã hóa giá trị của nó. Vì vậy, chúng ta có thể giảm số bit cần thiết để mã hóa một một ảnh khi giải tương quan. (13) Hình 2.6. Phân bố mức xám của hai thành phần khi xoay tương ứng với hình 2.5 Hầu hết các bức ảnh đều chứa sự tương quan ngẫu nhiên về dữ liệu do đó dẫn đến sự dư thừa dữ liệu. Phép biến đổi KL loại bỏ được sự dư thừa dữ liệu nhờ việc giải tương quan dữ liệu và giảm số chiều vì vậy ảnh có thể được lưu trữ hiệu quả hơn. Đó là ứng dụng của phép biển đổi KL. (14) III. PHÂN TÍCH THÀNH CÁC THÀNH PHẦN CHÍNH PCA TRONG TRÍCH CHỌN ĐẶC TRƯNG KHUÔN MẶT 1. Phép biến đổi PCA Tương tự như phép biến đổi KL, ta có không gian dữ liệu quan sát S = {x} gồm n vector dữ liệu mẫu N chiều. Dữ liệu thô tồn tại sự tương quan ngẫu nhiên giữa các thành phần do đó có sự dư thừa dữ liệu. Ý tưởng của phép biến đổi PCA là phân tích dữ liệu thành các thành phần không tương quan (gọi là các thành phần chính) để giảm độ dư thừa dữ liệu. Phép biến đổi PCA được định nghĩa như sau: uT x(2.1) Trong đó x = [ x1 x2 … xN ]T , u = [ u1 u2 … uN ]T – các thành phần thứ i và j không tương quan trong không gian mới. Và ma trận Ψ là ma trận phép biến đổi với kích thước N2 có dạng: NN N N T i T v...............1 1 11 (2.2) Ở đó vi là các vector riêng tương ứng với ma trận hiệp phương sai của các x quan sát được. T x x E C (2.3) Trong đó N ... 1 và n k ik i xn1 12. Phân tích thành chính PCA trong trích chọn đặc trưng khuôn mặt 2.1. Tính toán các vector riêng Giả sử chúng ta một tập luyện gồm M bức ảnh khuôn mặt I1, I2…IM. Các bức ảnh (15) Chúng ta biểu diễn mỗi bức ảnh Ii kích thước NxN bằng một vector có kích thước N2 chiều. (16) Bước 2: Tính i i với i = 1, 2…N2 Bước 3: Tính ma trân hiệp phương sai C của các vector quan sát: T M i T i i AAMC1 1(ma trận C có kích thước N2xN2) (2.4) Trong đó A 12...M (kích thước N2xM) Bước 4: Tìm các vector riêng ui của ma trận hiệp phương sai C hay của AAT. Tuy nhiên kích thước của ma trận này là N2xN2 quá lớn nên việc tìm vector riêng của ma trận này là không khả thi. Chúng ta xem xét ma trận ATA có kích thước MxM Tìm các vector riêng vi của ma trận ATA, ta có: i i i T vAvA(2.5) Quan hệ giữa ui và vi: i i i i i i T i i i T AvCAvAvAvAAvAvA(2.6) Từ đó suy ra ui = Avi Vì vậy, AAT và ATA có cùng các giá trị riêng và các vector riêng của chúng quan hệ với nhau theo ui = Avi Chú ý: - ATA có thể có M giá trị riêng và vector riêng. - AAT có thể có N2 giá trị riêng và vector riêng. - M giá trị riêng của ATA cùng với các vector riêng tương ứng với M giá trị riêng lớn nhất của AAT. (17) Sau đó tính M vector riêng ứng với M giá trị riêng lớn nhất của AAT: u i = Avi đồng thời chuẩn hóa vector ui sao cho ui 1 Bước 5: Chọn K vector riêng ứng với K giá trị riêng lớn nhất. 2.2. Biểu diễn khuôn mặt theo cơ sở tìm được Mỗi vector biểu diễn khuôn mặt (trừ đi vector trung bình) i trong tập luyện có thể được biểu diễn bằng một tổ hợp tuyến tính của K giá trị riêng tính ở trên. K j j j i meanwu1 ˆi=1,2…M. Trong đó i T j j uwchính là thành phần chính thứ j trong không gian mới. Và uj gọi là các ảnh riêng. Mỗi ảnh được chuẩn hóa trong tập luyện sẽ được biểu diễn trong cở sở này bởi vector: i K i i i www...2 1 trong đó i = 1, 2… M (18) IV. KẾT LUẬN Như vậy, em đã trình bày về các đặc trưng sinh trắc ảnh khuôn mặt, gồm các đặc trưng hình học và các đặc trưng về diện mạo khuôn mặt, cùng với đó là lý thuyết và ứng dụng của phép biến đổi KL và phân tích thành các thành phần chính PCA trong trích chọn đặc trưng khuôn mặt. Phép biến đổi KL ứng dụng trong bài toán giảm số chiều của không gian dữ liệu số lớn thành không gian có số chiều nhỏ hơn nhằm giảm giá tính toán và tăng hiệu quả của các kỹ thuật trong xử lý ảnh. Trong khi PCA giải tương tương quan dữ liệu và phân tích thành các thành phần chính ứng dụng trong trích chọn đặc trưng khuôn mặt. (19) TÀI LIỆU THAM KHẢO [1]. Bài giảng môn Xử lý ảnh – PGS.TS Nguyễn Hoàng Lan – Đại học Bách Khoa Hà Nội, 2010. [2]. Digital Image Processing – William K. Pratt- Wiley, 2007. [3]. Fundamental of Image Processing – Ian T. Young et al., 1998. [4]. Principal Component Analysis – I.T.Jolliffe – Springer, 2002. [5]. R.D. Dony “Karhunen-Lòeve Transform” - The Transform and Data Compression Handbook, 2001. |
