Estou procurando conselhos sobre como trabalhar com conjuntos de dados mundiais nos quais o GCS tem um intervalo positivo (0 a 360 graus). Os dados com os quais estou trabalhando são dados oceanográficos do NetCDF e possuem valores de coordenadas positivos, conforme mencionado. Quando exibido em um GCS WGS84 normal no ArcGIS 10 da ESRI, ele é compensado com outros dados que existem no intervalo de -180 a 180. Se eu o reprojetar, ocorrerão lacunas no meridiano de primeira linha (principalmente devido às suas origens, às vezes até desligado como 25 graus W). Meu pensamento é criar um sistema de projeção / coordenadas personalizado que possua os valores positivos e espero que outros conjuntos de dados mundiais sejam reprojetados rapidamente. Alguém vê algum problema com isso ou tem outras soluções? Espero que isso faça sentido. Agradeço antecipadamente qualquer ajuda.
** Atualização ** Esta é uma captura de tela de um conjunto de dados de países com gratícula na parte superior para mostrar o CRS normal do GCS. Basta importar os dados do NetCDF como estão, mas os valores não começam até 20,5 graus a leste do Prime Meridian.
GCS normal http://grafa.co/rnd/img/GCS_Normal.png
Se eu optar por usar o CRS do NetCDF (que é realmente o mesmo), ele reprojetará os dados mundiais em tempo real, se estiverem definidos. Observe que a gratícula está indefinida para que não seja reprojetada. Tipo de GCS reprojetado http://grafa.co/rnd/img/GCS_Other.png
Mas se eu tentar reprojetar tudo para um GCS normal com valores negativos, é como envolver um mapa plano dos dados em todo o mundo e desaparecer quando atingir o Prime Meridian. Tipo de GCS reprojetado http://grafa.co/rnd/img/reprojected.png
Agora, eu sei que não há valores de 0 a 20,5, conforme declarado nos metadados. Mas por que os valores do restante dos dados não podem ser exibidos na lacuna? Eu até tentei uma mudança nas ferramentas de varredura sem sucesso.
Respostas:
OK, acho que o principal problema aqui é que a ESRI usa o centro da grade para calcular a posição e, quando escrevi os dados, estou usando o valor inferior esquerdo da célula para calcular os dados. Como resolvi isso foi criar meu arquivo netCDF com uma origem da grade X do centro da primeira célula e calcular cada célula usando isso como ponto de partida. Deixe-me uma nota, no bate-papo, se quiser discutir mais isso, pois agora funciona.
fonte
Se necessário, você pode facilmente converter seus rasters para -180 180 com GDAL
http://eloiribeiro.wordpress.com/2009/07/08/coverter-imagens-0%c2%ba-360%c2%ba-a-180-180%c2%ba-com-gdal/
(em português, você precisará usar o Google Translate)
fonte
Este é um bug conhecido agora, pois já enfrentamos o mesmo problema. Com um pouco de folga e descontrole, o que acabou me enviando detalhes de panopólio, para ver os dados do netCDF sem a linha que faltava.
Curiosamente, isso não se limita apenas ao netCDF. Se você criar pontos, usando as grades dentro do netCDF e depois agregar os pontos para criar polígonos, ele também não criará um polígono sobre o meridiano principal.
O erro é NIM070286: a Camada de varredura do NetCDF não exibe dados localizados ao longo do meridiano principal
Solicitei uma correção rápida da ESRI Inc. Vou mantê-la como uma pergunta favorita e atualizá-la quando souber mais. Estou de férias em breve, então espero que isso aconteça nos próximos 10 dias.
fonte