|
1
|
#! /bin/sh
|
|
2
|
################################################################################
|
|
3
|
# script Name: mosaicNwHemiAster30ArcSecBL.sh
|
|
4
|
# This script constructs the CIGAR SRTM image mosaic for the Northweetern Hemisphere
|
|
5
|
# quadrant of the global mosaic
|
|
6
|
#
|
|
7
|
# Author: Rick Reeves, NCEAS June 9, 2011
|
|
8
|
#
|
|
9
|
# NOTE: We are using SRTM version 4.1 data distributed by CGIAR (6001 x 6001 tiles)
|
|
10
|
# Determining the -ts (row, column # cells):
|
|
11
|
# 1) Look at the first line of corresponding .vrt file: contains the spanning extent
|
|
12
|
# (# rows, columns) of all image tiles in the virtual mosaic, as well as the pixel
|
|
13
|
# resolution in 'native units' for the spatial reference frame (in our case,
|
|
14
|
# the geographic coordinate system and decimal degrees).
|
|
15
|
# Example for this file:
|
|
16
|
#
|
|
17
|
#Transform> -1.8000013888888890e+02, 2.7777777777778255e-04, 0.0000000000000000e+00,
|
|
18
|
# 8.3000138888888884e+01, 0.0000000000000000e+00, -2.77777777777782 55e-04</GeoTransform>
|
|
19
|
# in decimal cegrees: 1 arc-seconds, or approx 30 meters.
|
|
20
|
#<VRTDataset rasterXSize="648001" rasterYSize="82801">
|
|
21
|
#
|
|
22
|
# The extent of the mosaic is fixed. To change the cell size, adjust the X and Y resolution.
|
|
23
|
# In the present case, to set cellsize at 30 arcseconds, reduce the number of cells in
|
|
24
|
# X and Y by a factor of 30: from 21601 and 2761 to 64801 and 82801.
|
|
25
|
################################################################################
|
|
26
|
#
|
|
27
|
rm /home/reeves/active_work/OutProducts/WestHemi/NorthWest/NorthWestHemiAster_30ArcSecBL.img
|
|
28
|
gdalwarp -of HFA -ot Int16 -ts 21601 2761 -r bilinear -srcnodata -9999 -dstnodata -9999 \
|
|
29
|
/home/reeves/active_work/OutProducts/WestHemi/NorthWest/NorthWestHemiAster.vrt \
|
|
30
|
/home/reeves/active_work/OutProducts/WestHemi/NorthWest/NorthWestHemiAster_30ArcSecBL.img
|
|
31
|
date
|
