ヒストグラム（1D配列）のロケール最小値を探す（Python）

Question

私はレーダー画像を処理し、水分を検出するヒストグラムで極小値を見つけなければなりません。ヒストグラムは各エリアごとに少し異なりますので、ヒストグラムごとにローカルミニマムを自動的に見つける必要があります。

マイ入力アレイは、画像値（0.82154、0.012211、...）の1Dアレイです。私はnumpyとmatplotlibでヒストグラムを作成する方法を知っていますが、画像に表示されているロケールの最小値を決定するために何をすべきか分かりません。私はpython scipyライブラリを使用します。

ヒストグラムを滑らかにして、最小値を簡単に判別できるようにする必要があります。データを円滑にするために何を使用するか教えてください。このような何か：

Answer 1

あなたはnumpy.convolve()でnumpyのでデータをスムーズにすることができますが、次の機能を使用することができます。

：

import numpy 

def smooth(x,window_len=11,window='hanning'): 
    if x.ndim != 1: 
     raise ValueError, "smooth only accepts 1 dimension arrays." 

    if x.size < window_len: 
     raise ValueError, "Input vector needs to be bigger than window size." 


    if window_len<3: 
     return x 


    if not window in ['flat', 'hanning', 'hamming', 'bartlett', 'blackman']: 
     raise ValueError, "Window is on of 'flat', 'hanning', 'hamming', 'bartlett', 'blackman'" 


    s=numpy.r_[x[window_len-1:0:-1],x,x[-2:-window_len-1:-1]] 
    #print(len(s)) 
    if window == 'flat': #moving average 
     w=numpy.ones(window_len,'d') 
    else: 
     w=eval('numpy.'+window+'(window_len)') 

    y=numpy.convolve(w/w.sum(),s,mode='valid') 
    return y 

はまた、scipyのダウンロードドキュメントをご覧ください

• http://scipy-cookbook.readthedocs.io/items/SignalSmooth.html

あなたは彼らの隣人よりも1次元配列a小さいのすべてのエントリを探しているなら、あなたはscipyのダウンロードで

numpy.r_[True, a[1:] < a[:-1]] & numpy.r_[a[:-1] < a[1:], True] 

を試すことができますが> = 0.11あなたは以下のように使用することができます。

import numpy as np 
from scipy.signal import argrelextrema 

x = np.random.random(12) 

# for local minima 
argrelextrema(x, np.less)