Package graphicsのジェネリック関数を使用して SpatialPoints オブジェクトをプロットしようとしています。プロットの範囲は、最も外側のデータポイントのサイズと正確に一致する必要があります。以下は、デルタ y に対するデルタ x の比率 (最も外側の値) が 1 である簡単な例です。
library(sp)
xydata <- as.matrix(cbind(c(1,1,2,2,1.5),c(1,2,2,1,1.5)))
xy_sp = SpatialPoints(xydata)
bbox(xy_sp)
# min and max for x and y coordinates are identical, i.e. apect ratio = 1
png("output.png", width = 500, height = 500, units = "px", type = "cairo" )
par(mar=c(0,0,0,0), xaxs = "i")
plot(xy_sp, xlim = c(1,2), ylim = c(1,2), pch = "X")
# the following lines should NOT be visible
abline(a = 1, b = 0); abline(a = 2, b = 0); abline(v = 1);abline(v = 2)
dev.off()
私の理解では、上記のスクリプトの結果のプロットは、1 つのポイント (座標 1.5,1.5) のみが表示され、4 つの「ablines」はすべて非表示になるグラフィックを生成する必要があります。そうでないのはなぜですか？
