-
Notifications
You must be signed in to change notification settings - Fork 0
/
blue_whale_map.R
95 lines (81 loc) · 2.8 KB
/
blue_whale_map.R
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
92
93
94
95
pacman::p_load(
dplyr,
tidyverse,
ggplot2,
tidyr,
readr,
lubridate,
ggmap,
scales,
ggmap,
RColorBrewer
)
map_data <- read.csv("C:\\Users\\Kilia\\OneDrive\\Documents\\RStudio Projects\\blue_whale.csv")
# exploring in R; useful for reading list of colnames etc.
head(map_data)
colnames(map_data)
View(map_data)
which(
is.na(map_data$location_long),
arr.ind = TRUE
)
class(map_data$tag_local_identifier) = 'integer' #Setting tag_local_identifier as integer and string to learn more about ggmap functionality.
## MAP
?register_google #Query to set up google API
##ggmap call the map within a range specified by paramters.
mapped_whale <- get_map(
location = c(
lon = mean(map_data$location_long), #Setting the origin (x=0,y=0) as the average (x,y) coordinates of the data set
lat = mean(map_data$location_lat)
),
zoom = 4,
maptype = "satellite",
scale = 2,
extent = "device",
color = 'bw'
)
##Trying contour map
ggmap(mapped_whale)+
stat_density2d(
data = map_data,
aes(
x = location_long,
y = location_lat,
fill = ..level..
),
bins = 30,
h = 3,
n = 100,
alpha = 1,
geom = 'polygon',
size = 3
)+
scale_fill_gradientn(
colors = rev(brewer.pal(10, 'Spectral')),
labels = comma
)+
geom_point(
data = map_data,
aes(
x = location_long,
y = location_lat,
color = "Whale"
),
fill = 'red',
shape = 21,
alpha = .2
)+
guides(
fill = FALSE,
color = guide_legend(
title = "GPS Tag",
override.aes = list(size = 2)
)
)+
theme(
legend.key = element_rect(fill = 'white'),
axis.title.y = element_blank(),
axis.title.x = element_blank()
)+
ggtitle(" Azores Blue Whale Spatial Distribution, 04/2009-07/2016") +
xlab("Long") + ylab("Lat")