为何在R与QGIS中无法正确生成目标GPS轨迹缓冲区？

处理动物GPS轨迹数据集时，需在轨迹起点生成半径50km的缓冲区，但CRS始终出错：

• QGIS中，项目CRS为EPSG:32632，轨迹数据是EPSG:4326的点，但生成的缓冲区出现在赤道位置
• R中生成的缓冲区位置完全错误，导致裁剪后的codoni_clipped为空，绘图时缓冲区与轨迹错位（红色轨迹，右侧小紫点为缓冲区，形状尺寸正确但位置偏移）

数据结构 #

> str(codoni)
'data.frame':   149366 obs. of  8 variables:
$ id      : chr  "N34246@096AFE12" "N34246@096AFE12" "N34246@096AFE12" "N34246@096AFE12" ...
$ DateTime: POSIXct, format: "2022-02-23 03:36:00" "2022-02-23 05:36:00" "2022-02-23 07:36:00" "2022-02-23 09:36:00" ...
$ lat     : num  45.5 45.5 45.5 45.5 45.5 ...
$ lon     : num  12.6 12.6 12.6 12.6 12.6 ...
$ Altitude: int  0 0 0 6 0 9 0 82 0 0 ...
$ Sex     : int  1 1 1 1 1 1 1 1 1 1 ...
$ Age     : int  6 6 6 6 6 6 6 6 6 6 ...
$ Tag     : chr  "GPS-GSM_TECHNOSMART" "GPS-GSM_TECHNOSMART" "GPS-GSM_TECHNOSMART" "GPS-GSM_TECHNOSMART" ...

原始R代码 #

cavallino <- c(45.501518, 12.565928) 
cavallino_point <- st_sfc(st_point(cavallino))
st_crs(cavallino_point) <- 4326
# buffer around cavallino in EPSG:4326
buffer_4326 <- st_buffer(cavallino_point, dist = 50000)
st_crs(buffer_4326) <- st_crs(4326)
buffer_32632 <- st_transform(buffer_4326, 32632)
st_geometry(buffer_32632)
st_crs(buffer_32632)
# clip tracks to the extent of the buffer
codoni_sf <- st_as_sf(codoni, coords = c("lon", "lat"), crs = 4326) %>% st_transform(32632)
st_crs(codoni_sf)
codoni_clipped <- st_intersection(codoni_sf, buffer_32632)
st_geometry(codoni_clipped)
st_crs(codoni_clipped)
# plot tracks and buffer together
final_plot <- ggplot() +
geom_sf(data = codoni_sf, color = "red", size = 0.7) +
geom_sf(data = buffer_32632, fill = "blue", alpha = 0.3) +
labs(title = "Tracks and Buffer (50km radius)",
x = "Longitude",
y = "Latitude") +
theme_minimal()
print(final_plot)

样本数据输出 #

> dput(head(codoni, 5)) 
structure(list(id = c("N34246@096AFE12", "N34246@096AFE12", "N34246@096AFE12", "N34246@096AFE12", "N34246@096AFE12"), DateTime = structure(c(1645583760, 1645590960, 1645598160, 1645605360, 1645612560), class = c("POSIXct", "POSIXt"), tzone = ""), lat = c(45.504295, 45.50255, 45.50186, 45.5025867, 45.5005833), lon = c(12.5623333, 12.566, 12.56615, 12.56315, 12.5618667), Altitude = c(0L, 0L, 0L, 6L, 0L), Sex = c(1L, 1L, 1L, 1L, 1L), Age = c(6L, 6L, 6L, 6L, 6L), Tag = c("GPS-GSM_TECHNOSMART", "GPS-GSM_TECHNOSMART", "GPS-GSM_TECHNOSMART", "GPS-GSM_TECHNOSMART", "GPS-GSM_TECHNOSMART")), row.names = c(NA, 5L), class = "data.frame")

核心问题 #

1. 坐标顺序错误：st_point()要求坐标顺序为经度在前，纬度在后（符合EPSG:4326规范），但你传入的是c(纬度, 经度)，直接导致点的位置完全偏移。
2. 冗余设置缓冲区CRS：st_buffer()会继承输入对象的CRS，手动设置st_crs(buffer_4326) <- st_crs(4326)属于无效操作。
3. 地理坐标系下创建缓冲区：EPSG:4326是地理坐标系（单位为度），直接用dist=50000会被解析为50000度，完全不符合需求——正确做法是先转换到投影坐标系（如EPSG:32632，单位为米）再创建缓冲区。

修正后的R代码 #

# 正确设置坐标顺序：经度在前，纬度在后
cavallino <- c(12.565928, 45.501518) 
cavallino_point <- st_sfc(st_point(cavallino), crs = 4326)

# 先转换到投影坐标系EPSG:32632，再创建50km缓冲区
buffer_32632 <- cavallino_point %>%
  st_transform(32632) %>%
  st_buffer(dist = 50000)

# 处理轨迹数据：转换到投影坐标系
codoni_sf <- st_as_sf(codoni, coords = c("lon", "lat"), crs = 4326) %>% 
  st_transform(32632)

# 裁剪轨迹
codoni_clipped <- st_intersection(codoni_sf, buffer_32632)

# 绘图（修正坐标轴标签为投影坐标系单位）
final_plot <- ggplot() +
  geom_sf(data = codoni_sf, color = "red", size = 0.7) +
  geom_sf(data = buffer_32632, fill = "blue", alpha = 0.3) +
  labs(title = "Tracks and Buffer (50km radius)",
       x = "Easting (m)",
       y = "Northing (m)") +
  theme_minimal()
print(final_plot)

QGIS问题解决思路 #

QGIS中缓冲区出现在赤道，同样是坐标顺序或CRS匹配错误：

• 确保轨迹数据的CRS正确识别为EPSG:4326
• 创建缓冲区前，先将轨迹起点转换为EPSG:32632（投影坐标系，米为单位），再生成缓冲区——直接在EPSG:4326下创建50km缓冲区会因单位为度出现异常位置。

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