Re: 오늘의 마지막 문제. 오늘 배운 내용으로 시각화 하고 싶은 주제를 자유롭게 시각화 하시오 !

[19기_이동근]

# 데이터 불러오기
ty <- read.csv("ty_all.csv", header = T)
ty

# 지도 중심점 용 태풍 데이터와 태풍 6개 태풍별 경도 위도 뽑아서 데이터 프레임 생성
tya <- ty[ty$EXP_TIME == 0, c("TY_NAME", "ELON", "NLAT")]

ty1 <- tya[tya$TY_NAME == "쁘라삐룬", c("TY_NAME", "ELON", "NLAT")]
ty2 <-tya[tya$TY_NAME == "레끼마", c("TY_NAME", "ELON", "NLAT")]
ty3 <-tya[tya$TY_NAME == "마이삭", c("TY_NAME", "ELON", "NLAT")]
ty4 <-tya[tya$TY_NAME == "찬투", c("TY_NAME", "ELON", "NLAT")]
ty5 <-tya[tya$TY_NAME == "힌남노", c("TY_NAME", "ELON", "NLAT")]
ty6 <-tya[tya$TY_NAME == "카눈", c("TY_NAME", "ELON", "NLAT")]

# 중심점 계산
center <- c(mean(tya$ELON), mean(tya$NLAT+2))

# 구글 지도 가져오기
kor <- get_map(location = center, zoom = 5, maptype = "roadmap")
ggmap(kor)

# 지도 시각화
kor.map <- ggmap(kor) + 
  geom_point(data = ty1, aes(x = ELON, y = NLAT), size = 2, alpha = 0.7, col = "red") +
  geom_point(data = ty2, aes(x = ELON, y = NLAT), size = 2, alpha = 0.7, col = "orange") +
  geom_point(data = ty3, aes(x = ELON, y = NLAT), size = 2, alpha = 0.7, col = "yellow") +
  geom_point(data = ty4, aes(x = ELON, y = NLAT), size = 2, alpha = 0.7, col = "green") +
  geom_point(data = ty5, aes(x = ELON, y = NLAT), size = 2, alpha = 0.7, col = "blue") +
  geom_point(data = ty6, aes(x = ELON, y = NLAT), size = 2, alpha = 0.7, col = "navy")

# 지도 출력
print(kor.map)