library(sp)
r1 = cbind(c(180114, 180553, 181127, 181477, 181294, 181007, 180409,
180162, 180114), c(332349, 332057, 332342, 333250, 333558, 333676,
332618, 332413, 332349))
r2 = cbind(c(180042, 180545, 180553, 180314, 179955, 179142, 179437,
179524, 179979, 180042), c(332373, 332026, 331426, 330889, 330683,
331133, 331623, 332152, 332357, 332373))
r3 = cbind(c(179110, 179907, 180433, 180712, 180752, 180329, 179875,
179668, 179572, 179269, 178879, 178600, 178544, 179046, 179110),
c(331086, 330620, 330494, 330265, 330075, 330233, 330336, 330004,
329783, 329665, 329720, 329933, 330478, 331062, 331086))
r4 = cbind(c(180304, 180403,179632,179420,180304),
c(332791, 333204, 333635, 333058, 332791))
sr1=Polygons(list(Polygon(r1)),"r1")
sr2=Polygons(list(Polygon(r2)),"r2")
sr3=Polygons(list(Polygon(r3)),"r3")
sr4=Polygons(list(Polygon(r4)),"r4")
sr=SpatialPolygons(list(sr1,sr2,sr3,sr4))
srdf=SpatialPolygonsDataFrame(sr, data.frame(cbind(1:4,5:2),
row.names=c("r1","r2","r3","r4")))
srdf
data(meuse)
coordinates(meuse) = ~x+y
plot(meuse)
polygon(r1)
polygon(r2)
polygon(r3)
polygon(r4)
meuse[,1:4]
meuse
sr
over(sr, meuse[,1:4])
r1
plot(meuse)
polygon(r1)
polygon(r2)
polygon(r3)
polygon(r4)
over(sr, geometry(meuse))
over(sr, geometry(meuse), returnList = TRUE)
library(GeoFIS)
data(conductivity_2014)
zoning <- NewZoning(conductivity_2014)
zonable_data <- zoning$zonable_data()
data(conductivity_border)
zoning$border <- conductivity_border
zoning$perform_zoning()
map5 <- zoning$map(5)
map5
over(sr, meuse, returnList = TRUE)
over(map5, zonable_data, returnList = TRUE)
over(map5, geometry(zonable_data), returnList = TRUE)
len(zonable_data)
zonable_data.len
length(zonable_data)
library(RColorBrewer)
palette <- brewer.pal(3, "Blues")
breaks <- c(30, 50, 70, 100)
colors <- palette[findInterval(map5$conduct_mean, vec = breaks)]
par(mar = c(0, 0, 2, 0))
plot(map5, col = colors, main = "Map 5 zones")
legend("bottomright", legend = levels(cut(map5$conduct_mean, breaks)),
fill = palette, title = "Conductivity mean")
o=over(map5, geometry(zonable_data), returnList = TRUE)
o
class(o)
o[[0]]
o[0]
o[1]
length(o[1])
length(o[[1]])
length(o[[1]]) + length(o[[2]]) + length(o[[3]]) + length(o[[4]]) + length(o[[5]])
o=over(map5, zonable_data, returnList = TRUE)
o
CRS
CRS("EPSG:40400")
CRS("EPSG:2154")
CRS("EPSG:40400")
over(map5, zonable_data, returnList = TRUE)
over(zonable_data, map5, returnList = TRUE)
over(zonable_data, map5)
