相关疑难解决方法(0)

在以下代码中,可以使用这些属性导入部分LAS文件版本1.1:

名单.项目格式大小必需

1. X长4个字节*
2. Y长4个字节*
3. Z长4个字节*
4. 强度无符号短2个字节
5. 返回3位(位0,1,2)3位*
6. 返回数(给定脉冲)3位(位3,4,5)3位*
7. 扫描方向标志1位(位6)1位*
8. 飞行边缘线1位(位7)1位*
9. 分类unsigned char 1字节*
10. 扫描角度等级(-90到+90) - 左侧无符号字符1个字节*
11. 用户数据unsigned char 1个字节
12. 点源ID无符号短2字节*
13. GPS时间双8字节*

并返回x,y,z和强度(上面列表的第1到4行):

allbytes <- matrix(readBin(con, "raw", n = pointDataRecordLength * numberPointRecords, size = 1, endian = "little"),
                   ncol= pointDataRecordLength, nrow = numberPointRecords, byrow = TRUE)    
close(con)
mm <- matrix(readBin(t(allbytes[,1:(3*4)]), "integer", size = 4, n = 3 * numberPointRecords, endian = "little"), ncol = 3, byrow = TRUE)

mm[,1] <- mm[ ,1] * xyzScaleOffset[1,1] + xyzScaleOffset[1, 2]
mm[,2] <- mm[ ,2] …