Apa sistem koordinat raster netCDF ini?

10

Saya memperoleh file raster netCDF, tetapi saya tidak bisa mendapatkan metadata apa pun untuk mendapatkan nama sistem koordinat raster yang telah dibangun. Raster itu sendiri tidak memiliki sistem koordinat yang tertanam. Saya pikir itu hanya WGS84, dan kelihatannya seperti itu sekilas, tetapi dengan investigasi lebih lanjut di ArcMap, saya melihat bahwa itu adalah sistem yang agak tidak umum. Berikut ini tampilannya: masukkan deskripsi gambar di sini

Oranye raster adalah raster normal di WGS84 yang saya masukkan di sini untuk tujuan perbandingan. Yang ungu adalah raster dengan sistem koordinat yang tidak diketahui. Apakah Anda punya petunjuk apa ini?

Beberapa pembaruan: Inilah raster netCDF: https://www.dropbox.com/s/nottbl9yt6dwss6/sic_average_nclimate.nc?dl=0 Saya juga bisa mendapatkan beberapa metadata dari penyedia gambar:

netcdf sic_average_nclimate {
dimensions:
    nlon = 361 ;
    nlat = 90 ;
    nseas = 4 ;
variables:
    float SIC_Change(nlat, nlon) ;
        SIC_Change:Title = "Gridded Multi-Model Ensemble Mean Annual Mean Change in Ice Concentration 21C-20C" ;
    float SIC_Season_Change(nseas, nlat, nlon) ;
        SIC_Season_Change:Title = "Gridded Multi-Model Ensemble Mean Seasonal Mean Change in Ice Concentration 21C-20C" ;
    float SIC_Change_STD(nlat, nlon) ;
        SIC_Change_STD:Title = "Gridded Multi-Model Standard Deviation of the Annual Mean Change in Ice Concentration 21C-20C" ;
    float SIC_Season_Change_STD(nseas, nlat, nlon) ;
        SIC_Season_Change_STD:Title = "Gridded Multi-Model Standard Deviation of the Seasonal Mean Change in Ice Concentration 21C-20C" ;
    float LAT(nlat) ;
        LAT:Title = "Latitude" ;
    float LON(nlon) ;
        LON:Title = "Longitude" ;

// global attributes:
        :Title = "Ice Concentration metrics for Model subset as in Figure 1 of NCLIMATE paper" ;

Mereka menunjukkan batas panjang, tetapi tampaknya tidak ada informasi mengenai sistem koordinat.

multigoodverse
sumber
1
Antartika biasanya dalam Proyeksi Stereografi Kutub sehingga semuanya bergerak ke utara dari kutub selatan. nsidc.org/data/polar_stereo/ps_grids.html
Mapperz
1
Aneh, terlihat seperti 0-360, tetapi mengapa ia merencanakannya di belahan bumi utara?
mkennedy
1
Jika ekstensi yang dilaporkan di header NetCDF terlihat seperti lat / lon, tetapi lat: 0 hingga?, Lon: 0 hingga 360. Bisakah Anda meretasnya ke lat: -90 ke (? - 90) dan lon: -180 ke + 180?
mkennedy
@Mapperz & mkennedy, saya mencoba menetapkannya Proyeksi Stereografis Polar, tetapi tidak berhasil. Saya mungkin bisa memantulkannya secara vertikal ke bawah, tapi sisi kanan setengahnya masih di luar garis bujur 180E. Saya membuka netCDF dengan notepad, tetapi tidak dapat menemukan informasi yang terkait dengan lat, panjang. Saya lampirkan raster dalam pertanyaan jika Anda ingin membukanya di sisi Anda.
multigoodverse
Tautan itu tidak berfungsi - dropbox mengatakan "pemilik tidak memiliki akses yang diotorisasi".
BradHards

Jawaban:

6

Anda dapat menggeser data ke posisi yang benar menggunakan perintah GDAL ini:

gdal_translate -a_srs EPSG:4326 -a_ullr 0 0 360 -90 NETCDF:"sic_average_nclimate.nc":SIC_Change change360.tif
gdalwarp -t_srs WGS84 change360.tif change180.tif  -wo SOURCE_EXTRA=1000 --config CENTER_LONG 0

(dengan sedikit bantuan oleh Frank Warmerdam: Cara memproyeksikan ulang raster dari 0 360 ke -180 180 dengan memotong 180 meridian )

dan hasilnya terlihat seperti ini dalam proyeksi stereografi:

masukkan deskripsi gambar di sini

AndreJ
sumber
1

(Ini bukan jawaban, melainkan saran proses yang terlalu panjang untuk dijadikan komentar)

Anda dapat melakukan geo-referensi di QGIS menggunakan alat GeoReferencer karena Anda tampaknya dapat mengidentifikasi poin-poin umum pada kedua raster. Ini akan menghasilkan skrip GDAL dengan satu set poin GCP, dan parameter transformasi. Anda dapat menambahkan skrip ini dengan contoh poin-poin ini ke pertanyaan Anda, karena dapat menjelaskan apa perbedaan antara kedua proyeksi tersebut.

Skenario kasus terburuk, Anda akan memiliki perkiraan kasar dari nilai-nilai berdasarkan raster yang diposisikan ulang yang mungkin membantu Anda dengan tujuan Anda.

djq
sumber
+1 dan ya, georeferensi mungkin menjadi poin terakhir, tetapi prosesnya akan menghasilkan beberapa kesalahan karena sulit untuk mencocokkan poin kontrol. Selain itu, raster netCDF harus dibelah karena tidak memiliki meridian utama yang sama. Misalnya, ekor Antartika di sebelah kanan, sedangkan raster oranye menunjukkannya di sebelah kiri.
multigoodverse
0

Sudut kiri bawah berada pada (0,0) koordinat, tetapi bagi saya (0,0) adalah koordinat sudut kanan atas dalam kasus Anda dan saya tidak tahu sistem koordinat geografis dengan asalnya di -90.

Saat Anda bekerja di ArcGIS, Anda dapat menggunakan bilah georeferensi untuk memindahkan data Anda di lokasi yang tepat. Pindahkan kanan atas ke (0,0) dan kiri bawah ke (-360, -90) dan ini akan menghasilkan transformasi yang tepat tetapi Anda masih akan memiliki masalah dengan dateline.

masukkan deskripsi gambar di sini

Dengan demikian Anda dapat mencoba tautan ini yang akan membantu Anda membuat sistem koordinat geografis khusus. Singkatnya, dalam sistem referensi Anda, buat WGS 84 kustom dengan meridian utama pada -180, lalu transformasi ubahsuaian menggunakan rotasi bujur.

masukkan deskripsi gambar di sini

Georeferensi dalam sistem baru ini akan terdiri dari pengaturan asal pada (-180, -90). Satu titik kontrol sudah cukup karena Anda hanya perlu terjemahan. Pilih pixel kiri bawah dengan hati-hati, klik kanan dan tambahkan nilai XY yang tepat). Maka seharusnya terlihat seperti ini

masukkan deskripsi gambar di sini

Sebagai komentar terakhir, tidak ada kepastian bahwa sistem koordinat didasarkan pada WGS 84. Tetapi sudah pasti sistem koordinat geografis dalam derajat dan resolusi raster Anda mungkin lebih kasar daripada potensi kesalahan karena pemilihan datum yang salah.

radouxju
sumber