On 09/04/2015 03:29 PM, Ignacio Martinez wrote: > I'm trying to create a map using leaflet. I'm basically following this > <https://github.com/rstudio/leaflet/issues/169>but my data is a bit > different this time around. Additionally, I'm using dplyr instead of > data.table. > > This is the code i'm trying to run: > > counts <- > structure( > list( > STUSPS = c( > "CA", "NC", "TX", "FL", "VA", "OH", > "NY", "GA", "IL", "WA", "CO", "AZ", "MD", "LA", "SC", "KS", > "DC", > "TN", "MA", "MI", "MN", "NJ", "WI", "PA", "AL", "KY", "OK", > "MO", > "ID", "MS", "NM", "IN", "NV", "AR", "OR", "AK", "UT", "IA", > "MT", > "HI", "NE", "CT", "WV", "WY", "SD", "VT", "ND", "ME", "RI", > "NH", > "DE", "PR", "GU", "VI", "MP", "AS" > ), count = c( > 36511L, 27513L, > 25861L, 21099L, 19415L, 17012L, 15860L, 14362L, 13923L, 13753L, > 11671L, 10540L, 9777L, 8923L, 8355L, 8219L, 8142L, 8076L, 7985L, > 7770L, 7662L, 7531L, 7273L, 7212L, 7042L, 6708L, 6674L, 6429L, > 6239L, 5580L, 5395L, 5172L, 5013L, 4933L, 4911L, 4797L, 4694L, > 4458L, 3873L, 3757L, 3608L, 3111L, 2604L, 2217L, 2156L, 2143L, > 2030L, 1544L, 1417L, 1168L, 772L, 531L, 148L, 63L, 7L, 2L > ) > ), .Names = c("STUSPS", > "count"), class = c("tbl_df", "tbl", "data.frame"), > row.names = c(NA,-56L) > ) > > > x = c("leaflet", "rgdal", "maptools", "mapproj", "rgeos", "dplyr") > lapply(x, library, character.only = TRUE) > # From https://www.census.gov/geo/maps-data/data/cbf/cbf_state.html > states <- readOGR(dsn = "./cb_2014_us_state_5m.shp", > layer = "cb_2014_us_state_5m", verbose = FALSE) > # Make a copy of the SPDF attribute table, and then work normally, as > with any data.frame/data.table object > states.df <- states@data > # Create an explicit attribute to keep polygons IDs (useful to > "re-attach" the table to the polygons later) > states.df <- states.df %>% mutate(rn=row.names(states)) > > # join them > states.df <- full_join(x = states.df, y = counts, by="STUSPS") %>% > filter(!(STUSPS%in%c("PR", "na", "MH", "FM", "PW"))) %>% na.omit() %>% > mutate(popup=paste0('<b>',NAME,'</b>','<br>', > '<b>count:</b> ', > prettyNum(count,big.mark=",",scientific=FALSE)))
at this stage, you have:
> dim(states@data)
[1] 56 9
> dim(states.df)
[1] 55 12
>
> # Re-attach the attribute table to the SPDF
> states@data <- states.df
now, you created an invalid object:
> length(geometry(states))
[1] 56
> dim(states@data)
[1] 55 12
this confirms once more what has been said on this list so often, that
instead of using constructor functions such as
SpatialPolygonsDataFrame(), assigning slots directly is dangerous and
better left to those who know what they do (and check sanity):
> states <- SpatialPolygonsDataFrame(geometry(states), states.df)
Error in SpatialPolygonsDataFrame(geometry(states), states.df) :
Object length mismatch:
geometry(states) has 56 Polygons objects, but states.df has 55 rows
> # Make sure polygons IDs and data.frame row.names match
> states <- spChFIDs(states, states$rn)
>
> # Create map as usual...
>
> The error I get is:
>
> Error in spChFIDs(SP, x) : lengths differ
>
>
> Thanks for the help!
>
>
> Ignacio
>
> [[alternative HTML version deleted]]
>
> _______________________________________________
> R-sig-Geo mailing list
> R-sig-Geo@r-project.org
> https://stat.ethz.ch/mailman/listinfo/r-sig-geo
>
--
Edzer Pebesma
Institute for Geoinformatics (ifgi), University of Münster,
Heisenbergstraße 2, 48149 Münster, Germany; +49 251 83 33081
Journal of Statistical Software: http://www.jstatsoft.org/
Computers & Geosciences: http://elsevier.com/locate/cageo/
Spatial Statistics Society http://www.spatialstatistics.info
signature.asc
Description: OpenPGP digital signature
_______________________________________________ R-sig-Geo mailing list R-sig-Geo@r-project.org https://stat.ethz.ch/mailman/listinfo/r-sig-geo
