如何使用R将经纬度数据框写入PostGIS数据库？

当然可以实现，以下是两种常用方法，结合你提供的数据库连接代码进行扩展：

方法一：用sf包（推荐，高效处理空间数据） #

sf包专门用于空间数据处理，能直接生成PostGIS兼容的几何列，操作更简便。

1. 准备依赖包与数据 #

library(RPostgres)
library(sf)

# 示例含经纬度的数据框
df <- data.frame(
  id = 1:3,
  lon = c(-70.60, -70.61, -70.62),
  lat = c(41.55, 41.56, 41.57),
  value = c("A", "B", "C")
)

2. 转换为sf空间对象 #

将经纬度列转换为PostGIS支持的Point类型几何列，指定坐标系为WGS84（EPSG:4326）：

sf_df <- st_as_sf(df, coords = c("lon", "lat"), crs = 4326)

3. 连接PostGIS数据库 #

注意：你提供的代码中dbname = "mytable"应为数据库名称，不是表名，请修正为实际值：

con <- dbConnect(RPostgres::Postgres(),
                 host = "xxx", 
                 dbname = "your_database_name", # 替换为实际数据库名
                 port = 5432,
                 user = "xx",
                 password = "xx")

4. 写入数据库 #

用st_write直接将sf对象写入PostGIS，支持覆盖已有表：

# 将数据写入名为target_table的表，OVERWRITE=TRUE表示覆盖已有表（可选）
st_write(sf_df, con, "target_table", layer_options = "OVERWRITE=TRUE")

方法二：纯SQL语句插入（无需sf包） #

如果不想引入额外包，可以手动通过SQL创建几何列并插入数据，适合小数据集。

1. 创建目标表（若不存在） #

先在PostGIS中创建包含几何列的表：

dbExecute(con, "
  CREATE TABLE IF NOT EXISTS target_table (
    id INTEGER,
    value TEXT,
    geometry GEOMETRY(Point, 4326)
  )
")

2. 插入数据 #

使用参数化查询避免SQL注入，循环插入每条记录：

for (i in 1:nrow(df)) {
  dbExecute(con, "
    INSERT INTO target_table (id, value, geometry)
    VALUES ($1, $2, ST_SetSRID(ST_MakePoint($3, $4), 4326))
  ", params = list(df$id[i], df$value[i], df$lon[i], df$lat[i]))
}

验证写入结果 #

可以用你提供的查询代码验证数据是否成功写入：

lon <- -70.60
lat <- 41.55

dat <- dbGetQuery(
  con, 
  paste0("SELECT * FROM target_table
         WHERE 
         ST_INTERSECTS(target_table.geometry, ST_SETSRID(ST_MAKEPOINT(", lon, ",", lat, "), 4326))")
)

print(dat)

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