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:
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.
sumber
Jawaban:
Anda dapat menggeser data ke posisi yang benar menggunakan perintah GDAL ini:
(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:
sumber
(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.
sumber
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.
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.
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
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.
sumber