LocalUnwrap {RSEIS} | R Documentation |
unwrap the phase spectrum so it does not wrap around
LocalUnwrap(p, cutoff = cutoff)
p |
phase spectrum |
cutoff |
cut off angle = pi |
Unwrapped spectrum
Algorithm minimizes the incremental phase variation by constraining it to the range [-pi,pi]
Jonathan M. Lees<jonathan.lees.edu>
x = 1:512 amp = sin(1*2*pi*x/16) + sin(2*2*pi*x/16) + sin(3*2*pi*x/16) spc = fft(amp) plot(Mod(spc), type='l') angle = Arg(spc) plot(angle, type='l') unang = LocalUnwrap(angle, cutoff =pi ) plot(unang, type='l')