空間分析　第四章

第四章點事件的空間型態

讀取資料與初步整理

library(sf);library(aspace)
setwd("D:/SA_Book/Lab_Data_Ch4")
TW.town=st_read(dsn="Taiwan_town.shp")
mazhou=st_read(dsn="Taiwan_temple_mazhou.shp")

問題4-1：計算中心位置、離散程度及方向性

請問台灣本島媽祖廟宇的基本空間型態為何？請以中心位置、離散程度與方向性來描述。

# 計算與繪製平均中心
# st_coordinates可以匯出點資料的x,y座標
mazhou.xy=st_coordinates(mazhou)
mean.pt=mean_centre(points = mazhou.xy)

##   id CENTRE.x CENTRE.y
## 1  1 215074.5  2615632

par(mar=c(1,1,1,1))
plot(st_geometry(TW.town),border="grey60") #底圖
points(mazhou.xy,pch=20,cex=0.5,col='chartreuse3') #所有廟宇點位
points(mean.pt$CENTRE.x,mean.pt$CENTRE.y,pch=19,col='red',cex=2)
# 計算與繪製標準圓
calc_sdd(points = mazhou.xy)
plot_sdd(plotnew=F, plotcentre=F,plotpoints=F,sdd.col='red',titletxt="")
# 計算與繪製標準橢圓
calc_sde(points = mazhou.xy)
plot_sde(plotnew=F, plotcentre=F,plotpoints=F,sde.col='blue',titletxt="")

問題4-2：計算中心地理位置

若台灣宗教信仰的民間團體欲舉辦一場全台的媽祖繞境活動，各地信眾會從當地的媽祖廟出發。假設每間媽祖廟的信眾人數相同，且只考慮直線距離之下，若要選擇一間合適的媽祖廟作為會師終點，使各地信眾走到終點的距離總和最小，則應該選擇何座媽祖廟宇作為會師終點？

# 計算與繪製地理中心物件
CF.pt=CF(points = data.frame(mazhou.xy))

##   id     CF.x    CF.y
## 1  1 199423.9 2610921

par(mar=c(1,1,1,1))
plot(st_geometry(TW.town),border="grey60") #底圖
points(mazhou.xy,pch=20,cex=0.5,col='chartreuse3') #所有廟宇點位
points(CF.pt$CF.x,CF.pt$CF.y,pch=17,col='red',cex=2)

問題4-3：計算中位數中心點

承問題4-2所找到的會師終點，由於該媽祖廟的室內腹地不夠大，因此，主辦單位需再找一個新地點搭建活動舞台，考量的選址條件為希望各地信眾從出發到舞台的距離總和可以最小。請問，主辦單位應選擇何處？

# 計算與繪製中位數中心點
median.pt=median_centre(points = mazhou.xy)

##   id median.x median.y
## 1  1 198975.8  2616313

par(mar=c(1,1,1,1))
plot(st_geometry(TW.town),border="grey60") #底圖
points(mazhou.xy,pch=20,cex=0.5,col='chartreuse3') #所有廟宇點位
points(median.pt$median.x,median.pt$median.y,pch=15,col='red',cex=2)

第四章 點事件的空間型態

第四章點事件的空間型態