|
1
|
# Script for clipping each USGS GTOPO30 coverages to extent: N59 to N82
|
|
2
|
# degree lat. (to match extent of ASTER GDEM tiles, for the use in
|
|
3
|
# Check Part 2) The clipped data will be saved at
|
|
4
|
# "/data/project/organisms/DEM/Yuni/Data/GTOPO30/"
|
|
5
|
#
|
|
6
|
# Note: if I knew how to apply gdal_translate in R, the script could be
|
|
7
|
# less redundant.
|
|
8
|
#
|
|
9
|
# December 9th 2011
|
|
10
|
# Yuina Nunokawa
|
|
11
|
|
|
12
|
|
|
13
|
|
|
14
|
|
|
15
|
|
|
16
|
|
|
17
|
|
|
18
|
export GTOPO="/data/project/organisms/DEM/usgsGTOPO30"
|
|
19
|
export CLIP="/data/project/organisms/DEM/Yuni/Data/GTOPO30/"
|
|
20
|
|
|
21
|
gdal_translate -projwin 20 81.99 59.99 59 $GTOPO/e020n90/E020N90.DEM $CLIP/E020N82_Clipped.dem
|
|
22
|
|
|
23
|
gdal_translate -projwin 60 81.99 99.99 59 $GTOPO/e060n90/E060N90.DEM $CLIP/E060N82_Clipped.dem
|
|
24
|
|
|
25
|
gdal_translate -projwin 100 81.99 139.99 59 $GTOPO/e100n90/E100N90.DEM $CLIP/E100N82_Clipped.dem
|
|
26
|
|
|
27
|
gdal_translate -projwin 140 81.99 179.99 59 $GTOPO/e140n90/E140N90.DEM $CLIP/E140N82_Clipped.dem
|
|
28
|
|
|
29
|
gdal_translate -projwin -20 81.99 19.99 59 $GTOPO/w020n90/W020N90.DEM $CLIP/W020N82_Clipped.dem
|
|
30
|
|
|
31
|
gdal_translate -projwin -60 81.99 -20.01 59 $GTOPO/w060n90/W060N90.DEM $CLIP/W060N82_Clipped.dem
|
|
32
|
|
|
33
|
gdal_translate -projwin -100 81.99 -60.01 59 $GTOPO/w100n90/W100N90.DEM $CLIP/W100N82_Clipped.dem
|
|
34
|
|
|
35
|
gdal_translate -projwin -140 81.99 -100.01 59 $GTOPO/w140n90/W140N90.DEM $CLIP/W140N82_Clipped.dem
|
|
36
|
|
|
37
|
gdal_translate -projwin -180 81.99 -140.01 59 $GTOPO/w180n90/W180N90.DEM $CLIP/W180N82_Clipped.dem
|
|
38
|
|
|
39
|
|
