位相 unwrapping

位相の unwrapping について調べていました。

SAR ではなく、単なる信号の再現（逆フーリエ）でこの話題が出てきました。

調べてみると、ありますね。

位相アンラッピング(Phase unwrapping)
http://retrofocus28.blogspot.com/2013/12/phase-unwrapping_26.html

def mywrap(y):
    for i in range(1, len(y)):
        div = y[i] - y[i-1]
        if div < -pi:
            y[i:] += 2*pi
        elif div > pi:
            y[i:] -= 2*pi
        else:
            pass
    
    return y

numpy なら1行です。

np.unwrap() 

簡単ですね。
理想的な条件なら再現できるようですが、実際のデータでは難しいようです。