Similar to 2dimensional wavelet decomposition, for a given irrelularlattice field represented by spatial polygons dataframe, the method decompose the field into different scales and a global trend component by EEMD method. The scale components are also called also called intrinsic mode functions (IMFs), which represent different scale information in the spatial field.
1 2 
spPolysDf 
a SpatialPolygonsDataFrame object. 
valueField 
a field name that stores value. 
nMaxIMF 
maximum number of components to be decomposed. 
tolSift 
sift tolerence, a small number. 
neemd 
number of EEMD iterations, a large number can make a stable result. 
wnsd 
standard deviation of added noise; it is a ratio to the standard deviation of above data. 
fmodel 
surface fitting function ("thinplate", "gaussian", "cubic", "multiquadric"). 
fig 
whether plot decomposed results. 
A SpatialPolygonsDataFrame with original value, decomposed imfs and global trend.
Hu, M.G. and J.F. Wang, et al. A hierarchicalscale decomposition method for irregular lattice field. Computers & Geosciences, submitted.
Huang, N. E. and Z. Shen, et al. The empirical mode decomposition and the Hilbert spectrum for nonlinear and nonstationary time series analysis. Proceedings of The Royal Society A  Mathematical, Physical & Engineering Sciences, 1998, 454(1971): 903995.
1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23  ## Not run:
library(rgdal)
# polygon data
mydata < system.file("extdata/simu.shp", package = "ibeemd")
layer < basename(mydata)
layer < substr(layer, 1, nchar(layer)4)
mydataDf < readOGR(dsn=mydata, layer=layer)
#spplot(mydataDf)
rslt < iBEEMD(
spPolysDf = mydataDf,
valueField = "value",
nMaxIMF = 10,
tolSift = 0.05,
neemd = 500,
wnsd = 0.05,
fmodel = "thinplate",
fig = TRUE)
## End(Not run)
#spplot(rslt)

Questions? Problems? Suggestions? Tweet to @rdrrHQ or email at ian@mutexlabs.com.
All documentation is copyright its authors; we didn't write any of that.