如何用R为非结构化海岸空间点生成最优拟合LineString

问题说明 #

现有一组无序的沿海岸线分布坐标，需要将其拟合为可在leaflet中展示的空间LineString。直接通过sfheaders转换为LineString时，因点的顺序混乱导致结果不符合预期，希望实现类似smooth.spline生成的平滑拟合线，并转为空间数据格式。

用户提供的坐标数据：

x <- structure(list(anid = c(1, 1, 1, 1, 1, 1, 1,1,1,1,1,1,1, 1,1,1,1,1,1,1)
                  , longitude = c(-96.67,   -96.7,  -96.72, -96.5,  -96.97, -96.72, -96.55, -97.07, -96.48, -96.39, -96.45, -96.44, -96.98, -97.02, -96.78, -96.57, -96.36, -97.07, -96.53, -96.45)
                  , latitude = c(28.13, 28.11,  28.09,  28.27,  27.89,  28.09,  28.23,  27.81,  28.28,  28.36,  28.31,  28.32,  27.88,  27.85,  28.05,  28.21,  28.39,  27.8,   28.24,  28.31))
             , class = "data.frame", row.names = c(NA,   -20L))

用户已完成的可视化代码：

library(sf)
library(leaflet)
TX_coord <- st_as_sf(x, coords = c(2:3))
st_crs(TX_coord) <- 4326
leaflet(TX_coord) %>% addTiles() %>% addCircleMarkers()

直接转换LineString的问题代码：

library(sfheaders)
lsTX <- sf_linestring( obj = x
                      , x = "longitude"
                      , y = "latitude"
                      , linestring_id = "anid"
                      , keep = T
)
leaflet(lsTX) %>% addTiles() %>% addPolylines(weight=3)

期望的拟合效果（基于smooth.spline）：

sfit <- smooth.spline(x[,2] ~ x[,3], spar=0.80)
plot(x[,2], x[,3], pch=20, main="best fit through coordinates", 
     xlab="longitude", ylab="latitude", cex=0.50)
lines(sfit$y, sfit$x, lwd=1.25, col="blue")

实现步骤 #

1. 生成平滑拟合的坐标序列 #

先对原始点按纬度排序，确保拟合序列连续，再利用smooth.spline生成有序的拟合经纬度点。

2. 转换为空间LineString对象 #

将拟合得到的坐标序列转换为sf格式的LineString，并设置正确的坐标参考系（WGS84，EPSG:4326）。

3. 在leaflet中展示拟合结果 #

将原始点和拟合线同时加载到leaflet中，验证效果。

完整代码示例 #

# 加载所需包
library(sf)
library(leaflet)
library(stats)

# 1. 生成平滑拟合点
# 按纬度排序，保证拟合序列连续
x_sorted <- x[order(x$latitude), ]
# 执行平滑样条拟合，spar控制平滑程度
sfit <- smooth.spline(x_sorted$longitude ~ x_sorted$latitude, spar = 0.8)
# 提取拟合后的经纬度数据
fit_coords <- data.frame(
  longitude = sfit$y,
  latitude = sfit$x
)

# 2. 转换为sf LineString
fit_line <- st_sfc(
  st_linestring(as.matrix(fit_coords)),
  crs = 4326
)
fit_line_sf <- st_sf(geometry = fit_line)

# 3. 可视化原始点和拟合线
leaflet() %>%
  addTiles() %>%
  addCircleMarkers(data = TX_coord, radius = 3, color = "red") %>%
  addPolylines(data = fit_line_sf, weight = 3, color = "blue")

说明 #

• 按纬度（或经度）排序原始点，是为了避免smooth.spline因点无序生成跳跃的拟合线。
• spar参数可按需调整：值越大拟合线越平滑，值越小越贴近原始点。
• 转换后的fit_line_sf是标准sf空间对象，可直接用于leaflet或其他GIS工具展示。

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