Di mana offset raster ini diperkenalkan di kotak alat pemrosesan QGIS?

8

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:

reklasifikasi parameter

. . . dan tempat input raster terletak dengan benar dibandingkan dengan kisi hex ini:

masukan raster

. . . raster keluaran ditabrak tenggara oleh dua sel hex, atau sekitar 0,2 derajat:

masukkan deskripsi gambar di sini

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
Bill Morris
sumber
1
Anda harus menguji alat SAGA asli menggunakan GUI atau CLI dan melihat apakah Anda mendapatkan hasil yang sama.
Giovanni Manghi
1
Terlalu banyak jalan buntu kesalahan dari instal dokumen ini di komputer saya - Saya tidak tertarik menghapus konfigurasi GDAL / PROJ saya saat ini hanya untuk menguji SAGA secara asli.
Bill Morris
1
Jika Anda menggunakan SAGA di QGIS maka SAGA sudah terpasang di mesin Anda.
Giovanni Manghi
Setelah lebih dari 2 jam meneliti dokumen, saya tidak tahu bagaimana CARA MELUNCURKAN GUI DI OSX. Dokumentasi mengingatkan R dalam semua cara yang buruk. Tidak ada variasi pada "saga" (-gis, -gui, dll.) Menghasilkan apa pun di baris perintah, dan kompilasi dari sumber mengembalikan kesalahan yang sama sekali tidak membantu. Berterima kasihlah kepada surga yang berbelas kasih (devs inti, yaitu) bahwa kekacauan ini sekarang dapat diakses oleh manusia di sidebar QGIS. Saya sangat berharap bahwa SAGA bukan yang mengubah asal output raster saya, karena saya tidak akan pernah bisa men-debug itu.
Bill Morris
Apakah Anda dapat meluncurkan saga dengan saga_guiperintah?
SaultDon

Jawaban:

1

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?

Marc Pfister
sumber
Seluruh sistem saya mendapatkan GDAL 1.10.1; untuk alasan apa pun saya tidak bisa mendapatkan konsol python QGIS untuk melaporkan versi yang digunakan di sana
Bill Morris
terlihat jauh lebih dari 1/2 piksel pada gambar.
radouxju
Ya, tetapi ukuran piksel tidak aktif juga. Inilah gdalinfo untuk tiffs input dan output
Bill Morris
Ah, melihat gdalinfo sepertinya pergeserannya adalah .95 dan -0.5731265. Pergeseran 1/2 piksel di atas mungkin merupakan area yang tepat untuk konversi titik. Raster asli Anda memiliki beberapa batasan aneh, mungkin mereka mengacaukan konversi.
Marc Pfister
Poin yang bagus. Aneh bahwa itu tergantung dari ujung dunia. Raster adalah hasil konversi dari vektor; biarkan saya melihat apakah saya dapat mengulangi konversi itu dengan batas yang benar.
Bill Morris