问题标签 [spatstat]
For questions regarding programming in ECMAScript (JavaScript/JS) and its various dialects/implementations (excluding ActionScript). Note JavaScript is NOT the same as Java! Please include all relevant tags on your question; e.g., [node.js], [jquery], [json], [reactjs], [angular], [ember.js], [vue.js], [typescript], [svelte], etc.
spatstat - 从包 spatstat 将矩阵转换为 im 对象 - 行是相反的
第一次来这里的用户。我正在将矩阵从包 spatstat 转换为 im 对象。但是,矩阵 mat[1,1] 的左上角值写入图像的左下角,例如 im 函数从上到下读取矩阵的行,但从底部开始将它们写入 im 中。试图以相反的顺序将 im() 的 yrow 参数指定为 seq() 但结果仍然相同。建议如何解决这个问题?
r - 将图像转换为具有R中坐标和像素值列的数据框
我想将遥感图像转换data.frame
为 xy 坐标的两列和像素值的第三列。这是 的默认数据输入
ppmlasso
。
原始图像作为没有先验坐标的矩阵导入到 R 中。最小和最大坐标和分辨率在文本文件中给出,因此我使用它们im
为spatstat
. 这是一个简单的例子,可以说seq
十进制坐标,例如图像的分辨率为 10 m。
到目前为止,RS 数据的两种不同表示形式 -matrix
和im
. 现在我想要第三种表示形式 -data.frame
两列用于坐标对,第三列用于相应的 RS 值。直觉上我是这样做的:
去的原因raster
是使用coordinates
调用。然而结果是错误的 - 坐标对和矩阵中的值不匹配,因为coordinates
从顶行开始沿整行调用顺序,而as.vector
按列调用值。
建议?
r - R-spatstat im 方形像素?(用于导出到 Arc 栅格)
我在将密度 im 从 spatstat 导出为 ArcGIS 可读的文件格式时遇到问题。这是我的代码
一切正常。但是当我尝试导出到 ascii 光栅文件时
我得到这个问题
我使用的数据是针对加利福尼亚州的,因此纵横比不是 1。那么......我如何使密度 im 具有方形像素?
r - 安装 spatstaton R studio
我正在尝试在 R Studio 中上传 R 包“spatstat”,并收到以下警告消息。任何想法我做错了什么?
install.packages("spatst") 将软件包安装到“C:/Users/Mallarie/Documents/R/win-library/3.1”(因为未指定“lib”) install.packages 中的警告:软件包“spatst”不可用(作为 R 版本 3.1.3 的二进制包)
谢谢!
process - 点过程中标记和协变量之间的关系是什么
我对点过程中的标记和协变量感到困惑。我正在尝试通过使用 spatstat 在 R 中创建一个带有少量协变量的标记点模式模型,但我不确定标记和协变量之间的关系。有人可以帮我吗?
谢谢。
---- update 我有一个关于给定位置的定居点人口的给定点模式,并且很少有诸如土壤肥力、年降雨量等协变量。我想描述定居点的分布。
python - rpy2 convert numpy array into ppp (point pattern) spatstat
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 - 在 R 的同一图中用两个样本模拟雷区
我试图通过在同一个地块中绘制两个泊松分布的样本来模拟雷区,一个比另一个具有更高的强度和更小的面积。这是雷区,另一个只是被视为点的噪音(石头、洞、金属)。我无法让 R 在轴上绘制具有相同单位的点。无论我做什么,这些点都会覆盖整个情节,即使我只希望 X 点覆盖四分之一的情节。我的 R 代码如下:
如果你能帮助我,请告诉我。
r - eps 在密度.ppp() 中如何使用?
我正在使用 density.ppp() 来获取一系列坐标的强度函数。数据是投影的,以米为单位。我试图估计点数/公里,所以在查看文档后,我使用 eps = 1000 将输出图像的分辨率设置为 1km2。
但是 eps 是否仅指估计值的位置?或者强度函数估计的面积单位?
我是否仍需要将值重新缩放 1000x1000 才能获得 km2?
python - Spatstat:给定一个二维点列表,如何将它们连接成一个多边形,并进一步使其成为研究区域?
请允许我从一个最简单的任务开始这个问题:如果我有四个点是一个矩形的顶点,存储在一个 4x2 矩阵中,我怎样才能把它变成一个矩形窗口?(请不要使用任何特定于绘制矩形的特殊命令,因为矩形被升起只是为了表示一般类的规则几何对象)
让事情变得更复杂,假设我有一个 nx2 矩阵,我怎样才能连接所有的 n 个点,使它成为一个多边形?请注意,对象不一定是凸的。我认为主要困难在于,R 怎么知道哪个点应该与哪个点连接?
我问的原因是我正在对一条鱼进行一些图像处理,我设法通过在 python 中使用 opencv 找到轮廓来获得鱼的身体线条,并将其输出为 nx2 csv 文件。当我将 csv 文件读入 R 并尝试使用 sp 包中的 SpatialPolygnos 将其转换为多边形时,发生了一些非常意外的行为;中间某处似乎有一个中断,多边形被切成两半,即多边形的边界没有连接。无论如何我可以解决这个问题吗?
谢谢你。
编辑:有人好心地指出,这可能是另一个问题的重复:在 R 中绘制多边形。但是,该问题的解决方案依赖于绘制的形状是凸的,因此按天使排序是有意义的;然而,这里的形状不一定是凸的,它不会起作用。
r - 使用 Matlab 和 R 比较 Ripley 的 K 函数的估计值
我正在使用以下 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 值吗?
谢谢