RでGIS 4

いよいよ課題。 コドラート法とやらをやってみる。

ポイントが重複してるのは受け付けないんだそーだ。
LibreOfficeで座標の重複している行を削除して試行。

> StrokePoints <- read.table("tanigawa.csv", sep=";", header=TRUE)

> X <- ppp(StrokePoints$X, StrokePoints$Y, c(138.9,139), c(36.8,36.9))
> plot(X)

> quadratcount(X, nx=5, ny=5)
               x
y               [138.9,138.92] (138.92,138.94] (138.94,138.96] (138.96,138.98] (138.98,139]
  (36.88,36.9]               2               0               0               0            0
  (36.86,36.88]              0               0               0               0            0
  (36.84,36.86]              0              18               4               1            0
  (36.82,36.84]              0               7               6               0            0
  [36.8,36.82]               0               0               1               0            0
> plot(X)
> plot(quadratcount(X, nx=5, ny=5), cex=1.2, add=TRUE)

> plot(X, cex=0.2, main="")

> plot(quadrat.test(X, nx=5, ny=5),  col=1, cex=1.5, add=TRUE)

 警告メッセージ：
Some expected counts are small; chi^2 approximation may be inaccurate

なんで警告でるのかわからんけど。ここが僕の限界みたいです。