SRTM4 çizimi

Aşağıdaki adresten SRTMv4 indirilebilir.

http://srtm.csi.cgiar.org/SELECTION/inputCoord.asp

Bu veri seti SRTMv3'ün 5x5 gridlenmiş hali. Verideki boşluklar resmi SRTM versiyon 2'den farklı yöntemlerle doldurulmuş.

Bu siteden veri indirmek için istediğiniz çerçeveleri fare ile harita üzerinde seçin. File format olarak Arc Info ascci'yi seçin ve aynı hizada sağda "Click here to begin search" düğmesine basın. Açılan sayfada seçtiğiniz gridleri indirin.

Örneğin ben kuzeybatı Anadolu için iki çerçeve seçtim. Bunlar srtm_42_04.zip srtm_43_04.zip olarak geliyor. Bunları açıp içindeki srtm_42_04.asc ve srtm_43_04.asc dosyalarını bulunduğunuz klasör içine atın.

Şimdi bunları GMT grd'ye çevirelim

 xyz2grd -E srtm_42_04.asc -Gsrtm_42_04.grd
 xyz2grd -E srtm_43_04.asc -Gsrtm_43_04.grd

Bunu plot etmek için SRTM çizimi bölümüne bakınız.

Bu gridlerin GMT ile kolayca mozaik hale getirilmesi zor. Elle gridlerin kordinatların değiştirilmesi gerekiyor. srtm_42_04.asc'i text editörü ile açıp bakarsak başlık kısmında aşağıdaki satırları görebilirsiniz.

ncols         6001
nrows         6001
xllcorner     24.999583987092
yllcorner     39.999583575447
cellsize      0.00083333333333333
NODATA_value  -9999

Bunu şu şekilde düzeltelim.

ncols         6001
nrows         6001
xllcorner     25
yllcorner     40
cellsize      0.00083333333333333
NODATA_value  -9999

6001x0.00083333333333333=5.000833333 yapar. Dolayısıyla bu gridin doğu kenarı 25+5.000833333=30.000833333 olması gerekir. İşte bu yüzden sağdaki devamı olan ikinci gridin başlangıç koordinatını 30.00083333333333333 olarak değiştirmemiz gerekir.

ncols 6001 nrows 6001 xllcorner 30.00083333333333333 yllcorner 40 cellsize 0.00083333333333333 NODATA_value -9999

Daha sonra ikisini birleştirmek için aşağıdaki komutu verebilirsiniz

grdpaste    42.grd  43.grd  -Gmoz.grd

Bu mantığı altalta (yani kuzey güney) olan gridlerin birleştirilmesindede kullanmak gerekiyor. Ancak oldukça çetrefilli olduğundan başka yöntemlerle birleştirilmesi daha mantıklı bir yaklaşım olur.

--Ziyadin 14:55, 8 Nisan 2009 (CEST)