Problemas longitudinales con datos de NetCDF

8

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.

geografa
fuente
geografa: si define una proyección personalizada, ¿no podría establecer el documento en otro crs (uno de los otros conjuntos de datos mundiales) y no exigir que "sobre la marcha". Pero el netcdf lo haría?
Brad Nesom
Por favor, publique una copia de NetCDF en alguna parte. Puede haber una manera de hacer lo que quiere, pero sería más fácil describirlo con un ejemplo real. Esencialmente, podría ser posible hacer algo como en la parte "Tying it All Together" de esta presentación: download.osgeo.org/gdal/presentations/… . O podría ser posible simplemente crear un nuevo archivo mundial con rangos de coordenadas estándar y sobrecargar la georreferenciación incorporada en NetCDF. De cualquier manera, tener una muestra de datos ayudaría.
lagerratrobe
Me di cuenta de que PUEDO reproyectar la trama. Solo tenía que darle el desplazamiento correcto. Publicaré una respuesta más sucinta en los próximos dos días. Para obtener datos, consulte data-portal.ecmwf.int/data/d/edit/personal/temporary/netcdf29 .
geografa
Me interesaría establecer el desplazamiento correcto aquí. Una cosa a tener en cuenta es que, creo, 0 y 360 son lo mismo.
Peludo

Respuestas:

1

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.

Peludo
fuente
0

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.

Peludo
fuente