|
1
|
#! /bin/sh
|
|
2
|
# script uses gdalwarp to create mosaic of the South Western Hemisphere
|
|
3
|
# (to 60 degrees North Latitude) with 3ArcSec resolution
|
|
4
|
# using CGIAR/SRTM tiles adjacent to 60 Degrees North Latitude.
|
|
5
|
#
|
|
6
|
# 2011/06/02: -te (georef extent) and
|
|
7
|
# -ts (number of rows, cols to cover w/3 arcsecond pixels) combination
|
|
8
|
# to get exact 3 arc second pixels.
|
|
9
|
# To get # pixels (ts parameter), multiply image extent, in degrees,
|
|
10
|
# by 3600 (number of arc seconds per degree)
|
|
11
|
# Author: Rick Reeves, NCEAS May 11, 2011
|
|
12
|
# Note: This script accesses the (no longer used) 'even 5 degree' 6000 column version
|
|
13
|
# of the CGIAR SRTM data set.
|
|
14
|
## NOTE: Computing the -ts parameter:
|
|
15
|
# 1) Calculate difference between bounding latitude (longitude) pairs,
|
|
16
|
# 2) multiply by 120 (3600 arc seconds per degree / 30 arc second resolution)
|
|
17
|
# NOTE: This script uses the SRTM version 4.1 data distributed by CGIAR (6001 x 6001 tiles)
|
|
18
|
# <VRTDataset rasterXSize="126001" rasterYSize="66001"> from VRT native SRTM res is 3 arc sec
|
|
19
|
date
|
|
20
|
rm /home/reeves/active_work/OutProducts/EastHemi/NorthEastEast/NorthEastEastHemiCgiar_30ArcSecBL.img
|
|
21
|
gdalwarp -of HFA -ot Int16 -ts 12601 6601 -r bilinear -srcnodata -9999 -dstnodata -9999 \
|
|
22
|
/home/reeves/active_work/OutProducts/EastHemi/NorthEastEast/NorthEastEastHemiCgiarSrtm.vrt \
|
|
23
|
/home/reeves/active_work/OutProducts/EastHemi/NorthEastEast/NorthEastEastHemiCgiar_30ArcSecBL.img
|
|
24
|
date
|
