Fixing maps library data for Pacific centred (0°-360° longitude) display

by

You could use the fact that internally, a map object returned by the map() function can be recalculated and used again in the map() function. I’d create a list with individual polygons, check which ones have very different longitude values, and rearrange those ones. I gave an example of this approach in the function below*, which allows something like :

plot.map("world", center=180, col="white",bg="gray",
   fill=TRUE,ylim=c(-60,90),mar=c(0,0,0,0))

to get

Corrected map center 180

If I were you, I’d shift everything a bit more, like in :

plot.map("world", center=200, col="white",bg="gray",
   fill=TRUE,ylim=c(-60,90),mar=c(0,0,0,0))

Corrected map center 200

The function :

plot.map<- function(database,center,...){
    Obj <- map(database,...,plot=F)
    coord <- cbind(Obj[[1]],Obj[[2]])

    # split up the coordinates
    id <- rle(!is.na(coord[,1]))
    id <- matrix(c(1,cumsum(id$lengths)),ncol=2,byrow=T)
    polygons <- apply(id,1,function(i){coord[i[1]:i[2],]})

    # split up polygons that differ too much
    polygons <- lapply(polygons,function(x){
        x[,1] <- x[,1] + center
        x[,1] <- ifelse(x[,1]>180,x[,1]-360,x[,1])
        if(sum(diff(x[,1])>300,na.rm=T) >0){
          id <- x[,1] < 0
          x <- rbind(x[id,],c(NA,NA),x[!id,])
       }
       x
    })
    # reconstruct the object
    polygons <- do.call(rbind,polygons)
    Obj[[1]] <- polygons[,1]
    Obj[[2]] <- polygons[,2]

    map(Obj,...)
}

*Note that this function only takes positive center values. It’s easily adapted to allow for center values in both directions, but I didn’t bother anymore as that’s trivial.

More Related Contents:

  1. rgdal package installation
  2. Plotting bar charts on map using ggplot2?
  3. How to put a geom_sf produced map on top of a ggmap produced raster
  4. How to GeoCode a simple address using Data Science Toolbox
  5. Choropleth map in ggplot with polygons that have holes
  6. Plot coordinates on map
  7. Labeling center of map polygons in R ggplot
  8. Maps, ggplot2, fill by state is missing certain areas on the map
  9. Convert latitude and longitude coordinates to country name in R
  10. Sum by two variables
  11. Drop unused factor levels in a subsetted data frame
  12. Use dynamic name for new column/variable in `dplyr`
  13. Changing fonts in ggplot2
  14. geom_rect and alpha – does this work with hard coded values?
  15. Using loops with knitr to produce multiple pdf reports… need a little help to get me over the hump
  16. Argument is of length zero in if statement
  17. How do I name the “row names” column in r
  18. How to change ‘Maximum upload size exceeded’ restriction in Shiny and save user file inputs?
  19. Dynamically add column names to data.table when aggregating
  20. Detect text language in R
  21. Finding the index inside a vector satisfying a condition
  22. finding close match from data frame 1 in data fame 2
  23. Showing string in formula and not as variable in lm fit
  24. Combine several data frames in the global environment by row (rbind)
  25. aggregate() puts multiple output columns in a matrix instead
  26. select columns based on multiple strings with dplyr contains()
  27. Select along one of n dimensions in array
  28. Calculate number of days between two dates in r
  29. Importing csv file into R – numeric values read as characters
  30. Read csv from specific row

Leave a Comment