(TIỂU LUẬN) báo cáo bài tập lớn đại số TUYẾN TÍNH giới thiệu phép biến đổi haar viết chương trình sử dụng phép biến đổi haar để nén dữ liệu

BỘ GIÁO DỤC VÀ ĐÀO TẠO TRƯỜNG ĐẠI HỌC BÁCH KHOA TP HCM ĐẠI HỌC QUỐC GIA TP HỒ CHÍ MINH BÁO CÁO BÀI TẬP LỚN ĐẠI SỐ TUYẾN TÍNH ĐỀ TÀI 23 Giới thiệu phép biến đổi Haar Viết chương trình sử dụng phép biến đổi Haar để nén liệu MT1007_L19_N09 BỘ GIÁO DỤC VÀ ĐÀO TẠO TRƯỜNG ĐẠI HỌC BÁCH KHOA TP HCM ĐẠI HỌC QUỐC GIA TP HỒ CHÍ MINH BÁO CÁO BÀI TẬP LỚN ĐẠI SỐ TUYẾN TÍNH ĐỀ TÀI 23 Giới thiệu phép biến đổi Haar Viết chương trình sử dụng phép biến đổi Haar để nén liệu NHÓM 09 – ĐSTT(MT1007)_L19 STT Lê Hoàn Lê Gia H Phạm H Trần Ho Đỗ Trần Văn Phú MỤC LỤC DANH SÁCH THÀNH VIÊN …………………………… MỤC LỤC …………………………………………………… LỜI MỞ ĐẦU ………………………………………………… CHƯƠNG I: CƠ SỞ LÝ THUYẾT I AI GIỚI THIỆU CHUNG …………………………… ………… Sơ lược thuật toán haar ……………………………………… Tác dụng ….…………………………………………………… CÁCH GIẢI THUẬT TOÁN ……………………… ……… Cách giải thuật toán ……… ………………………………… Nén liệu …….………………………………………… Nén không liệu ……….………………………………… 4 10 CHƯƠNG II: MATLAB I AI BI GIỚI THIỆU MATLAB ………………………………… … 15 CÁC LỆNH MATLAB ĐƯỢC SỬ DỤNG …………………… 15 ĐOẠN CODE MATLAB ……………………………………… 16 CHƯƠNG III: TỔNG KẾT I AI BI KIẾN THỨC …………………………………………………… 17 KĨ NĂNG ……………………………………………………… 17 TÀI LIỆU THAM KHẢO ……………………………………… 17 LỜI CẢM ƠN …………………………… …………… 18 LỜI MỞ ĐẦU Cuộc sống phát triển nhu cầu thông tin người phong phú, dẫn đến phát triển mạnh mẽ khoa học kĩ thuật, loại hình thơng tin vơ tuyến, hình thức xử lí tín hiệu, đặc biệt cơng nghệ xử lí ảnh Xử lí ảnh ngành khoa học tương đối mẻ so với nhiều ngành khoa học khác tập trung nghiên cứu phát triển ứng dụng thực tiễn nhiều ngành, lĩnh vực khác Trong “nén ảnh” phần xử lí ảnh có ứng dụng to lớn truyền thông lưu trữ, có nhiều phương pháp nén ảnh đời không ngừng cải tiến để ngày hoàn thiện đem lại hiệu nén cao cho chất lượng ảnh tốt Nén ảnh kĩ thuật mã hóa ảnh số hóa nhằm giảm số lượng bit liệu cần thiết để biểu diễn ảnh Vấn đề đặt yêu cầu ngày cao việc xử lí tín hiệu để đảm bảo vừa nén liệu, tiết kiệm dung lượng đường truyền tín hiệu, vừa đảm bảo trừ nhiễu tín hiệu có khả khơi phục lại tín hiệu với chất lượng tốt Có nhiều phương pháp xử lí tín hiệu với nhiều thuật toán, biến đổi toán học nghiên cứu Trong số đó, biến đổi Haar xem phép biến đổi mới, có nhiều tiềm năng, phát triển mạnh mẽ với ưu điểm vượt trội so với phép biến đổi truyền thống Haar cho phép phân tích tín hiệu miền thời gian tần số Do đó, biến đổi Haar ứng dụng rộng rãi nhiều lĩnh vực, từ y sinh tới công nghệ xử lí ảnh Trong báo cáo này, chúng em xin phép giới thiệu về: “Giới thiệu phép biến đổi Haar viết phương trình sử dụng phép biến đổi Haar để nén liệu” CHƯƠNG I: CƠ SỞ LÝ THUYẾT I GIỚI THIỆU CHUNG Sơ lược thuật toán Haar Biến đổi Haar chức biến đổi lâu đời nhất, đề xuất vào năm 1910 nhà toán học người Hungary, Alfréd Haar Nó tìm thấy hiệu ứng dụng nén tín hiệu hình ảnh kỹ thuật điện máy tính cung cấp cách tiếp cận đơn giản hiệu mặt tính tốn để phân tích khía cạnh địa phương tín hiệu Hàm biến đổi nén Haar cách đơn giản để thể hai phương pháp nén liệu nén không liệu Nó phụ thuộc vào trung bình giá trị khác ma trận hình ảnh tạo ma trận thưa gần thưa Một ma trận thưa lưu trữ cách hiệu quả, dẫn tới việc làm nhỏ kích cỡ ma trận Tác dụng Việc sử dụng máy tính cho cơng việc khác tăng lên chóng mặt Với đời máy ảnh kĩ thuật số, ứng dụng phổ biến lưu trữ, thao tác truyền hình ảnh kỹ thuật số Tuy nhiên, tệp bao gồm hình ảnh lớn nhanh chóng chiếm dung lượng nhớ ổ cứng máy tính Một hình ảnh tỷ lệ xám 256 x 256 pixel có 65.636 phần tử để lưu trữ hay hình ảnh màu 640 x 480 điển hình có gần triệu Kích thước tệp khiến trình tải xuống từ Internet trở nên lâu chậm Biến đổi Haar wavelet cung cấp phương tiện mà nén hình ảnh để chiếm khơng gian lưu trữ hơn, truyền điện tử nhanh mức độ chi tiết tăng dần AI CÁCH GIẢI THUẬT TOÁN Cách giải thuật toán Ở đây, tập trung vào thước xám (Gray Scale hiểu theo tiếng Việt có nghĩa thước xám hay cịn gọi “xám chuẩn”, “thẻ xám”, dụng cụ giúp kiểm tra, đánh giá độ bền màu sản phẩm vải, sản phẩm nhuộm.); nhiên, ảnh theo dạng rbg kiểm sốt cách nén màu khác theo lớp riêng biệt khác Phương pháp để bắt đầu với ảnh A, ta xem ma trận cỡ m ×n , với giá trị từ đến 255 Trong matlab, ma trận với 8-bit giá trị nguyên Chúng ta chia nhỏ ảnh thành × Ví dụ : Dưới hình ảnh thang độ xám 512 × 512 pixel trụ bay Nhà thờ Đức Bà Paris: Và ma trận × từ ảnh trên:  88      A=         Chúng ta tập trung vào r = (88 88 89 90 hàng đầu: Quá trình chuyển đổi diễn bước, với bước nhóm theo cặp: [88,88]; [89,90]; [96,97] [92,94]; Tiếp theo, thay cột r1 giá trị nửa tổng cột cuối nửa hiệu cặp Hàng ký hiệu r1h1 : r1h1 = ( 8889.5 r1 93 96.5 −0.5 −1 −0.5) mục gọi hệ số xấp xỉ mục cuối gọi hệ số chi tiết Tiếp theo, nhóm cột hàng [88,89.5] nửa tổng cặp số cột sau nửa hiệu cặp số Chúng ta giữ nguyên cột cuối r1h1 Và biểu thị hàng thứ hai r1h1h2 : r1h1h2 = ( 88.75 94.75 −0.75 −1.75 −0.5 −1 −0.5) Cuối cùng, nhóm mục r1h1h2 lại với [88.75,94.75] thay cột r1h1h2 nửa tổng cột thứ hai r1h1h2 nửa hiệu cặp số Ta giữ nguyên cột cuối r1h1h2 Ta ký hiệu hàng cuối r1h1h2 h3 r1h1h2 h3 = ( 91.75 −3 −0.75 −1.75 −0.5 −1 −0.5 ) Ta lặp lại quy trình cho hàng cịn lại A Sau đó, lặp lại quy trình tương tự cho cột A, nhóm hàng theo cách với cột Ma trận kết là:    96 −2.4375   −1.125 −2.6875   −0.6875    −0.1875   −0.875  −1.25 Lưu ý ma trận kết có số mục nhập hầu hết mục nhập lại gần Điều kết khác biệt thực tế pixel liền kề hình ảnh nói chung khơng khác nhiều   Bây giờ, ta thảo luận cách thực trình cách sử dụng phép nhân ma trận đây, ma trận biến đổi Haar wavelet Nếu ta có:  0.5     0.5 H =           AH1 tương đương với bước áp dụng cho tất hàng A AH −0.5 0 0.5 0.5 0.5 AH HH Nếu ta đặt H tích ma trận này, thì: H = HH H Để áp dụng quy trình cho cột, cần nhân A bên trái với Ma trận kết là: H       AH = H T    96 −2.43 −1.1 −2.68 −0.68 −0.18  −0.8    −1.2  Phần quan trọng trình phép biến đổi nghịch đảo Đặc biệt, H −1 A= ma trận khả nghịch Nếu B = H AH ) nên (HT nén A Hơn nữa, H khả nghịch nén khơng liệu phép BH Nén liệu Phép biến đổi Haar biểu diễn ma trận nén liệu để giữ lại chất lượng hình ảnh Đầu tiên, tỷ lệ nén hình ảnh tỷ lệ phần tử khác ảnh gốc với phần tử khác ảnh nén Gọi A ma trận dạng × H T AH wavelet H phạm vi từ đến ε biến thành số Ma trận có nhiều giá trị thể hình ảnh nén nhiều Khi giải nén hình ảnh cách sử dụng biến đổi Haar wavelet ngược, thu kết với hình ảnh gần với gốc Ví dụ, với ma trận A lấy từ ví dụ, chọn ε = 0.25 , có kết với ma trận:       96 0 −2.4375 −0.3125   −1.125 −0.375 −2.6875 0 0 0 −0.6875     −0.875   −1.25  Ma trận đại diện cho tỉ lệ nén 27 −0.3125 48 0 0 =1.7 (phần tử số 48 tức số 96 chia đôi phần mẫu số 27 số khác ma trận trên) Các ảnh ví dụ hình ảnh ban đầu, nén theo tỷ lệ khác nhau: Ảnh gốc Tỷ lệ nén 10:1 Tỷ lệ nén 30:1 Tỷ lệ nén 50:1 Nén khơng liệu Ví dụ: Xét phép nén liệu dựa biến đổi Haar Giả sử ta có đoạn liệu X0 = (154 150 156 152 160 160 152 156)T Ta chia liệu thành cặp tìm nửa tổng a nửa hiệu d cặp Khi ta dãy X1 = (152 154 160 154 2 −2)T với số nửa tổng số nửa hiệu cặp Có thể dùng ma trận để mơ tả q trình sau: 10  0.5      0  A=  0.5        Các cột A tạo nên họ trực giao Nếu ta chia cột cho độ dài họ trực chuẩn ma trận tương ứng ma trận trực giao H sau:                      −      Khi nghịch đảo H HT Ta dùng ma trận H để nén liệu thay cho A Ta có Y =HX Giải nén liệu ta được, quan trọng bảo tồn khoảng cách góc Thực vậy, ta có: HX = (HX)T (HX) = XTHTHX = XTX = X Tức độ dài vectơ X độ dài vectơ HX Ngồi góc hai vectơ u v thỏa cosθ = (u , v) u.v Qua phép biến đổi H, góc hai vectơ Hu Hv thỏa 11 cos α = Hay góc hai vectơ u v với góc hai ảnh Hu Phép biến đổi trực giao H bảo toàn góc khoảng cách nên khơng làm thay đổi hình dạng ảnh trình nén Tiếp tục trình nén sau Vector Y=H Với H1 ma trận H ví dụ trước, có phần tử sau nhỏ cị Giữ bốn số sau dãy lại Đối với bốn số đầu, ta chia làm hai cặp, tìm nửa tổng nửa Tiếp tục trình trên, giữ số cuối Y2 hiệu chúng, ta có Y3 = H 3Y2 = ( 310 lại Hai số đầu thay nửa tổng nửa −4 12 Tóm lại ta có Y3 = H H H1 X ⇔ Y3 = QX0 Véctơ Y3 Giải nén dãy liệu − − Dưới ảnh minh họa cho phép biến đổi nén không liệu: 13 Ảnh gốc Ảnh sau lần biến đổi Ảnh sau lần biến đổi Ảnh sau lần biến đổi 14 CHƯƠNG II: MATLAB I GIỚI THIỆU MATLAB MATLAB môi trường tính tốn số lập trình, thiết kế cơng ty MathWorks MATLAB cho phép tính tốn số với ma trận, vẽ đồ thị hàm số hay biểu đồ thơng tin, thực thuật tốn, tạo giao diện người dùng liên kết với chương trình máy tính viết nhiều ngơn ngữ lập trình khác MATLAB giúp đơn giản hóa việc giải tốn tính tốn kĩ thuật so với ngơn ngữ lập trình truyền thống C, C++, Fortran MATLAB sử dụng nhiều lĩnh vực, bao gồm xử lý tín hiệu ảnh, truyền thơng, thiết kế điều khiển tự động, đo lường kiểm tra, phân tích mơ hình tài chính, tính tốn sinh học, bên cạnh đó, MATLAB cịn sử dụng để phát triển thuật tốn Với hàng triệu kĩ sư nhà khoa học làm việc môi trường công nghiệp môi trường hàn lâm, MATLAB ngơn ngữ tính tốn khoa học Thông qua báo cáo chúng em tìm hiểu ứng dụng MATLAB Đại số tuyến tính, cụ thể tìm hiểu phép biến đổi Haar ứng dụng vào việc nén liệu, đồng thời viết chương trình MATLAB sử dụng phép biến đổi Haar để nén hình ảnh 15 II CÁC LỆNH MATLAB ĐƯỢC SỬ DỤNG Lệnh clc clear all close all imread('filename') rgb2gray(image) imresize(A,[NUMROWS NUMCOLS]) imagesc( ) im2double( ) imwrite(A,'filename') colormap gray axis off axis square title('text') length(X) norm( ) zeros(N) subplot(2,1,n) 16 BI ĐOẠN CODE MATLAB clc clear all close all IM = imread('anh.jpg'); coverIM = rgb2gray(IM); s=size(coverIM); coverIM = imresize(im2double(coverIM),[s(1,1)/2 s(1,1)/2]); subplot(2,1,1); figure(1) imagesc(coverIM) colormap gray title('Hinh anh ban dau') axis off axis square dim=length(coverIM); m=dim/2; At=zeros(dim); for i=1:m At(i,2*i-1)=1/2; At(i,2*i)=1/2; At(i+m,2*i-1)=1/2; At(i+m,2*i)=-1/2; end for i=1:dim H1(:,i)=At(:,i)/(norm(At(:,i))); end Y=H1*coverIM; figure(1) subplot(2,1,2); imagesc(Y) colormap gray title('Hinh anh sau bien doi Haar') axis off axis square imwrite(Y,'anh_sau.jpg') 17 CHƯƠNG III: TỔNG KẾT I KIẾN THỨC Dựa vào kiến thức lí thuyết học để thực phép biến đổi Haar - Tìm hiểu sử dụng phần mềm Matlab viết chương trình sử dụng phép biển dổi Haar để nén liệu - II KĨ NĂNG Sau thời gian làm việc chung, nhóm đạt kết sau: + Nhóm hoạt động ổn định, có tổ chức, tôn trọng lẫn nhau, thường xuyên tương tác, trao đổi qua lại thành viên + Thành viên có trách nhiệm, tham gia tất bàn luận hồn thành nhiệm vụ phân cơng + Mọi cơng việc rõ ràng nên nhóm hoạt động cách tối ưu đạt hiệu công việc cao + Các thành viên nâng cao kĩ làm việc nhóm khả giao tiếp, thảo luận với nhau, đưa quan điểm cá nhân,… + Giữa thành viên có thống nhất, người có tính kỷ luật, tn theo nội quy nhóm BI TÀI LIỆU THAM KHẢO Đặng Văn Vinh, Giáo trình Đại Số Tuyến Tính, NXB Đại học Quốc Gia TP Hồ Chí Minh, 2020 David R Bull and Fan Zhang, Intelligent Image and Video Compression Comunicating Pictures, Second Edition Greg Ames, Image Compression, 07/12/2002 https://web.archive.org/web/20110125080404/http://online.redwoods.cc.ca.us/instruct/dar nold/laproj/Fall2002/ames/paper.pdf Help Center for MATLAB, Simulink and other MathWorks products 18 LỜI CẢM ƠN Trước tiên với tình cảm chân thành nhất, chúng em xin bày tỏ lòng biết ơn đến tất cá nhân tổ chức tạo điều kiện hỗ trợ, giúp đỡ chúng em suốt trình học tập nghiên cứu đề tài Trong thời gian từ bắt đầu học tập trường đến nay, chúng em nhận nhiều quan tâm, giúp đỡ quý Thầy Cô bạn bè Với lòng biết ơn sâu sắc nhất, chúng em xin gửi lời cảm ơn đến thầy Bùi Anh Tuấn truyền đạt vốn kiến thức quý báu cho chúng em thời gian học mơn Đại số tuyến tính Nhờ có lời hướng dẫn, giảng dạy Thầy nên đề tài tập lớn chúng em hồn thành Bản báo cáo thực khoảng thời gian ngắn nên hạn chế chúng em cịn nhiều bỡ ngỡ, khơng tránh khỏi thiếu sót, mong nhận ý kiến đóng góp q báu q Thầy Cơ bạn để kiến thức chúng em cải thiện đồng thời có điều kiện bổ sung, nâng cao kĩ Chúng em xin chân thành cảm ơn! 19 ... TRƯỜNG ĐẠI HỌC BÁCH KHOA TP HCM ĐẠI HỌC QUỐC GIA TP HỒ CHÍ MINH BÁO CÁO BÀI TẬP LỚN ĐẠI SỐ TUYẾN TÍNH ĐỀ TÀI 23 Giới thiệu phép biến đổi Haar Viết chương trình sử dụng phép biến đổi Haar để nén liệu. .. em tìm hiểu ứng dụng MATLAB Đại số tuyến tính, cụ thể tìm hiểu phép biến đổi Haar ứng dụng vào việc nén liệu, đồng thời viết chương trình MATLAB sử dụng phép biến đổi Haar để nén hình ảnh 15... biến đổi Haar viết phương trình sử dụng phép biến đổi Haar để nén liệu? ?? CHƯƠNG I: CƠ SỞ LÝ THUYẾT I GIỚI THIỆU CHUNG Sơ lược thuật toán Haar Biến đổi Haar chức biến đổi lâu đời nhất, đề xuất vào