세미프로젝트04 시각화

7. 지도모형으로 데이터 시각화

필요한 패키지 불러오기(순서 중요)

install.packages("ggmap")
install.packages("ggplot2")
install.packages("raster")
install.packages("rgeos")
install.packages("maptools")
install.packages("rgdal")

library(ggmap)
library(ggplot2)
library(raster)
library(rgeos)
library(maptools)
library(rgdal)

전국 노인인구 비율 시각화 그래프
data가져오기

setwd('C:/ITWILL/2_Rwork/00/SemiProject')
P3 = total
P2 = read_excel("korea_id.xlsx")
P3 = P3[,c(1,16)]

join

P3 = P3 %>% left_join(P2)

결측치 조회,제거

sum(is.na(P3))
P3[is.na(P3)] = 0

시각화 준비

map = shapefile("C:/ITWILL/2_Rwork/00/SemiProject/TL_SCCO_SIG.shp") #지리 정보 데이터셋

[error] 삭제했던 다른 형식 파일 복원 후 로드 가능

좌표계 변환

map = spTransform(map, CRSobj = CRS('+proj=longlat +ellps=WGS84 +datum=WGS84 +no_defs'))

 

map을 df로 변환. SIG_CD칼럼이 ID로 변환

new_map = fortify(map, region = 'SIG_CD')
new_map$id = as.numeric(new_map$id)

JOIN

P_merge = merge(new_map, P3, by='id')

빈도그래프 그리기

ggplot() + geom_polygon(data = P_merge, aes(x=long, y=lat, group=group), fill = 'white', color='black')
ggplot() + geom_polygon(data = P_merge, aes(x=long, y=lat, group=group, fill = pop_old_be))

plot = ggplot() + geom_polygon(data = P_merge, aes(x=long, y=lat, group=group, fill = pop_old_be))

plot + scale_fill_gradient(low = "#ffffe5", high = "#ffb825", space = "Lab", guide = "colourbar", labs(fill="인구 비율")) + theme_bw() + labs(title = "노인인구 비율") + theme(panel.grid.major.x = element_blank(), panel.grid.minor.x = element_blank(), panel.grid.major.y = element_blank(), panel.grid.minor.y = element_blank(), plot.title = element_text(face = "bold", size = 18, hjust = 0.5))

전국 비노인인구 비율 시각화 그래프

setwd('C:/ITWILL/2_Rwork/00/SemiProject')
P3 = total
P2 = read_excel("korea_id.xlsx")
P3 = P3[,c(1,15)]

join

P3 = P3 %>% left_join(P2)

결측치 조회,제거

sum(is.na(P3))
P3[is.na(P3)] = 0

JOIN

P_merge = merge(new_map, P3, by='id')

빈도그래프 그리기

ggplot() + geom_polygon(data = P_merge, aes(x=long, y=lat, group=group), fill = 'white', color='black')
ggplot() + geom_polygon(data = P_merge, aes(x=long, y=lat, group=group, fill = pop_be))

plot = ggplot() + geom_polygon(data = P_merge, aes(x=long, y=lat, group=group, fill = pop_be))

plot + scale_fill_gradient(low = "#ffffe5", high = "#ffb825", space = "Lab", guide = "colourbar", labs(fill="인구 비율")) + theme_bw() + labs(title = "비노인인구 비율") + theme(panel.grid.major.x = element_blank(), panel.grid.minor.x = element_blank(), panel.grid.major.y = element_blank(), panel.grid.minor.y = element_blank(), plot.title = element_text(face = "bold", size = 18, hjust = 0.5))