问题标签 [spatstat]

第一次来这里的用户。我正在将矩阵从包 spatstat 转换为 im 对象。但是，矩阵 mat[1,1] 的左上角值写入图像的左下角，例如 im 函数从上到下读取矩阵的行，但从底部开始将它们写入 im 中。试图以相反的顺序将 im() 的 yrow 参数指定为 seq() 但结果仍然相同。建议如何解决这个问题？

我想将遥感图像转换data.frame为 xy 坐标的两列和像素值的第三列。这是 的默认数据输入 ppmlasso。

原始图像作为没有先验坐标的矩阵导入到 R 中。最小和最大坐标和分辨率在文本文件中给出，因此我使用它们im为spatstat. 这是一个简单的例子，可以说seq十进制坐标，例如图像的分辨率为 10 m。

到目前为止，RS 数据的两种不同表示形式 -matrix和im. 现在我想要第三种表示形式 -data.frame两列用于坐标对，第三列用于相应的 RS 值。直觉上我是这样做的：

去的原因raster是使用coordinates调用。然而结果是错误的 - 坐标对和矩阵中的值不匹配，因为coordinates从顶行开始沿整行调用顺序，而as.vector按列调用值。

建议？

我在将密度 im 从 spatstat 导出为 ArcGIS 可读的文件格式时遇到问题。这是我的代码

一切正常。但是当我尝试导出到 ascii 光栅文件时

我得到这个问题

我使用的数据是针对加利福尼亚州的，因此纵横比不是 1。那么......我如何使密度 im 具有方形像素？

我正在尝试在 R Studio 中上传 R 包“spatstat”，并收到以下警告消息。任何想法我做错了什么？

install.packages("spatst") 将软件包安装到“C:/Users/Mallarie/Documents/R/win-library/3.1”（因为未指定“lib”） install.packages 中的警告：软件包“spatst”不可用（作为 R 版本 3.1.3 的二进制包）

谢谢！

我对点过程中的标记和协变量感到困惑。我正在尝试通过使用 spatstat 在 R 中创建一个带有少量协变量的标记点模式模型，但我不确定标记和协变量之间的关系。有人可以帮我吗？

谢谢。

---- update 我有一个关于给定位置的定居点人口的给定点模式，并且很少有诸如土壤肥力、年降雨量等协变量。我想描述定居点的分布。

I have a numpy array as given, corresponding to x,y coordinates. I am searching for a way to convert this into ppp point structure to perform spatial analysis using spatstat tool in R via rpy2.

Ideally I would like to convert this into ppp using rpy2. I also tried to write numpy into csv and read that into R. However it still doesn't recognize it as ppp structure.

我试图通过在同一个地块中绘制两个泊松分布的样本来模拟雷区，一个比另一个具有更高的强度和更小的面积。这是雷区，另一个只是被视为点的噪音（石头、洞、金属）。我无法让 R 在轴上绘制具有相同单位的点。无论我做什么，这些点都会覆盖整个情节，即使我只希望 X 点覆盖四分之一的情节。我的 R 代码如下：

如果你能帮助我，请告诉我。

我正在使用 density.ppp() 来获取一系列坐标的强度函数。数据是投影的，以米为单位。我试图估计点数/公里，所以在查看文档后，我使用 eps = 1000 将输出图像的分辨率设置为 1km2。

但是 eps 是否仅指估计值的位置？或者强度函数估计的面积单位？

我是否仍需要将值重新缩放 1000x1000 才能获得 km2？

请允许我从一个最简单的任务开始这个问题：如果我有四个点是一个矩形的顶点，存储在一个 4x2 矩阵中，我怎样才能把它变成一个矩形窗口？（请不要使用任何特定于绘制矩形的特殊命令，因为矩形被升起只是为了表示一般类的规则几何对象）

让事情变得更复杂，假设我有一个 nx2 矩阵，我怎样才能连接所有的 n 个点，使它成为一个多边形？请注意，对象不一定是凸的。我认为主要困难在于，R 怎么知道哪个点应该与哪个点连接？

我问的原因是我正在对一条鱼进行一些图像处理，我设法通过在 python 中使用 opencv 找到轮廓来获得鱼的身体线条，并将其输出为 nx2 csv 文件。当我将 csv 文件读入 R 并尝试使用 sp 包中的 SpatialPolygnos 将其转换为多边形时，发生了一些非常意外的行为；中间某处似乎有一个中断，多边形被切成两半，即多边形的边界没有连接。无论如何我可以解决这个问题吗？

谢谢你。

编辑：有人好心地指出，这可能是另一个问题的重复：在 R 中绘制多边形。但是，该问题的解决方案依赖于绘制的形状是凸的，因此按天使排序是有意义的；然而，这里的形状不一定是凸的，它不会起作用。

我正在使用以下 Matlab 代码来估计 Ripley 的 K 函数。

“RipleysK”功能可在以下网址找到：http: //www.colorado.edu/geography/class_homepages/geog_4023_s07/labs/lab10/RipleysK.m

相比之下，我使用的是以下 R 脚本。

Matlab 为指定的 r 值（即 dist）估计 K，而 R 脚本不估计（估计直到 r = 12.5）。

有人可以评论吗？哪一个是正确的？我们可以在 R 脚本中指定 r 值吗？

谢谢