Luận văn thạc sĩ ứng dụng xử lý ảnh đếm các đối tượng có ảnh chạm nhau

Trang phụ bìa ĐẠI HỌC THÁI NGUYÊN TRƯỜNG ĐẠI HỌC CÔNG NGHỆ THÔNG TIN VÀ TRUYỀN THÔNG NGUYỄN THỊ THANH XUÂN ỨNG DỤNG XỬ LÝ ẢNH ĐẾM CÁC ĐỐI TƯỢNG CÓ ẢNH CHẠM NHAU LUẬN VĂN THẠC SĨ KHOA HỌC MÁY TÍNH Thái Nguyên năm 2020 i LỜI CAM ĐOAN Tên là: Nguyễn Thị Thanh Xuân Sinh ngày: 10/10/1977 Học viên lớp cao học CHK17A - Trường Đại học Công nghệ Thông tin Truyền thông - Đại học Thái Nguyên Hiện công tác tại: Trường THPT Đồng Hỷ - tỉnh Thái Nguyên Xin cam đoan: Đề tài “Ứng dụng xử lý ảnh đếm đối tượng có ảnh chạm nhau” TS Phạm Đức Long hướng dẫn cơng trình nghiên cứu riêng tơi Tất tài liệu tham khảo có nguồn gốc, xuất xứ rõ ràng Tác giả xin cam đoan tất nội dung luận văn nội dung đề cương yêu cầu thầy giáo hướng dẫn Nếu sai tơi hồn tồn chịu trách nhiệm trước hội đồng khoa học trước pháp luật Thái Nguyên, ngày 01 tháng 10 năm 2020 Tác giả luận văn Nguyễn Thị Thanh Xuân ii LỜI CẢM ƠN Sau thời gian nghiên cứu làm việc nghiêm túc, động viên, giúp đỡ hướng dẫn tận tình Thầy giáo hướng dẫn TS Phạm Đức Long, luận văn với Đề tài “Ứng dụng xử lý ảnh đếm đối tượng có ảnh chạm nhau” Tơi xin bày tỏ lòng biết ơn sâu sắc đến: Thầy giáo hướng dẫn TS Phạm Đức Long tận tình dẫn, giúp đỡ tơi hồn thành luận văn Phịng đào tạo Sau đại học Trường Đại học Công nghệ Thông tin Truyền thơng giúp đỡ tơi q trình học tập thực luận văn Trường THPT Đồng Hỷ - tỉnh Thái Nguyên nơi công tác tạo điều kiện tối đa cho thực khóa học Tơi xin chân thành cảm ơn bạn bè, đồng nghiệp gia đình động viên, khích lệ, tạo điều kiện giúp đỡ suốt trình học tập, thực hồn thành luận văn Thái Nguyên, ngày 01 tháng 10 năm 2020 Tác giả luận văn Nguyễn Thị Thanh Xuân LỜI CAM ĐOAN LỜI CẢM ƠN MỤC LỤC DANH MỤC CÁC BẢNG DANH MỤC CÁC HÌNH VẼ DANH MỤC CÁC KÝ HIỆU, CÁC CHỮ VIẾT TẮT ĐẶT VẤN ĐỀ CHƯƠNG TỔNG QUAN ĐẾM ĐỐI TƯỢNG VÀ NHU CẦU ĐẾM BẰNG XỬ LÝ ẢNH 1.1 Đếm đối tượng thực tế 1.1.1 Đếm thủ công mắt 1.1.2 Đếm qua thiết bị sensor 1.1.3 Đếm qua ảnh 1.2 Đặc điểm đếm đối tượng qua ảnh 1.2.1 Ảnh số 1.2.2 Những thuận lợi áp dụng đếm xử lý ảnh .8 1.2.3 Khó khăn đếm đối tượng qua ảnh KẾT LUẬN CHƯƠNG 10 CHƯƠNG CÁC PHƯƠNG PHÁP ĐẾM BẰNG XỬ LÝ ẢNH 11 2.1 Đếm đối tượng tách rời 11 2.1.1 Thuật toán kinh điểm đếm đối tượng tách rời 11 2.1.2 Đếm sai có đối tượng chạm 12 2.2 Đếm đối tượng có ảnh dính dùng hình thái học kinh điển 14 2.2.1 Hình thái học kinh điển 14 2.2.2 Hạn chế đếm đối tượng dính 16 2.3 Phối hợp Biến đổi Watershed biến đổi khoảng cách 17 2.3.1 Biến đổi khoảng cách (DT - Distance Transform) 17 iv 2.3.2 Phân vùng Watershed (WS-Watershed Segmentation) .19 2.3.3 Phối hợp biến đổi Watershed biến đổi khoảng cách .27 2.4 Sử dụng mạng nơ ron 31 2.5 Đếm đối tượng có ảnh dính dùng hình thái học định hướng 35 2.5.1 Khảo sát Hình thái học gradients kinh điển 35 2.5.2 Hình thái học có định hướng 37 2.6 Xây dựng thuật toán đếm đối tượng có ảnh dính dùng hình thái học định hướng 41 KẾT LUẬN CHƯƠNG 42 CHƯƠNG ỨNG DỤNG XỬ LÝ ẢNH ĐẾM CÂY THÉP QUA ẢNH ĐẦU BÓ 43 3.1 Cơ sở chọn thuật toán nhu cầu thực tế 43 3.2 Đếm thép qua ảnh đầu bó 51 3.3 Kết thực nghiệm 52 3.4 Nhận xét đánh giá, so sánh 55 KẾT LUẬN CHƯƠNG 56 KẾT LUẬN 56 HƯỚNG PHÁT TRIỂN 57 TÀI LIỆU THAM KHẢO 58 PHỤ LỤC I viv DANH MỤC CÁC BẢNG Bảng 2.1 Một số ví dụ biến đổi khoảng cách 18 Bảng 3.1 Quy cách thép cường lực (thép có vằn) 45 Bảng 3.2 Quy cách thép tròn trơn 46 Bảng 3.3 Độ xác kết đếm qua ảnh đầu bó thép (%) 55 v vi DANH MỤC CÁC HÌNH VẼ Hình 1.1 Đếm sơ khai trừu tượng với số Hình 1.2 Đếm thủ cơng đầu bó thép dùng sơn mắt thường .2 Hình 1.3 Đếm sensor Hình 1.4 Đếm để ước tính sản lượng Hình 1.5 Đếm hồng cầu qua kính hiển vi Hình 1.6 Đếm tơ đường cao tốc Hình 1.7 Nguyên lý máy ảnh số Hình 1.8 Nguyên lý sensor CCD Trong hình vi điểm ảnh .6 Hình 1.9 Sensor CCD sensor CMOS Hình 1.10 Sensor CMOS Hình 1.11 Ảnh thiếu sáng chất lượng kém, bị không tin Hình 1.12 Ảnh hồng cầu bị chạm, chồng lên .9 Hình 1.13 Ảnh nhiễu khơng nhiễu Hình 2.1 Tạo điểm ảnh biên cho ảnh I 11 Hình 2.2 Minh họa phương pháp đếm kinh điển 12 Hình 2.3 a) ảnh nguyên b) ảnh nhị phân 13 Hình 2.4 Thực phép co phép dãn nở ảnh nhị phân 13 Hình 2.5 a) Ảnh nguyên b) Sau dùng phép co c) Sau dùng phép dãn nở 14 Hình 2.6 Ảnh nhị phân nguyên sau thực phép hình thái học 15 Hình 2.7 Hình thái học ảnh xám 16 Hình 2.8 Thực hình thái học ảnh xám 16 Hình 2.9 Dùng phép co để tách đối tượng 17 Hình 2.10 Do phép co chưa đủ số lần thực 17 Hình 2.11 Tác dụng ba kiểu biến đổi khoảng cách 19 Hình 2.12 Nguyên lý biến đổi lưu vực sông (watershed transform) nơi giá trị cường độ xác định đồi lưu vực Đối với mục đích phân vùng, lưu vực bị ngập để kết hợp vùng tương ứng 19 vii Hình 2.13 Địa hình chưa bị ngập nước 20 Hình 2.14 Bắt đầu cho ngập thung lũng 20 Hình 2.15 Vị trí Watershed 20 Hình 2.16 Xây dựng tường ngăn cách vị trí Watershed 20 Hình 2.17 Tồn địa hình bị nhấn chìm 21 Hình 2.18 Ví dụ phân vùng ảnh Watershed 21 Hình 2.19 Watershed transform lưới ô vuông, cho kết nối khác .26 Hình 2.20 Phân vùng mức nhiễu áp dụng cho lát cắt trục hình ảnh CT Ảnh cuối: lưu vực riêng lẻ (individual basins) hợp để tạo thành vùng lớn liên tiếp 27 Hình 2.21 Thuật toán phối hợp biến đổi Watershed biến đổi khoảng cách .29 Hình 2.22 Phân vùng đồng xu dính 30 Hình 2.23 Tách hạt cà phê dính 31 Hình 2.24 Trước tiên, mạng CNN-DC thực phát điểm tâm thép mạng CNN (deep CNN), sau áp dụng thuật tốn Phân cụm khoảng cách để có thông tin số lượng thép vị trí tâm thép 32 Hình 2.25 Minh họa q trình phân cụm khoảng cách: (1) Tính khoảng cách 33 Hình 2.26 Kiến trúc CNN 34 Hình 2.27 Quá trình sử dụng phương pháp dùng mạng nơ ron chập cuộn CNN cho đêm định tâm thép 34 Hình 2.28 Ví dụ hình thái học gradient với ảnh đường 36 Hình 2.29 Hình thái học gradient với tín hiệu 1D có tác dụng cải thiện tín hiệu hai vùng 36 Hình 2.30 Ví dụ hình thái học gradient với ảnh mạch 37 Hình 2.31 Quy ước hướng ảnh với thực hình thái học có định hướng 37 Hình 2.32 Hình thái học có hướng 38 Hình 2.33 Các vị trí hướng thực Hình thái học định hướng 39 Hình 2.34 Góc  thay đổi qua trục thể điểm lõm (a) lồi (b) biên.…39 viii Hình 2.35 Điểm A đổi chiều gradient theo hướng Y-Y 40 Hình 2.36 Tìm tọa độ điểm thực hình thái học định hướng 40 Hình 2.37 Các mặt nạ hình thái học theo hướng a) X-X, b) Y-Y c) Y-X d) X-Y 40 Hình 3.1 Thép xây dựng a) dạng b) dạng cuộn 44 Hình 3.2 Thép trịn trơn thép vằn 45 Hình 3.3 Phơi thép cán 48 Hình 3.4 Đầu thép tròn thép cường lực 48 Hình 3.5 Dây chuyền cán cấu đếm DANIEL (Ytalia) nhà máy Thái Trung – Thái Nguyên 49 Hình 3.6 Hệ thống đếm thép Camera của dây chuyền cán Thái Trung (Thái Ngun) sử dụng cơng nghệ DANIEL 49 Hình 3.7 Dây chuyền cán NasteelVina sử dụng đếm thủ công 49 Hình 3.8 a) Hệ thống đếm thép băng tải SINOM GROUP LTD Sai số 0.7 % [11] b) ảnh chương trình xử lý 50 Hình 3.9 Đưa số đếm vào chương trình mã vạch để quản lý .50 Hình 3.10 Tách đối tượng dùng mẫu 52 Hình 3.11 Mẫu ảnh dùng thực nghiệm 53 Hình 3.12 Các bước thực tách đối tượng ảnh a) ảnh nguyên b) ảnh xám c) ảnh nhị phân d), e), f) bước trung gian g) tìm biên để xác định hưởng thực hình thái học h) kết tách đối tượng 54 Hình 3.13 Kết đếm số đầu bó 55 Hình 3.14 Trường hợp đầu thép bị thụt vào 55 vi ix DANH MỤC CÁC KÝ HIỆU, CÁC CHỮ VIẾT TẮT AI- Artificial Intelligence Trí tuệ nhân tạo CNN - Convolutional Neural Networks Mạng nơ ron chập cuộn (phần mềm) CNN- Cellular Neural Netwwork Mạng nơ ron tế bào (phần cứng) Closing Toán tử đóng Closing Dilation Phép dãn nở DT - Distance Transform Biến đổi khoảng cách Deep learning Học sâu Erosion Phép co Morphological Operator Tốn tử hình thái học Morphological Gradient Hình thái học Gradient Opening Tốn tử mở Opening OM - Oriented Morphology Hình thái học định hướng Sensor CCD (Charge Coupled Device) Cảm biến CCD Sensor Complementary Metal-Oxide-Semiconductor (CMOS) Cảm biến bán dẫn ô xit kim loại Watershed Segmentation Biến đổi Watershed 52 Hình 3.10 Tách đối tượng dùng mẫu Trong phần thực nghiệm phương pháp đếm số thép qua ảnh đầu bó trình bày phần 2.6 áp dụng để thực xử lý số ảnh đầu bó thép 3.3 Kết thực nghiệm Mẫu ảnh: Với mục đích thử nghiệm để đưa vào ứng dụng thực tế, mẫu ảnh để xử lý thu thập từ thực tế nhà máy cán thép NasteelViNa cán thép Thái Trung Thái Nguyên Ảnh đầu bó thép D10-D16 trơn D10 – D16 cường lực chụp máy ảnh dân dụng thơng thường ống kính 6.2 đến 18.6mm Độ phân giải 10.1 Mpixels Ánh sáng chụp kho nhà máy ánh sáng tự nhiên khơng dùng đèn flash 53 Hình 3.11 Mẫu ảnh dùng thực nghiệm Thực chương trình Ngơn ngữ sử dụng lập trình thực thuật toán 2.6 Matlab phiên 2017  Đọc file ảnh chuyển sang ảnh xám, ảnh nhị phân a) c) b) d) 54 e) f) g) h) Hình 3.12 Các bước thực tách đối tượng ảnh a) ảnh nguyên b) ảnh xám c) ảnh nhị phân d), e), f) bước trung gian g) tìm biên để xác định hưởng thực hình thái học h) kết tách đối tượng Kết quả: Số đếm ảnh N=108 (hình 3.13) n =108 55 Hình 3.13 Kết đếm số đầu bó Bảng 3.3 Độ xác kết đếm qua ảnh đầu bó thép (%) Loại thép 10 vằn 12 vằn > 14  Trong trường hợp với ảnh thép 10 vằn 12 vằn mà đầu bó thép khơng khơng đạt kết 100% Lý giải thích hình 3.14 Với sai số cần có biện pháp hỗ trợ khí “vỗ” để tạo bó thép đầu Việc thực giống cầm bó đũa dỗ đầu xuống bó đũa phía đầu dỗ xuống đầu đũa tiếp xúc với mặt phẳng có bó đũa đầu Khi chụp ảnh đầu bó đũa phía dỗ xuống khơng có ảnh đầu đũa bị thụt vào  Với bó đầu kết đạt 100% Đầu thép bị thụt vào Hình 3.14 Trường hợp đầu thép bị thụt vào 3.4 Nhận xét đánh giá, so sánh 56 Trường hợp đặc biệt: trường hợp hai mũi tên (hình 2.32 trang 39 mục 2.5.2) Khi cần kết hợp watershed xử lý xử lý theo hình thái học định hướng lại có xu hướng kht sâu vào đối tượng mà khơng có hiệu trường hợp khác Hướng phát triển: Đã có nghiên cứu cứng hóa thuật toán xử lý ảnh Do việc nghiên cứu hồn thiện phương pháp đề xuất sau thực thuật tốn FPGA để có thiết bị cầm tay sử dụng chuyên dụng cho đếm thép bó loại ảnh đầu bó tương tự hoàn toàn khả thi KẾT LUẬN CHƯƠNG Chương thực việc tiến hành thực nghiệm để tách thép ảnh đầu bó thép để thực q trình đếm Với loại thép trịn tính đồng chúng nên ảnh đầu bó thực tách đếm tương đối dễ dàng sử dụng hình thái học kinh điển Với ảnh đầu bó thép vằn phương pháp hình thái học có định hướng hữu ích tách đối tượng ảnh 57 KẾT LUẬN Phương pháp đếm số đối tượng qua ảnh áp dụng nhiều sống nhiều phương pháp để thực việc đếm Khó khăn lớn việc đếm qua ảnh số trường hợp ảnh đối tượng dính với chồng lấn lên Luận văn khảo sát trường hợp ảnh đối tượng dính Để tách đối tượng dính phương pháp hình thái học co ảnh áp dụng Tuy nhiên nhiều trường hợp khơng thể áp dụng phương pháp thực pháp co ảnh qua lần thực diện tích đối tượng (hoặc nhóm đối tượng) bị giảm có giảm đến chưa tách Trong trường hợp thực hình thái học phương pháp hữu ích Trong thời gian thực luận văn em hệ thống lại kiến thức xử lý ảnh nghiên cứu thêm số nghiên cứu giả hình thái học gradient, watersheed Thực nghiệm thực với trường hợp cụ thể công nghiệp đếm số thép qua ảnh đầu bó thép Em xin cảm ơn thày, giáo Khoa Công nghệ thông tin trường đại học Công nghệ Thông tin Truyền thông Đại học Thái Nguyên tận tình trình truyền đạt kiến thức Cảm ơn thày giáo hướng dẫn TS Phạm Đức Long giải thích, giúp đỡ em q trình thực luận văn HƯỚNG PHÁT TRIỂN Ứng dụng kỹ thuật xử lý ảnh sống nói chung sản xuất cơng nghiệp nói chung hướng nghiên cứu ứng dụng có nhiều tiếm Những phát triển em tập trung vào:  Xây dựng thuật tốn khắc phục trường hợp hình 2.32  Xây dựng mạch phát lỗi ví dụ mạch phá lỗi mạch linh kiện dây chuyền sản xuất  Nghiên cứu cứng hóa chương trình xử lý ảnh; chương trình xử lý ảnh tách rời khỏi máy tính hoạt động ổn định hơn, đáp ứng môi trường khắc nghiệt công nghiệp, chiến tranh 58 Qua trình thực luận văn có nhiều cố gắng lực cịn có hạn chế lĩnh vực nghiên cứu yêu cầu kiến thức rộng nên luận văn em cịn có khiếm khuyết Em mong nhận bảo, giúp đỡ thày cô giáo Em xin chân thành cảm ơn Thái Nguyên Ngày 01 tháng 10 năm 2020 HỌC VIÊN NGUYỄN THỊ THANH XUÂN 59 TÀI LIỆU THAM KHẢO [1] Amruta Pandit, Jyoti Rangole, Literature Review on Object Counting using Image Processing Techniques, International Journal of Advanced Research in Electrical, Electronics and Instrumentation Engineering, Vol 3, Issue 4, April 2014 [2] Sanaullah Khan, Aamir Khan, Faisal Saleh Khattak, Arslan Naseem, An Accurate and Cost Effective Approach to Blood Cell Count, International Journal of Computer Applications (0975 – 8887) Volume 50 – No.1, July 2012 [3] Jeydeo K.Dharpure, Madhukar B.Potdar, Counting Objects using Homogeneous Connected Components, International Jourpplications , Volume 63 – No.21, February 2013 [4] Pham Duc Long, Construcing hardwave to Count Number of Steel Bar in the Steel Bunches by Image Processing, TNU Journal of Science and Technology, 62018 [5] Jos B.T.M Roerdink, Arnold Meijster, The Watershed Transform: Definitions, Algorithms and Parallelization Strategies, IOS Press, Fundamenta Informaticae 41, 187–228, 2001 [6] Amandeep Kaur, Aayushi, Image Segmentation Using Watershed Transform, International Journal of Soft Computing and Engineering (IJSCE) ISSN: 22312307, Volume-4, Issue-1, March 2014 [7] Zhun Fan, Jiewei Lu, Benzhang Qiu, Tao Jiang, Kang An, Alex Noel Josephraj, and Chuliang Wei, Automated Steel Bar Counting and Center Localization with Convolutional Neural Networks, Published in ArXiv 2019 [8] Serge Beucher, P Soille, J.F Rivest, Morphological gradients, Proc SPIE "Image Science and Technology", San Jose, California, Feb 1992 [9] Pham Duc Long, Oriented Morphology and Application to Count the Number of Steelbars in the Bunches by Image Processing, TNU Journal of Science and Technology, 6- 2020 [10] LO Chua, L Yang, Cellular Neral Network Theory, IEEE Transactions on circuits and systems, 1988 60 [11] Zuoxian Nie, Mao-Hsiung Hung, Jing Huang, A Novel Algorithm of Rebar Counting on Conveyor Belt Based on Machine Vision, Journal of Information Hiding and Multimedia Signal Processing Volume 7, Number 2, March 2016 [12] Weiyan Hou, Zhengwei Duan, Xiaodan Liu, “A Template Covering based Algorithm to Count the Bundled Steel Bars,” 2011 4th International Congress on Image and Signal Processing, 2014, pp 1813-1816 I PHỤ LỤC Một số mã nguồn chương trình Hiình thái học gradient %% Dilation, erosion, and the morphological gradient % The morphological operator dilation acts like a local maximum operator % to Erosion acts like a local minimum operator You can use them together % compute something called the _morphological gradient_ % The basic form of grayscale image dilation computes, for each image % is pixel, the maximum value of its neighboring pixels The neighborhood % defined by the _structuring element_ For example, this structuring % element: % se1 = strel([1 1]) %% % defines a neighborhood consisting of the pixel itself, together with its % left and right neighbors Whereas this structuring element: se2 = strel([1; 1; 1]) %% % defines a neighborhood consisting of the pixel itself, together with the % pixels above and below it % % Let's try dilation on everyone's favorite sample MATLAB matrix, the magic % square: II m5 = magic(5) %% d1 = imdilate(m5, se1) %% % or |d1(3,3)|, which is 20, is the maximum of |[m5(3,2), m5(3,3), m5(3,4)]|, % |[6 13 20]| % Dilating with |se2| computes the 3-pixel local maximum vertically: % d2 = imdilate(m5, se2) %% % |d2(3,3)| is the maximum of |[m5(2,3), m5(3,3), m5(4,3)]| %% Erosion % Grayscale image erosion computes the minimum of each pixel's % neighborhood e1 = imerode(m5, se1) %% % |e1(3,3)| is the minimum of |[m5(3,2), m5(3,3), m5(3,4)]| %% Morphological gradient % Dilation and erosion are often used in combination to produce a desired % image processing effect One simple combination is the _morphological III % gradient_ P Soille, in section 3.8 of the second edition of % _Morphological Image Analysis: Principles and Applications_, talks about % three kinds of basic morphological gradients: % * dilated_image - eroded_image % * original_image - eroded_image % * dilated_image - original_image % Soille calls the first one the basic morphological gradient, and you % compute it this way using MATLAB and the Image Processing Toolbox: % Read in circuit board image, crop it so it's not too big for the blog % page, and convert it to grayscale: % % rgb = imread('board.tif'); I = rgb2gray(rgb(1:256,1:256,:)); se = strel(ones(3,3)); basic_gradient = imdilate(I, se) - imerode(I, se); subplot(1,2,1), imshow(I), imcredit('Image courtesy of Alexander V Panasyuk') subplot(1,2,2), imshow(basic_gradient, []) %% % The second form is called the _half-gradient by erosion_ or _internal % gradient_ internal_gradient = I - imerode(I, se); IV subplot(1,2,2), imshow(internal_gradient, []) %% % "The internal gradient enhances internal boundaries of objects brighter % than their background and external boundaries of objects darker than % their background For binary images, the internal gradient generates a % mask of the internal boundaries of the foreground image objects." % [Soille, page 86] % The third form is called the _half-gradient by dilation_ or _external % gradient_: % external_gradient = imdilate(I, se) - I; %% Direction gradients % By using line segments as structuring elements, you can compute % directional gradients seh = strel([1 1]); sev = strel([1;1;1]); horizontal_gradient = imdilate(I,seh) - imerode(I,seh); vertical_gradient = imdilate(I,sev) - imerode(I,sev); subplot(1,2,1) imshow(horizontal_gradient, []), title('Horizontal gradient') subplot(1,2,2) V imshow(vertical_gradient, []), title('Vertical gradient') %% % Try this on your own grayscale images You can experiment with larger % structuring elements, too, such as |strel('disk',7)| % % I'll be covering ways to combine erosion, dilation, and other fundamental % morphological operators in the near future %% % _Steve Eddins_ % _Copyright 2006 The MathWorks, Inc._ /////////////////////////// %read a file image RGB = imread('anh_toan_bo1.png'); figure(1),imshow(RGB) %chuyen sang gray I = rgb2gray(RGB); figure(2),imshow(I) BW = im2bw(I, 0.46); figure(3),imshow(BW); se = strel('disk',12); closeBW = imopen(BW,se); holes = imfill(closeBW, [1 1]); closeBW(~holes) = 1; figure (6), imshow(closeBW) % Thuc hien DT BW1=bwdist(BW,'euclidean'); %BW1=bwdist(BW); VI figure(4),imshow(BW1); BW2=bwdist(BW1,'cityblock'); figure(5),imshow(BW2); BW3 = imread('cac_truong_hop.png'); BW4 = edge(BW3,'canny'); figure (7), imshow(BW4) /////////////////////////////////////////// Đếm thép đầu bó %RGB = imread('dongxu.png');% RGB = imread('dau_bo_chua_eosion.png'); %figure(1),imshow(RGB) %chuyen sang gray I = rgb2gray(RGB); figure(1),imshow(I) BW = im2bw(I, 0.46);%BW = im2bw(I, 0.47); %BW=imread('dau_bo_chua_eosion.png'); figure(2),imshow(BW); se = strel('disk',12); closeBW = imopen(BW,se); holes = imfill(closeBW, [1 1]); closeBW(~holes) = 1; figure (3), imshow(closeBW) erodedBW = imerode(closeBW,se); figure (8), imshow(closeBW) L = bwlabel(BW,8); [L,n] = bwlabel(closeBW) %title(['Number =',n]); ... nghiệm đếm ảnh sản phẩm phần mềm Matlab Do em lựa chọn đề tài ? ?Ứng dụng xử lý ảnh đếm đối tượng có ảnh chạm nhau? ?? cho đề tài luận văn tốt nghiệp 2 CHƯƠNG TỔNG QUAN ĐẾM ĐỐI TƯỢNG VÀ NHU CẦU ĐẾM... Những thuận lợi áp dụng đếm xử lý ảnh .8 1.2.3 Khó khăn đếm đối tượng qua ảnh KẾT LUẬN CHƯƠNG 10 CHƯƠNG CÁC PHƯƠNG PHÁP ĐẾM BẰNG XỬ LÝ ẢNH 11 2.1 Đếm đối tượng tách rời ... thể xét 50 pixel) 5) Đếm N1 đối tượng có S  Scb đối tượng 6) Xóa N1 đối tượng có S  Scb đối tượng 7) Lặp với nhóm đối tượng cịn lại ảnh có S  Scb: { 7.1 Tìm biên nhóm có có S  Scb 7.2 Xác định