The Marie-Galante centroids has been extracted from OpenStreetMap.
Map data copyrighted OpenStreetMap contributors and available from 'http://www.openstreetmap.org'
x |
longitude |
y |
latitude |
1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 | library(gwfa)
data("mariegalante")
test=gwfa(points=mariegalante,q=0,radius=(20*2^((0:6)/2)),
bandwith=1600,sample_size=500,cell_size=2000)
test=test[test$count>100,]#select the cells with at least 100 points.
#estimate the fractal dimension on the 7 radius
X=cbind(rep(1,length(test@radius)),log2(test@radius))
fit_frac_dim=(do.call(cbind,test[,4:10]))%*%t(solve(t(X)%*%X)%*%t(X))
test$dimfrac=fit_frac_dim[,2]
#create spatial polygon dataframe
shp=grid_to_spdf(test,"2970")
## Not run:
library(cartography)
choroLayer(spdf=shp,nclass=5,var="dimfrac",method="fisher-jenks")
## End(Not run)
|
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.