1 |
7526fb1c
|
Jim Regetz
|
# 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
|