AskOverflow.Dev

AskOverflow.Dev Logo AskOverflow.Dev Logo

AskOverflow.Dev Navigation

  • Início
  • system&network
  • Ubuntu
  • Unix
  • DBA
  • Computer
  • Coding
  • LangChain

Mobile menu

Close
  • Início
  • system&network
    • Recentes
    • Highest score
    • tags
  • Ubuntu
    • Recentes
    • Highest score
    • tags
  • Unix
    • Recentes
    • tags
  • DBA
    • Recentes
    • tags
  • Computer
    • Recentes
    • tags
  • Coding
    • Recentes
    • tags
Início / coding / Perguntas / 78368476
Accepted
M. Beausoleil
M. Beausoleil
Asked: 2024-04-23 03:24:11 +0800 CST2024-04-23 03:24:11 +0800 CST 2024-04-23 03:24:11 +0800 CST

Por que e como consertar o Mapview não mostra todos os pontos e o st_buffer mescla algumas áreas R?

  • 772

Estou construindo um mapa com o pacote mapview. Porém, quando uso st_sample e gero pontos, alguns não são mapeados. Além disso, st_buffer mescla algumas áreas (onde deveriam estar separadas).

library(sf)
library(mapview)
library(magrittr)
set.seed(145)

chemins.32188 = st_zm(chemins2) %>% st_transform(crs = 32188)

rd.pts.sentiers <- st_sample(x = chemins.32188, size = 10)

buff.pt = st_buffer(x = rd.pts.sentiers, 
                    dist = units::set_units(x = 100, 
                                           value = "m"))
mapviewOptions(fgb = FALSE)
mapview(chemins.32188) + 
  mapview(rd.pts.sentiers) +
  mapview(buff.pt)

insira a descrição da imagem aqui

Quando uso plot, funciona... mas o buffer ainda está mesclado.

plot(chemins.32188$geom) 
plot(buff.pt, add= TRUE, col = scales::alpha("blue", .5))
plot(rd.pts.sentiers, add= TRUE, pch = 19, col = "red")

insira a descrição da imagem aqui

Estes são os dados

