R语言空间数据匹配：将站点数据与最近点土壤信息关联

可以使用sp和rgdal包中的函数来进行空间连接和距离计算。具体步骤如下：

1. 将表1和表2转换成SpatialPointsDataFrame对象，其中经度和纬度作为坐标，站点名和点名作为属性。

library(sp)
library(rgdal)

# 将表1转换成SpatialPointsDataFrame对象
coordinates(df1) <- c('lon1', 'lat1')
proj4string(df1) <- CRS('+proj=longlat +datum=WGS84')
df1_sp <- SpatialPointsDataFrame(df1, data = df1)

# 将表2转换成SpatialPointsDataFrame对象
coordinates(df2) <- c('lon2', 'lat2')
proj4string(df2) <- CRS('+proj=longlat +datum=WGS84')
df2_sp <- SpatialPointsDataFrame(df2, data = df2)

2. 使用sp::over函数将每个站点匹配到最近的点，并将点的SLLL、SDUL、SSAT、SBDM信息添加到表1中。

# 使用over函数将每个站点匹配到最近的点
df1_points <- over(df1_sp, df2_sp)

# 将点的SLLL、SDUL、SSAT、SBDM信息添加到表1中
df1$SLLL <- df1_points$SLLL
df1$SDUL <- df1_points$SDUL
df1$SSAT <- df1_points$SSAT
df1$SBDM <- df1_points$SBDM

3. 可以使用sp::spDists函数计算每个站点和最近点之间的距离，以便进行质量控制。

# 计算每个站点和最近点之间的距离
df1_dist <- spDists(df1_sp, df2_sp, longlat = TRUE)

# 将距离添加到表1中
df1$dist <- apply(df1_dist, 1, min)

注意：在计算距离时，需要注意坐标系的选择和单位的转换。在本例中，使用了WGS84坐标系和经纬度单位。