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
|