chemins2 = structure(list(Name = "Path Measure", description = NA_character_, 
               timestamp = structure(NA_real_, class = c("POSIXct", "POSIXt"
               )), begin = structure(NA_real_, class = c("POSIXct", "POSIXt"
               )), end = structure(NA_real_, class = c("POSIXct", "POSIXt"
               )), altitudeMode = NA_character_, tessellate = -1L, extrude = 0L, 
               visibility = -1L, drawOrder = NA_integer_, icon = NA_character_, 
               sentier = NA_character_, layer = "path Parc — Path Measure", 
               path = "~/Desktop/path Parc.kml|layername=Path Measure|geometrytype=LineString25D|uniqueGeometryType=yes", 
               geom = structure(list(structure(list(structure(c(299434.130314458, 
                                                                300531.43655987, 299434.130314458, 5046816.92873926, 5047257.53966348, 
                                                                5046816.92873926), dim = 3:2), structure(c(299434.130314458, 
                                                                                                           300401.157828111, 299434.130314458, 5046816.92873926, 5046249.12770288, 
                                                                                                           5046816.92873926), dim = 3:2)), class = c("XY", "MULTILINESTRING", 
                                                                                                                                                     "sfg"))), n_empty = 0L, crs = structure(list(input = "EPSG:32188", 
                                                                                                                                                                                                  wkt = "PROJCRS[\"NAD83 / MTM zone 8\",\n    BASEGEOGCRS[\"NAD83\",\n        DATUM[\"North American Datum 1983\",\n            ELLIPSOID[\"GRS 1980\",6378137,298.257222101,\n                LENGTHUNIT[\"metre\",1]]],\n        PRIMEM[\"Greenwich\",0,\n            ANGLEUNIT[\"degree\",0.0174532925199433]],\n        ID[\"EPSG\",4269]],\n    CONVERSION[\"MTM zone 8\",\n        METHOD[\"Transverse Mercator\",\n            ID[\"EPSG\",9807]],\n        PARAMETER[\"Latitude of natural origin\",0,\n            ANGLEUNIT[\"degree\",0.0174532925199433],\n            ID[\"EPSG\",8801]],\n        PARAMETER[\"Longitude of natural origin\",-73.5,\n            ANGLEUNIT[\"degree\",0.0174532925199433],\n            ID[\"EPSG\",8802]],\n        PARAMETER[\"Scale factor at natural origin\",0.9999,\n            SCALEUNIT[\"unity\",1],\n            ID[\"EPSG\",8805]],\n        PARAMETER[\"False easting\",304800,\n            LENGTHUNIT[\"metre\",1],\n            ID[\"EPSG\",8806]],\n        PARAMETER[\"False northing\",0,\n            LENGTHUNIT[\"metre\",1],\n            ID[\"EPSG\",8807]]],\n    CS[Cartesian,2],\n        AXIS[\"easting (E(X))\",east,\n            ORDER[1],\n            LENGTHUNIT[\"metre\",1]],\n        AXIS[\"northing (N(Y))\",north,\n            ORDER[2],\n            LENGTHUNIT[\"metre\",1]],\n    USAGE[\n        SCOPE[\"Engineering survey, topographic mapping.\"],\n        AREA[\"Canada - Quebec between 75°W and 72°W.; Canada - Ontario - east of 75°W.\"],\n        BBOX[44.98,-75,62.53,-72]],\n    ID[\"EPSG\",32188]]"), class = "crs"), class = c("sfc_MULTILINESTRING", 
                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                      "sfc"), precision = 0, bbox = structure(c(xmin = 299434.130314458, 
                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                ymin = 5046249.12770288, xmax = 300531.43655987, ymax = 5047257.53966348
                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                      ), class = "bbox"))), row.names = 1L, sf_column = "geom", agr = structure(c(Name = NA_integer_, 
                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                  description = NA_integer_, timestamp = NA_integer_, begin = NA_integer_, 
                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                  end = NA_integer_, altitudeMode = NA_integer_, tessellate = NA_integer_, 
                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                  extrude = NA_integer_, visibility = NA_integer_, drawOrder = NA_integer_, 
                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                  icon = NA_integer_, sentier = NA_integer_, layer = NA_integer_, 
                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                  path = NA_integer_), class = "factor", levels = c("constant", 
                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                                    "aggregate", "identity")), class = c("sf", "data.frame"))
  • 1 1 respostas
  • 16 Views

