CHẾ TẠO THIẾT BỊ ĐẾM SẢN PHẨM BẰNG PHƯƠNG PHÁP XỬ LÝ ẢNH

BỘ GIÁO DỤC VÀ ĐÀO TẠO TRƯỜNG ĐẠI HỌC NÔNG LÂM THÀNH PHỐ HỒ CHÍ MINH KHĨA LUẬN TỐT NGHIỆP CHẾ TẠO THIẾT BỊ ĐẾM SẢN PHẨM BẰNG PHƯƠNG PHÁP XỬ LÝ ẢNH Họ tên sinh viên: TRẦN VỸ ANH Ngành: ĐIỀU KHIỂN TỰ ĐỘNG Niên khóa: 2007-2011 Tháng 05/2011 CHẾ TẠO THIẾT BỊ ĐẾM SẢN PHẨM BẰNG PHƯƠNG PHÁP XỬ LÝ ẢNH Tác giả TRẦN VỸ ANH Luận văn đệ trình để hồn tất u cầu cấp Kỹ sư Điều khiển tự động Ngành Điều khiển tự động Giáo viên hướng dẫn: Th.S Đặng Phi Vân Hài Tháng 05 năm 2011 i LỜI CẢM ƠN Em xin trân trọng cám ơn tất quý Thầy/Cô trường Đại học Nơng Lâm TP.Hồ Chí Minh, q Thầy/Cơ Khoa Cơ Khí trang bị cho em kiến thức quý báu giúp đỡ em suốt trình học tập trường Em xin chân thành cám ơn Thầy/Cô môn Điều Khiển Tự Động, cám ơn thầy Nguyễn Văn Bạn giúp đỡ em nhiệt tình học tập thời gian thực Đề tài Em xin bày tỏ biết ơn chân thành cô Đặng Phi Vân Hài tận tình hướng dẫn em suốt trình làm Luận văn tốt nghiệp Đặc biệt, em xin cám ơn Thầy/Cô Hội Đồng dành thời gian nhận xét góp ý để Luận văn em hoàn thiện Cuối cùng, xin kính gửi đến Cha Mẹ lòng biết ơn sâu sắc cơng ơn sinh thành dưỡng dục cho có ngày hơm Đồng thời, cám ơn người thân bạn bè tập thể lớp DH07TD động viên, ủng hộ tạo cho tơi điều kiện thuận lợi suốt q trình hoàn thành Luận văn TP.HCM, ngày 20 tháng 05 năm 2011 Sinh viên thực TRẦN VỸ ANH ii TÓM TẮT Trong khoảng thời gian 10 năm trở lại đây, việc ứng dụng xử lý ảnh vào sống phát triển, đặc biệt lĩnh vực tự động hóa sản xuất, như: Đếm số lượng sản phẩm, phân loại sản phẩm, nhận diện bảng số xe vào cổng, theo dõi phát triển môi trường sống loại thủy sản nuôi trồng thủy sản, v.v… Từ đó, với gợi ý thầy Lê Văn Bạn giúp đỡ cô Đặng Phi Vân Hài thầy môn Điều Khiển Tự Động, em thực đề tài với mong muốn: Đưa lý thuyết xử lý ảnh áp dụng vào quy trình đếm số lượng sản phẩm, nhằm góp phần tự động hóa sản xuất, tăng suất lao động Cụ thể, em nghiên cứu chế tạo thiết bị đếm số lượng cá giống với quan sát camera xử lý phần mềm Matlab Với mục đích cung cấp thiết bị cho hộ bán cá giống để phân phối cá cho hộ nuôi trồng thủy sản Do thời gian thực đề tài hạn chế, chưa sâu sát thực tế sản xuất nên đề tài chắn khơng tránh khỏi thiếu xót Em mong nhận đóng góp ý kiến quý Thầy Cô bạn bè để đề tài em hoàn thiện Em xin chân thành cám ơn iii MỤC LỤC Trang Trang tựa i LỜI CẢM ƠN ii TÓM TẮT iii MỤC LỤC .iv DANH SÁCH CÁC HÌNH vi CHƯƠNG I: MỞ ĐẦU 1.1 Tổng quan 1.2 Mục đích yêu cầu nghiên cứu 1.3 Thời gian địa điểm nghiên cứu đề tài 1.3.1 Địa điểm tiến hành đề tài 1.3.2 Phân bố thời gian tiến hành đề tài 1.4 Đối tượng thiết bị nghiên cứu 1.4.1 Đối tượng nghiên cứu 1.4.2 Thiết bị nghiên cứu 1.5 Phương pháp nghiên cứu CHƯƠNG II: TỔNG QUAN 2.1 Tổng quan xử lý ảnh 2.1.1 xử lý ảnh gì? 2.1.2 loại ảnh không gian màu 2.2 Tổng quan Matlab 2.2.1 Matlab gì? 2.2.2 Các cửa sổ làm việc Matlab 10 2.2.3 Các kiểu liệu MATLAB .10 2.2.4 Hệ thống lệnh MATLAB 11 2.2.4.1 Lệnh hệ thống Matlab .11 2.2.4.2 Lệnh số học .11 2.2.4.3 Lệnh cấu trúc 12 iv 2.2.5 Ứng dụng MATLAB 13 2.2.6 Ứng dụng Matlab xử lý ảnh .13 2.2.6.1 Image Acquisition Toolbox 13 a Thu nhận ảnh tĩnh vào Matlab 14 b.Thu nhận ảnh động vào Matlab 16 2.2.6.2 Image processing toolbox 16 2.3 Cấu tạo webcamera .19 CHƯƠNG II: KẾT QUẢ VÀ THẢO LUẬN 20 3.1 Phần cứng thiết bị 20 3.2 Phần mềm xử lý 22 3.2.1 Lưu đồ xử lý ảnh 22 3.2.2 Cấu trúc lệnh xử lý sử dụng .23 3.2.2.1 Các lệnh thu nhận ảnh 23 3.2.2.2 Các lệnh xử lý ảnh 25 a Chuyển ảnh RGB sang ảnh xám 25 b Lấy cạnh (edge) 26 c Xóa pixel nhiễu ảnh 27 d làm giãn ảnh (Dilate) 27 e làm co ảnh (erode) 28 f Đếm số lượng cá 28 3.2.3 Thiết kế giao diện GUI .29 CHƯƠNG IV: KẾT LUẬN VÀ ĐỀ NGHỊ 30 4.1 Kết đạt 30 4.2 Hướng phát triển đề tài 30 TÀI LIỆU THAM KHẢO 31 Phụ lục 1: Chương trình xử lý ảnh Trên Matlab 32 Phụ lục 2: Bảng khảo nghiệm số lượng cá đếm 33 Phụ lục 3: Matlab help 35 Phụ lục 4: Một số định dạng ảnh mà Matlab hỗ trợ .37 v DANH SÁCH CÁC HÌNH Hình 2.1 Q trình xử lý ảnh Hình 2.2 Sơ đồ tổng quát hệ thống xử lý ảnh Hình 2.3 Pixel ảnh (nguồn internet) Hình 2.4 Mơ hình ảnh RGB (nguồn internet) Hình 2.5 Mơ hình màu HSV (nguồn internet) Hình 2.6 Mơ hình màu L*a*b (nguồn internet) Hình 2.7 Biểu tượng Matlab Hình 2.8 Giao diện làm việc Matlab 10 Hình 2.9 Sơ đồ kiểu liệu 11 Hình 2.10 minh họa lệnh imread 14 Hình 2.11 minh họa lệnh imshow .15 Hinh 2.12 Ảnh chụp từ camera 16 Hình 2.13 Giãn ảnh .18 Hình 2.14 Co ảnh 18 Hình 3.1 Cấu tạo thiết bị đếm cá giống 20 Hình 3.2 Bản vẽ chi tiết thiết bị đếm cá giống 21 Hình 3.3 Lưu đồ xử lý ảnh 22 Hình 3.4 thông tin thiết bị camera 23 Hình 3.5 Ảnh động hiển thị Matlab .24 Hình 3.6 ảnh sau chụp 25 Hình 3.7 Chuyển ảnh RGB sang ảnh xám 26 Hình 3.8 lấy cạnh từ ảnh xám .26 Hình 3.9 xóa pixel nhiễu sau lấy cạnh 27 Hình 3.10 Làm giãn ảnh .27 Hình 3.11 kết sau co ảnh 28 Hình 3.12 kết sau dung lệnh đếm số lượng .28 Hình 3.13 Giao diện làm việc cho người sử dụng 29 vi CHƯƠNG I MỞ ĐẦU 1.1 Tổng quan Khi gia nhập WTO nơng nghiệp Việt Nam có nhiều hội để phát triển Và nói ngành thủy sản đóng vai trò quan trọng kinh tế quốc dân nước ta Có đặc điểm nguồn lợi thủy sản mang tính tái tạo, tái sinh Nhưng người khai thác qua khả tái sinh nguồn lợi bị cạn kiệt Trên thực tế sản lượng thủy sản mà người khai thác ngày bị suy giảm Nếu người khơng tiến hành giải pháp khác nguy cạn kiệt nguồn lợi điều dể dàng nhận thấy Vì vậy, ni trồng thủy sản vừa nhằm mục đích phục vụ nhu cầu nước đồng thời xuất nói giải pháp hữu hiệu giai đoạn phát triển hội nhập Theo số liệu tổng cục Thống kê Bộ Thủy sản, sản lượng thủy sản Việt Nam trước năm 1985 có sa sút, yếu từ năm 1985 trở lại liên tục phát triển qua năm Năm 1986, tong sản lượng thủy sản đạt 830.500 tấn; năm 1990 1.019.000 tấn; đến năm 2000 đạt 2.003.000 tấn; năm 2004 đạt 3.073.600 Trong khai thác hải sản đạt tương ứng 487.000 tấn; 709.000 tấn; 1.280.000 tấn; 1.230.500 nuôi trồng thủy sản 243.500 tấn; 310.000 tấn; 723.000 tấn; 1.150.000 Trong năm gần nghành thủy sản nước ta liên tục phát triển trở thành bốn mặt hàng xuất chủ lực Việt Nam Năm 2010 xuất thủy sản đạt 4,95 tỷ USD, chiếm 6,9% tổng kim nghạch xuất nước ta Xuất thủy sản Việt Nam trở thành tượng thị trường xuất thủy sản giới, sản lượng nuôi trồng đứng thứ 3, sau Trung Quốc Ấn Độ Hiện nay, việc sản xuất ươm giống thủy sản thực sở tư nhân, trại sản xuất cá tơm giống, có hỗ trợ viện, trung tâm nghiên cứu ngày chuyên nghiệp gắn với vùng ni Đã hình thành cụm ni trồng – chế biến – xuất thủy sản – phụ trợ tỉnh Đồng sơng Cửu Long Tính chun nghiệp động doanh nghiệp phân đoạn nâng lên Nhiều doanh nghiệp làm ăn bản, sử dụng quy trình cơng nghệ sản xuất đại, quy mơ lớn, áp dụng mơ hình quản trị tiên tiến Đến sản phẩm thủy sản Việt Nam có mặt 80 nước vùng lãnh thổ Trong Mỹ Nhật Bản thị trường lớn Với thống kê cho ta thấy, nghành thủy sản Việt Nam phát triển mạnh mẽ Từ đó, em tiến hành thực đề tài: “ chế tạo thiết bị đếm số lượng cá giống phương pháp xử lý ảnh “ hướng dẫn cô Đặng Phi Vân Hài, với mong muốn làm tăng suất lao động, giảm sai sót người đếm thủ cơng, góp phần tự động hóa sản xuất từ khâu cung cấp giống 1.2 Mục đích yêu cầu nghiên cứu Mục đích yêu cầu nghiên cứu nghiên cứu chế tạo thành công thiết bị đếm cá giống phương pháp xử lý ảnh ứng dụng phần mềm Matlab 1.3 Thời gian địa điểm nghiên cứu đề tài 1.3.1 Địa điểm tiến hành đề tài Đề tài tiến hành vào ngày 28/03/2011 – 22/06/2011 Trường Đại Học Nông Lâm Thành Phố Hồ Chí Minh Khoa Cơ Khí Trường Đại Học Bách Khoa Thành Phố Hồ Chí Minh 1.3.2 Phân bố thời gian tiến hành đề tài Tìm hiểu đề tài nghiên cứu lí thuyết : từ ngày 28/03 đến 10/04 Thiết kế chế tạo thiết bị đếm cá giống : 10/04 đến 20/04 Tiến hành chụp hình, xử lý ảnh : 20/04 đến 15/05 Viết báo cáo chỉnh sửa: 15/05 đến 31/05 1.4 Đối tượng thiết bị nghiên cứu 1.4.1 Đối tượng nghiên cứu - Nghiên cứu chế tạo thiết bị đếm cá giống phương pháp xử lý ảnh - Phần mềm Matlab 7.8 (R2009a) 1.4.2 Thiết bị nghiên cứu Web camera Colovis Máy tính Nguồn sáng Phần mềm Matlab 7.8 1.5 Phương pháp nghiên cứu Vận dụng phương pháp xử lý ảnh phần mềm Matlab, công cụ xử lý ảnh phần mềm Matlab có khả thu nhận ảnh từ camera, xử lý ảnh, đếm số lượng, cho kết 3.2.2 Cấu trúc lệnh xử lý sử dụng 3.2.2.1 Các lệnh thu nhận ảnh Trước viết lệnh thu nhận ảnh ta tiến hành kiểm tra thiết bị kết nối với máy tính hay chưa Lệnh imaqwinfo cung cấp thông tin loại adaptor thiết bị cài đặt máy tính Cấu trúc: a = imaqhwinfo  + Sau biết tên adaptor ta tiến hành kiểm tra xem có thiết bị (camera) kết nối, số deviceID thiết bị cách dung lệnh imaqhwinfo(‘adaptor name’) Cấu trúc: a = imaqhwinfo(‘winvideo’)  % winvideo tên adaptor + Lệnh imaqwinfo (‘tên adaptor’,deviceID) cung cấp thơng tin định dạng (format) tín hiệu vào Cấu trúc: a = imaqhwinfo(‘tên adaptor’, deviceID) Hình 3.4 mơ tả lệnh Hình 3.4 thơng tin thiết bị camera Sau kết nối camera vào Matlab, kiểm tra thông tin camera Ta bắt đầu công đoạn chụp ảnh thu nhận vào Matlab, xem hình chụp, lệnh sau: 23 Trước tiên ta tạo luồng liệu lệnh videoinput Cấu trúc: Vid = videoinput(‘winvideo’,1,’rgb24_320x240’);  Preview(vid) Trong đó: Tên thiết bị Matlab quy định winvideo DeviceID: số thiết bị kết nối Định dạng ảnh thu nhận lại ta mặc định rgb24_320x240 Preview(obj): xem trực tiếp ảnh động trình lấy ảnh Hình 4.5 Kết lấy ảnh sau sử dụng lệnh videoinput Hình 3.5 Ảnh động hiển thị Matlab Ta dùng lệnh getsnapshot chụp hình sử dụng hàm imshow để xem hình chụp Cấu trúc: hinh = getsnapshot(vid); Imshow(hinh) 24 Hình 3.6 ảnh sau chụp Chụp hình xong ta lưu hình lại lệnh imwrite Cấu trúc: imwrite(hinh,'ca.jpg'); % ‘ca.jpg’ filename định dạnh hình lưu 3.2.2.2 Các lệnh xử lý ảnh a Chuyển ảnh RGB sang ảnh xám Trong trình xử lí ảnh, ảnh RGB ảnh gồm ma trận chiều R, G, B gồm lớp red, green blue cấu thành, tùy theo tỷ lệ màu R, G, B mà chúng tạo nên ảnh màu RGB Vì tính tốn trực tiếp ảnh RGB phức tạp nên cần phải chuyển sang ảnh xám để dễ dàng tính tốn Hình 4.6 kết chuyển ảnh RGB sang ảnh xám Cấu trúc: Xam = rgb2gray (hinh) % dòng lệnh chuyển ảnh RGB sang ảnh xám imshow(Xam) % hiển thị ảnh xám 25 Hình 3.7 Chuyển ảnh RGB sang ảnh xám b Lấy cạnh (edge) Biên ảnh đường biên giới hay đường bao xung quanh ảnh, mà thể rõ tính chất vật lý ảnh, đặc biệt thay đổi cường độ ảnh Có nhiều phương pháp phát biên như: phương pháp Gradient, phương pháp Laplace, phương pháp zero-cross, v.v… Tuy nhiên, phạm vi đề tài ta chọn phương pháp Gradient để sử dụng (sử dụng mặt nạ Sobel) Cấu trúc : bw = edge(Xam,'sobel',0.025); Trong đó: - ‘sobel’: Là phương pháp lấy biên gradient theo kiểu Sobel - 0.025: giá trị ngưỡng gradient(do thực nghiệm) Hình 3.8 lấy cạnh từ ảnh xám 26 c Xóa pixel nhiễu ảnh xn = bwareaopen(bw,3); % xóa pixel có giá trị nhỏ ( giá trị thực nghiệm) Hình 3.9 xóa pixel nhiễu sau lấy cạnh d làm giãn ảnh (Dilate) t = bwmorph(xn,'dilate'); Hình 3.10 Làm giãn ảnh 27 e làm co ảnh (erode) se = strel('disk',1); d = imerode(t,se); với strel phần tử cấu trúc dạng đĩa trước co ảnh sau co ảnh Hình 3.11 kết sau co ảnh f Đếm số lượng cá dem = bwconncomp(t); Hình 3.12 kết sau dung lệnh đếm số lượng 28 3.2.3 Thiết kế giao diện GUI Hình 3.13 Giao diện làm việc cho người sử dụng 29 CHƯƠNG IV KẾT LUẬN VÀ ĐỀ NGHỊ 4.1 Kết đạt Trong trình nghiên cứu tìm hiểu thực hiện, đề tài đạt số kết sau: ● Chế tạo thiết bị đếm cá giống cho kết nhanh xác so với đếm thủ công truyền thống người đếm ● Đối với nguồn sáng: Để tăng chất lượng hình ảnh cho việc xử lý dễ dàng ta bố trí thêm nguồn sáng nằm sau thiết bị nhằm tăng độ tương phản cá với ảnh ● Đối với hình ảnh xử lý: Ảnh xử lý ảnh 2D ảnh chụp ta chủ động từ màu ánh sáng nên ảnh chụp để xử lý bị nhiễu có độ xác cao 4.2 Hướng phát triển đề tài ● Chất lượng hình ảnh phụ thuộc nhiều vào độ phân giải camera Trong phạm vi đề tài sử dụng camera có độ phân giải 640x480, chọn camera có độ phân giải cao cho độ xác cao ● Để hình ảnh cho kết nhanh chóng xác ngồi việc chọn mua webcam phù hợp việc bố trí nguồn sáng quan trọng (tránh bị trường hợp phản xạ ánh sáng) ● Đề tài nghiên cứu cho vấn đề đếm cá giống, ta phát triển thêm cho việc đếm tơm loại hải sản có kích thước nhỏ ● Số lượng cá mà thiết bị đếm sau lần không 100 độ phân giải camera thấp kích thước thiết bị hạn chế Nếu thiết bị đếm có kích thước lớn camera có độ phân giải cao ta đếm nhiều cá với độ xác cao 30 TÀI LIỆU THAM KHẢO I Sách tham khảo Lim S J., Xử lý ảnh (Nguyễn Văn Ngọ dịch, Nguyễn Viết Kính hiệu đính), Đại học Quốc Gia (Khoa Công Nghệ), Hà Nội Nguyễn Quang Hoan, 2006, Xử lý ảnh , Học Viện Cơng Nghệ Bưu Chính Viễn Thơng (Chương trình đào tạo đại học từ xa), Hà Nội Lương Mạnh Bá, Nguyễn Thanh Thủy, 1999, nhập môn xử lý ảnh số, nhà xuất Khoa Học Kỹ Thuật, Hà Nội Nguyễn Đức Thành, MATLAB ứng dụng điều khiển, Nxb Đại học Quốc Gia Tp Hồ Chí Minh Matlab help phần mềm Matlab 7.8.0 (R2009a) Đỗ Năng Tồn, Phạm Việt Bình, 2007, Giáo Trình Xử lý ảnh, Đại học Thái Nguyên, Khoa công nghệ thông tin II Internet http://tailieu.vn/xem-tai-lieu/matlab-cho-x u-ly-anh.260003.html http://www.ngohaibac.net/category/programming-my_work/matlab/ http://www.ebook.edu.vn/?page=1.39&view=8395 http://www.deeforum.net/forum/index.php?topic=1628.0 31 Phụ lục Chương trình xử lý ảnh Trên Matlab Chụp lưu ảnh vid = videoinput('winvideo',1,'rgb24_320x240') preview(vid) h = getsnapshot(vid); imwrite(h,'ca.jpg'); imshow (h) Đọc ảnh vào matlab xử lý i = imread('D:\de tai\hinh_vy\canho2\ca.jpg'); i1 = rgb2gray(i); bw = edge(i1,'sobel',0.025); xn = bwareaopen(bw,3); imshow(xn) t = bwmorph(xn,'dilate'); se = strel('disk',1); d = imerode(t,se); imshow(d) dem = bwconncomp(d, 4) 32 Phụ lục Bảng khảo nghiệm số lượng cá đếm Số lần Kết đếm đếm thiết bị Thực tế Ghi 53 53 Chính xác 68 68 Chính xác 53 55 Do dính 71 71 Chính xác 70 70 Chính xác 74 78 Do dính 61 64 Do dính 69 71 Camera không ổn định 71 75 Do dính 10 85 89 Do dính Trường hợp 1: Do hai dính 33 Trường hợp 2: Do camera khơng ổn định - Gọi A độ xác qua lần đo, ta có: A = (kqđ x 100)/tt Với: kqđ: Là kết đo máy tt: Là số lượng cá thực tế - Độ xác thiết bị qua 10 lần đo là: Atb = (A1 + A2 + … + A10) / 10 = 97,38 % Với: Atb độ xác thiết bị A1, A2,…,A10 độ xác qua lần đo 34 Phụ lục Matlab help Phần tử cấu trúc - Phần tử cấu trúc tạo Matlab dòng lệnh sau SE = strel ( shape, parameters ) - Phần tử cấu trúc SE loại hình dạng (shape) định rõ ( hình vng, hình thoi, hình elip,…) Tùy thuộc vào hình dạng phần tử tạo mà có thơng số thêm vào tương ứng Bảng 1: Các phần tử cấu trúc dạng mặt phẳng 'arbitrary' 'pair' 'diamond' 'periodicline' 'disk' 'rectangle' 'line' 'square' 'octagon' Cách tạo phần tử cấu trúc: ● SE = strel('diamond', R) % R bán kính phần tử hình thoi ● SE = strel('disk', R, N) % R bán kính phần tử hình đĩa, N có giá trị 0, 4, 6, Nếu N khơng chọn giá trị mặc định 35 ● SE = strel('line', LEN, DEG) % LEN DEG chiều dài góc tương ứng ● SE = strel('octagon', R) ● SE = strel('rectangle', MN) ● SE = strel('square', W) % R bán kính phần tử đa giác % MN vector phần tử % W chiều rộng phần tử hình vng 36 Phụ lục Một số định dạng ảnh mà Matlab hỗ trợ ĐỊNH DẠNG ẢNH TÊN ĐẦY ĐỦ (FORMAT) (FULL NAME) ‘bmp’ Windows Bitmap (BMP) ‘cur’ Windows Cursor Resources (CUR) ‘gif’ GraphicsInterchange Format (GIF) ‘hdf’ Hierarchical Data Format (HDF4) ‘ico’ Windows Icon resources (ICO) ‘jpg’ ‘jpef’ Joint Photographic Experts Group (JPEG) ‘pcx’ Windows Paintbrush (PCX) ‘pgm’ Portable Graymap (PGM) ‘png’ Portable Network Graphics (PNG) ‘pnm’ Portable Anymap (PNM) ‘ppm’ Portable Pixmap (PPM) ‘ras’ ‘tif’ ‘tiff’ ‘xwd’ Sun Raster (RAS) Tagged Image File Format (TIFF) X Windows Dump (XWD) 37 ... thiện Em xin chân thành cám ơn iii MỤC LỤC Trang Trang tựa i LỜI CẢM ƠN ii TÓM TẮT iii MỤC LỤC .iv DANH SÁCH CÁC HÌNH vi CHƯƠNG... biến – xuất thủy sản – phụ trợ tỉnh Đồng sơng Cửu Long Tính chun nghiệp động doanh nghiệp phân đoạn nâng lên Nhiều doanh nghiệp làm ăn bản, sử dụng quy trình cơng nghệ sản xuất đại, quy mơ lớn,... Imshow(‘filename.fmt’) imshow(I) % I tên biến gán cho ma trận ảnh Ví dụ: i = imread('D:de taihinh _vy canho2c1.jpg'); >> imshow (i) Hình 2.11 minh họa lệnh imshow 15 b.Thu nhận ảnh động vào Matlab