Estoy buscando asesoramiento para trabajar con conjuntos de datos mundiales donde el GCS tiene un rango positivo (0 grados a 360 grados). Los datos con los que estoy trabajando son datos oceanográficos de NetCDF y tienen valores de coordenadas positivos como se mencionó. Cuando se muestra en un GCS WGS84 normal en ArcGIS 10 de ESRI, se compensa con otros datos que existen en el rango de -180 a 180. Si lo vuelvo a proyectar, se producen huecos alrededor del meridiano principal (principalmente debido a sus orígenes, a veces tan lejos) apagado como 25 grados W). Mi idea es crear un sistema de proyección / coordenadas personalizado que tenga valores positivos y espero que otros conjuntos de datos mundiales se vuelvan a proyectar sobre la marcha. ¿Alguien ve algún problema con esto o tiene otras soluciones? Espero que esto tenga sentido. Agradezco cualquier ayuda por adelantado.
** Actualización ** Esta es una captura de pantalla del conjunto de datos de un país con una retícula en la parte superior para mostrar el CRS normal de GCS. Solo importando los datos de NetCDF tal como están, los valores no comienzan hasta 20.5 grados al este del primer meridiano.
GCS normal http://grafa.co/rnd/img/GCS_Normal.png
Si elijo usar el CRS de NetCDF (que en realidad es el mismo), reproyectará los datos del mundo sobre la marcha si está definido. Tenga en cuenta que la retícula no está definida, por lo que no se reproyecta. Una especie de GCS reproyectado http://grafa.co/rnd/img/GCS_Other.png
Pero luego, si trato de reproyectar todo a un GCS normal con valores negativos, es como envolver un mapa plano de los datos en todo el mundo y luego desaparece cuando llega al Meridiano principal. Una especie de GCS reproyectado http://grafa.co/rnd/img/reprojected.png
Ahora, sé que no hay valores de 0-20.5 como se indica en los metadatos. Pero, ¿por qué los valores del resto de los datos no se muestran en la brecha? Incluso probé un cambio en las herramientas ráster en vano.
Respuestas:
OK, creo que el problema principal aquí es que ESRI usa el centro de la cuadrícula para calcular la posición, y cuando he escrito los datos, estoy usando el valor inferior izquierdo de la celda para calcular los datos. La forma en que resolví esto fue creando mi archivo netCDF con un origen de cuadrícula X del centro de la primera celda, luego calculando cada celda usando eso como punto de partida. Envíenme una nota, en el chat, si desea discutir esto más a fondo, ya que ahora funciona.
fuente
Si es necesario, puede convertir fácilmente sus rásteres a -180 180 con GDAL
http://eloiribeiro.wordpress.com/2009/07/08/coverter-imagens-0%c2%ba-360%c2%ba-a-180-180%c2%ba-com-gdal/
(en portugués, deberá utilizar el Traductor de Google)
fuente
Este es un error conocido ahora, ya que hemos experimentado el mismo problema. Con un poco de tooing y fromming, que terminó conmigo enviándoles detalles de panopolia, para ver los datos de netCDF sin la línea faltante.
Curiosamente, esto no se limita únicamente al netCDF. Si crea puntos, usando las cuadrículas dentro del netCDF, y luego agrega los puntos para construir polígonos, no construirá un polígono sobre el meridiano principal también, serán separados.
El error es NIM070286: NetCDF Raster Layer no muestra datos ubicados a lo largo del meridiano principal
He solicitado un parche de ESRI Inc. Lo mantendré como una pregunta favorita y lo actualizaré todo una vez que tenga más información. Estoy de vacaciones pronto, así que espero que esto ocurra en los próximos 10 días.
fuente