Bagaimana cara membuat 4096x4096 png dari kumpulan data 1440x720 menggunakan GDAL?

8

Saya punya waktu mencoba melakukan sesuatu yang tampaknya cukup mudah. Saya mencoba untuk mengambil file netcdf dari ftp://eclipse.ncdc.noaa.gov/pub/OI-daily-v2/NetCDF/2011/AVHRR-AMSR/ dan mengubahnya menjadi 4096x4096 png sehingga saya dapat memasang ubin itu untuk digunakan dengan peta google. Saya telah menggunakan perintah berikut tetapi sepertinya selalu membuat area kosong di perbatasan timur dan barat:

/usr/local/bin/gdal_translate -a_srs EPSG:4326 -outsize 4096 4096 -a_ullr 0 90 360 -90 NETCDF:"input.nc":sst -of netCDF output.nc
/usr/local/bin/gdaldem color-relief output.nc -b 1 -alpha colorscale.txt -of PNG output.png
/usr/bin/convert output.png -transparent #FFFFFF output.png

Apakah saya hanya mendekati ini dengan cara yang salah? Apakah saya menggunakan alat yang salah?

FWIW - gdal_translate adalah v1.8.0

Terima kasih sebelumnya.

UPDATE: "Jawaban" di bawah dari Dija - kolega saya - sebenarnya adalah tanggapan atas komentar pada pertanyaan itu sendiri.

UPDATE Secara tidak sengaja menambahkan jawaban, bukan komentar. Pindah ke sini.

Saya pikir saya akan menimbang di sini. Saya telah mengerjakan masalah ini untuk sementara waktu sekarang. Saya sudah mencoba banyak iterasi proses. Saya mengubah gambar ini ke 4096x4096 karena sebenarnya menciptakan gambar yang lebih baik daripada jika tidak. Artinya, jika saya membiarkan gambar dalam ukuran aslinya, 1440x720 saya mendapatkan celah yang lebih besar di meridian utama (untuk kumpulan data ini) dan pada meridian ke-180 untuk set data lainnya. Saya mengubah gambar output menjadi kelipatan 256 (agar lebih cocok dengan google) dan mendapatkan berbagai hasil.

Jika saya skalakan ke 1536x768, celah akan muncul di sebelah kanan 0 meridian. Jika saya skala ke 1280x640 itu ke kiri. Jika saya biarkan sendiri, itu cukup banyak membelahnya.

Saya telah mengambil gambar dasar yang dimuntahkan gdaldem, memotongnya menjadi dua dan mencocokkan ujungnya ... mereka berbaris sempurna. Ini membuat saya berpikir bahwa ada sesuatu yang aneh terjadi dengan gdal2tiles.py

Saya menulis tiler saya sendiri bahwa saya sedang bersiap-siap untuk melepaskan gambar kami untuk melihat apakah ada celah.

Saya akan mengumpulkan beberapa cuplikan layar untuk menunjukkan apa yang saya bicarakan dengan bagaimana kesenjangannya serta melaporkan apa yang dilakukan tiler saya.

Terima kasih atas masukan yang Anda miliki. Itu dihargai.

PEMBARUAN 2 Saya menjalankan maptiler yang saya tulis terhadap raster asli dan berfungsi dengan baik tanpa ada celah. gdal2tiles.py tampaknya menjadi masalah. Ini memperkenalkan tepi transparan ke ubin terakhir yang menyentuh meridian ke-180.

malonso
sumber
Sesuatu harus memberi ketika Anda mengubah gambar dengan rasio aspek 1: 2 menjadi satu dengan rasio aspek 1: 1: baik terjadi kliping, Anda mendapatkan batas kosong untuk mengisi ruang kosong, atau Anda memerlukan peregangan 2: 1.
whuber
Saya telah menjalankan proses yang sama pada salah satu file dari situs FTP yang Anda berikan, dan saya memiliki gambar yang sesuai dengan batas-batasnya dengan benar. Bisakah Anda memberikan tangkapan layar dari efek yang Anda dapatkan? Apakah efeknya terjadi setelah langkah gdaldem, atau langkah konversi?
MerseyViking
Ah ya, seperti kata whuber, rasio aspek Anda salah jika Anda menggunakan parameter yang Anda berikan dalam teks, tetapi parameter baris perintah, yang merupakan parameter yang saya gunakan, memiliki rasio aspek yang benar.
MerseyViking
Karena itu, saya baru saja mencoba membuatnya sebagai 4096 x 4096 dan QGIS memuatnya tanpa masalah. Selain itu, memuat kedua versi ke GIMP menunjukkan peta yang mengarah ke kanan dan ke kiri.
MerseyViking
@whuper - Tentu saja, saya sepenuhnya setuju. Saya menyadari bahwa mencoba membuat dataset 1: 2 menjadi gambar 1: 1, Anda hampir membuat data tidak valid. Karena itu, harapan saya adalah untuk mendapatkan output yang pada dasarnya adalah peregangan 2: 1, daripada ruang kosong. Saya pada belas kasihan dari persyaratan bisnis yang ditetapkan di depan saya :( Oh hal-hal yang akan kita lakukan untuk gambar yang bagus. @MerseyViking - Menarik, biarkan saya kembali kepada Anda tentang hal itu. dan jadi saya perlu mendapatkan masukan mereka. Terima kasih atas umpan baliknya.
malonso

Jawaban:

2

gdal_grid memiliki parameter yang memungkinkan Anda menentukan koordinat luas raster keluaran Anda.

gdal_grid -txe [xcoord1] [xcoord2] -tye [ycoord1] [ycoord2]

di mana koordinat kiri bawah adalah (xcoord1, ycoord1) dan kanan atas adalah (xcoord2, ycoord2)

GDAL akan menginterpolasi data Anda sesuai dengan luaran output Anda

WxMacMan
sumber