小编psy*_*sky的帖子

如何将激光雷达格式 las 转换为 data.frame?

激光雷达数据只是 3D 坐标,通常采用las文件格式。\xd0\xa1内容示例

\n
library(rgdal)\nlibrary(raster)\nlibrary(tmap)\nlibrary(tmaptools)\nlibrary(lidR)\nlibrary(RStoolbox)\nlas_cat <- readLAScatalog("C:/1/078-638.las")\nsummary(las_cat)\nopt_chunk_size(las_cat) <- 500\nplot(las_cat, chunk_pattern = TRUE)\nlas_cat\n#> class       : LAScatalog (v1.2 format 1)\n#> extent      : 637999, 638240.5, 6077999, 6079999 (xmin, xmax, ymin, ymax)\n#> coord. ref. : NA \n#> area        : 483081.1 units\xc2\xb2\n#> points      : 3.68 million points\n#> density     : 7.6 points/units\xc2\xb2\n#> density      : 5.6 pulses/units\xc2\xb2\n#> num. files  : 1\n
Run Code Online (Sandbox Code Playgroud)\n

有没有办法获得data.frameR 的典型点坐标?\n作为示例,我们可以使用来自此http://data.wvgis.wvu.edu/elevation/的数据。\n此外,有没有办法获得角度反射,data.frame激光雷达文件中的反射数量?

\n

r geospatial lidar las lidr

5
推荐指数
1
解决办法
1054
查看次数

如何在 R 中创建运行给定时间的空循环

我需要创建一个运行给定时间(例如 2 小时)的空循环。该循环只是无用地运行,无论它做什么,重要的是它加载 R 执行正好 2 小时。

例如,让我们有某种脚本

model=lm(Sepal.Length~Sepal.Width,data=iris)
Run Code Online (Sandbox Code Playgroud)

在这一行之后有一个空循环,它执行了整整 2 个小时的操作

for i....
Run Code Online (Sandbox Code Playgroud)

空循环经过 2 小时完成后,继续执行后续行

summary(model)
predict(model,iris)
Run Code Online (Sandbox Code Playgroud)

(无论哪一行,重要的是在代码的某个位置循环浪费了 2 个小时)

怎样才能做到呢?感谢您的帮助。

r

5
推荐指数
1
解决办法
74
查看次数

如何在 R 中为给定函数绘制图表

假设有这个函数sqrt(x^2)+0.9*sqrt(3.3-x^2)*sin(30*pi*x)这个函数生成一个形状的图heart

在此输入图像描述

有没有办法使用 ggplot2 重现此函数来获得红色heart作为输出?感谢您的帮助。

r ggplot2

5
推荐指数
3
解决办法
609
查看次数

如何比较多个矩阵并计算 R 中的百分比差异

假设这是我的数据示例

nrow<-4
ncol<-5
m1 <- matrix(rbinom(nrow*ncol,1,.5),nrow,ncol)
m2 <- matrix(rbinom(nrow*ncol,1,.5),nrow,ncol)
m3 <- matrix(rbinom(nrow*ncol,1,.5),nrow,ncol)
Run Code Online (Sandbox Code Playgroud)

我需要依次比较 3 个矩阵。根据这个原则。例如

米1

[,1] [,2] [,3] [,4] [,5]
[1,] 1 0 1 0 1
[2,] 1 1 1 0 1
[3,] 1 0 0 1 1
[4,] 0 0 0 1 0
Run Code Online (Sandbox Code Playgroud)

平方米

[,1] [,2] [,3] [,4] [,5]
[1,] 1 0 0 0 1
[2,] 0 1 1 0 0
[3,] 1 0 0 0 0
[4,] 0 1 0 1 0
Run Code Online (Sandbox Code Playgroud)

。现在计算矩阵每列中的匹配数。取两个矩阵的第一列

    m1 m2

    1 1 …
Run Code Online (Sandbox Code Playgroud)

r matrix

4
推荐指数
1
解决办法
67
查看次数

如何将一个值乘以 R 中的随机值范围

假设我有一些值,那么它dat$x1= c(23L, 45L)。如何将这些值中的每一个随机乘以 -2 到 2 范围内的任何值?比如23*1,5或者23*-0,2等等。有什么简单的方法可以做到呢?谢谢。

r

1
推荐指数
1
解决办法
178
查看次数

如何在R中选择一个数据集中但不在另一个数据集中的类别

有2个数据集

s=structure(list(var1 = c("a", "f", "k", "tt", "ee"), var2 = c("b", 
"g", "l", "qq", "rr"), var3 = c("c", "h", "m", "ff", "cc"), var4 = c("d", 
"i", "n", "gg", "vv"), var5 = c("e", "j", "o", "aa", "xx"), metric_var = c(100L, 
200L, 300L, 567L, 789L)), class = "data.frame", row.names = c(NA, 
-5L))
Run Code Online (Sandbox Code Playgroud)

medagger=structure(list(var1 = c("a", "z", "w", "f", "k"), var2 = c("b", 
"u", "e", "g", "l"), var3 = c("c", "p", "r", "h", "m"), var4 = c("d", 
"q", "q", "i", "n"), var5 …
Run Code Online (Sandbox Code Playgroud)

r dplyr

0
推荐指数
1
解决办法
660
查看次数

标签 统计

r ×6

dplyr ×1

geospatial ×1

ggplot2 ×1

las ×1

lidar ×1

lidr ×1

matrix ×1