Geotiff 이미지와 해당 Lidar 데이터 (x, y, z)를 UTM 좌표로 지정했습니다. Lidar 데이터를 이미지의 RGB 값과 병합해야합니다.
즉, 결국에는 Geotiff 이미지의 해당 RGB 값으로 코딩 된 LiDAR 구름 색의 각 점을 (3D) 플롯해야합니다.
QGIS를 사용하여 Lidar 데이터를 shapefile로 변환했습니다. 다음에 어떻게해야합니까?
R에서는 plot3D
함수를 시도했지만 작동하지 않았습니다. 텍스트 doc , shapefile 및 tif 이미지를 첨부하고 있습니다
편집하다:
아래와 같이 다음 프로그램을 수행했습니다.
require(raster)
require(maptools) # to take shape files
#require(car) # for scatter3D
require(plot3Drgl)
##setwd("C:\\Users\\Bibin Wilson\\Documents\\R")
##source('Lidar.r')
data = read.csv("C:\\Users\\Bibin Wilson\\Desktop\\Lidar\\lidardata.csv")
#nr = nrow(data)
nc = ncol(data)
nr = 500
require(rgdal)
X = readGDAL("C:\\Users\\Bibin Wilson\\Desktop\\Lidar\\image.tif")
topx = 4.968622208855732e+05;
topy = 5.419739403811632e+06;
final = matrix(nrow = nr, ncol = nc+2)
for(i in 1:nr) {
x = data[i,1]
y = data[i,2]
rr = round((topy-y)/0.0833)
cc = abs(round((x-topx)/0.0833))
if(rr == 0) {
rr = 1
}
if(cc == 0) {
cc = 1
}
final[i,1] = x
final[i,2] = y
final[i,3] = data[i,3]
final[i,4] = rr
final[i,5] = cc
}
for(i in 1:nr) {
x = final[i,1]
y = final[i,2]
z = final[i,3]
rr = final[i,4]
cc = final[i,5]
if(rr <= 5086 && cc<=3265) {
r = X[rr,cc,1]/255
g = X[rr,cc,2]/255
b = X[rr,cc,3]/255
c = cbind(r,g,b)
scatter3D(x,y,z,2,c)
}
}
그러나 그래프를 그리는 동안 다음 오류가 표시됩니다.
오류가
[.data.frame
(X @ 데이터, I, J, ..., 드롭 = FALSE) 미사용 인수 (1)
편집하다:
아래와 같이 RGB가없는 3D 모델을 얻었습니다.