경기도 입시 학원 분포 │R ggmap 활용

경기도 지도 위에 입시 학원 분포 정도를 그려보겠습니다. 지난 포스팅에 이어서 R의 ggmap 패키지를 이용합니다. 위치 정보 데이터는 경기지역경제 포탈의 소상공인 정보(2020년 10월 기준)에서 갖고 왔고, R Script는 최대한 친절하게 쓴 R로 지도에 점찍고, 선 긋고 색칠하기를 참고했습니다. 개인적으로는 몇 년 만에 ggmap을 다뤄 보는데요. 좀 더 이것저것 해봐야겠지만, API 키를 등록하는 것 말고는 큰 차이가 없는 것 같습니다. 실습 파일은 포스팅 하단에 업로드해둘 테니 참고 바랍니다. 

 

ggmap 과 API key

ggmap 패키지 호출 후에 구글맵에서 부여받은 API 키를 먼저 등록해야 합니다. 

R ggmap 활용을 위한 구글맵 API 키 발급 방법

library(ggmap)
register_google(key='부여받은 API 키 입력')

 

한국 지도 출력

한국 지도를 호출해보겠습니다. zoom의 숫자를 조정하여 적절한 크기를 찾습니다. 

map<-get_map(location = '한국', zoom=6) 
ggmap(map)

R ggmap 한국지도

 

경기도 지도 출력

경기도 지역만 호출해 보겠습니다. 경도값과 위도 값 호출 방법은 하단에 업로드한 실습 파일에 있으니, 참고 바랍니다. 

ggd <- ggmap(get_googlemap(center = c(lon=127.190292, lat=37.567167),zoom=9,maptype='roadmap',color='bw'))
ggd

R ggmap 경기도 지도 출력

 

경기도 입시 학원 위치 데이터 호출

경기지역경제포탈에서 내려받은 소상공인 정보 파일에서 입시학원만 추출하여 따로 저장하였습니다. 아래 업로드한 파일을 내려받으면 됩니다. 

a_ggd <-read.csv("C:/test/academy.csv", header=T)
head(loc_ggd)

academy.csv

1.43MB

 

 

경기도 입시 학원 위치 출력

경기도 지도에 입시 학원 위치를 표시힙니다. 수원, 성남, 부천 등에 입시 학원이 많이 소재하는 것으로 보입니다. 

ggd + geom_point(data=a_ggd, aes(x=경도, y =위도), size=0.5, alpha=0.3, colour = "steelblue")

R ggmap 경기도 입시 학원 위치 출력

 

경기도 입시 학원 밀도 출력

입시 학원 소재지를 밀도 기준 히트맵으로 구현하면, 좀 더 분명해 보입니다. 즉 입시 학원이 많이 밀집한 지역일수록 분홍색으로 표시됩니다. 

p <- ggd + stat_density_2d(data=a_ggd, aes(x=경도, y=위도, fill=.. level.., alpha=.. level..), geom='polygon', size=7, bins=30)
p + scale_fill_gradient(low='light sky blue', high='pink', guide=F) + scale_alpha(range=c(0.02, 0.8), guide=F)

R ggmap 경기도 입시 학원 밀도 출력

 

실습파일

아래 파일을 내려받은 후, 구글맵에서 부여받은 API key를 입력하면 위와 같은 결과를 구현할 수 있습니다. 

R ggmap_B.R

0.00MB

아시아 지도 출력 │R ggmap 활용 │중앙 아시아 │홍콩