Я получил растровый файл 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" ;
Они показывают границу в продольном направлении, но, по-видимому, не имеют никакой информации относительно системы координат.
источник
Ответы:
Вы можете переместить данные в правильную позицию, используя эти команды GDAL:
(с небольшой помощью Фрэнка Вармердама: как перепроектировать растр от 0 360 до -180 180 с обрезкой 180 меридианов )
и результат выглядит следующим образом в стереографической проекции:
источник
(Это не ответ, а скорее предложение процесса, которое слишком длинное, чтобы быть комментарием)
Вы можете геопривязать его в QGIS, используя инструмент GeoReferencer, так как вы, похоже, сможете определить общие точки на обоих растрах. Это создаст сценарий GDAL с набором точек GCP и параметрами преобразования. Вы можете добавить этот сценарий с образцом этих моментов к вашему вопросу, поскольку он может объяснить разницу между двумя проекциями.
В худшем случае у вас будет приблизительная оценка значений на основе измененного растра, который может помочь вам в достижении вашей цели.
источник
Нижний левый угол находится в координатах (0,0), но для меня (0,0) это координата верхнего правого угла в вашем случае, и я не знаю какой-либо географической системы координат с источником в -90.
Работая в ArcGIS, вы можете использовать панель инструментов с географической привязкой для перемещения данных в нужное место. Переместите верхний правый угол (0,0), а нижний левый - (-360, -90), и это должно привести к точному преобразованию, но у вас все равно будет проблема с датой.
Таким образом, вы можете попробовать эту ссылку, которая должна помочь вам создать собственную географическую систему координат. Чтобы подвести итог, в вашей системе отсчета создайте пользовательский WGS 84 с начальным меридианом -180, затем пользовательское преобразование с использованием вращения по долготе.
Геопривязка в этой новой системе будет заключаться в установке начала координат в (-180, -90). Достаточно одной контрольной точки, потому что вам нужен только перевод. Тщательно выберите нижний левый пиксель, щелкните правой кнопкой мыши и добавьте точные значения XY). Тогда это должно выглядеть так
И последнее замечание: нет уверенности в том, что система координат основана на WGS 84. Но это, безусловно, географическая система координат в градусах, и разрешение вашего растра, вероятно, более грубое, чем потенциальная ошибка из-за неправильного выбора данных.
источник