无法安装R包rgdal及相关功能使用受阻求助

rgdal已从CRAN正式移除并停止维护，无需再尝试安装。推荐使用sf（替代空间数据读写、矢量操作）和terra（替代原raster包的栅格功能）完成需求，以下是完整实现代码：

安装依赖包 #

install.packages(c("sf", "terra", "dplyr"))

完整实现代码 #

# 加载核心包
library(sf)
library(terra)
library(dplyr)

# 1. 读取GPS点位数据（替换为你的CSV文件路径）
# 假设CSV包含lon(经度)、lat(纬度)列
gps_raw <- read.csv("your_gps_data.csv") %>%
  # 转换为sf空间对象，初始坐标系为WGS84(EPSG:4326)
  st_as_sf(coords = c("lon", "lat"), crs = 4326)

# 2. 转换为投影坐标系（必须用平面坐标系计算核密度，示例用中国东部UTM区EPSG:32650，按需替换）
gps_proj <- st_transform(gps_raw, crs = 32650)

# 3. 创建栅格模板（控制核密度结果的分辨率，这里设为100米，可调整）
raster_template <- rast(ext(gps_proj), resolution = 100, crs = st_crs(gps_proj)$wkt)

# 4. 计算核密度（adjust参数控制平滑度，数值越大结果越平滑）
kd_result <- density(rast(gps_proj), adjust = 2, template = raster_template)

# 5. 提取核密度等值线（自定义阈值，示例为0.001到0.01，步长0.002）
kd_contours <- as.contour(kd_result, levels = seq(0.001, 0.01, by = 0.002)) %>%
  st_as_sf()

# 6. 保存结果
# 保存核密度栅格（TIFF格式）
writeRaster(kd_result, "kernel_density.tif", overwrite = TRUE)
# 保存等值线矢量文件（支持Shapefile和GeoJSON）
st_write(kd_contours, "kd_contours.shp", overwrite = TRUE)
st_write(kd_contours, "kd_contours.geojson", overwrite = TRUE)

关键参数说明 #

• 投影坐标系：核密度计算基于平面距离，必须将经纬度数据转换为投影坐标系（如UTM），需根据你的GPS数据所在区域选择对应EPSG代码。
• 栅格分辨率：resolution = 100代表每个栅格单元边长为100米，数值越小结果精度越高，但计算速度越慢。
• 平滑参数adjust：调整核函数的带宽，数值越大，密度分布越平滑，可根据点位分布密度微调。
• 等值线阈值：levels参数可自定义需要提取的密度值，按需修改序列范围和步长。

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