如何在R语言中使用sf包将ASC文件转换为栅格？

ASC是ASCII格式的栅格文件，直接用read.table()读取会丢失关键的空间元数据（比如行列数、坐标范围、分辨率等），没法直接转成有效的空间对象。以下是正确的处理流程：

步骤说明 #

1. 先用栅格专用包读取ASC文件，完整保留空间信息
2. 按需转换为sf格式的空间对象（sf主打矢量数据，若仅处理栅格，直接用栅格包操作更高效）

示例代码 #

# 安装并加载所需包
install.packages(c("terra", "sf"))
library(terra)
library(sf)

# 读取ASC栅格文件（terra包的rast函数会自动解析文件头的空间元数据）
raster_data <- rast("Myfile.asc")

# 将栅格转换为sf矢量对象（每个栅格单元转为独立的多边形要素）
sf_object <- st_as_sf(raster_data)

# 保存sf空间对象（示例保存为GeoPackage格式，也可选择shp等格式）
st_write(sf_object, "output_sf.gpkg", driver = "GPKG")

关键提示 #

• 不要用read.table()读取ASC栅格：它只会提取文件里的数值部分，完全丢失空间定位信息，得到的纯数值表格无法直接转为带空间属性的sf对象。
• 如果不需要转矢量格式，直接用terra或raster包操作栅格即可，这两个包专门针对栅格数据设计，处理效率更高。

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