Яка система координат це растр netCDF?


10

Я отримав растровий файл netCDF, але не зміг отримати метадані, щоб отримати ім'я системи координат, на якій був побудований растр. Сам растр не має вбудованої системи координат. Я думав, що це буде просто WGS84, і на перший погляд це виглядало так, але при подальшому дослідженні в ArcMap я побачив, що це досить поширена система. Ось як це відображається: введіть тут опис зображення

Помаранчевий растр - це звичайний растр у WGS84, який я вставлю сюди для порівняння. Фіолетовий - це растр із невідомою системою координат. Чи є у вас поняття, що це може бути?

Деякі оновлення. Ось растр netCDF: https://www.dropbox.com/s/nottbl9yt6dwss6/sic_average_nclimate.nc?dl=0 Мені також вдалося отримати метадані від постачальника зображень:

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" ;

Вони показують межу широкої ширини, але, мабуть, не мають жодної інформації щодо системи координат.


1
Антарктида зазвичай знаходиться в Полярній стереографічній проекції, тому все йде на північ від південного полюса. nsidc.org/data/polar_stereo/ps_grids.html
Mapperz

1
Дивно, схоже на 0-360, але чому це побудовано в північній півкулі?
mkennedy

1
Якщо згадані розширення у заголовку NetCDF виглядають як lat / lon, але вони є lat: 0 to?, Lon: 0 to 360. Чи можете ви зламати їх у lat: -90 to (? - 90) та lon: -180 to + 180?
mkennedy

@Mapperz & mkennedy, я спробував призначити йому Полярну стереографічну проекцію, але це не вийшло. Я, мабуть, можу відобразити його вертикально вниз, але його права бічна половина все ще буде понад 180E довготи. Я відкрив netCDF із блокнотом, але не міг знайти жодної інформації, пов’язаної з lat, довго. Я додаю растр у питанні, якщо вам цікаво відкрити його на боці.
багатоповерховий

Це посилання не працює - у вікні "Dropbox" написано "власник не має права доступу"
BradHards

Відповіді:


6

Ви можете перемістити дані у правильне положення за допомогою цих команд GDAL:

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

(за невеликої допомоги Френка Вармердама: Як відремонтувати растр від 0 360 до -180 180 з вирізанням 180 меридіана )

і результат виглядає таким чином у стереографічній проекції:

введіть тут опис зображення


1

(Це не відповідь, а швидше пропозиція щодо процесу, яка занадто довга, щоб коментувати)

Ви можете зробити геопосилання на QGIS за допомогою інструмента GeoReferencer, оскільки, здається, зможете визначити загальні точки на обох растрах. Це створило б сценарій GDAL з набором точок GCP та параметрами перетворення. Ви можете додати цей скрипт із зразком цих пунктів до свого питання, оскільки це може пояснити, у чому різниця між двома прогнозами.

Найгірший випадок, ви мали б приблизну оцінку значень на основі репозиціонованого растру, який може допомогти вам досягти вашої мети.


+ 1 і так, геореференція може бути останньою кулею, але процес призведе до деяких помилок, оскільки складно зіставити контрольні точки. Крім того, растр netCDF повинен був би бути розділений, оскільки він не має однакового меридіана. Наприклад, хвіст Антарктиди знаходиться праворуч, а помаранчевий растр - зліва.
multigoodverse

0

У нижньому лівому куті розміщено координати (0,0), але для мене (0,0) - це координата верхнього правого кута у вашому випадку, і я не знаю жодної географічної системи координат з початком в -90.

Під час роботи в ArcGIS ви можете використовувати панель інструментів геореференції для переміщення своїх даних у потрібному місці. Перемістіть верхній правий кут до (0,0), а нижній лівий - (-360, -90), і це має призвести до точного перетворення, але у вас все-таки виникнуть проблеми з дателіном.

введіть тут опис зображення

Таким чином, ви можете спробувати це посилання, яке допоможе вам створити власну географічну систему координат. Підсумовуючи це, у вашій системі відліку створіть спеціальний WGS 84 з основним меридіаном у -180, а потім власну трансформацію за допомогою обертання довготи.

введіть тут опис зображення

Геореференція у цій новій системі полягала б у встановленні джерела (-180, -90). Одного контрольного пункту достатньо, тому що вам просто потрібен переклад. Уважно виберіть нижній лівий піксель, клацніть правою кнопкою миші та додайте точні значення XY). Тоді це повинно виглядати так

введіть тут опис зображення

На завершення, немає впевненості, що система координат базується на WGS 84. Але це точно географічна система координат у градусах, і роздільна здатність вашого растру, ймовірно, грубіша, ніж потенційна помилка через неправильний вибір дат.

Використовуючи наш веб-сайт, ви визнаєте, що прочитали та зрозуміли наші Політику щодо файлів cookie та Політику конфіденційності.
Licensed under cc by-sa 3.0 with attribution required.