तेजी से 2 डी औसत फ़िल्टरिंग के लिए scipy.weave.inline का उपयोग

Question

मेरे पास 2 डी माध्य फ़िल्टर (3x3 विंडो) में एक बाधा है जिसका उपयोग मैं छवियों के एक बहुत बड़े सेट पर करता हूं, और मैं इसे आजमाने और अनुकूलित करना चाहता हूं। मैंने scipy.ndimage median_filter, साथ ही PIL, scipy.signal और scikits-image का परीक्षण किया है। हालांकि, एसओ में ब्राउजिंग में मैंने सीखा है कि सी में एक तेज़ ओ (एन) औसत फिल्टर है (कॉन्सटेंट टाइम में मेडियन फ़िल्टरिंग Rolling median algorithm in C देखें), और मुझे आश्चर्य हुआ कि क्या मैं इसे scipy.weave.inline का उपयोग करके पायथन में कार्यान्वित कर सकता हूं? वैकल्पिक मार्ग पर कोई सुझाव?

Answer 1

मुझे अंतर्निहित एल्गोरिदम नहीं पता है, लेकिन scikits-image में rolling median filter है।

अन्यथा, मैं इसे Cython (सी/पायथन पिजिन भाषा) में लिखने की अनुशंसा करता हूं। Numpy arrays के साथ काम करने के लिए convolution example/tutorial जांचना सुनिश्चित करें।

Answer 2

यदि आपका अभी भी रुचि मैं कोशिश करता हूँ numpy के आकृति बदलें और मंझला:

a= some big array 
a.reshape(N,3,3) #N being specific to your array 
[numpy.median(m) for m in a] 

मैं कैसे इस तराजू अपने testet तरीकों की तुलना में पता नहीं है, लेकिन आप के लिए जकड़ना सकता है अगर आप सी के साथ अनुकूलित करना चाहते हैं सूची समझ में पाश ...

Answer 3

इस प्रयास करें: Rolling median in C - Turlach implementation

http://ideone.com/8VVEa

उपयोग:

Mediator* m = MediatorNew(9); 
for (...) 
{ 
     MediatorInsert(m, value); 
     median = MediatorMedian(m); 
} 

मेरा मानना ​​है कि इस आर algo के रूप में ही है, लेकिन क्लीनर (आश्चर्यजनक तो, वास्तव में)।

आप या तो इसे लपेट सकते हैं, या इसे बंद कर सकते हैं और नुंबा (या साइथन) का उपयोग कर सकते हैं। मुझे लगता है कि मैं सिथॉन पर नुंबा की सिफारिश करता हूं, अगर कुछ और नहीं क्योंकि सादा पुराना पायथन कोड है।

मैं scikits को यह जोड़ने का सुझाव है, अगर यह तेजी से पहले से ही scikits में एक :)