Как Reorder Формального класс Пространственных Полигоны

Question

Я ищу способ повторного заказа набор Формального класс пространственных многоугольников с помощью

Я использую данные переписи США (Limited в Техас) и хочу создать 33 многоугольники из разных уездных комбинаций.

library(tmap) 
library(maptools) 
library(ggplot2) 
library(rgeos) 
library(sp) 
library(mapdata) 
library(rgdal) 
library(raster) 

# Download the map of texas and get the LMAs boundaries 
# Download shape 
f <- tempfile() 
download.file("http://www2.census.gov/geo/tiger/GENZ2010/gz_2010_us_050_00_20m.zip", destfile = f) 
unzip(f, exdir = ".") 
US <- read_shape("gz_2010_us_050_00_20m.shp") 

# Select only Texas 
Texas <- US[(US$STATE %in% c("48")),] 


# Load the LMA append data 
LMAs = read.table('LMA append data.csv',header=T, sep=',') 

# Append LMA data to Texas shape 
Texas$FIPS <- paste0(Texas$STATE, Texas$COUNTY) 
Texas <- append_data(Texas, LMAs, key.shp = "FIPS", key.data = "FIPS") 
Texas <- Texas[order(Texas$LMA),] 

# Create shape object with LMAs polygons 
Texas_LMA <- unionSpatialPolygons(Texas, IDs=Texas$LMA) 

Я попытался преобразовать Texas_LMA в SpatialPolygonsDataFrame с

# Create shape object with LMAs polygons 
Texas_LMA <- unionSpatialPolygons(Texas, IDs=Texas$LMA) 
spp <- SpatialPolygonsDataFrame(Texas_LMA,data=matrix(1:33,nrow=33,ncol=1)) 

Но это не работает для меня.

Answer 1

Ваш вопрос не очень ясен. Но я думаю, что это то, что вы после:

library(raster) 
f <- tempfile() download.file("http://www2.census.gov/geo/tiger/GENZ2010/gz_2010_us_050_00_20m.zip", destfile = f) 
unzip(f, exdir = ".") 
US <- shapefile("gz_2010_us_050_00_20m.shp") 
Texas <- US[(US$STATE %in% c("48")),] 
LMAs = read.csv('LMA append data.csv') 
Texas$FIPS <- paste0(Texas$STATE, Texas$COUNTY) 

Вы не предоставили данные LMA, поэтому угадать немного отсюда:

Texas <- merge(Texas, LMAs, by="FIPS") 
Texas_LMA <- aggregate(Texas, by='LMA')