ขจัดเสียงรบกวนจากการถ่ายภาพรังสีทันตกรรม


9

ฉันกำลังทำงานในโครงการของการใช้รูปแบบการใช้งานรูปร่าง เพื่อหาฟันในภาพรังสีของฟัน สำหรับผู้ที่คุ้นเคยกับเทคนิคฉันกำลังพยายามที่จะสุ่มตัวอย่างตามเวกเตอร์ปกติสำหรับแต่ละจุดสังเกต กระดาษแนะนำให้ใช้อนุพันธ์ของพิกเซลตัวอย่าง: "เพื่อลดผลกระทบของการเปลี่ยนแปลงความเข้มทั่วโลกเราจะสุ่มตัวอย่างอนุพันธ์ตามโปรไฟล์แทนที่จะเป็นค่าระดับสีเทาแน่นอน"

ดังนั้นปัญหาของฉันคือวิธีกรองภาพรังสีด้วยวิธีที่ดีที่สุดในการเตรียมพวกเขาสำหรับการใช้ผู้ประกอบการอนุพันธ์ ขณะนี้ฉันกำลังใช้การผสมผสานของตัวกรองมัธยฐานเพื่อลบสิ่งที่ฉันคิดว่าส่วนใหญ่เป็นเสียงควอนตัม (mottle) มันตามมาด้วยตัวกรองทวิภาคี จากนั้นฉันก็ใช้โอเปอเรเตอร์ Scharrเพื่อคำนวณการไล่ระดับสีจริงซึ่งควรจะสุ่มตัวอย่าง

ผลลัพธ์จะแสดงด้านล่าง: ผล

ภาพแรกแสดงข้อมูลต้นฉบับ ในภาพที่สองและสามข้อมูลที่ถูกกรองจะถูกนำเสนอก่อนเป็นขนาดของสเปกตรัมหลังจาก FFT และจากนั้นเป็นข้อมูลภาพที่ถูกกรอง ภาพที่สี่แสดงผลลัพธ์ของการใช้ตัวดำเนินการ Scharr กับภาพที่สาม

คำถามของฉันคือ:

  • มีวิธีการที่รู้จักกันดีในการลดเสียงรบกวนในภาพรังสีด้วยวิธีทันตกรรมที่จะแตกต่างจากวิธีการของฉันหรือไม่?
  • อะไรทำให้เกิดลักษณะ "ควัน" ของขอบและพื้นที่ "แบน" (ไม่ใช่ขอบ)? มันเป็นสัญญาณรบกวนที่เหลืออยู่ในภาพที่ถูกกรองหรือมันมีอยู่ในตัวดำเนินการไล่ระดับสี ถ้าเป็นเสียงตัวกรองใดที่เหมาะสมที่สุดที่จะใช้ ค่ามัธยฐานตัวกรองดีในการลบ blobs ที่มีเสียงดังขนาดเล็ก แต่เคอร์เนลขนาดใหญ่ทำให้ขอบเบลอมากเกินไป ดังนั้นตัวกรองทวิภาคีจะใช้ในการกรอง blobs ขนาดใหญ่ขึ้นและทำให้สีเท่ากันทั่วบริเวณโดยไม่ทำอันตรายต่อขอบ แต่ไม่สามารถกรองโครงสร้างควันนี้ได้
  • มีตัวเลือกที่ดีกว่าตัวดำเนินการ Scharr เพื่อสร้างการไล่ระดับสีในกรณีนี้หรือไม่?
  • โบนัส: สิ่งนี้จะได้รับการพิจารณาว่าเป็นอินพุตที่ดีสำหรับ Active Shape Model หรือไม่? ฉันยังไม่ทราบว่ามันแข็งแกร่งแค่ไหน

1
คุณสามารถลองใช้วิธีการกรองแบบเลื่อนได้ เกี่ยวกับภูมิภาคที่มีควันไม่มากที่คุณสามารถทำได้ Scharr ไม่เป็นไร Canny จะดีกว่าถ้าคุณกำลังมองหาขอบที่แท้จริง
Rosa Gronchi

ฉันสามารถตอบคำถามข้อที่ 1 ก่อนอื่นคุณต้องระบุประเภทของเสียงที่มีผลต่อภาพทันตกรรม จากนั้นลองค้นหาวิธีที่มีชื่อเสียงที่สามารถกำจัดเสียงรบกวนประเภทนั้นได้
แมกซ์เวล

คำตอบ:


3

เท่าที่ฉันเข้าใจโดยการมาของภาพคุณหมายถึงการแยกขอบ ฉันขอแนะนำให้กรองภาพด้วยตัวกรอง Gaussian ที่ค่อนข้างใหญ่ หากต้นทุนการคำนวณของภาพที่ได้มานั้นไม่สำคัญกับงานของคุณฉันขอแนะนำให้ใช้ตัวตรวจจับขอบแสนรู้ มันไวต่อเสียงรบกวนน้อยกว่าและไม่หลอกเสียงรบกวนและค้นหาขอบที่อ่อนแอพร้อมกับขอบที่แข็งแรง คำสั่ง Matlab สำหรับสิ่งนั้น:

   [MinThresh MaxThresh]=[-0.3 0.5];
   EDGE_No_SMOKE=edge(im,'canny',[MinThresh MaxThresh]);

และผลลัพธ์คือ (ฉันรู้ว่ามันอาจไม่ใช่ผลลัพธ์ที่คุณค้นหาอย่างไรก็ตามการเล่นกับตัวแปร Threshold และขนาดตัวกรองจะทำให้คุณได้ผลลัพธ์ที่ต้องการ):

ป้อนคำอธิบายรูปภาพที่นี่

โปรดทราบว่าคุณไม่เห็นเอฟเฟกต์ควันอีกต่อไป นอกจากนี้เกี่ยวกับขอบผิดเหล่านั้นคุณสามารถลบออกได้โดยใช้เทคนิคการเปิดและปิดภาพ

โดยการใช้ไซต์ของเรา หมายความว่าคุณได้อ่านและทำความเข้าใจนโยบายคุกกี้และนโยบายความเป็นส่วนตัวของเราแล้ว
Licensed under cc by-sa 3.0 with attribution required.