spl0 = readOGR(pd.pcs$lake)
sp.lake = gSimplify(spl0, bm.para$tol.wb)