Apakah ada standar untuk koordinat piksel dalam raster yang di-georeferensi?

Kami bekerja dengan kumpulan data raster yang berbeda, dan koordinat piksel disimpan di tengah piksel atau di satu sudut. Informasi ini biasanya disediakan dalam metadata, tetapi kadang-kadang kita perlu menebak dan posting ini menggambarkan bahwa itu tidak selalu mudah ditemukan.

Karena kami juga menghasilkan beberapa layer, saya bertanya-tanya apakah ada aturan standar dan / atau praktik terbaik untuk menetapkan koordinat untuk piksel, idealnya dengan referensi atau dengan pro / kontra.

Ini agak berantakan. Anda harus membaca setidaknya dokumen GeoTIFF ini

http://www.remotesensing.org/geotiff/spec/geotiff2.5.html#2.5.2.2 dan beberapa pertimbangan GDAL

http://trac.osgeo.org/gdal/wiki/rfc33_gtiff_pixelispoint .

Sebagai aturan praktis semua raster (foto udara, satelit) menggunakan pixel-is-area dan data pengukuran seperti DEM menggunakan pixel-is-point. Tetapi bagaimana itu mencerminkan georeferensi dan luasan gambar sayangnya tidak begitu sederhana. Untuk GDAL titik jangkar piksel selalu di sudut kiri atas piksel - jika Anda menyimpan Ground Control Point (GPC) baru, koordinat yang diberikan selalu terikat ke sudut kiri atas piksel yang sesuai dan itu penting jika piksel didefinisikan sebagai area atau titik. GDAL juga melaporkan luasan yang sama untuk gambar pixel-is-area dan pixel-is-point. Terserah pengguna atau perangkat lunak yang menggunakan GDAL untuk memperlakukan mereka secara berbeda jika perlu.

Aturan OGC saat ini untuk georeferensi cakupan adalah dengan menggunakan kisi yang ditentukan dalam spesifikasi GML v.3.2.1 yang panjangnya 437 halaman. Itu membuat setengah piksel perbedaan dalam lokasi asal untuk pixel-is-point dan pixel-is-area karena dalam kasus pixel area asal akan diatur ke tengah pixel. Kutipan berikut berasal dari GML dalam standar JPEG2000:

Catatan: GMLJP2 mengikuti definisi kisi-kisi di GML 3.2.1 [OGC 07-036] klausa 19.2.2: “Ketika titik kisi digunakan untuk mewakili ruang sampel (misalnya piksel gambar), titik kisi mewakili pusat dari ruang sampel (lihat ISO 19123: 2005, 8.2.2) ”. Ini sesuai dengan nilai pixelInCell dari ImageCRS yang diatur ke CellCenter sebagaimana ditentukan dalam ISO 19111. Ini dapat diartikan sebagai asal RectifiedGrid adalah titik tengah piksel sudut.

Untuk membuat hal-hal tidak terlihat terlalu sederhana, beberapa program seperti ERDAS bayangkan (AFAIK) mengikuti internal skema pixel-is-point serta MapServer hingga versi 7.0 yang akan datang, dengan mencoba mengubah on-the-fly menjadi skema pixel-is-area. Membaca tentang ini: http://mapserver.org/fr/development/rfc/ms-rfc-107.html .

Untuk interoperabilitas terbaik, saya akan merekomendasikan untuk menggunakan pixel-is-area dan GeoTIFFs untuk semua raster alami (gambar udara dan satelit, peta topografi). Ini adalah kombinasi yang terkenal dan dengan demikian dapat diandalkan. Jika Anda bermain dengan DEM dan data pengukuran / sensor lainnya dan setengah pergeseran piksel dapat bermakna, sadari bahwa mereka mungkin dilabeli sebagai titik-pixel dan program yang berbeda mungkin atau mungkin tidak memperkenalkan pergeseran setengah piksel.

Keduanya umum, dan keduanya tidak dapat dianggap sepenuhnya standar.

Untuk GeoTIFF, keduanya dimungkinkan - lihat Bagian 2.5.2.2 untuk GTRasterTypeGeoKeyyang menjelaskan metode interpretasi. The GeoTIFF FAQ menyarankan menggunakan default ( PixelIsArea) nilai yang tag untuk kompatibilitas dengan versi GDAL.

The Dunia Format file menggunakan pusat pixel.

Memotong gambar dari NITF (menggunakan ekstensi ICHIPB) secara eksplisit mengkodekan posisi piksel pertama dalam gambar terkelupas (subset) sebagai koordinat (0,5. 0,5), dan memetakan yang kembali ke koordinat tertentu pada gambar induk (superset) misalnya ( 2.5, 1.5).

Saya tidak tahu apakah ada konvensi, tetapi koordinat pada gambar biasanya merujuk ke pusat pixel. Ini dapat berbeda dalam hal raster yang dihasilkan dari, misalnya, model komputer.

Jika Anda tidak yakin, Anda dapat menggunakan GMT , Alat Pemetaan Generik, untuk menguji Anda raster dan mengonversinya ke beberapa format yang sesuai, karena itu secara eksplisit menawarkan opsi untuk menentukan pendaftaran kisi atau node. Lihat bagian tutorial ini, dan juga manual untuk fungsi grdedit .