SRTM4 çizimi

From GMT Türkiye Wiki
(Sürümler arası farklar)
Jump to: navigation, search
38. satır: 38. satır:
 
olması gerekir. İşte bu yüzden sağdaki devamı olan ikinci gridin başlangıç koordinatını 30.00083333333333333 olarak değiştirmemiz gerekir.  
 
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
+
ncols        6001
nrows        6001
+
nrows        6001
xllcorner    30.00083333333333333
+
xllcorner    30.00083333333333333
yllcorner    40
+
yllcorner    40
cellsize      0.00083333333333333
+
cellsize      0.00083333333333333
NODATA_value  -9999
+
NODATA_value  -9999
 
   
 
   
 
Daha sonra ikisini birleştirmek için aşağıdaki komutu verebilirsiniz
 
Daha sonra ikisini birleştirmek için aşağıdaki komutu verebilirsiniz

Sayfanın 16:14, 8 Nisan 2009 tarihindeki hali

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)

Personal tools