Saya menggunakan toolbox pemrosesan (luar biasa!) Di QGIS 2.0 untuk mengklasifikasi ulang beberapa nilai NoData dalam raster dari -9999 menjadi 0, tetapi pada titik tertentu dalam proses tersebut ada offset sekitar 0,2 derajat yang diperkenalkan ke lapisan output.
Saya menggunakan alat "Reklasifikasi (Sederhana)" (saya percaya itu bagian dari SAGA), dengan beberapa parameter yang cukup mendasar:
. . . dan tempat input raster terletak dengan benar dibandingkan dengan kisi hex ini:
. . . raster keluaran ditabrak tenggara oleh dua sel hex, atau sekitar 0,2 derajat:
Asal input adalah -180, 83.623. . . tetapi asal output adalah -179.05,83.05. Biasanya saya akan menganggap ini adalah hasil dari konversi WGS84 yang salah, tetapi input, output dan proyek semuanya EPSG: 4326.
Dari mana asal offset ini berasal, dan apakah layak mengajukan tiket SAGA atau hanya sesuatu yang saya abaikan?
Perbarui Log proses menunjukkan transformasi dilakukan pada asal input:
Transformation x' = -179.999989 + x * 0.100000 + y * 0.000000
Transformation y' = 83.623125 + x * 0.000000 + y * -0.100000
dan kemudian beberapa baris kemudian parameter baru muncul, sekarang salah:
Grid system: 0.1; 3600x 1735y; -179.949989x -89.826875y
saga_gui
perintah?Jawaban:
Offset Anda adalah 1/2 piksel (0,1 derajat).
Garis bujur:
-179.999989 - -179.949989 = 0,05
Latitude kurang jelas karena satu kotak direferensikan dari kiri atas, yang lain dari kiri bawah:
83.623125 - (1735 * 0.1) - -89.826875 = 0,05
Saya menduga ini berkaitan dengan masalah RasterPixelIsPoint atau RasterPixelIsArea:
http://www.remotesensing.org/geotiff/faq.html#PixelIsPoint
Apakah Anda menggunakan gdal 1.7 atau 1.8?
sumber