Pendek:
Saya ingin mengubah resolusi raster dan menghaluskan jalan abu-abu seperti yang ditunjukkan pada gambar di bawah. Preferensi adalah menggunakan GDAL, PIL atau Numpy.
Deskripsi:
Saya mengarahkan poin ke raster dengan resolusi keluaran 20 meter dengan High Performance Geostatistic Library . Saya tidak ingin mengubah resolusi output karena waktu interpolasi meningkat secara eksponensial.
Dengan resolusi ini gambar output jelek (pixelated dan alias). Saya tidak tahu apakah itu benar secara konseptual tetapi saya ingin gambar menjadi lebih halus seperti pada contoh di bawah ini. Ini seperti 'reinterpolasi' gambar menjadi resolusi yang lebih baik. Saya menggunakan python sehingga preferensi saya adalah GDAL, Python Imaging Library atau Numpy. Jawabannya bisa teoretis, seperti menunjukkan nama algoritma atau konsep operasi semacam ini.
Sumber:
Tujuan:
Hasil EDIT dengan gdalwarp cubic spline:
Jawaban:
1) Cara yang sulit: Dengan sedikit pengkodean (relatif) mudah untuk menerapkan interpolasi bilinear untuk mencapai resampling yang layak.
2) Cara mudah: gunakan GDAL seperti yang dijelaskan dalam posting GISSE sebelumnya , tetapi secara terbalik (mengurangi ukuran piksel).
sumber
Gunakan GDALReprojectImage , yang diekspos dengan Python:
Untuk interpolasi yang halus, gunakan metode bilinear atau kubik. Fungsi ini canggung, karena tidak mengambil argumen kata kunci, sehingga Anda perlu menemukan posisi:
Mungkin bagian yang sulit adalah pengaturan
dst_ds
, yang perlu memiliki geotransform serupasrc_ds
, tetapi dengan ukuran sel yang dimodifikasi.sumber
GRA_CubicSpline
(coba ini dulu) vs.GRA_Cubic
(konvolusi).Untuk memperlancar variasi, Anda memerlukan filter low-pass. Anda bisa menulis sendiri menggunakan GDAL, atau ada yang menggunakan GRASS. Saya belum mencobanya, tapi ini panduan http://wiki.awf.forst.uni-goettingen.de/wiki/index.php/Exercise_31
Anda mungkin ingin mengambil sampel raster Anda terlebih dahulu sebelum menerapkan filter low-pass untuk mendapatkan output resolusi yang lebih baik.
sumber
Anda dapat menggunakan filter peringkat / median dengan jari-jari = 5, yaitu ukuran kernel ukuran = 11, (untuk setiap saluran rgb).
sumber