Decision-Making Information

ProximityOne   information resources & solutions   (888) DMI-SOLN   (888) 364-7656 Put data to work more effectively. Certificate in Data Analytics Visual analysis with ProximityOne tools click graphic for info; hover to pause   visually analyze clients/markets site analysis using 1 mile radius $median household income patterns by census tract - Houston 113th Congressional Districts US Asian Indian population 2010 geocoded students and school McKinney ISD, TX geocoded students and school with tax parcels & streets high school attendance zones with schools by type %high school graduates by census tract - Puerto Rico Austin, TX MSA counties & places 10K+ population markers Appalachia counties (green) & coalfield counties (orange) China provinces percent urban & cities (markers) by state plan Honolulu census tracts (red) & census blocks Central Park area NYC Kansas City Metro & Counties Home Depot locations (markers) World Cities; focus on Spain Data Analytics Blog Mapping Statistical Data Support & Technical Assistance help using these resources

#the following 4 commands only need to be run once
# ... in the event you run the remainder of the command multiple times
install.packages('leafletR')
install.packages('rgdal')
install.packages('rgeos')
install.packages('sp')

library(leafletR)
library(rgdal) #for reading/writing geo files
library(rgeos) #for simplification
library(sp)

# large file; may take a few minutes to download
url<-"http://www2.census.gov/geo/tiger/TIGER2010DP1/County_2010Census_DP1.zip"
downloaddir<-"c:/Leaflet"
destname<-"tiger.zip"
download.file(url, destname)
unzip(destname, exdir=downloaddir, junkpaths=TRUE)

######

filename<-list.files(downloaddir, pattern=".shp", full.names=FALSE)
filename<-gsub(".shp", "", filename)

# may take 1 minute to read
dat<-readOGR(downloaddir, filename)

# Create a subset of New York counties; 48 is the FIPS code for NY;
# change 48 to 20 to make Kansas map.
subdat<-dat[substring(dat$GEOID10, 1, 2) == "48",]

subdat<-spTransform(subdat, CRS("+init=epsg:4326"))

# assign pseudo names to field names
names(subdat)[names(subdat) == "DP0010001"]<-"Population"
names(subdat)[names(subdat) == "NAMELSAD10"]<-"County"

# save required items (only 2 needed from broader set)
subdat_data<-subdat@data[,c("County", "Population")]

subdat<-gSimplify(subdat,tol=0.01, topologyPreserve=TRUE)
subdat<-SpatialPolygonsDataFrame(subdat, data=subdat_data)

# Start LeafletR operations

# Create GeoJSON
leafdat<-paste(downloaddir, "/", filename, ".geojson", sep="")

writeOGR(subdat, leafdat, layer="", driver="GeoJSON")

cuts<-round(quantile(subdat$Population, probs = seq(0, 1, 0.20), na.rm = FALSE), 0)
cuts[1]<-0
# Fields to use in the map click popup
popup<-c("County", "Population")

# Gradulated legend based on an attribute (population field)
sty<-styleGrad(prop="Population", breaks=cuts, right=FALSE, style.par="col",
# style.val=rev(heat.colors(5)), leg="Population (2010)", lwd=1)
style.val=c("#eff3ff","#bdd7e7","#6baed6","#3182bd","#08519c"),
leg="Population 2010",
fill.alpha=0.7, rad=8,
alpha=1,lwd=2,col="#006400")

# Create map and load into browser
map<-leaflet(data=leafdat, dest=downloaddir, style=sty,
title="index", base.map="osm",
incl.data=TRUE, popup=popup)

# view HTML source code and optionally copy/use in your own web page
browseURL(map)

# ... in the event you run the remainder of the command multiple times
install.packages('leafletR')
install.packages('rgdal')
install.packages('rgeos')
install.packages('sp')

# Create a subset of New York counties; 48 is the FIPS code for NY;
# change 48 to 20 to make Kansas map.
subdat<-dat[substring(dat$GEOID10, 1, 2) == "48",]