Databases
StackOverflow2013
Postsdbo
PORelocating Alaska and Hawaii on thematic map of the USA with ggplot2
 

Note that there are some explanatory texts on larger screens.

plurals
  1. PORelocating Alaska and Hawaii on thematic map of the USA with ggplot2
    primarykey
    data
    text
    <p>I am trying to create a thematic map showing all 50 US states, but I am having trouble relocating Alaska and Hawaii in a reliable way. I have a couple of ideas but none of them work well. I will demonstrate them now.</p> <p>First we need to import the data; using the data in the <code>maps</code> package is not enough because it does not include Hawaii and Alaska.</p> <pre><code>setwd(tempdir()) download.file("https://dl.dropbox.com/s/wl0z5rpygtowqbf/states_21basic.zip?dl=1", "usmapdata.zip", method = "curl") # This is a mirror of http://www.arcgis.com/home/item.html? # id=f7f805eb65eb4ab787a0a3e1116ca7e5 unzip("usmapdata.zip") require(rgdal) all_states &lt;- readOGR("states_21basic/", "states") require(ggplot2); require(maptools); require(rgeos); require(mapproj); all_states &lt;- fortify(all_states, region = "STATE_NAME") </code></pre> <p>Now we define some plot aesthetics:</p> <pre><code>p &lt;- ggplot() + geom_polygon( aes(x=long, y=lat, group = group, fill = as.numeric(as.factor(id))), colour="white", size = 0.25 ) + coord_map(projection="azequalarea") + scale_fill_gradient(limits = c(1,50)) </code></pre> <p>Now we remove all background etc so they don't clash when we overlap the non-contiguous states:</p> <pre><code>p &lt;- p + 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(), panel.background=element_blank(), panel.border=element_blank(), panel.grid.major=element_blank(), panel.grid.minor=element_blank(), plot.background=element_blank()) </code></pre> <h2>Using viewports</h2> <p>The first idea I had was to use viewports:</p> <pre><code>AK &lt;- p %+% subset(all_states, id == "Alaska") + theme(legend.position = "none") HI &lt;- p %+% subset(all_states, id == "Hawaii") + theme(legend.position = "none") contiguous &lt;- p %+% subset(all_states, id != "Alaska" &amp; id != "Hawaii") grid.newpage() vp &lt;- viewport(width = 1, height = 1) print(contiguous, vp = vp) subvp1 &lt;- viewport(width = 0.25, height = 0.25, x = 0.18, y = 0.33) print(AK, vp = subvp1) subvp2 &lt;- viewport(width = 0.12, height = 0.12, x = 0.32, y = 0.27) print(HI, vp = subvp2) </code></pre> <p><img src="https://i.stack.imgur.com/UytXg.png" alt="enter image description here"></p> <p>This looks nice but it is not satisfactory because it is very sensitive to slight changes in the figure, for example resizing or changes in the size and shape of the legend.</p> <h2>Manually moving Alaska and Hawaii</h2> <pre><code>all_states_AKHImoved &lt;- within(all_states, { lat[id == "Alaska"] &lt;- lat[id == "Alaska"] - 45 long[id == "Alaska"] &lt;- long[id == "Alaska"] + 40 lat[id == "Hawaii"] &lt;- lat[id == "Hawaii"] + 0 long[id == "Hawaii"] &lt;- long[id == "Hawaii"] + 70 }) p %+% all_states_AKHImoved </code></pre> <p><img src="https://i.stack.imgur.com/XwvZS.png" alt="enter image description here"></p> <p>This is not satisfactory because Alaska is usually not to scale on most US maps so it looks very big. Also relocating Alaska and Hawaii changes the distortion introduced by the map projection.</p> <h2>Question</h2> <p>Does anyone have any better approaches? </p>
    singulars
    1. PO
      singulars
      1. PTAnswer
    1. This table or related slice is empty.
    1. PTQuestion
    1. USorizon
    1. USorizon
    plurals
    1. This table or related slice is empty.
    1. This table or related slice is empty.
    1. This table or related slice is empty.
    1. PO
      singulars
      1. PTAnswer
    2. PO
      singulars
      1. PTAnswer
    1. VO
      singulars
      1. PORelocating Alaska and Hawaii on thematic map of the USA with ggplot2
      1. This table or related slice is empty.
      1. VTUpMod
    2. VO
      singulars
      1. PORelocating Alaska and Hawaii on thematic map of the USA with ggplot2
      1. This table or related slice is empty.
      1. VTUpMod
    3. VO
      singulars
      1. PORelocating Alaska and Hawaii on thematic map of the USA with ggplot2
      1. This table or related slice is empty.
      1. VTUpMod
    1. COmissing require(mapproj), and ggplot(data =all_states)
      singulars
      1. PORelocating Alaska and Hawaii on thematic map of the USA with ggplot2
      1. USagstudy
    2. COThank you, I have added ```require(mapproj)```. I do not think ```ggplot(data = all_states)``` is required because the data is added later using the ```%+%``` operator. This permits the construction of maps showing just Alaska, Hawaii or the 48 contiguous states.
      singulars
      1. PORelocating Alaska and Hawaii on thematic map of the USA with ggplot2
      1. USorizon
 

Querying!

 
Guidance
A row detail

Detail views are divided into sections. All the information in the data section comes from columns in the selected row. The other sections display data from other, related rows.

Related data can be related in a to-one or a to-many fashion. Captions of data related in a to-many fashion link to a list view showing a filtered view of the table.

Try moving around until you find a non-empty to-many entry and click on the label to get to one. You can move back to the root by clicking on the database name in the header.

SQuiL has stopped working due to an internal error.

If you are curious you may find further information in the browser console, which is accessible through the devtools (F12).

Reload