I have 3-D geospatial data of Geopotential anomalies that has 8658 time points, and 31x31 spatial points ( 8658time x 31 Lat x 31 Lon). I need to filter the data to retain only high frequency anomalies, using LANCZOS FILTER. The errors on execution point out that Lanczos filtering is done on 2-D time series (any way past this issue?). Now I have hustled through a variety of permutations; closest one being converting that data into (31x31 = 961) 2-D arrays of 8658x1 elements (i.e. separate time-series for each spatial point), applied the filter and then concatenated the resulting arrays into the original 8658x31x31. I am totally unsure about the authenticity of this. In a nutshell I need to apply high pass filter on 3-D data of dimensions [time x latitude x longitude]. Any thoughts?
