Qual sistema de coordenadas é essa varredura do netCDF?

10

Obtive um arquivo de varredura do netCDF, mas não consegui nenhum metadado para obter o nome do sistema de coordenadas em que a varredura foi criada. A varredura em si não possui nenhum sistema de coordenadas incorporado. Eu pensei que seria apenas um WGS84, e parecia que, à primeira vista, mas com uma investigação mais aprofundada no ArcMap, vi que era um sistema pouco comum. Aqui está como ele é exibido: insira a descrição da imagem aqui

A varredura laranja é uma varredura normal no WGS84 que eu inseri aqui para fins de comparação. O roxo é o raster com o sistema de coordenadas desconhecido. Você tem alguma idéia do que possa ser isso?

Algumas atualizações: Aqui está a varredura do netCDF: https://www.dropbox.com/s/nottbl9yt6dwss6/sic_average_nclimate.nc?dl=0 Eu também consegui obter alguns metadados do provedor de imagem:

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

Eles mostram o limite a longo prazo, mas aparentemente nenhuma informação sobre o sistema de coordenadas.

multigoodverse
fonte
11
A Antártica geralmente está na Projeção Estereográfica Polar, então tudo está indo para o norte a partir do pólo sul. nsidc.org/data/polar_stereo/ps_grids.html
Mapperz
11
Estranho, parece 0-360, mas por que ele está tramando no hemisfério norte?
mkennedy
11
Se as extensões relatadas no cabeçalho do NetCDF forem lat / lon, mas forem lat: 0 a?, Lon: 0 a 360. Você pode cortá-las para lat: -90 a (? - 90) e lon: -180 a + 180?
mkennedy
@ Mapperz & mkennedy, tentei atribuir uma projeção estereográfica polar, mas não funcionou. Provavelmente posso espelhá-lo verticalmente para baixo, mas a metade direita ainda estará além das 180E de longitude. Abri o netCDF com o bloco de notas, mas não consegui encontrar nenhuma informação relacionada a lat, long. Estou anexando a varredura na pergunta, se você estiver curioso para abri-la do seu lado.
multigoodverse
Esse link não funciona - o dropbox diz "o proprietário não autorizou o acesso".
BradHards

Respostas:

6

Você pode mudar os dados para a posição correta usando estes comandos 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

(com uma pequena ajuda de Frank Warmerdam: Como reprojetar raster de 0 360 a -180 180 com corte de 180 meridianos )

e o resultado fica assim em uma projeção estereográfica:

insira a descrição da imagem aqui

AndreJ
fonte
1

(Esta não é uma resposta, mas uma sugestão de processo que é muito longa para ser um comentário)

Você pode fazer uma referência geográfica no QGIS usando a ferramenta GeoReferencer, pois parece ser capaz de identificar pontos comuns em ambas as rasters. Isso geraria um script GDAL com um conjunto de pontos GCP e parâmetros de transformação. Você pode adicionar esse script com uma amostra desses pontos à sua pergunta, pois pode explicar qual é a diferença entre as duas projeções.

No pior cenário, você teria uma estimativa aproximada dos valores com base na varredura reposicionada, que pode ajudá-lo com seu objetivo.

djq
fonte
+ 1 e sim, o georreferenciamento pode ser o último marcador, mas o processo geraria alguns erros porque é difícil corresponder aos pontos de controle. Além disso, a varredura do netCDF precisaria ser dividida porque não possui o mesmo meridiano principal. Por exemplo, a cauda da Antártica está à direita, enquanto a varredura laranja mostra à esquerda.
multigoodverse
0

O canto inferior esquerdo está em (0,0) coordenadas, mas para mim (0,0) é a coordenada do canto superior direito no seu caso e não conheço nenhum sistema de coordenadas geográficas com origem em -90.

Enquanto trabalha no ArcGIS, você pode usar a barra de ferramentas de georreferenciamento para mover seus dados no local certo. Mova o canto superior direito para (0,0) e o canto inferior esquerdo para (-360, -90) e isso deve gerar uma transformação exata, mas você ainda terá um problema com a linha de dados.

insira a descrição da imagem aqui

Assim, você pode tentar este link, que deve ajudá-lo a criar um sistema de coordenadas geográficas personalizado. Para resumir, em seu sistema de referência, crie um WGS 84 personalizado com o meridiano principal em -180 e, em seguida, uma transformação personalizada usando rotação de longitude.

insira a descrição da imagem aqui

O georreferenciamento nesse novo sistema consistiria em definir a origem em (-180, -90). Um ponto de controle é suficiente porque você só precisa de uma tradução. Selecione cuidadosamente o pixel inferior esquerdo, clique com o botão direito e adicione os valores XY exatos). Então deve ficar assim

insira a descrição da imagem aqui

Como observação final, não há certeza de que o sistema de coordenadas seja baseado no WGS 84. Mas é com certeza um sistema de coordenadas geográficas em graus e a resolução da sua varredura é provavelmente mais grossa do que o erro potencial devido a uma seleção de dados incorreta.

radouxju
fonte