setwd('Desktop/DROUGHT/website/data')

library(ggmap)

var <- read.csv('dly_state_extremes.csv', header=TRUE, sep=",")
attach(var)

cols = c('firebrick1','dodgerblue','rosybrown1','blue','green4','lightblue1')
months = c('January','February','March','April','May','June','July','August','September','October','November','December','Annual')
lab = c('High Maximum Temperature','Low Maximum Temperature','High Minimum Temperature','Low Minimum Temperature','High Precipitation','High Snowfall')

lon = array(data=NA, dim=6); lat = lon;

i=1
for (i in 1:13) {
  lon[1] = maxt.hi.lon[i]
  lon[2] = maxt.lo.lon[i]
  lon[3] = mint.hi.lon[i]
  lon[4] = mint.lo.lon[i]
  lon[5] = pcpn.hi.lon[i]
  lon[6] = snow.hi.lon[i]
  
  lat[1] = maxt.hi.lat[i]
  lat[2] = maxt.lo.lat[i]
  lat[3] = mint.hi.lat[i]
  lat[4] = mint.lo.lat[i]
  lat[5] = pcpn.hi.lat[i]
  lat[6] = snow.hi.lat[i]
  
  record.loc <- data.frame(lon,lat)  
  
  outpic = paste('recs',i,'.png',sep='')
  png(file=outpic)
  
  map <- get_map(location = c(lon=-105.5,lat=39), zoom=7)
  mapPoints <- ggmap(map, extent="device") + geom_point(data=record.loc, aes(x=lon, y=lat, size=1.5), show.legend=FALSE)
  mapPoints
  
  dev.off()
  
outfile <- 'extremes_image_map.html'
writeout <- file(outfile, "w")


scale=1.
iwidth = 475. * scale; iheight = 353. * scale;
xoffset = 5. * scale; yoffset = 64. * scale;

j=1
for (j in 1:6) {
##<area alt="" title="High Maximum Temperature" href="#" shape="circle" data-tooltip="sticky1" coords="398,323,10" />
  
  x = ((lon[j] + 109.) * (iwidth / 7.)) + xoffset
  y = ((41. - lat[j]) * (iheight / 4.)) + yoffset
  
  x = round(x); y = round(y);

  out1 = '<area title="'
  out2 = '" href="#" shape="circle" data-tooltip="sticky'
  out3 = '" coords="'
  out4 = ','
  out5 = ',10" />'
  
  sout = paste(out1,lab[j],out2,j,out3,x,out4,y,out5,sep="")
  write(sout, file=writeout, append=TRUE)
}

close(writeout)

  
}