Project

General

Profile

Download (5.89 KB) Statistics
| Branch: | Revision:
1 77ed568c Jim Regetz
library(rgdal)
2
aster.dir <- "DEM/asterGdem"
3
4
#Check that path=correct:
5
a<- GDALinfo("DEM/asterGdem/ASTGTM_N59E069_dem.tif", 
6
    silent=TRUE,returnRAT=TRUE)
7
8
# for a given file, return TRUE if file name matches the internally
9
# specified lower left corner, otherwise return string with those
10
# coordinates
11
check <- function(tilename, path=".", silent=TRUE) {
12
    # build expected filename
13
    origin <- round(GDALinfo(file.path(path, tilename),
14
        silent=silent)[c("ll.x", "ll.y")])
15
    ly <- origin["ll.y"]
16
    y <- sprintf("%s%02d", if (ly>=0) "N" else "S", abs(ly))
17
    lx <- origin["ll.x"]
18
    x <- sprintf("%s%03d", if (lx>=0) "E" else "W", abs(lx))
19
    expected.name <- paste("ASTGTM_", y, x, "_dem.tif", sep="")
20
    # compare to actual filename
21
    if (tilename==expected.name) {
22
        TRUE
23
    } else {
24
        paste(y, x)
25
    }
26
}
27
28
aster.tiles <- list.files(aster.dir, pattern="^ASTGTM.*_dem.tif$") 
29
titlecheck <- sapply(aster.tiles, check, path=aster.dir)
30
31
BadAster_df<- data.frame(expected=titlecheck[titlecheck!="TRUE"])
32
33
# ASTGTM_N59E069_dem.tif  N63 E109
34
#ASTGTM_N63E113_dem.tif  N69 E107
35
# ASTGTM_N63E117_dem.tif	N69 E113
36
#ASTGTM_N64E098_dem.tif	N70 E117
37
#ASTGTM_N65E104_dem.tif	N73 E084
38
#	ASTGTM_N65E111_dem.tif	N73 E098
39
#	ASTGTM_N65E117_dem.tif	N66 E130
40
41
KnownBadcheck <- function(badtilename, path="DEM/asterGdem/incorrectTilesJuly2011", silent=TRUE) {
42
    # build expected filename
43
    origin <- round(GDALinfo(file.path(path, badtilename),
44
        silent=silent)[c("ll.x", "ll.y")])
45
    ly <- origin["ll.y"]
46
    y <- sprintf("%s%02d", if (ly>=0) "N" else "S", abs(ly))
47
    lx <- origin["ll.x"]
48
    x <- sprintf("%s%03d", if (lx>=0) "E" else "W", abs(lx))
49
    Badexpected.name <- paste("ASTGTM_", y, x, "_dem.tif", sep="")
50
    # compare to actual filename
51
    if (badtilename==Badexpected.name) {
52
        TRUE
53
    } else {
54
        paste(y, x)
55
    }
56
}
57
58
bad.tiles <- list.files("DEM/asterGdem/incorrectTilesJuly2011", pattern="^ASTGTM.*_dem.tif$") 
59
badtitlecheck <- sapply(bad.tiles, KnownBadcheck, path="DEM/asterGdem/incorrectTilesJuly2011")
60
61
BadTitle_df<- data.frame(expected=badtitlecheck[badtitlecheck!="TRUE"])
62
63
#----------------------------------------------------------------------------------------------------------------------
64
#Check num.tif files
65
b<- GDALinfo("DEM/asterGdem/ASTGTM_N49E000_num.tif", 
66
silent=TRUE,returnRAT=TRUE)
67
b
68
69
Numcheck <- function(Numtilename, path=".", silent=TRUE) {
70
    # build expected filename
71
    origin <- round(GDALinfo(file.path(path, Numtilename),
72
        silent=silent)[c("ll.x", "ll.y")])
73
    ly <- origin["ll.y"]
74
    y <- sprintf("%s%02d", if (ly>=0) "N" else "S", abs(ly))
75
    lx <- origin["ll.x"]
76
    x <- sprintf("%s%03d", if (lx>=0) "E" else "W", abs(lx))
77
    Numexpected.name <- paste("ASTGTM_", y, x, "_num.tif", sep="")
78
    # compare to actual filename
79
    if (Numtilename==Numexpected.name) {
80
        TRUE
81
    } else {
82
        paste(y, x)
83
    }
84
}
85
86
Num.tiles <- list.files("DEM/asterGdem/NewTiles", pattern="^ASTGTM.*_num.tif$") 
87
Numtitlecheck <- sapply(Num.tiles, Numcheck, path="DEM/asterGdem/NewTiles")
88
89
NumTitle_df<- data.frame(expected=Numtitlecheck[Numtitlecheck!="TRUE"])
90
91
#Located one bad num.tif file
92
N49E007_NumRead<- GDALinfo("DEM/asterGdem/ASTGTM_N49E007_num.tif", 
93
silent=TRUE,returnRAT=TRUE)
94
95
GDALinfo("DEM/asterGdem/ASTGTM_N59E069_dem.tif", 
96
silent=TRUE,returnRAT=TRUE)
97
98
#--------------------------------------------------------------------------------
99
#Check new tiles uploaded to replace faulty tiles
100
101
Newcheck <- function(Newtilename, path=".", silent=TRUE) {
102
    # build expected filename
103
    origin <- round(GDALinfo(file.path(path, Newtilename),
104
        silent=silent)[c("ll.x", "ll.y")])
105
    ly <- origin["ll.y"]
106
    y <- sprintf("%s%02d", if (ly>=0) "N" else "S", abs(ly))
107
    lx <- origin["ll.x"]
108
    x <- sprintf("%s%03d", if (lx>=0) "E" else "W", abs(lx))
109
    Newexpected.name <- paste("ASTGTM_", y, x, "_dem.tif", sep="")
110
    # compare to actual filename
111
    if (Newtilename==Newexpected.name) {
112
        TRUE
113
    } else {
114
        paste(y, x)
115
    }
116
}
117
118
New.tiles <- list.files("DEM/asterGdem/NewTiles", pattern="^ASTGTM.*_dem.tif$") 
119
Newtitlecheck <- sapply(New.tiles, Newcheck, path="DEM/asterGdem/NewTiles")
120
121
NewAster_df<- data.frame(Newexpected=Newtitlecheck[Newtitlecheck!="TRUE"])
122
123
#Yeah, no mistakes in first 4!
124
125
g<- list.files("DEM/asterGdem/N59to81_E60to99", pattern="^ASTGTM.*_dem.tif$")
126
h<- list.files("DEM/asterGdem/N59to81_E60to99/Non_DEM's", pattern="^ASTGTM.*_num.tif$")
127
128
GDALinfo("DEM/asterGdem/N59to81_E60to99/ASTGTM_N64E098_dem.tif", 
129
silent=TRUE,returnRAT=TRUE)
130
131
#--------------------------------------------------------------------------
132
#new folder check
133
nfcheck <- function(nftilename, path=".", silent=TRUE) {
134
    # build expected filename
135
    origin <- round(GDALinfo(file.path(path, nftilename),
136
        silent=silent)[c("ll.x", "ll.y")])
137
    ly <- origin["ll.y"]
138
    y <- sprintf("%s%02d", if (ly>=0) "N" else "S", abs(ly))
139
    lx <- origin["ll.x"]
140
    x <- sprintf("%s%03d", if (lx>=0) "E" else "W", abs(lx))
141
    nfexpected.name <- paste("ASTGTM_", y, x, "_dem.tif", sep="")
142
    # compare to actual filename
143
    if (nftilename==nfexpected.name) {
144
        TRUE
145
    } else {
146
        paste(y, x)
147
    }
148
}
149
150
nf.tiles <- list.files("DEM/asterGdem/N59to81_E60to99", pattern="^ASTGTM.*_dem.tif$") 
151
nftitlecheck <- sapply(nf.tiles, nfcheck, path="DEM/asterGdem/N59to81_E60to99")
152
153
nfAster_df<- data.frame(nfexpected=nftitlecheck[nftitlecheck!="TRUE"])
154
155
156
Numcheck <- function(Numtilename, path=".", silent=TRUE) {
157
    # build expected filename
158
    origin <- round(GDALinfo(file.path(path, Numtilename),
159
        silent=silent)[c("ll.x", "ll.y")])
160
    ly <- origin["ll.y"]
161
    y <- sprintf("%s%02d", if (ly>=0) "N" else "S", abs(ly))
162
    lx <- origin["ll.x"]
163
    x <- sprintf("%s%03d", if (lx>=0) "E" else "W", abs(lx))
164
    Numexpected.name <- paste("ASTGTM_", y, x, "_num.tif", sep="")
165
    # compare to actual filename
166
    if (Numtilename==Numexpected.name) {
167
        TRUE
168
    } else {
169
        paste(y, x)
170
    }
171
}