Início / coding / Perguntas / 79575071
Accepted
erc
Asked: 2025-04-15 19:55:24 +0800 CST

Crie elipses com base no dataframe

  • 772

Tenho um dataframe com coordenadas e raios para uma elipse em cada linha.

dat <- data.frame(id = c("a", "b", "c"),
                  x = c(1, 2, 3),
                  y = c(1, 2, 3),
                  sx = c(.3, .5, .7),
                  sy = c(.2, .4, .6))
> dat
  id x y  sx  sy
1  a 1 1 0.3 0.2
2  b 2 2 0.5 0.4
3  c 3 3 0.7 0.6

Como posso criar um objeto espacial com todas as elipses?

Tentei o seguinte, o que resultou em um erro.

library(dplyr)
library(sfdep)
    dat %>%
      rowwise() %>%
      ellipse(x = x, y = y, sx = sx, sy = sy, n = 10, rotation = 0)

Error in ellipse(., x = x, y = y, sx = sx, sy = sy, n = 10, rotation = 0) : 
  unused argument (.)

2 respostas

  1. Obviamente, sfdep::ellipse()não tem um argumento de dados onde você tenta passar .. Operações por linha são ineficientes de qualquer maneira, então, em vez de se preocupar com isso, você pode usá-lo Map()para passar as colunas como vetores.

    > res <- do.call('Map', c(list(f=sfdep::ellipse), dat[-1])) |> setNames(dat$id)

    Dá

    > lapply(res, head, 3)
$a
            x        y
[1,] 1.300000 1.000000
[2,] 1.299408 1.012558
[3,] 1.297634 1.025067

$b
            x        y
[1,] 2.500000 2.000000
[2,] 2.499013 2.025116
[3,] 2.496057 2.050133

$c
            x        y
[1,] 3.700000 3.000000
[2,] 3.698619 3.037674
[3,] 3.694480 3.075200

    Trace-o

    rng <- asplit(matrixStats::colRanges(do.call(rbind, res)), 1)
par(mar=c(4.5, 4, 2, 2))
with(rng, plot(x, y, type='n'))
clr <- hcl.colors(length(res))
mapply(lines, res, col=clr)
legend('topleft', title='foo', legend=names(res), lty=1, col=clr)

    insira a descrição da imagem aqui

    • 3
  2. Best Answer

    Tente transformar data.frameem sfobjeto de pontos primeiro, depois você pode usar st_ellipse(). Com rowwise()você ainda precisa usar mutate().

    library(sfdep)
library(dplyr)
library(sf)
#> Linking to GEOS 3.13.1, GDAL 3.10.2, PROJ 9.5.1; sf_use_s2() is TRUE

dat <- data.frame(id = c("a", "b", "c"),
                  x = c(1, 2, 3),
                  y = c(1, 2, 3),
                  sx = c(.3, .5, .7),
                  sy = c(.2, .4, .6))

dat_sf <- 
  dat |> 
  st_as_sf(coords = c("x","y"), remove = FALSE) |> 
  rowwise() |> 
  mutate(geometry = st_ellipse(geometry, sx, sy, rotation = 0, n = 10)) |> 
  ungroup()

dat_sf
#> Simple feature collection with 3 features and 5 fields
#> Geometry type: LINESTRING
#> Dimension:     XY
#> Bounding box:  xmin: 0.7 ymin: 0.8097887 xmax: 3.7 ymax: 3.570634
#> CRS:           NA
#> # A tibble: 3 × 6
#>   id        x     y    sx    sy                                         geometry
#>   <chr> <dbl> <dbl> <dbl> <dbl>                                     <LINESTRING>
#> 1 a         1     1   0.3   0.2 (1.3 1, 1.242705 1.117557, 1.092705 1.190211, 0…
#> 2 b         2     2   0.5   0.4 (2.5 2, 2.404508 2.235114, 2.154508 2.380423, 1…
#> 3 c         3     3   0.7   0.6 (3.7 3, 3.566312 3.352671, 3.216312 3.570634, 2…

plot(dat_sf[,"id"], axes = TRUE)

    • 3

relate perguntas