dla_mateusza.R

library(sp)
library(rgdal)

## rgdal: version: 1.2-16, (SVN revision 701)
##  Geospatial Data Abstraction Library extensions to R successfully loaded
##  Loaded GDAL runtime: GDAL 2.2.1, released 2017/06/23
##  Path to GDAL shared files: /usr/share/gdal/2.2
##  GDAL binary built with GEOS: TRUE 
##  Loaded PROJ.4 runtime: Rel. 4.9.2, 08 September 2015, [PJ_VERSION: 492]
##  Path to PROJ.4 shared files: (autodetected)
##  Linking to sp version: 1.2-5

library(maptools)

## Checking rgeos availability: TRUE

wgs84 <-  CRS('+init=epsg:4326')

ramka <- data.frame(lon=rnorm(1000, 20, 2),lat=rnorm(1000, 50, 2), zjawiska=round((runif(1000)*5))) # jakies wspolrzedne

coordinates(ramka) <- ~lon+lat
proj4string(ramka) <- wgs84


library(raster)
X=seq(from=-10, to=40, by=0.25)
Y=rev(seq(from=35, to=75, by=0.25))

r <- raster(ncol=length(X), nrow=length(Y), xmn=X[1],
            xmx=X[length(X)], ymn=Y[length(Y)], ymx=Y[1]) ## stworzenie obiektu rastrowego
res(r) <- 0.25
proj4string(r) <- CRS('+init=epsg:4326')

tatam <- rasterize(ramka, r, field='zjawiska', fun=sum)
plot(tatam)
library(mapdata)

## Loading required package: maps

map('world',add=T)