SpatialPointsDataFrame'i SpatialLinesDataFrame ürününe dönüştürün R

Question

Hurricane parçalarını çizmek için HURDAT veri kümesinde çalışıyorum. Şu anda yıl 2004.

> str(cluster.2004.sdf) 
Formal class 'SpatialPointsDataFrame' [package "sp"] with 5 slots 
    ..@ data  :'data.frame': 2693 obs. of 4 variables: 
    .. ..$ Sid  : int [1:2693] 1331 1331 1331 1331 1331 1331 1331 1331 1331 1331 ... 
    .. ..$ clusterid: num [1:2693] 2 2 2 2 2 2 2 2 2 2 ... 
    .. ..$ name  : Factor w/ 269 levels "","ABBY  ",..: 6 6 6 6 6 6 6 6 6 6 ... 
    .. ..$ WmaxS : num [1:2693] 78.9 82.8 80.9 70.9 76.9 ... 
    ..@ coords.nrs : num(0) 
    ..@ coords  : num [1:2693, 1:2] 754377 612852 684956 991386 819565 ... 
    .. ..- attr(*, "dimnames")=List of 2 
    .. .. ..$ : NULL 
    .. .. ..$ : chr [1:2] "lon" "lat" 
    ..@ bbox  : num [1:2, 1:2] -3195788 1362537 4495870 9082812 
    .. ..- attr(*, "dimnames")=List of 2 
    .. .. ..$ : chr [1:2] "lon" "lat" 
    .. .. ..$ : chr [1:2] "min" "max" 
    ..@ proj4string:Formal class 'CRS' [package "sp"] with 1 slots 
    .. .. ..@ projargs: chr "+proj=lcc +lat_1=60 +lat_2=30 +lon_0=-60 +ellps=WGS84" 

    > summary(cluster.2004.sdf) 
Object of class SpatialPointsDataFrame 
Coordinates: 
     min  max 
lon -3195788 4495870 
lat 1362537 9082812 
Is projected: TRUE 
proj4string : 
[+proj=lcc +lat_1=60 +lat_2=30 +lon_0=-60 +ellps=WGS84] 
Number of points: 2693 
Data attributes: 
     Sid   clusterid    name   WmaxS  
Min. :1331 Min. :1.000 IVAN  :517 Min. : 14.83 
1st Qu.:1334 1st Qu.:2.000 FRANCES :403 1st Qu.: 31.35 
Median :1337 Median :3.000 JEANNE :379 Median : 50.04 
Mean :1337 Mean :2.898 KARL  :283 Mean : 61.66 
3rd Qu.:1339 3rd Qu.:4.000 DANIELLE :271 3rd Qu.: 90.40 
Max. :1341 Max. :4.000 BONNIE :253 Max. :142.52 
           (Other) :587 

Her fırtına "Sid" etiketli eşsiz bir fırtına kimliği referansı vardır için şuna benzer R bir SpatialPointsDataFrame nesnesi ürettiler. SpatialPointsDataFrame öğesini "Sid" ile gruplamak ve tüm noktaları bir Satıra dönüştürmek istiyorum.

Plyr paketinden ddply ile devam ettim, ama ne yaptığımı bilmiyorum. Bunu, veri çerçevesindeki her satırın bir döngüsünde döngü yaparak ve bir listeye koordinat ekleyerek ve ardından sp listesinden Satırlar işlevini kullanarak bu listeyi dönüştürerek yapabilirim. Bununla birlikte, daha fazla R dönüşümünü tercih ederim. Teşekkür Richard

Answer 1

mdsumner çözümüyle ilgili sorun çıkan data.frame her satır için bir satır olması gerektiğini, ancak onun kodunda her bir nokta için bir satır vardır. düzeltilmiş kod şöyle olacaktır:

## example data 
d <- data.frame(x=runif(7), y=runif(7), id = c(rep("a", 3), rep("b", 4))) 

library(sp)  
coordinates(d) <- ~x+y 

## list of Lines per id, each with one Line in a list 
x <- lapply(split(d, d$id), function(x) Lines(list(Line(coordinates(x))), x$id[1L])) 

# the corrected part goes here: 
lines <- SpatialLines(x) 
data <- data.frame(id = unique(d$id)) 
rownames(data) <- data$id 
l <- SpatialLinesDataFrame(lines, data) 

Yani temelde sorun kimliğiyle gruplandırılmış hatlar için bir data.frame (her hat için bir satır) oluşturmak zorunda olmasıdır. Yukarıdaki durumda, id'dan başka bir veri olmadığı zaman oldukça kolaydır. sonuçta bir liste ile anlaşma

data <- sqldf(' 
select id, max(something), sum(something_else) 
from d 
group by id 
') 

Answer 2

## example data 
d <- data.frame(x=runif(7), y=runif(7), id = c(rep("a", 3), rep("b", 4))) 
##split(d, d$id) 

library(sp)  
coordinates(d) <- ~x+y 

## list of Lines per id, each with one Line in a list 
x <- lapply(split(d, d$id), function(x) Lines(list(Line(coordinates(x))), x$id[1L])) 

## or one Lines in a list, with all Line objects 
## x <- list(Lines(lapply(split(d, d$id), function(x) Line(coordinates(x))), paste(unique(d$id), collapse = "_"))) 

## etc. 
SpatialLines(x, CRS(as.character(NA))) 

## need to be careful here, assuming one Lines per original row 
## and we trash the original rownames . . . 
SpatialLinesDataFrame(SpatialLines(x, CRS(as.character(NA))), d[,"id", drop = FALSE], match.ID = FALSE)