-
Notifications
You must be signed in to change notification settings - Fork 0
/
compactness.R
52 lines (42 loc) · 1.59 KB
/
compactness.R
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
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
library(geosphere)
library(rgeos)
#Compactness
PolsbyPopperFunction<-function(spatialObject){
perimeterValue<-perimeter(spatialObject); #in meters
areaValue<-areaPolygon(spatialObject); #in square meters
#bboxES<-spatialObject@bbox
#distanceVertical<-distm(c(bboxES[1,2],bboxES[2,1]),c(bboxES[1,2],bboxES[2,2]),fun = distHaversine)
#distanceHorizontal<-distm(c(bboxES[1,1],bboxES[2,1]),c(bboxES[1,2],bboxES[2,1]),fun = distHaversine)
compactnessValue<-4*pi*areaValue/(perimeterValue*perimeterValue)
return(compactnessValue)
}
AreaOverConvexHullFunction<-function(spatialObject){
areaValue<-areaPolygon(spatialObject); #in square meters
ch<-gConvexHull(spatialObject)
compactnessValue<-areaValue/area(ch)
return(compactnessValue)
}
#if typeOfCompact==1 PolsbyPopperFunction ----if typeOfCompact==2 AreaOverConvexHullFunction
TotalCompactnessofaResult<-function(ASolution,typeCom){
WholeSpatial<-readOGR(dsn="data5",layer="ASIAASIA")
totalCompactScore=0;
if(typeCom==1){
for(i in 1:length(ASolution)){
for(j in 1:2){
mergedShape<-aggregate(rbind(WholeSpatial[ASolution[[i]],]))
totalCompactScore<-PolsbyPopperFunction(mergedShape)+totalCompactScore;
}
}
}else{
for(i in 1:length(ASolution)){
for(j in 1:2){
mergedShape<-aggregate(rbind(WholeSpatial[ASolution[[i]],]))
totalCompactScore<-AreaOverConvexHullFunction(mergedShape)+totalCompactScore;
}
}
}
return(totalCompactScore)
}
ch = gConvexHull(Europe)
plot(Europe,col='blue',border="blue");
plot(ch,add=TRUE)