空間分析　 W4課堂練習

助教 杜承軒 2021.03.14

library(sf);library(tmap)
setwd("D:/1092SA/Data")
FF=st_read("Tpe_Fastfood.shp", options="ENCODING=BIG5")
TP=st_read("Taipei_Vill.shp", options="ENCODING=BIG5")

DAAN=TP[TP$TOWN=="大安區",]
FF.DAAN=FF[DAAN,] #在大安區的速食店
xtabs(~STORE,FF.DAAN)

## STORE
## KFC MIC 
##   2  11

FF.TP=st_intersection(FF,TP)
xtabs(~TOWN+STORE,FF.TP)

##         STORE
## TOWN     KFC MIC
##   士林區   2   8
##   大同區   1   3
##   大安區   2  11
##   中山區   4   9
##   中正區   2   8
##   內湖區   2   5
##   文山區   1   7
##   北投區   1   5
##   松山區   2   8
##   信義區   1   8
##   南港區   0   3
##   萬華區   2   3

M=FF[FF$STORE=="MIC",]
M.1k.sep=st_buffer(M,1000)
M.1k=st_union(M.1k.sep)
tm_shape(M.1k)+tm_polygons()+tm_shape(M)+tm_dots()