激光雷达数据只是 3D 坐标,通常采用las文件格式。\xd0\xa1内容示例
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\nRun Code Online (Sandbox Code Playgroud)\n有没有办法获得data.frameR 的典型点坐标?\n作为示例,我们可以使用来自此http://data.wvgis.wvu.edu/elevation/的数据。\n此外,有没有办法获得角度反射,data.frame激光雷达文件中的反射数量?
我需要创建一个运行给定时间(例如 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 个小时)
怎样才能做到呢?感谢您的帮助。
假设有这个函数sqrt(x^2)+0.9*sqrt(3.3-x^2)*sin(30*pi*x)这个函数生成一个形状的图heart
有没有办法使用 ggplot2 重现此函数来获得红色heart作为输出?感谢您的帮助。
假设这是我的数据示例
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) 假设我有一些值,那么它dat$x1= c(23L, 45L)。如何将这些值中的每一个随机乘以 -2 到 2 范围内的任何值?比如23*1,5或者23*-0,2等等。有什么简单的方法可以做到呢?谢谢。
有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)