R Shiny结合Leaflet、geosphere实现点位距离过滤时出现pointsToMatrix错误

报错来自两处代码逻辑的不匹配：

• 最关键的写法错误：写经销商筛选代码时，漏了data.frame行索引必须的逗号。R中对data.frame执行单参数方括号访问时，仅传逻辑向量会默认按列维度匹配，而非按行筛选。对比城市筛选代码data[data$City == input$City, ]是带逗号的，所以城市功能正常；但经销商的筛选代码Ddata <- Ddata[Ddata$Dealer == input$Dealer]没有逗号，最终得到的Ddata是纯向量结构，不存在经纬度列，后续构造距离计算的输入矩阵时得到空值，传入distHaversine后格式检查不通过，就抛出了.pointsToMatrix的报错。
• 距离计算的输入格式冗余：原代码用cbind拼接单列经纬度的写法容易被R自动压缩成非预期的向量结构，导致格式不匹配。

library(shiny)
library(DT)
library(dplyr)
library(leaflet)
library(geosphere)

# 内置示例数据
City <- c("Boston","Boston", "Boston", "Lowell","Lowell", "Lowell","Worcestor", "Worcestor","Worcestor","Springfield","Springfield","Springfield")
lat <- c(42.35, 42.355, 42.345, 42.63,42.625,42.635,42.27,42.265,42.275, 42.1,42.105,42.095)
lng <- c(-71.05,-71.045,-71.055,-71.316,-71.315,-71.317,-71.79,-71.785,-71.795,-72.6,-72.595,-72.605)
MassLocations <- data.frame(City, lat, lng)

Dealer <- c("West","Central", "East")
dealer_lat <- c(42.1, 42.0, 42.2)
dealer_lng <- c(-72.5,-71.8, -71.1)
MassDealers <- data.frame(Dealer, lat = dealer_lat, lng = dealer_lng)

ui <- fluidPage(
  titlePanel("Mass mpg by location"),
  fluidRow(
    column(4,
      selectInput("Dealer", "Dealer:", c("All", unique(as.character(MassDealers$Dealer)))),
      sliderInput("Distance", "Distance from Dealer:", min = 1, max = 100, value = 100),
      selectInput("City", "City:", c("All", unique(as.character(MassLocations$City))))
    )
  ),
  leafletOutput("map01", height = "500px"),
  DT::dataTableOutput("table")
)

server <- function(input, output) {
  dataShow <- reactive({
    data <- MassLocations
    Ddata <- MassDealers
    if (input$Dealer != "All") {
      # 核心修复：行筛选补逗号，明确按行匹配经销商
      Ddata <- Ddata[Ddata$Dealer == input$Dealer, ]
      # 手动构造2列矩阵，避免自动格式压缩
      dealer_point <- matrix(c(Ddata$lng, Ddata$lat), ncol = 2)
      # 修正单位换算：1英里=1609.34米，替换原近似值1600
      data$Distance <- (distHaversine(dealer_point, cbind(data$lng, data$lat))) / 1609.34
      # 距离判断包含阈值边界，符合常规交互认知
      data <- data[data$Distance <= input$Distance, ]
    }
    # 城市筛选逻辑保持原有正确写法
    if (input$City != "All") {
      data <- data[data$City == input$City, ]
    }
    # 传递选中的经销商点位给地图渲染模块
    selected_dealer <- Ddata[Ddata$Dealer == input$Dealer, ]
    list(area_data = data, dealer_info = selected_dealer)
  })

  output$table <- DT::renderDataTable(DT::datatable(dataShow()$area_data))

  output$map01 <- renderLeaflet({
    render_data <- dataShow()
    qMap <- leaflet(data = render_data$area_data) %>%
      addTiles() %>%
      addCircles(radius = 300, color = "red", label = ~City)
    # 选中经销商时，添加蓝色参考标记并自动调整视图
    if (input$Dealer != "All" && nrow(render_data$dealer_info) > 0) {
      qMap <- qMap %>% 
        addCircleMarkers(
          data = render_data$dealer_info,
          radius = 10,
          color = "blue",
          label = ~Dealer
        ) %>% 
        setView(lng = render_data$dealer_info$lng, lat = render_data$dealer_info$lat, zoom = 10)
    }
    qMap
  })
}

shinyApp(ui = ui, server = server)

• 补全经销商筛选的行索引逗号，解决根本的格式报错
• 距离计算环节手动构造标准2列矩阵，避免R自动压缩单行结构导致的格式不匹配
• 米转英里的换算系数替换为标准值，计算结果更精准
• 地图新增选中经销商的蓝色醒目标记，自动定位到参考点位置，避免用户找不到距离计算的基准点
• 距离筛选判断从严格小于调整为小于等于，匹配滑块设置阈值的用户预期

内容的提问来源于stack exchange，提问作者MikeHealth21