How to use R and datasets, to make GeoLocated Maps

Philly Crime

I assure I commited neither. I got the dataset here.

The Code

library("ggmap") #Get Philadelphia Crime Data, this takes a while
philCrimeData <- read.csv("~/Downloads/crime.csv",stringsAsFactors = FALSE)
# Get Philly Map, Zoom in slightly
philCrimeMap <- get_map(location = "philadelphia",maptype = "roadmap", zoom = 12)
#Create Map
philCrimeMapGG <- ggmap(philCrimeMap) #Create My Address
me <- data.frame( Lon = c(-75.112021), Lat = c(39.991058), Crime = "Me", stringsAsFactors = FALSE
) #Data Needs to be arranged so that the legend prints properly
#Create rapes table
rapes <- philCrimeData[which( philCrimeData$Text_General_Code == 'Rape' ),]
rapes <- rapes[,c("Lat","Lon")]
#Create Crime Col
rapes$Crime <- "Rapes"
#Create Murders Table
murders <- philCrimeData[grep(pattern = "Homicide",philCrimeData$Text_General_Code),]
murders <- murders[,c("Lat","Lon")]
#Create Crime Col
murders$Crime <- "Homicides"
#Bind Tables
Crimes <- rbind(rapes, murders)
#Add Dots
pngf <- philCrimeMapGG + #Make Crime dots, There's a lot So make them Transparent geom_point(show.legend = TRUE,data = Crimes, aes(x = Lon, y = Lat, color = Crime),size = 0.5, alpha = 0.2) + #Make An Outline for the Me dot, 1 unit bigger than the me dot, Don't want this to show up in the legend, Not Transparent geom_point(show.legend = NA,data = me, aes(x = Lon, y = Lat), color = "black", size = 3) + #Make the Me dot, size 2 compared to our .5 dot Crime dots, Not Transparent geom_point(show.legend = TRUE,data = me, aes(x = Lon, y = Lat, color = Crime), size = 2)
#The Legend is created Automagically, By utilizing the Crimes and Me Tables data Structure
#The colors for Murders and Rapes are automagically assigned, Except the Me dot which is specified
#Add Labels
final <- pngf + xlab(label = "Longitude") + ylab(label = "Latitude") + ggtitle("Philly",subtitle = "Criminal Activity") + theme(plot.title = element_text(hjust = 0.5), plot.subtitle = element_text( hjust = 0.5)) #Save as PNG
png(filename = 'Philly.png', width=1000, height=625)

The code is probably better highlighted on GitHub.

Adult Business Map

I’ve never been to any of these. The dataset can be found here.

Getting My Lat and Long

First I got a google maps api key, and I saved it to “~/.gmaps”. Then I wrote this Script:

#!/bin/bash - #===============================================================================
# # USAGE: ./ 1600 Pennsylvania Ave NW, Washington, DC 20006
# REQUIREMENTS: gsed for Mac, if your on linux change gsed to sed
# BUGS: ---
# NOTES: One Shot, Get JSON for Lat and Lon for an address
# AUTHOR: Bsdpunk, Dusty Carver
# ORGANIZATION: # CREATED: 12/08/2017 19:25
#=============================================================================== set -o nounset # Treat unset variables as an error
#Put Your API Key in ~/.gmaps
#Just Copy it, and Paste it into vim as the above file
KEY=$(cat ~/.gmaps)
for i in $ARGARR;
#Remove Commas
APPEND=$(printf $APPEND | gsed 's/,//g')
#Get JSON, and reduce it
curl -sKL "$APPEND&key$KEY" | egrep 'lat|lng' | head -n2 

Creating The Map

Now It’s just like the other one. Just as simple, and I tried to comment the hell out of my code.

library("ggmap") #Get Nashville Adult Data
nashAdultData <- read.csv("~/data/adult2.csv",stringsAsFactors = FALSE)
# Get Nashville Map, zoom in 2 additional levels
nashAdultMap <- get_map(location = "nashville",maptype = "roadmap", zoom = 12)
#Create Map
nashAdultMapGG <- ggmap(nashAdultMap) #Create My Address
me <- data.frame( Long = c(-86.74), Lat = c(36.22), AdultBusiness = "Me", stringsAsFactors = FALSE
#Create Adult Frame
adultGPS <- nashAdultData[,c("Lat","Long")]
#Add type
adultGPS$AdultBusiness <- "Adult Business"
#Combine tables
AdultCommerce <- rbind(me, adultGPS)
#Save Map, with additions to variables
final <- nashAdultMapGG + #Add black dots for erryone, do not show in legend geom_point(show.legend = FALSE,data = AdultCommerce, aes(x = Long, y = Lat),size = 3, color = "black") + #Create different Color dots based on type(AdultBusiness or Not), do show in legend geom_point(show.legend = TRUE,data = AdultCommerce, aes(x = Long, y = Lat, color = AdultBusiness),size = 2) + #Label X axis Longitude xlab(label = "Longitude") + #Label Y axis Longitude ylab(label = "Latitude") + #Add Title and Subtitle ggtitle("Registered Adult Businesses",subtitle = "Nashville") + #Center Title and Subtitle theme(plot.title = element_text(hjust = 0.5), plot.subtitle = element_text( hjust = 0.5)) #Save as PNG
png(filename = 'Nashville.png', width=1000, height=800)

Cause /r/data_irl is awesome.

Cause maybe someone on /r/data_irl will see this, and realize making shit in R is ridic easy. And we get more data, IRL.

Guess what you didn’t make a map, you just made a graph. That thing is just decimal numbers mapped out.

I am @bsdpunk[@bsdpunk], and I reccomend buying this R book cause machine learning is awesome.

Related posts

Leave a Reply

Be the First to Comment!

Notify of