ปรับให้เรียบ / reinterpolating แรสเตอร์ด้วย GDAL?

สั้น:
ฉันต้องการเปลี่ยนความละเอียดของแรสเตอร์และปรับทางลาดสีเทาตามที่แสดงในภาพตะโกน การตั้งค่าคือการใช้ GDAL, PIL หรือ Numpy

คำอธิบาย:
ฉัน Kriging จุดลงไปในเชิงภาพที่มีความละเอียดการส่งออกของ 20 เมตรกับGeostatistical ห้องสมุดที่มีประสิทธิภาพสูง ฉันไม่ต้องการเปลี่ยนความละเอียดเอาต์พุตเนื่องจากเวลาการแก้ไขเพิ่มขึ้นเป็นทวีคูณ
ด้วยความละเอียดนี้ภาพที่ส่งออกจะน่าเกลียด (pixelated และ aliased) ฉันไม่รู้ว่ามันถูกต้องในแนวความคิดหรือไม่ แต่ฉันต้องการภาพที่นุ่มนวลเหมือนในตัวอย่างร้อง มันเหมือนกับการ 'ปรับกำลังภาพ' ใหม่ให้เป็นภาพที่มีความละเอียดมากกว่า ฉันใช้ python ดังนั้นการตั้งค่าของฉันคือ GDAL, Python Imaging Library หรือ Numpy คำตอบอาจเป็นเชิงทฤษฎีเช่นชี้ให้เห็นชื่ออัลกอริทึมหรือแนวคิดของการดำเนินการชนิดนี้

ที่มา:

ปลายทาง:

แก้ไขผลลัพธ์ด้วย gdalwarp cubic spline:

1) วิธีที่ยาก: ด้วยการเข้ารหัสเล็กน้อย (ค่อนข้าง) ง่ายต่อการใช้การแก้ไขแบบสองทางเพื่อให้ได้การ resampling ที่เหมาะสม

2) วิธีที่ง่าย: ใช้ GDAL ตามที่อธิบายไว้ในโพสต์ GISSE ก่อนหน้านี้แต่กลับด้าน (ลดขนาดพิกเซล)

ใช้GDALReprojectImageซึ่งแสดงใน Python:

from osgeo import gdal
help(gdal.ReprojectImage)

สำหรับการแก้ไขที่ราบรื่นให้ใช้วิธี bilinear หรือลูกบาศก์ ฟังก์ชั่นนี้น่าอึดอัดใจเนื่องจากไม่มีข้อโต้แย้งคำหลักดังนั้นคุณต้องค้นหาตำแหน่ง:

gdal.ReprojectImage(src_ds, dst_ds, None, None, gdal.GRA_Bilinear)

อาจเป็นเรื่องยากที่จะตั้งค่าส่วนdst_dsที่จำเป็นต้องมี geotransform คล้ายกับsrc_dsแต่มีการปรับขนาดเซลล์

เพื่อให้การเปลี่ยนแปลงราบรื่นขึ้นคุณต้องใช้ตัวกรองความถี่ต่ำ คุณสามารถเขียนของคุณเองโดยใช้ GDAL หรือมีคนใช้ GRASS ฉันยังไม่ได้ลอง แต่นี่เป็นคำแนะนำhttp://wiki.awf.forst.uni-goettingen.de/wiki/index.php/Exercise_31

คุณอาจต้องการตัวอย่างแรสเตอร์ของคุณก่อนที่จะใช้ฟิลเตอร์ low-pass เพื่อให้ได้ความละเอียดที่ดีขึ้น

คุณสามารถใช้ตัวกรองอันดับ / มัธยฐานกับรัศมี = 5 เช่นขนาดเคอร์เนล = 11, (สำหรับแต่ละช่อง rgb)