Project

General

Profile

Download (2.65 KB) Statistics
| Branch: | Revision:
1 a87c0b1d Jim Regetz
# Quick exploration and file listing-based comparison of the CGIAR SRTM
2
# tiles available in the following places:
3
#
4
#   - those downloaded by Ming
5
#       jupiter.nceas.ucsb.edu:~organisms/SRTM_90m_ASCII_v4.1/
6
#
7
#   - those downloaded by Rick
8
#       jupiter.nceas.ucsb.edu:~organisms/CgiarSrtmAll/5_5x5_ascii/
9
#
10
#   - those currently available at CGIAR (ASCII format)
11
#       ftp://srtm.csi.cgiar.org/SRTM_v41/SRTM_Data_ArcASCII/
12
#
13
# Note that in addition to the comparisons made below, I also noticed
14
# that a few timestamps on zips Ming downloaded are newer than the
15
# timestamps on the CGIAR ftp site. It appears that these are bad zips
16
# on the CGIAR site; in subsequent email exchange with Ming, he said he
17
# believes he replaced them by downloading the corresponding GeoTIFFs
18
# from CGIAR and converting them to ASCII.
19
#
20
# Jim Regetz
21
# NCEAS
22
# Created on 21-May-2011
23
24
# list of tiles Ming downloaded onto jupiter
25
cmd <- "cd ~organisms/SRTM_90m_ASCII_v4.1 && ls -l srtm_*zip"
26
files <- system(cmd, intern=TRUE)
27
jup.v4.1 <- with(read.table(textConnection(files)), data.frame(tile=V8,
28
    size=V5, date=as.Date(V6, format="%Y-%m-%d"), time=V7))
29
30
# list of tiles Rick downloaded onto jupiter
31
cmd <- "cd ~organisms/CgiarSrtmAll/5_5x5_ascii && ls -l srtm_*zip"
32
files <- system(cmd, intern=TRUE)
33
jup.5_5x5 <- with(read.table(textConnection(files)), data.frame(tile=V8,
34
    size=V5, date=as.Date(V6, format="%Y-%m-%d"), time=V7))
35
36
# current list of tiles at CGIAR ftp site
37
cmd <- "curl ftp://srtm.csi.cgiar.org/SRTM_v41/SRTM_Data_ArcASCII/"
38
files <- system(cmd, intern=TRUE)
39
cgiar <- with(read.table(textConnection(files)), data.frame(tile=V4,
40
    size=V3, date=as.Date(V1, format="%m-%d-%y"), time=V2))
41
42
# list CGIAR tiles not found in jup.5_5x5
43
setdiff(cgiar$tile, jup.5_5x5$tile)
44
##  "srtm_32_05.zip" "srtm_68_03.zip"
45
46
# list jup.5_5x5 tiles not found in CGIAR (none!)
47
setdiff(jup.5_5x5$tile, cgiar$tile)
48
##  character(0)
49
50
# list CGIAR tiles not found in jup.v4.1
51
setdiff(cgiar$tile, jup.v4.1$tile)
52
##  "srtm_65_01.zip" "srtm_65_05.zip" "srtm_65_06.zip" "srtm_65_07.zip"
53
##  "srtm_65_08.zip" "srtm_65_10.zip" "srtm_65_11.zip" "srtm_65_13.zip"
54
##  "srtm_65_15.zip" "srtm_65_16.zip"
55
56
# list jup.v4.1 tiles not found in CGIAR (none!)
57
setdiff(jup.v4.1$tile, cgiar$tile)
58
##  character(0)
59
60
# show that all common tiles between CGIAR and jup.v4.1 have the same file
61
# size
62
common <- merge(cgiar, jup.v4.1, by="tile")
63
identical(common$size.x, common$size.y)
64
##  [1] TRUE
65
66
# but note that three files in jup.v4.1 have different (newer, as it
67
# turns out) timestamps
68
identical(common$date.x, common$date.y)
69
##  [1] FALSE
70
common$tile[common$date.x!=common$date.y]
71
##  "srtm_19_06.zip" "srtm_21_14.zip" "srtm_42_14.zip"