วิธีการตรวจจับจุดสูงสุด

อะไรคืออัลกอริธึมการตรวจจับสูงสุดที่มีอยู่? ฉันมีข้อมูลที่มีเสียงดังและฉันต้องการที่จะใช้การตรวจจับสูงสุดสำหรับข้อมูลนี้ ข้อมูลกลับด้านจริง ๆ แล้วฉันพยายามหาจุดล่าง

นี่คือภาพรวมของข้อมูลใน Excel ฉันชอบที่จะตรวจสอบพื้นทั้งสอง ฉันคิดถึงการส่งข้อมูลผ่านตัวกรอง low-pass และจากนั้นทำค่าเฉลี่ยเคลื่อนที่ซึ่งฉันกำหนดยอดเขาและภายในค่าเฉลี่ยเคลื่อนที่ฉันทำการค้นหาอีกครั้ง ฉันมีพื้นหลังเป็นศูนย์ DSP; นี่เป็นเพียงแนวทางสามัญสำนึก ฉันต้องการได้ยินสิ่งที่ผู้เชี่ยวชาญแนะนำ

อย่างไรก็ตาม, มีหลายวิธีที่คุณสามารถใช้ที่นี่ วิธีการหนึ่งที่ผมอยากจะแนะนำคือการใช้แบบไม่ต่อเนื่องแปลงเวฟเล็ต (DWT) และโดยเฉพาะอย่างยิ่งดูที่Daubechies เวฟ ฉันจะเลือกพูด Daub-14 / Daub-Tetra

โดยทั่วไปสิ่งที่คุณต้องทำจริงๆคือ 'แนวโน้ม' สัญญาณของคุณจากนั้นก็ทำการเลือกขั้นต่ำหรือสูงสุด นั่นจะกำจัดค่าผิดปกติของคุณ การแปลงเวฟเล็ต daub-14 / daub-tetra สามารถช่วยคุณทำสิ่งนี้ได้และสิ่งนี้จะช่วยได้โดยเฉพาะอย่างยิ่งเนื่องจากคุณไม่ทราบลักษณะของสัญญาณของคุณ (โดยใช้ daub-14 คุณสามารถแสดงสัญญาณพหุนามอย่างถูกต้องในระดับ 14/2 = 7 และดูเหมือนว่าคุณไม่ต้องการมากกว่านั้น)

การคำนวณเวฟเล็ตนี้จะแปลงพลังงานของคุณให้เป็นพลังงาน สิ่งบ่งชี้เหล่านั้นแสดงถึงน้ำหนักของเวกเตอร์พื้นฐาน ส่วนที่เหลือของน้ำหนัก (นึกคิด) จะใกล้ศูนย์ เมื่อคุณมีสัญญาณรบกวนในสัญญาณของคุณ (เช่นเดียวกับคุณ) น้ำหนักเหล่านั้นที่ปกติเป็นศูนย์จะมีน้ำหนักบางส่วนในตอนนี้ แต่คุณสามารถลดสัญญาณเหล่านั้นออกและ 'denoise' สัญญาณของคุณ เมื่อเสร็จแล้วคุณสามารถทำการตรวจหา max / min อย่างง่ายได้

มีรายละเอียดเพิ่มเติมที่เกี่ยวข้องคุณอาจส่งอีเมลถึงฉันหากคุณต้องการหารือเกี่ยวกับวิธีการใช้งาน ฉันเคยทำสิ่งนี้มาก่อน

แก้ไข: นี่คือภาพบางส่วนที่แสดงให้เห็น Daub-Tetra Denoiser:

ฉันยังห่างไกลจากการเป็นผู้เชี่ยวชาญ แต่นี่คือสิ่งที่ฉันจะทำ:

คุณดูเหมือนจะมีสัญญาณที่แตกต่างกันอย่างช้าๆซ้อนทับกับความผันผวน จุดสูงสุดที่คุณค้นหานั้นมีความผันผวนมากขึ้นดังนั้นฉันจะตรวจพบโดยใช้สิ่งนั้น

1. ให้เป็นสัญญาณดิบของคุณ ใช้ค่าเฉลี่ยเคลื่อนที่มากกว่าจำนวนตัวอย่างที่เหมาะสมเพื่อสร้างผู้ให้บริการที่ราบรื่น$X(T)$$Y(t)$

2. หากคุณประสบความสำเร็จจะต้องมีค่าเฉลี่ยเป็นศูนย์ คำนวณส่วนเบี่ยงเบนมาตรฐานของXY$X-Y$$\sigma$$X-Y$

