Project

General

Profile

Download (1.32 KB) Statistics
| Branch: | Revision:
1 674bbf95 Rick Reeves
#! /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
date
19
rm /home/reeves/active_work/OutProducts/EastHemi/SouthEastEast/SouthEastEastHemiCgiar_30ArcSecBL.img
20
gdalwarp -of HFA -ot Int16 -te 64.999583 -50.0004165 180.0004166 5.0004166 -ts 13800 6600  -r bilinear -srcnodata -9999 -dstnodata -9999 \
21
/home/reeves/active_work/OutProducts/EastHemi/SouthEastEast/SouthEastEastHemiCgiarSrtm.vrt \
22
/home/reeves/active_work/OutProducts/EastHemi/SouthEastEast/SouthEastEastHemiCgiar_30ArcSecBL.img
23
date