Hello dear R-Users
I have a little problem, and i don't know what is happening there. I have a
txt file containing for 480 points xCoord,yCoord,pm10, and time. For those
of you who dont know, pm10 are little air pollutants. What i wand to do is a
simple plot of these pm10 values with theis x and y coordinate. I use
following code fragment :

library(sp)
library(lattice)
library(rgdal)

walk1 <-
read.table("/home/hacke/Diplom/Unidaten/ComputerData/Walk1.txt",sep=',',
header = TRUE)
walk2 <-
read.table("/home/hacke/Diplom/Unidaten/ComputerData/Walk2.txt",sep=',',
header = TRUE)
walk3 <-
read.table("/home/hacke/Diplom/Unidaten/ComputerData/Walk3.txt",sep=',',
header = TRUE)

coordinates(walk1) <- ~xCoord+yCoord
pic <- spplot((walk1,zcol = "pm10",col.regions =
c("green","lightgreen","yellow","orange","red" ),main = "Concentration of
PM10 Computed Walk1",xlab = "Longitude",ylab = "Latitude"))
plot(pic)

As I am not doing this the first time, everything should be fine. But the
plot is a straight line from the lower left to the upper right corner.
I am 100 % sure not all of the coordinates are on a straight line, prove is
here :

> x <-
matrix(c(walk1$xCoord[1],walk1$xCoord[2],walk1$xCoord[50],walk1$yCoord[1],walk1$yCoord[2],walk1$yCoord[50],1,1,1),nrow
=3,ncol=3)
+ x
+ det(x)
        [,1]    [,2] [,3]
[1,] 3405540 5758268    1
[2,] 3405790 5758518    1
[3,] 3407790 5760518    1
[1] 6.043561e-07

There are also only about 100 points plotted, not all the 480.

Coordinates are GK3 coordinates ( i think this is only common in Germany,
but should be okay in my opinion ) In other plots where all worked fine i
think I used UTM coordinates. Anybody knows perhaps what is going wrong?

Best wishes
Jan Hackenberg

	[[alternative HTML version deleted]]

