1
|
# R script for creating tables of negative values and its frequency for each
|
2
|
# for each nine region. For region 6 and 7, the list (table) were too long
|
3
|
# to show at once, resulting in separating first part of the table and
|
4
|
# the latter part of the table to print.
|
5
|
#
|
6
|
#
|
7
|
# 9-Dec- 2011
|
8
|
# Yuina Nunokawa
|
9
|
|
10
|
|
11
|
|
12
|
|
13
|
|
14
|
library(raster)
|
15
|
|
16
|
v1 <- values(raster("/data/project/organisms/DEM/Yuni/Data/aster2/aster2_e020e059_82N.tif"))
|
17
|
t1 <- as.data.frame(table(v1[v1 < 0]))
|
18
|
|
19
|
v2 <- values(raster("/data/project/organisms/DEM/Yuni/Data/aster2/aster2_e060e099_82N.tif"))
|
20
|
t2 <- as.data.frame(table(v2[v2 < 0]))
|
21
|
|
22
|
v3 <- values(raster("/data/project/organisms/DEM/Yuni/Data/aster2/aster2_e100e139_82N.tif"))
|
23
|
t3 <- as.data.frame(table(v3[v3 < 0]))
|
24
|
|
25
|
v4 <- values(raster("/data/project/organisms/DEM/Yuni/Data/aster2/aster2_e140e179_82N.tif"))
|
26
|
t4 <- as.data.frame(table(v4[v4 < 0]))
|
27
|
|
28
|
v5 <- values(raster("/data/project/organisms/DEM/Yuni/Data/aster2/aster2_w020e019_82N.tif"))
|
29
|
t5 <- as.data.frame(table(v5[v5 < 0]))
|
30
|
|
31
|
v6 <- values(raster("/data/project/organisms/DEM/Yuni/Data/aster2/aster2_w060w021_82N.tif"))
|
32
|
t6 <- as.data.frame(table(v6[v6 < 0]))
|
33
|
|
34
|
|
35
|
v7 <- values(raster("/data/project/organisms/DEM/Yuni/Data/aster2/aster2_w100w061_82N.tif"))
|
36
|
t7 <- as.data.frame(table(v7[v7 < 0]))
|
37
|
|
38
|
v8 <- values(raster("/data/project/organisms/DEM/Yuni/Data/aster2/aster2_w140w101_82N.tif"))
|
39
|
t8 <- as.data.frame(table(v8[v8 < 0]))
|
40
|
|
41
|
v9 <- values(raster("/data/project/organisms/DEM/Yuni/Data/aster2/aster2_w180w141_82N.tif"))
|
42
|
t9 <- as.data.frame(table(v9[v9 < 0]))
|
43
|
|
44
|
|
45
|
|
46
|
# print the table
|
47
|
|
48
|
t1
|
49
|
t2
|
50
|
t3
|
51
|
t4
|
52
|
t5
|
53
|
t6[1:100]
|
54
|
t6
|
55
|
head(t7, 41)
|
56
|
t7
|
57
|
t8
|
58
|
t9
|