2016-07-15 111 views
1

我想從ImageJ中的一個粒度分析中繪製一個pSize對象,其大小爲sizeSize imageSizeX,image objectsY粒子$ X,粒子$ Y翻轉的Java座標(y0 = top-剩下)。 我有一種工作液(例如用定製):在翻轉的Java座標中繪製

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) 

其結果如下圖(其省略了傳說):

enter image description here

不過,我需要創建一個空使用翻轉的座標繪製命令(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)):

enter image description here

有沒有辦法翻轉在沒有第一個定義繪圖命令的情況下在spatstat中繪製座標?

回答

0

這是plot.ppp中的一個錯誤。我已經在spatstat目前的開發版本,版本1.46-1.010可固定從the spatstat repository on GitHub

spatstat上CRAN當前版本(1.46-1)的Windows,plot.owin情節方法並不確認xlim, ylim,你想,而plot.ppp纔不是。所以如果X是點模式,你可以做

W <- Window(X) 
plot(W, xlim=rev(W$xrange), main="The title") 
plot(X, add=TRUE, ...) 
+0

好。感謝您的快速修復和答覆! – Marcel

+0

只是在窗口版本中翻轉y座標的小錯字: plot(W,ylim = rev(W $ yrange),main =「The title」) – Marcel