1 respostas

  • Voted
  1. Best Answer
    margusl
    2024-04-23T06:05:27+08:002024-04-23T06:05:27+08:00

    No seu reprex você tem MULTILINESTRING2 LINESTRINGs, st_sample()retorna 2 MULTIPOINTs e com essa semente específica, o primeiro está vazio (ou seja, todos os pontos são do mesmo sub LINESTRING):

    set.seed(145)
    chemins.32188 = st_zm(chemins2) %>% st_transform(crs = 32188)
    st_as_text(chemins.32188$geom)
    #> [1] "MULTILINESTRING ((299434.1 5046817, 300531.4 5047258, 299434.1 5046817), (299434.1 5046817, 300401.2 5046249, 299434.1 5046817))"
    
    rd.pts.sentiers <- st_sample(x = chemins.32188, size = 10)
    rd.pts.sentiers
    #> Geometry set for 2 features  (with 1 geometry empty)
    #> Geometry type: MULTIPOINT
    #> Dimension:     XY
    #> Bounding box:  xmin: 299516.7 ymin: 5046255 xmax: 300390.5 ymax: 5046768
    #> Projected CRS: NAD83 / MTM zone 8
    #> MULTIPOINT EMPTY
    #> MULTIPOINT ((300390.5 5046255), (299516.7 50467...
    

    Aparentemente mapview(ou leaflet) não está muito satisfeito com geometrias vazias (ponto) e você também deve receber o seguinte aviso de mapview():

    Warning message:
    In validateCoords(lng, lat, funcName) :
      Data contains 1 rows with either missing or invalid lat/lon values and will be ignored
    

    E quando você cria um buffer de MULTIPOINT rd.pts.sentiers, cada recurso será um MULTIPOLYGON, portanto, sobreposições em um único recurso causarão uma mesclagem de polígono.

    Para corrigir, você pode querer verificar amostras em busca de geometrias vazias e converter amostras em POINTs antes de criar buffers:

    library(sf)
    #> Linking to GEOS 3.11.2, GDAL 3.6.2, PROJ 9.2.0; sf_use_s2() is TRUE
    library(mapview)
    library(magrittr)
    
    set.seed(145)
    chemins.32188 <- st_zm(chemins2) %>% st_transform(crs = 32188)
    rd.pts.sentiers <- st_sample(x = chemins.32188, size = 10)
    rd.pts.sentiers
    #> Geometry set for 2 features  (with 1 geometry empty)
    #> Geometry type: MULTIPOINT
    #> Dimension:     XY
    #> Bounding box:  xmin: 299516.7 ymin: 5046255 xmax: 300390.5 ymax: 5046768
    #> Projected CRS: NAD83 / MTM zone 8
    #> MULTIPOINT EMPTY
    #> MULTIPOINT ((300390.5 5046255), (299516.7 50467...
    
    # keep only non-empty samples, cast from MULTIPOINT to POINT
    rd.pts.sentiers <- 
      rd.pts.sentiers[!st_is_empty(rd.pts.sentiers)] %>% 
      st_cast("POINT")
    rd.pts.sentiers
    #> Geometry set for 10 features 
    #> Geometry type: POINT
    #> Dimension:     XY
    #> Bounding box:  xmin: 299516.7 ymin: 5046255 xmax: 300390.5 ymax: 5046768
    #> Projected CRS: NAD83 / MTM zone 8
    #> First 5 geometries:
    #> POINT (300390.5 5046255)
    #> POINT (299516.7 5046768)
    #> POINT (300002.8 5046483)
    #> POINT (299737 5046639)
    #> POINT (300193.1 5046371)
    
    # get 10 POLYGON features instead of 2(1) MULTIPOLYGONs
    buff.pt = st_buffer(x = rd.pts.sentiers, 
                        dist = units::set_units(x = 100, 
                                                value = "m"))
    buff.pt
    #> Geometry set for 10 features 
    #> Geometry type: POLYGON
    #> Dimension:     XY
    #> Bounding box:  xmin: 299416.7 ymin: 5046155 xmax: 300490.5 ymax: 5046868
    #> Projected CRS: NAD83 / MTM zone 8
    #> First 5 geometries:
    #> POLYGON ((300490.5 5046255, 300490.3 5046250, 3...
    #> POLYGON ((299616.7 5046768, 299616.5 5046763, 2...
    #> POLYGON ((300102.8 5046483, 300102.7 5046478, 3...
    #> POLYGON ((299837 5046639, 299836.9 5046634, 299...
    #> POLYGON ((300293.1 5046371, 300293 5046366, 300...
    
    mapviewOptions(fgb = FALSE)
    mapview(chemins.32188) + 
      mapview(rd.pts.sentiers) +
      mapview(buff.pt)
    

    Criado em 23/04/2024 com reprex v2.1.0

    • 1

relate perguntas

  • Adicionar número de série para atividade de cópia ao blob

  • A fonte dinâmica do empacotador duplica artefatos

  • Selecione linhas por grupo com 1s consecutivos

  • Lista de chamada de API de gráfico subscritoSkus estados Privilégios insuficientes enquanto os privilégios são concedidos

  • Função para criar DFs separados com base no valor da coluna

Sidebar

