Jakim układem współrzędnych jest ten raster netCDF?

10

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: wprowadź opis zdjęcia tutaj

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.

multigoodverse
źródło
1
Antarktyda jest zwykle w polarnej projekcji stereograficznej, więc wszystko biegnie na północ od bieguna południowego. nsidc.org/data/polar_stereo/ps_grids.html
Mapperz
1
Dziwne, wygląda na 0-360, ale dlaczego knuje na półkuli północnej?
mkennedy
1
Jeśli zgłaszane zakresy w nagłówku NetCDF wyglądają jak lat / lon, ale są lat: 0 do?, Lon: 0 do 360. Czy możesz je zhakować do lat: -90 do (? - 90) i lon: -180 do + 180?
mkennedy
@Mapperz & mkennedy, próbowałem przypisać mu projekcję stereograficzną Polar, ale to nie zadziałało. Prawdopodobnie mogę to zrobić pionowo w dół, ale jego prawa połowa będzie nadal znajdować się poza długością 180E. Otworzyłem netCDF za pomocą notatnika, ale nie mogłem znaleźć żadnych informacji związanych z lat, długo. Dołączam raster do pytania, jeśli chcesz go otworzyć po swojej stronie.
multigoodverse
Ten link nie działa - Dropbox mówi „właściciel nie autoryzował dostępu”.
BradHards,

Odpowiedzi:

6

Możesz przesunąć dane do właściwej pozycji za pomocą następujących poleceń 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

(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:

wprowadź opis zdjęcia tutaj

AndreJ
źródło
1

(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.

djq
źródło
+ 1 i tak, georeferencje mogą być ostatnią kulą, ale proces przyniósłby pewne błędy, ponieważ trudno jest dopasować punkty kontrolne. Ponadto raster netCDF musiałby zostać podzielony, ponieważ nie ma tego samego głównego południka. Na przykład ogon Antarktydy znajduje się po prawej stronie, a pomarańczowy raster pokazuje go po lewej.
multigoodverse
0

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.

wprowadź opis zdjęcia tutaj

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.

wprowadź opis zdjęcia tutaj

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ć

wprowadź opis zdjęcia tutaj

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.

radouxju
źródło