title: Plotting density points on a map in R
body: <p>I've been trying to plot the number of tourists from different regions in Canada on a map in R using ggplot2 and openstreet map, but I seem to be missing a step as my points all end up in the bottom right corner of the map and my map shrinks in size.</p>
  
  <p>Here are some of the data I used in the dataset map.tourists. </p>
  
  <pre><code>    id  Nb.Touristes    Nb.Nuitees
  1001    939.9513    1879.903
  1004    1273.4336   2546.867
  1006    776.5203    3882.602
  1010    3118.4872   18598.194
  1102    921.7354    3971.677
  1103    622.8770    1245.754
  </code></pre>
  
  <p>And here is the code I have until now. The data with the coordinates is in the Statistic Canada file that I download in the code below.</p>
  
  <pre><code>download.file("http://www12.statcan.gc.ca/census-recensement/2011/geo/bound-limit/files-fichiers/gcd_000b11a_e.zip", destfile="gcd_000b11a_e.zip")
  unzip("gcd_000b11a_e.zip")
  library(maptools)
  canada&lt;-readShapeSpatial("gcd_000b11a_e")
  
  library(GISTools)
  CDCenters &lt;- coordinates(canada) 
  CDCenters &lt;- SpatialPointsDataFrame(coords=canada, data=canada@data, 
                                  proj4string=CRS("+proj=longlat +ellps=clrk66"))
  CDCenters=data.frame(CDCenters, row.names=NULL , id=CDCenters$CDUID)
  
  canada_map &lt;- merge(CDCenters, map.tourists, by="id")
  
  list &lt;- ls()
  list &lt;- list[-grep("canada_map", list)] 
  rm(list=list)
  rm(list)
  
  Sys.setenv(NOAWT=1)
  library(OpenStreetMap)
  library(rgdal)
  library(stringr)
  library(ggplot2)
  
  mp &lt;- openmap(c(71, -143), c(40, -50), zoom=4, type="osm",mergeTiles=TRUE)
  library(ggplot2)
  
  autoplot(mp) +
  geom_point(data=canada_map, alpha = I(8/10), aes(x=coords.x1,y=coords.x2, size=Nb.Touristes, color=Nb.Touristes)) +
  theme(axis.line=element_blank(), axis.text.x=element_blank(), axis.text.y=element_blank(), axis.ticks=element_blank(), axis.title.x=element_blank(), axis.title.y=element_blank()) + 
  scale_size_continuous(range= c(1, 25)) +  
  scale_colour_gradient(low="blue", high="red") +
  labs(title="Nombre de touristes à Montréal en 2010 selon la division de recensement d’origine")
  </code></pre>
  
  <p>I'd post the image that I get, but I don't have enough reputation yet !</p>
  
  <p>I get two legends, the map is concentrated in the top left corner and all the points seem to be in the bottom right corner...</p>
  
  <p>What can I do?</p>
  
  <p>Thank you!!</p>