Stats

  • Perguntas 205573
  • respostas 270741
  • best respostas 135370
  • utilizador 68524
  • Highest score
  • respostas
  • Marko Smith

    Vue 3: Erro na criação "Identificador esperado, mas encontrado 'import'" [duplicado]

    • 1 respostas
  • Marko Smith

    Por que esse código Java simples e pequeno roda 30x mais rápido em todas as JVMs Graal, mas não em nenhuma JVM Oracle?

    • 1 respostas
  • Marko Smith

    Qual é o propósito de `enum class` com um tipo subjacente especificado, mas sem enumeradores?

    • 1 respostas
  • Marko Smith

    Como faço para corrigir um erro MODULE_NOT_FOUND para um módulo que não importei manualmente?

    • 6 respostas
  • Marko Smith

    `(expression, lvalue) = rvalue` é uma atribuição válida em C ou C++? Por que alguns compiladores aceitam/rejeitam isso?

    • 3 respostas
  • Marko Smith

    Quando devo usar um std::inplace_vector em vez de um std::vector?

    • 3 respostas
  • Marko Smith

    Um programa vazio que não faz nada em C++ precisa de um heap de 204 KB, mas não em C

    • 1 respostas
  • Marko Smith

    PowerBI atualmente quebrado com BigQuery: problema de driver Simba com atualização do Windows

    • 2 respostas
  • Marko Smith

    AdMob: MobileAds.initialize() - "java.lang.Integer não pode ser convertido em java.lang.String" para alguns dispositivos

    • 1 respostas
  • Marko Smith

    Estou tentando fazer o jogo pacman usando apenas o módulo Turtle Random e Math

    • 1 respostas
  • Martin Hope
    Aleksandr Dubinsky Por que a correspondência de padrões com o switch no InetAddress falha com 'não cobre todos os valores de entrada possíveis'? 2024-12-23 06:56:21 +0800 CST
  • Martin Hope
    Phillip Borge Por que esse código Java simples e pequeno roda 30x mais rápido em todas as JVMs Graal, mas não em nenhuma JVM Oracle? 2024-12-12 20:46:46 +0800 CST
  • Martin Hope
    Oodini Qual é o propósito de `enum class` com um tipo subjacente especificado, mas sem enumeradores? 2024-12-12 06:27:11 +0800 CST
  • Martin Hope
    sleeptightAnsiC `(expression, lvalue) = rvalue` é uma atribuição válida em C ou C++? Por que alguns compiladores aceitam/rejeitam isso? 2024-11-09 07:18:53 +0800 CST
  • Martin Hope
    The Mad Gamer Quando devo usar um std::inplace_vector em vez de um std::vector? 2024-10-29 23:01:00 +0800 CST
  • Martin Hope
    Chad Feller O ponto e vírgula agora é opcional em condicionais bash com [[ .. ]] na versão 5.2? 2024-10-21 05:50:33 +0800 CST
  • Martin Hope
    Wrench Por que um traço duplo (--) faz com que esta cláusula MariaDB seja avaliada como verdadeira? 2024-05-05 13:37:20 +0800 CST
  • Martin Hope
    Waket Zheng Por que `dict(id=1, **{'id': 2})` às vezes gera `KeyError: 'id'` em vez de um TypeError? 2024-05-04 14:19:19 +0800 CST
  • Martin Hope
    user924 AdMob: MobileAds.initialize() - "java.lang.Integer não pode ser convertido em java.lang.String" para alguns dispositivos 2024-03-20 03:12:31 +0800 CST
  • Martin Hope
    MarkB Por que o GCC gera código que executa condicionalmente uma implementação SIMD? 2024-02-17 06:17:14 +0800 CST

Hot tag

python javascript c++ c# java typescript sql reactjs html

Explore

  • Início
  • Perguntas
    • Recentes
    • Highest score
  • tag
  • help

Footer

AskOverflow.Dev

About Us

  • About Us
  • Contact Us

Legal Stuff

  • Privacy Policy

Language

  • Pt
  • Server
  • Unix

© 2023 AskOverflow.DEV All Rights Reserve