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.
sumber
Jawaban:
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
sumber