Ứng dụng xử lý ảnh đếm các đối tượng có ảnh chạm nhau

Đếm các đối tượng trong ảnh là một bài toán được đặt ra nhiều trong thực tế như đếm các tế bào, đếm hồng cầu trong máu qua ảnh từ kính hiển vi điện tử, đếm các sản phẩm công nghiệp qua các camera công nghiệp… Khó khăn lớn nhất với việc đếm bằng xử lý ảnh là ảnh của các đối tượng cần đếm thường dính với nhau gây ra kết quả sai khác so với thực tế. Cho đến nay đã có nhiều phương pháp làm tăng độ chính xác của kết quả đếm bằng xử lý ảnh. Tuy nhiên kết quả của các phương pháp đó vẫn chưa cho độ chính xác 100%. Trong bài báo này chúng tôi đề xuất một phương pháp mới áp dụng hình thái học có định hướng để tách các đối tượng bị dính nhau trong ảnh để đưa độ chính xác kết quả của các quá trình đếm có khả năng đạt đến độ chính xác tuyệt đối trong một số ứng dụng trong công nghiệp hiện nay. Thực nghiệm được thực hiện trên ảnh của các bó thép cho số lượng cây thép có trong bó có thể đạt độ chính xác tuyệt đối. Kết quả này có thể được ứng dụng thay cho việc đếm bằng thủ công trong các nhà máy cán thép hiện nay. Nhu cầu đếm các đối tượng có rất nhiều trong thực tiễn. Đếm thủ công bằng mắt là một phương pháp tự nhiên, đã có lâu đời. Nhưng ngày nay việc đếm thủ công đã không còn đáp ứng được nữa do tốc độ xuất hiện và tồn tại các đối tượng cần đếm rất nhanh, do số lượng đối tượng cần đếm vô cùng lớn, do các đối tượng trong nhóm dính hoặc chồng lấn lên nhau Mục tiêu của luận văn là nghiên cứu Xây dựng thuật toán đếm đối tượng có ảnh dính nhau dùng hình thái học định hướng. Chúng ta thấy morphological gradient chỉ có tác dụng nổi biên ở các phía trong và ngoài biên cũng như nổi biên so với nền ảnh trong một số trường hợp và không tác dụng tách rời các đối tượng dính nhau. Với mục tiêu tách đối tượng tách rời các đối tượng dính nhau chúng tôi đưa ra một loại hình thái học mới thực hiện hình thái học dựa trên đặc điểm của biên ảnh tại những vị trí có biểu hiện đặc điểm của hai đối tượng “dính nhau”. Từ đó tiến hành xây dựng mô hình thực nghiệm đếm ảnh sản phẩm trên phần mềm Matlab. Do vậy em lựa chọn đề tài “Ứng dụng xử lý ảnh đếm các đối tượng có ảnh chạm nhau” cho đề tài luận văn tốt nghiệp.

Trang phụNGUN bìa ĐẠI HỌC THÁI 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 tơi 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 đỡ tơi suốt q trình học tập, thực hoà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 iii MỤC LỤC LỜI CAM ĐOAN i LỜI CẢM ƠN ii MỤC LỤC iii DANH MỤC CÁC BẢNG iv DANH MỤC CÁC HÌNH VẼ v DANH MỤC CÁC KÝ HIỆU, CÁC CHỮ VIẾT TẮT vi ĐẶT VẤN ĐỀ Error! Bookmark not defined 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 .5 1.2.1 Ảnh số 1.2.2 Những thuận lợi áp dụng đếm xử lý ảnh 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 vvi DANH MỤC CÁC HÌNH VẼ Hình 1.1 Đếm sơ khai trừu tượng với số .2 Hình 1.2 Đếm thủ cơng đầu bó thép dùng sơn mắt thường 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 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 Hình 1.13 Ảnh nhiễu khơng nhiễu .9 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 vng, 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 tố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 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 Nguyên) 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 % b) ảnh chương trình xử lý [11] 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 Tố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 Toá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 tố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) g) f) 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 Số ảnh Bằng đầu xử lý Không đầu 10 vằn 100 100% 98% 12 vằn 100 100% 99% > 14 100 100% 100%  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 khoét 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 tố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 á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, cô 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 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 % Erosion acts like a local minimum operator You can use them together to % compute something called the _morphological gradient_ % %% Dilation % The basic form of grayscale image dilation computes, for each image % pixel, the maximum value of its neighboring pixels The neighborhood is % 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) %% % |d1(3,3)|, which is 20, is the maximum of |[m5(3,2), m5(3,3), m5(3,4)]|, or % |[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]); ... ĐẾ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... 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... nhiên đếm xử lý ảnh có khó khăn chất phương pháp phụ thuộc ánh sáng có sai số đối tượng bị dính /chạm khơng có phương pháp loại trừ thích hợp 11 CHƯƠNG CÁC PHƯƠNG PHÁP ĐẾM BẰNG XỬ LÝ ẢNH 2.1 Đếm đối