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