Otrzymałem plik rastrowy netCDF, ale nie mogłem uzyskać żadnych metadanych, aby uzyskać nazwę układu współrzędnych, na którym raster został zbudowany. Sam raster nie ma wbudowanego żadnego układu współrzędnych. Myślałem, że będzie to po prostu WGS84 i na pierwszy rzut oka tak to wyglądało, ale po dalszych badaniach w ArcMap zauważyłem, że jest to raczej nietypowy system. Oto jak to wyświetla:
Pomarańczowy raster jest normalnym rastrem w WGS84, który tu wstawiłem w celach porównawczych. Fioletowy to raster z nieznanym układem współrzędnych. Czy masz jakieś pojęcie, co to może być?
Niektóre aktualizacje: Oto raster netCDF: https://www.dropbox.com/s/nottbl9yt6dwss6/sic_average_nclimate.nc?dl=0 Udało mi się również uzyskać metadane od dostawcy obrazu:
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" ;
Pokazują one długość, ale najwyraźniej nie zawierają żadnych informacji dotyczących układu współrzędnych.
źródło
Odpowiedzi:
Możesz przesunąć dane do właściwej pozycji za pomocą następujących poleceń GDAL:
(z niewielką pomocą Franka Warmerdama: Jak zmienić projekt rastra z 0 360 na -180 180 za pomocą cięcia 180 południków )
a wynik wygląda tak w projekcji stereograficznej:
źródło
(To nie jest odpowiedź, ale sugestia procesu, która jest zbyt długa, aby być komentarzem)
Możesz georeferencję w QGIS za pomocą narzędzia GeoReferencer, ponieważ wydaje się, że jesteś w stanie zidentyfikować wspólne punkty na obu rastrach. Wygeneruje to skrypt GDAL z zestawem punktów GCP i parametrami transformacji. Możesz dodać ten skrypt z próbką tych punktów do swojego pytania, ponieważ mogłoby to wyjaśnić różnicę między tymi dwiema projekcjami.
W najgorszym przypadku będziesz miał przybliżony szacunek wartości na podstawie zmienionego położenia rastra, który może ci pomóc w osiągnięciu celu.
źródło
Lewy dolny róg ma współrzędne (0,0), ale dla mnie (0,0) to współrzędna prawego górnego rogu w twoim przypadku i nie znam żadnego układu współrzędnych geograficznych z początkiem -90.
Podczas pracy w ArcGIS możesz użyć paska narzędzi do georeferencji, aby przenieść dane we właściwe miejsce. Przesuń prawy górny do (0,0) i dolny lewy do (-360, -90), co powinno dać dokładną transformację, ale nadal będziesz miał problem z linią danych.
Możesz zatem wypróbować ten link, który powinien pomóc w utworzeniu niestandardowego układu współrzędnych geograficznych. Podsumowując, w systemie odniesienia utwórz niestandardowy WGS 84 z głównym południkiem na -180, a następnie niestandardową transformację za pomocą obrotu długości geograficznej.
Georeferencje w tym nowym systemie polegałyby na ustawieniu źródła na (-180, -90). Jeden punkt kontrolny wystarczy, ponieważ potrzebujesz tylko tłumaczenia. Ostrożnie wybierz lewy dolny piksel, kliknij prawym przyciskiem myszy i dodaj dokładne wartości XY). To powinno tak wyglądać
Na koniec, nie ma pewności, że układ współrzędnych oparty jest na WGS 84. Jest to jednak na pewno układ współrzędnych geograficznych w stopniach, a rozdzielczość twojego rastra jest prawdopodobnie z grubsza niż potencjalny błąd z powodu niewłaściwego wyboru układu odniesienia.
źródło