Introduction to spatial analysis

library(sp)
library(gstat)
data(meuse)
?meuse
class(meuse)

[1] "data.frame"

names(meuse)

 [1] "x"       "y"       "cadmium" "copper"  "lead"    "zinc"    "elev"    "dist"   
 [9] "om"      "ffreq"   "soil"    "lime"    "landuse" "dist.m" 

Transform meuse to spatial data frame

#coordinates(meuse)=~x+y
head(coordinates(meuse))

       x      y
1 181072 333611
2 181025 333558
3 181165 333537
4 181298 333484
5 181307 333330
6 181390 333260

class(meuse)

[1] "SpatialPointsDataFrame"
attr(,"package")
[1] "sp"

summary(meuse)

Object of class SpatialPointsDataFrame
Coordinates:
     min    max
x 178605 181390
y 329714 333611
Is projected: NA 
proj4string : [NA]
Number of points: 155
Data attributes:
    cadmium           copper            lead            zinc             elev       
 Min.   : 0.200   Min.   : 14.00   Min.   : 37.0   Min.   : 113.0   Min.   : 5.180  
 1st Qu.: 0.800   1st Qu.: 23.00   1st Qu.: 72.5   1st Qu.: 198.0   1st Qu.: 7.546  
 Median : 2.100   Median : 31.00   Median :123.0   Median : 326.0   Median : 8.180  
 Mean   : 3.246   Mean   : 40.32   Mean   :153.4   Mean   : 469.7   Mean   : 8.165  
 3rd Qu.: 3.850   3rd Qu.: 49.50   3rd Qu.:207.0   3rd Qu.: 674.5   3rd Qu.: 8.955  
 Max.   :18.100   Max.   :128.00   Max.   :654.0   Max.   :1839.0   Max.   :10.520  
                                                                                    
      dist               om         ffreq  soil   lime       landuse       dist.m      
 Min.   :0.00000   Min.   : 1.000   1:84   1:97   0:111   W      :50   Min.   :  10.0  
 1st Qu.:0.07569   1st Qu.: 5.300   2:48   2:46   1: 44   Ah     :39   1st Qu.:  80.0  
 Median :0.21184   Median : 6.900   3:23   3:12           Am     :22   Median : 270.0  
 Mean   :0.24002   Mean   : 7.478                         Fw     :10   Mean   : 290.3  
 3rd Qu.:0.36407   3rd Qu.: 9.000                         Ab     : 8   3rd Qu.: 450.0  
 Max.   :0.88039   Max.   :17.000                         (Other):25   Max.   :1000.0  
                   NA's   :2                              NA's   : 1                   

plot

plot(meuse)
title("Location map")

plot(meuse,pch=2,col=3)

add features

image(meuse.grid,col="lightgrey")
plot(meuse.pol,col="grey",add=TRUE)
plot(meuse, pch=1, add=TRUE,cex=0.7)
plot(meuse, axes=T, cex=0.3)

plot(meuse.pol, add=T)
plot(meuse.grid,col="lightgrey",add=T)

title("Location Map")