-
Notifications
You must be signed in to change notification settings - Fork 0
/
Copy patheurocentric.r
36 lines (30 loc) · 1.14 KB
/
eurocentric.r
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
convert.to.eurocentric <- function(raster.in) {
rl <- crop(raster.in,
extent(c(xmin(raster.in),xmax(raster.in)-180,ymin(raster.in),
ymax(raster.in))))
# set extent from rr to rleft:
rr <- (crop(raster.in,
extent(c(xmax(raster.in)-180,xmax(raster.in),ymin(raster.in),
ymax(raster.in)))))
extent(rr) <- extent(c(xmin(rr)-360,xmax(rr)-360,ymin(raster.in),
ymax(raster.in)))
r.out = merge(rr,rl)
names(r.out) = names(raster.in)
return(r.out)
}
# zum einlesen:
# source("convert.to.eurocentric.R")
convert.to.pacificcentric <- function(raster.in) {
rl <- crop(raster.in,
extent(c(xmin(raster.in)+180,xmax(raster.in),ymin(raster.in),
ymax(raster.in))))
# set extent from rr to rleft:
rr <- (crop(raster.in,
extent(c(xmin(raster.in),xmax(raster.in)-180,ymin(raster.in),
ymax(raster.in)))))
extent(rr) <- extent(c(xmin(rr)+360,xmax(rr)+360,ymin(raster.in),
ymax(raster.in)))
r.out = merge(rr,rl)
names(r.out) = names(raster.in)
return(r.out)
}