Project

General

Profile

Download (1.43 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
# taset rasterXSize="392401" rasterYSize="79201">
19
#<VRTDataset rasterXSize="219601" rasterYSize="79201"> divide by 3.0 to get dimensions for 3 arc seconds resolution
20
date
21
rm /home/reeves/active_work/OutProducts/EastHemi/NorthWestEast/NorthWestEastHemiAster3ArcSec.img
22
gdalwarp -of HFA -ot Int16 -ts 73193 26400 -r bilinear -srcnodata -9999 -dstnodata -9999 \
23
/home/reeves/active_work/OutProducts/EastHemi/NorthWestEast/NorthWestEastHemiAster.vrt \
24
/home/reeves/active_work/OutProducts/EastHemi/NorthWestEast/NorthWestEastHemiAster3ArcSec.img
25
date