小编Mar*_*box的帖子

我希望有经验的人可以帮助我们如何从xyz数据中准备形状文件.虽然没有提供创建形状文件的前面步骤,但可以在这里看到彗星Churyumov-Gerasimenko 的精心准备的数据集的一个很好的例子.

我试图更好地理解如何将曲面应用于给定的XYZ坐标集.使用笛卡尔坐标是直接使用R包"rgl",但是环绕的形状看起来更难.我找到了R包geometry,它提供了QHULL函数的接口.我尝试使用它来计算Delaunay三角剖面,然后我可以用它绘制rgl.我无法弄清楚与该功能相关的一些选项delaunayn,可能控制这些方面的最大距离.我希望这里有人可能对xyz数据改进表面结构有一些想法.

使用"Stanford bunnny"数据集的示例:

library(onion)
library(rgl)
library(geometry)
data(bunny)

#XYZ point plot
open3d()
points3d(bunny, col=8, size=0.1)
#rgl.snapshot("3d_bunny_points.png")

#Facets following Delaunay triangulation
tc.bunny <- delaunayn(bunny)
open3d()
tetramesh(tc.bunny, bunny, alpha=0.25, col=8)
#rgl.snapshot("3d_bunny_facets.png")

这个答案让我相信Qhull的R实现可能存在问题.此外,我现在尝试了各种设置(例如delaunayn(bunny, options="Qt")),效果甚微.Qhull选项在此处列出

编辑:

这是一个额外的(更简单的)球体示例.即使在这里,小平面的计算并不总能找到最近的相邻顶点(如果你旋转球,你会看到一些小平面穿过内部).

library(rgl)
library(geometry)
set.seed(1)
n <- 10
rho <- 1
theta <- seq(0, 2*pi,, n) # azimuthal coordinate running from 0 to 2*pi 
phi <- seq(0, pi,, n) # polar …

这实际上应该非常简单,但我很难找到解决这个问题的方法.

我在R中有两个非常简单的数字向量.我只是想用它们绘制直方图.但是我希望它们在同一个图表上.棘手的部分是R默认重叠这两个直方图.我希望这些垃圾箱可以简单地并排放置,这样我就可以更好地直观地显示数据.

基本上这就是我想要做的

我对R和统计计算语言一般都是新手,所以如果你能回答我令人沮丧的问题,我将不胜感激.

鉴于对可重复科学的压力,我想知道我最近的工作是否保证在出版物中包含示例代码.我使用的数据集非常大,因此发布这些必要性是没有意义的 - 但是,我在R中应用的统计方法通常不为我的观众所知(尽管我认为它们应该是).

我在我的分析中使用经验正交函数分析(EOF)和广义加性模型(GAM).特别是GAM广泛用于生态学研究,但在物理科学中却不那么重要 - 我的工作跨越两个学科.

我肯定会参考我使用的R包,并且审阅者/读者自己寻找这些引用(和包含的示例)并不是很困难.所以,我的问题是,什么情况最适合在出版物中包含可重现的代码？

我正在尝试直接从URL地址加载.jpeg图像.我想知道是否有一种基本的方法来通过使用url连接来做到这一点.

我第一次尝试:

require(biOps)
con <- url("http://upload.wikimedia.org/wikipedia/commons/9/95/Apollonian_spheres.jpg")
pic <- readJpeg(con)
#Error in readJpeg(con) : Cannot open file.

另一个问题似乎是相同的,但对于.png文件.我试图适应.jpeg,但也有错误.

require(biOps)
require(RCurl)
myurl <- "http://upload.wikimedia.org/wikipedia/commons/9/95/Apollonian_spheres.jpg"
pic <-  readJpeg(getURLContent(myurl))
#Error in readJpeg(getURLContent(myurl)) : Cannot open file.

任何帮助将不胜感激!

有没有人知道一种方法来转动contourLines多边形的输出,以绘制为填充轮廓,如同filled.contours.是否有关于如何绘制多边形以便查看所有可用级别的顺序？以下是一段不起作用的示例代码段:

#typical plot
filled.contour(volcano, color.palette = terrain.colors)

#try
cont <- contourLines(volcano)
fun <- function(x) x$level
LEVS <- sort(unique(unlist(lapply(cont, fun))))
COLS <- terrain.colors(length(LEVS))
contour(volcano)
for(i in seq(cont)){
    COLNUM <- match(cont[[i]]$level, LEVS)
    polygon(cont[[i]], col=COLS[COLNUM], border="NA")
}
contour(volcano, add=TRUE)

我正在尝试计算一个类似Bezier的样条曲线,该曲线通过一系列xy坐标.一个例子就像cscvnMatlab中的函数的以下输出(示例链接):

我相信(不再维护)grid包用于执行此操作(grid.xspline功能？),但我无法安装包的存档版本,并且没有找到任何与我想要的完全一致的示例.

该bezier软件包看起来很有前景,但速度非常慢,我也无法做到这一点:

library(bezier)

set.seed(1)
n <- 10
x <- runif(n)
y <- runif(n)
p <- cbind(x,y)
xlim <- c(min(x) - 0.1*diff(range(x)), c(max(x) + 0.1*diff(range(x))))
ylim <- c(min(y) - 0.1*diff(range(y)), c(max(y) + 0.1*diff(range(y))))
plot(p, xlim=xlim, ylim=ylim)
text(p, labels=seq(n), pos=3)

bp <- pointsOnBezier(cbind(x,y), n=100)
lines(bp$points)
arrows(bp$points[nrow(bp$points)-1,1], bp$points[nrow(bp$points)-1,2],
  bp$points[nrow(bp$points),1], bp$points[nrow(bp$points),2]
)

如您所见,除了结束值之外,它不会通过任何点.

我非常感谢这里的一些指导!

有谁知道如何计算R中2个或更多个多边形之间的共同面积？我希望这种计算的输出是该重叠区域的新多边形的坐标.干杯

我正在尝试将以下日期转换为R中的时间类.出于某种原因,当未包含日值时,我将返回NAs.我希望能够按原样处理这个字符向量,而不必粘贴任意日值.谢谢你的帮助.

TS <- c("2004-12", "2005-01", "2005-02", "2005-03", "2005-04", "2005-05", 
"2005-06", "2005-07", "2005-08", "2005-09", "2005-10", "2005-11", 
"2005-12", "2006-01", "2006-02", "2006-03", "2006-04", "2006-05", 
"2006-06", "2006-07", "2006-08")
TSd <- paste(TS, "01", sep="-")

#doesn't work
as.Date(TS, format="%Y-%m")
as.POSIXlt(TS, format="%Y-%m")

#works
as.Date(TSd, format="%Y-%m-%d")
as.POSIXlt(TSd, format="%Y-%m-%d")

我试图在R包中包含一个(有些)大数据集.在Rstudio检查期间我一直收到警告说我可以通过压缩来节省空间:

* checking data for ASCII and uncompressed saves ... WARNING

  Note: significantly better compression could be obtained
        by using R CMD build --resave-data
          old_size new_size compress
  slp.rda    499Kb    310Kb    bzip2
  sst.rda    1.3Mb    977Kb       xz

我已经尝试添加-- resave-data到RStudio的"配置Buid工具"无效.

我看到我最近的包构建不再创建手册,即使我没有在构建设置中添加"--no-manual".之前我没有遇到过这个问题 - 我最近更新了devtools和roxygen2.

我错过了什么吗？

以下是我设置项目选项的方法: