ĐẠI HỌC QUỐC GIA THÀNH PHỐ HỒ CHÍ MINH TRƯỜNG ĐẠI HỌC KHOA HỌC TỰ NHIÊN KHOA VẬT LÝ-VẬT LÝ KỸ THUẬT BỘ MƠN VẬT LÝ HẠT NHÂN KHĨA LUẬN TỐT NGHIỆP ĐẠI HỌC NGHIÊN CỨU KỸ THUẬT DỰNG ẢNH CHỤP CẮT LỚP GAMMA CHO HỆ ĐO THÙNG THẢI PHÓNG XẠ CHỨA NGUỒN ĐIỂM SVTH: HUỲNH TĂNG HIỀN CBHD: TS HOÀNG THỊ KIỀU TRANG CBPB: TS TRẦN DUY TẬP TP HỒ CHÍ MINH−2014 GIỚI THIỆU ĐỀ TÀI • Tình hình chung  Dựng ảnh từ xạ ion hóa khám phá vào cuối năm 1890 [1]  Từ ảnh kỹ thuật số đời xuất phương thức chụp ảnh cắt lớp  Tái tạo xử lý ảnh cắt lớp ngày nhận nhiều quan tâm y học công nghiệp Việt Nam GIỚI THIỆU ĐỀ TÀI • Vấn đề: Dựng ảnh chụp ảnh cắt lớp gamma cho hệ đo thùng thải phóng xạ • Hướng giải quyết: -Tìm hiểu phương pháp tái tạo xử lý ảnh cắt lớp phù hợp với hệ đo thùng thải phóng xạ -Tái tạo xử lý ảnh cắt lớp cho hệ đo thùng thải phóng xạ chứa nguồn điểm NỘI DUNG Giới thiệu hệ đo thùng thải phóng xạ Dựng ảnh cho hệ đo thùng thải phóng xạ Kết luận kiến nghị Giới thiệu hệ đo thùng thải phóng xạ • Cấu hình hệ đo • Phương thức ghi nhận hình chiếu Cấu hình hệ đo Cấu hình hệ đo (tt.) Phương thức ghi nhận hình chiếu Dựng ảnh cắt lớp cho hệ đo thùng thải phóng xạ • • • • Tái tạo ảnh với nguồn đặt tâm Hồi phục ảnh Khảo sát giới hạn phát ảnh Tái tạo ảnh với nguồn điểm đặt vị trí khác tâm Tái tạo ảnh với nguồn đặt tâm Sinogram với thời gian đo 50 giây, góc quay 3° 10 Sự nhòe ảnh chiếu ngược 36 Hàm lọc Hann • Hàm lọc Hann: 𝑊 𝜔 = 𝜔 ∙ 𝑟𝑒𝑐𝑡(𝜔/2) ∙ 0.5 + 0.5 ∙ cos(𝜔) với hàm 𝑟𝑒𝑐𝑡 𝜔/2 định nghĩa sau 𝑘ℎ𝑖 𝜔 > 1 𝑟𝑒𝑐𝑡 𝜔/2 = 𝑘ℎ𝑖 𝜔 = 𝑘ℎ𝑖 𝜔 < 37 Phương pháp Lucy-Richardson Do nhiễu nguồn phóng xạ tuân theo phân bố Poisson nên ta có cơng thức: 𝑓 (𝐻∗𝑋) 𝐻∗ 𝑋 𝑃 𝑓 𝑋 = 𝑒 𝑓! Để hội xuất 𝑓 cực đại phải thỏa phương trình: 𝑑 𝐿𝑛 𝑃 𝑓 𝑋 𝑑𝑋 ⇒ 𝑋 𝑘:1 𝑘 = 𝑋 ∙ 𝑇 =0 ⇒ 𝐻 ∗ 𝑇 𝐻 ∗ 𝑓 =1 𝐻∗ 𝑋 𝑓 𝐻∗ 𝑋𝑘 38 PSF hệ đo Từ hàm làm khớp 𝑓 ta có: 𝑓 = −0.0046 ∙ 𝑥 + 0.9955 ∙ 𝑥 − 0.1 + 0.59 ∙ 𝑒 𝑥;10.58 2∙1.4439 ; • Dịch chuyển hàm tới gốc tọa độ 𝑓 = −0.0046 ∙ 𝑥 + 10.58 + 0.9955 ∙ (𝑥 + 10.58) − 0.1 + 0.59 ∙ 𝑒 ⇒ 𝑓 = −0.0046 ∙ 𝑥 + 0.0001615 ∙ 𝑥 + 0.4275 + 0.59 ∙ 𝑒 ⇒ 𝑓(𝑥, 𝑦) = −0.0046 ∙ (𝑥 +𝑦 ) + 0.0001615 ∙ 𝑥2 + 𝑦2 ; ; 𝑥 2∙1.4439 𝑥 2∙1.4439 + 0.4275 + 0.59 ∙ 𝑒 ; 𝑥 2:𝑦 2∙1.4439 • Tạo ma trận PSF từ hàm làm khớp hai chiều 39 Số liệu thực nghiệm • Số liệu vài hình chiếu nguồn tâm thùng 10 11 12 13 14 15 16 17 18 19 20 0° 4677 4573 5044 5179 5707 5782 6321 6397 7087 7554 7824 7008 6635 6693 6734 6256 6166 5428 5145 4712 30° 4534 4688 5055 5111 5687 5819 6339 6425 7051 7660 7764 6963 6523 6746 6576 6215 5853 5615 5115 4757 60° 4520 4681 4956 5324 5753 5658 6305 6351 7052 7510 7820 6967 6556 6690 6647 6214 5974 5482 5057 4675 90° 4658 4694 4972 5194 5721 5937 6297 6365 7114 7572 7607 7126 6643 6550 6551 6219 5940 5516 5162 4681 120° 4647 4661 4836 5101 5686 5841 6537 6469 7176 7654 7535 6976 6602 6689 6645 6347 6031 5507 5132 4747 150° 4552 4688 5041 5145 5613 5741 6331 6416 7122 7713 7767 6976 6588 6733 6704 6263 6150 5482 5133 4653 180° 4576 4829 4909 5091 5600 5796 6305 6476 6911 7496 7796 7056 6419 6840 6585 6222 5955 5579 5008 4610 210° 4565 4764 5042 5148 5616 5864 6369 6453 7050 7476 7623 6981 6691 6757 6670 6364 5991 5396 5144 4572 240° 4544 4778 4993 5068 5739 5860 6508 6384 7004 7577 7886 6957 6572 6669 6606 6167 5991 5701 5084 4687 270° 4546 4682 5017 5090 5624 5732 6287 6312 7147 7641 7780 6913 6380 6650 6720 6332 6084 5582 5080 4668 300° 4676 4814 5276 5062 5666 5901 6362 6472 7029 7581 7610 7078 6469 6585 6451 6228 6077 5696 4991 4670 330° 4477 4714 4923 5207 5539 5837 6279 6583 7021 7522 7617 6948 6470 6688 6605 6317 6047 5575 5183 4638 40 Số liệu thực nghiệm (tt.) • Số liệu hình chiếu nguồn vị trí khác tâm thùng 10 11 12 13 14 15 16 17 18 19 20 0° 22488 24165 26893 30377 34962 39634 45415 49731 54035 54526 54004 50481 45380 40170 35227 31504 28171 26755 22488 24165 30° 10395 10942 12366 13495 15449 17231 19771 22934 26439 31345 34893 37409 38617 38912 38108 35839 32794 30276 10395 10942 60° 4182 4342 4691 5148 5664 6372 7282 8673 9868 11369 13796 15802 16996 19068 19992 20700 20746 20678 4182 4342 90° 1650 1929 1954 2059 2237 2396 2748 3261 3636 4587 5313 5938 6258 6833 7627 8158 8788 9276 1650 1929 120° 996 957 1111 996 1120 1220 1369 1630 1817 1977 2323 2351 2524 2685 2866 2917 3156 3702 996 957 150° 913 916 976 985 882 914 923 997 1229 1167 1048 1246 1179 1251 1400 1422 1523 1410 913 916 180° 1049 1174 946 999 1048 912 1073 945 956 1062 1241 1160 866 933 860 917 735 854 1049 1174 210° 2079 2061 2150 2215 1842 1713 1558 1351 1225 1316 1025 1162 1014 968 919 877 807 707 2079 2061 240° 6196 5721 4967 4390 4033 3485 3148 2640 2481 2566 2270 2029 1808 1643 1483 1423 1373 1210 6196 5721 270° 14587 13983 13424 12420 11079 10305 9628 8691 7291 6324 5386 4737 3953 3525 3192 2843 2634 2493 14587 13983 300° 28796 30099 30698 30333 29075 26444 23693 20730 18438 15448 13648 11834 10834 9568 9043 8261 7937 7597 28796 30099 330° 31601 35721 40075 44849 48002 50125 49829 46509 42747 37798 32817 28586 25445 23222 21062 19391 17534 16130 31601 35721 41 Sự phân bố cường độ ảnh • Sự phân bố cường độ ảnh tái tạo có lọc 10 11 12 13 14 15 16 17 18 19 20 0,27 0,27 0,27 0,27 0,27 0,27 0,18 0,12 0,14 0,14 0,16 0,17 0,2 0,26 0,31 0,35 0,28 0,27 0,27 0,27 0,27 0,27 0,27 0,27 0,19 0,09 0,06 0,05 0,06 0,06 0,07 0,08 0,1 0,14 0,21 0,26 0,27 0,27 0,27 0,27 0,27 0,27 0,22 0,12 0,05 0,03 0,02 0,02 0,02 0,02 0,03 0,04 0,05 0,07 0,26 0,17 0,08 0,05 0,02 0,01 0,15 0,15 0,08 0,04 0,02 0,01 0,01 0,02 0,04 0,05 0,07 0,08 0,09 0,09 0,08 0,07 0,06 0,05 0,08 0,15 0,27 0,12 0,05 0,03 0,02 0,02 0,04 0,06 0,09 0,12 0,14 0,15 0,16 0,16 0,14 0,12 0,09 0,06 0,06 0,11 0,04 0,03 0,03 0,05 0,08 0,12 0,15 0,19 0,22 0,25 0,25 0,25 0,22 0,19 0,15 0,08 0,03 0,04 0,05 0,08 0,13 0,18 0,24 0,29 0,33 0,35 0,36 0,35 0,32 0,28 0,22 0,16 0,1 0,05 0,03 0,01 0 0,1 0,14 0,2 0,25 0,27 0,27 0,01 0,02 0,03 0,03 0,04 0,06 0,08 0,12 0,19 0,27 0,27 0,01 0,01 0,03 0,03 0,04 0,05 0,05 0,04 0,05 0,06 0,13 0,22 0,27 0,1 0,1 0,19 0,07 0,07 0,12 0,1 0,08 0,14 10 0,09 0,03 0,05 0,07 0,13 0,19 0,26 0,32 0,39 0,44 0,47 0,48 0,47 0,43 0,38 0,31 0,23 0,15 0,11 0,13 11 0,06 0,04 0,06 0,1 0,17 0,25 0,33 0,42 0,49 0,55 0,59 12 0,08 0,04 0,07 0,13 0,21 0,3 0,4 0,5 0,6 0,59 0,55 0,49 0,4 0,31 0,22 0,16 0,17 0,5 0,39 0,28 0,19 0,17 0,59 0,66 0,7 0,72 0,71 0,66 0,59 13 0,05 0,04 0,08 0,15 0,24 0,35 0,46 0,57 0,67 0,75 0,8 0,82 0,81 0,76 0,69 0,58 0,47 0,34 0,25 0,21 14 0,15 0,05 0,09 0,16 0,26 0,38 0,5 0,63 0,73 0,82 0,88 0,91 15 0,27 0,07 0,09 0,17 0,27 0,4 0,53 0,66 0,78 0,87 0,93 0,96 0,95 0,91 0,83 0,72 0,58 0,43 0,29 0,22 16 0,27 0,19 0,09 0,16 0,27 0,4 0,53 0,67 0,79 0,89 0,96 0,9 0,85 0,77 0,66 0,53 0,4 0,28 0,2 0,99 0,94 0,85 0,74 0,59 0,44 0,27 0,18 17 0,27 0,27 0,14 0,16 0,26 0,38 0,52 0,65 0,78 0,88 0,95 0,99 0,98 0,94 0,84 0,73 0,56 0,39 0,27 0,26 18 0,27 0,27 0,24 0,15 0,24 0,36 0,48 0,62 0,73 0,84 0,91 0,95 0,94 0,88 0,78 0,63 0,46 0,31 0,27 0,27 19 0,27 0,27 0,27 0,2 0,21 0,3 0,44 0,56 0,67 0,76 0,83 0,85 0,84 0,76 0,62 0,41 0,27 0,27 0,27 0,27 20 0,27 0,27 0,27 0,26 0,16 0,25 0,35 0,45 0,56 0,6 0,68 0,66 0,68 0,5 0,28 0,27 0,27 0,27 0,27 0,27 42 Sự phân bố cường độ ảnh (tt.) • Sự phân bố cường độ ảnh hồi phục từ ảnh tái tạo có lọc 10 11 12 13 14 15 16 17 18 19 20 0 0 0 0 0 0 0 0 0,01 0 0 0 0 0 0 0 0 0,04 0,04 0 10 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 11 0 0 0 0 0,02 0,01 0,01 0,03 0 0 0 0 12 0 0 0 0 0,04 0,03 0,01 0,04 0 0 0 0 13 0 0 0 0 0 0 0 0 0 0 14 0 0 0 0 0 0,02 0,24 0,03 0 0 15 0 0 0 0 0,09 0,12 0,2 0,02 0,03 0 0 16 0 0 0 0 0,19 0,74 0,24 0 0 0 0 17 0 0 0 0 0,08 0 0 0 0 0 18 19 20 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 43 Sự phân bố cường độ ảnh (tt.) • Sự phân bố cường độ ảnh tăng độ tương phản 10 11 12 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 10 0 0 0,01 0,02 0,05 11 0 0 0,02 0,05 0,01 0,21 0,31 0,37 0,39 0,37 0,31 0,21 12 0 0 0,01 0,04 13 0 0 0,02 0,07 0,17 0,34 0,52 0,66 0,75 0,78 0,76 0,68 0,55 0,36 0,18 0,07 0,02 0,01 14 0 0 0,02 0,09 0,24 0,43 0,63 0,78 0,86 15 0 0 0,03 0,11 0,27 16 0 0 0,02 0,1 0,28 0,15 0,72 0,88 0,96 17 0 0 0,04 0,1 0,25 0,49 0,71 0,86 0,95 0,9 18 0 0 0,01 0,07 0,21 0,42 0,61 19 0 0 0,01 0,05 0,15 0,32 0,52 0,68 0,78 0,82 0,81 0,68 0,46 0,16 20 0 0 13 14 15 16 18 19 20 0 0 0 0,01 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0,01 0,02 0,04 0,06 0,01 0,01 0,02 0,01 0,01 0,01 0,02 0,02 0,02 0,01 0,01 0,01 0,01 0,03 0,05 0,07 0,07 0,07 0,05 0,03 0,01 0,1 0,1 0,14 0,19 0,2 0,18 0,14 0,09 0,05 0,02 0,23 0,36 0,49 0,58 0,61 0,59 0,5 0,4 0,02 0,07 0,18 0,33 0,4 0,5 0,1 0,37 0,23 0,04 0,01 0,1 0,21 0,01 0,88 0,82 0,69 0,49 0,28 0,11 0,03 0,01 0,71 0,85 0,93 0,96 0,95 0,8 17 0,9 0,78 0,6 0,35 0,15 0,04 0,01 0,99 0,94 0,82 0,64 0,38 0,15 0,03 0,98 0,93 0,81 0,62 0,33 0,11 0 0 0,09 0,95 0,94 0,86 0,53 0,49 0,54 0,3 0,7 0,01 0,1 0,44 0,18 0,03 0 0 0 44 Sự phân bố cường độ ảnh (tt.) • Sự phân bố cường độ ảnh hồi phục từ ảnh tăng độ tương phản 10 11 12 13 14 15 16 17 18 19 20 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 10 0 0 0 0 0 0 0 0 0 0 11 0 0 0 0 0 0 0 0 0 0 12 0 0 0 0 0 0 0 0 0 0 13 0 0 0 0 0 0 0 0 0 0 14 0 0 0 0 0 0 0 0 0 0 15 0 0 0 0 0 0,15 0,22 0,12 0 0 0 16 0 0 0 0 0 0,13 0,18 0 0 0 17 0 0 0 0 0 0,07 0 0 0 0 18 0 0 0 0 0 0 0 0 0 0 19 0 0 0 0 0 0 0 0 0 0 20 0 0 0 0 0 0 0 0 0 0 45 Biến đổi Fourier rời rạc 46 47 Ống chuẩn trực 48 49 50 ... cắt lớp cho hệ đo thùng thải phóng xạ chứa nguồn điểm NỘI DUNG Giới thiệu hệ đo thùng thải phóng xạ Dựng ảnh cho hệ đo thùng thải phóng xạ Kết luận kiến nghị Giới thiệu hệ đo thùng thải phóng xạ. .. đề: Dựng ảnh chụp ảnh cắt lớp gamma cho hệ đo thùng thải phóng xạ • Hướng giải quyết: -Tìm hiểu phương pháp tái tạo xử lý ảnh cắt lớp phù hợp với hệ đo thùng thải phóng xạ -Tái tạo xử lý ảnh cắt. .. phục ảnh (tt.) Ảnh chưa hồi phục Ảnh hồi phục 18 Dựng ảnh cắt lớp cho hệ đo thùng thải phóng xạ • • • • Tái tạo ảnh với nguồn đặt tâm Hồi phục ảnh Khảo sát giới hạn phát ảnh Tái tạo ảnh với nguồn