2

我想基于一个公共列合并两个数据集。数据集 A 是一个 geoTIFF 图像，代表一个区域的 RGB 值。数据集 B 是具有相同区域 xyz 值的点云。

我想将图像中的 RGB 信息合并到 3d 数据中。我想使用两个数据集的 xy 坐标（它们在同一个坐标系中）。我写了一个脚本，灵感来自于 stackoverflow 中的代码片段，但我需要实现我的整个代码（源代码是1、2和3）。

问题是两个文件中的 xy 坐标具有不同的精度（十进制数）。数据集 A 有 0 到 2 个数字；数据集 B 有更多。我将数据集 B 的数字四舍五入为 2。现在，当数据集 A 的数字小于 2 时，我想用零填充，这样最终的合并有望起作用。

考虑到我的数据集有 >280000 行，一个简单的 if 语句会好吗？还是我应该去索引？无论如何，我在使用 R 方面相当新，所以我希望可能的海报能帮助我提供代码示例。下面是我的代码：

require(raster)
require(rgl)

setwd("C:/my/folder")

# Read tiff file
img <- stack("image.tif")

vals <- extract(img, 1:ncell(img))
coord <- xyFromCell(img, 1:ncell(img))
combine <- cbind(coord, vals)
remove(vals)
remove(coord)

# read POINTCLOUD and assign names
lidar <- read.table("lidardata.txt")
names(lidar) <- c("x","y","z")

decimalplaces <- function(x) {
  if ((x %% 1) != 0) {
    nchar(strsplit(sub('0+$', '', as.character(x)), ".", fixed=TRUE)[[1]][[2]])
  } else {
    return(0)
  }
}


# HERE I SHOULD PAD THE LIDAR VARIABLE WITH ZEROS IN DECIMAL POSITIONS WHEN THE DIGITS ARE LESS THAN 2!!!
lidar$xy <- do.call(paste0,lidar[,1:2])

combine$x <- round(combine$x, digits = 2)
combine$y <- round(combine$y, digits = 2)
combine$xy <- do.call(paste0,combine[1:2])

finaldata <- merge(combine,lidar,by = 'xy', all = FALSE)

编辑 1

正如@Heroka 所建议的，这也是激光雷达（数据集A）的外观示例，以及用零填充后应该如何。

激光雷达（原件）

x     y     z
12    9     87
11    23.4  100

激光雷达（已更改，并添加了“xy”列以加入）

x     y     z     xy
12.00 9.00  87    12.009.00
11.00 23.40 100   11.0023.40

编辑 2

counting <- sapply(lidar$x, decimalplaces) 我以某种方式设法检索了我的“激光雷达”变量（数据集 B）的所有 x 和 y 中的位数和 [0 1] 用于第二 (y) 列。我应该能够在我的 xy 数据集中找到值为 0 或 1 作为数字（不是 2）的每一行，并用 0 填充，就像上面的 LIDAR-altered 一样。

4

1 回答 1

1

我不明白为什么你需要用零填充。如果坐标属于类numeric并且两者都使用四舍五入round（这应该避免浮点精度问题），您应该能够仅通过它们合并。像这样的东西：

lidar$x <- round(lidar$x, 2)
lidar$y <- round(lidar$y, 2)
combine$x <- round(combine$x, digits = 2)
combine$y <- round(combine$y, digits = 2)

finaldata <- merge(combine, lidar, by = c("x", "y") , all = FALSE)
于 2016-01-11T15:12:10.270 回答