taszar.R

library(raster)

## Loading required package: sp

library(magic)

## Loading required package: abind
## 
## Attaching package: 'magic'
## 
## The following object is masked from 'package:raster':
## 
##     shift

setwd("/home/bartosz/Pulpit/taszar/")
pliki=dir(pattern="TTIn.grd")


mapa=as.matrix(raster(pliki[1])) # wczytujemy sobie jakis pierwszy z brzegu raster

wynik=array(NA,dim = c(dim(mapa),24))

#rozpoczecie petli
for (petla in 1:length(pliki)){

  mapa_tmp=as.matrix(raster(pliki[petla]))

  wynik[,,petla]=mapa_tmp

}

dim(wynik)

## [1] 101 113  24

maksiu=apply(wynik,1:2,max)# liczenie np. maxa
image(maksiu)

# wydaje mi sie, ze pakiet raster obraca pliki przy wczytywaniu
# dlatego tez chyba trzeba bedzie uzyc funkcji fliplr lub flipud z paczki matlab