3. คุณสามารถตรวจพบยอดเขาเมื่อโดยที่เป็น 2, 3, 4$|(X-Y)(t)| > \alpha \times \sigma$$\alpha$

สามารถตรวจพบการโอเวอร์โหลดหรือการอันเดอร์โดยการลบค่าสัมบูรณ์และใช้การทดสอบที่เหมาะสม นั่นคือสิ่งที่คุณกำลังมองหา?

การตรวจจับจุดสูงสุดมีแอปพลิเคชั่นค่อนข้างน้อยสำหรับสัญญาณ 1D หรือหลายมิติ นี่คือตัวอย่างเล็ก ๆ น้อย ๆ ที่แสดงว่าสัญญาณเหล่านี้มีความหลากหลายและการตีความจุดสูงสุดสามารถ:

• ข้อมูล 1D ของโปสเตอร์ต้นฉบับ

• การแปลงภาพเป็นจุดสูงสุดแต่ละจุดสูงสุดสอดคล้องกับบรรทัดในภาพต้นฉบับ

• ความสัมพันธ์ของภาพอัตโนมัติจุดสูงสุดแต่ละจุดสอดคล้องกับความถี่ที่เผยให้เห็น "รูปแบบเป็นระยะ";

• "ความสัมพันธ์ทั่วไป" ของภาพและเทมเพลตจุดสูงสุดแต่ละจุดสอดคล้องกับการเกิดขึ้นของเทมเพลตในภาพ (เราอาจสนใจที่จะตรวจสอบเฉพาะจุดที่ดีที่สุดหรือยอดเขาหลายแห่ง)

• ผลการกรองภาพสำหรับมุมของแฮร์ริสแต่ละจุดสูงสุดตรงกับมุมในภาพต้นฉบับ

เหล่านี้เป็นคำจำกัดความและเทคนิคการตรวจจับของยอดเขาที่ฉันพบ - แน่นอนว่ามีคนอื่นที่ฉันลืมหรือไม่ทราบและหวังว่าคำตอบอื่น ๆ จะครอบคลุมพวกเขา

เทคนิคการประมวลผลล่วงหน้ารวมถึงการทำให้ราบรื่นและ denoising @ คำตอบของ Mohammad นั้นเกี่ยวกับ wavelets และคุณสามารถดูการใช้งานต่าง ๆ ของพวกมันได้ในเอกสารของ Mathematica WaveletThreshold (ซึ่งฉันก็เอาตัวอย่างมาจากทาง)

จากนั้นคุณค้นหา maxima ขึ้นอยู่กับแอปพลิเคชันของคุณคุณต้องการเพียง maxima ทั่วโลก (เช่นการลงทะเบียนภาพ), maxima ท้องถิ่นเพียงไม่กี่ (เช่นการตรวจจับเส้น) หรือ maxima ในท้องถิ่นจำนวนมาก จากนั้นลบภูมิภาครอบ ๆ จุดสูงสุดที่เลือกเป็นต้นจนกว่าค่าสูงสุดที่เหลือจะต่ำกว่าขีด จำกัด หรือคุณสามารถค้นหา maxima ในพื้นที่ภายในขนาดพื้นที่ใกล้เคียงและเก็บเฉพาะค่าสูงสุดในพื้นที่ที่มีค่าสูงกว่าเกณฑ์ - บางคนแนะนำให้รักษา maxima ในพื้นที่ตามระยะทางจนถึงส่วนที่เหลือของ maxima ท้องถิ่น (เพิ่มเติม ดีกว่า) อาร์เซนอลยังมีการดำเนินงานทางสัณฐานวิทยา: การขยายสูงสุดและการเปลี่ยนรูปหมวกสามารถเหมาะสมได้

ดูผลลัพธ์ของสามเทคนิคเหล่านี้ในภาพที่ถูกกรองสำหรับมุมของ Harris:

นอกจากนี้บางแอปพลิเคชันพยายามค้นหาจุดสูงสุดที่ความละเอียดพิกเซลย่อย การแก้ไขซึ่งสามารถใช้เฉพาะแอปพลิเคชันมีประโยชน์

เท่าที่ฉันรู้ไม่มีกระสุนเงินและข้อมูลจะบอกว่าเทคนิคใดทำงานได้ดีที่สุด

มันจะดีมากหากมีคำตอบเพิ่มเติม มาจากสาขาวิชาอื่น

ผมคิดว่าวิธีการตรวจสอบยอดหนึ่งโดยทั่วไปคือเช่นนี้ที่เป็น refpeak(bottom)

for i=1,N   {
       if i=1   {  ref＝data(i)　｝
       else ｛ if　data(i)＜＝ref　｛ref ＝data(i) ｝｝
}