1

我想从一个灰度图像中的 ImageJ 中的粒子分析中绘制一个 ppp 对象,该图像的大小为 imageSizeX、imageSizeY 的点对象 Particles$X、Particles$Y 在翻转的 Java 坐标中(y0 = 左上角)。我有一种可行的解决方案(自定义示例):

library(spatstat)

X <- ppp(Particles$X, Particles$Y, c(0, imageSizeX), c(0, imageSizeY))
plot(x = 0, y = 0, xlim = c(0, imageSizeX), ylim = c(imageSizeY, 0), type = "n", main = "Density", 
    asp = 1, axes = F, xlab = "X", ylab = "Y")
plot(density(X), xlim = c(1, imageSizeX), ylim = c(imageSizeY, 0), add = T)
plot(X, axes = TRUE, xlim = c(1, imageSizeX), ylim = c(imageSizeY, 0), add = T)
axis(1)
axis(2, las = 2)

这导致以下情节(省略了图例):

在此处输入图像描述

但是我需要创建一个带有翻转坐标的空绘图命令(ylim = c(imageSizeY, 0)),然后必须添加 spatstat 图。

如果我尝试绘制:

library(spatstat)
X <- ppp(Particles$X, Particles$Y, c(0, imageSizeX), c(0, imageSizeY))
plot(density(X), xlim = c(1, imageSizeX), ylim = c(imageSizeY, 0))
plot(X, axes = TRUE, xlim = c(1, imageSizeX), ylim = c(imageSizeY, 0), add = T)
axis(1)
axis(2, las = 2)

坐标未翻转(ylim = c(imageSizeY,0)):

在此处输入图像描述

有没有办法在没有第一个定义绘图命令的情况下翻转 spatstat 中的绘图坐标?

4

1 回答 1

0

这是plot.ppp. 我已在 GitHub 上的 spatstat 存储库中的当前开发版本中修复spatstat1.46-1.010

spatstat在CRAN ( )的当前版本中1.46-1,windows 的 plot 方法确实可以按照您的意愿进行plot.owin识别,而不能。所以如果 X 是一个点模式,你可以做xlim, ylimplot.ppp

W <- Window(X)
plot(W, xlim=rev(W$xrange), main="The title")
plot(X, add=TRUE, ...)
于 2016-07-16T05:54:03.937 回答