8 votos

Cómo crear un 4096x4096 png de un 1440x720 conjunto de datos utilizando GDAL?

Estoy teniendo un oso de un tiempo tratando de hacer algo que parece bastante sencillo. Estoy tratando de tomar de un archivo netcdf de ftp://eclipse.ncdc.noaa.gov/pub/OI-daily-v2/NetCDF/2011/AVHRR-AMSR/ y convertirlo en un 4096x4096 png para que yo pueda mosaico para usar con google maps. He estado usando los siguientes comandos pero siempre parece crear zonas en blanco en el este y oeste de las fronteras:

/usr/local/bin/gdal_translate -a_srs EPSG:4326 -outsize 4096 4096 -a_ullr 0 90 360 -90 NETCDF:"input.nc":sst -of netCDF output.nc
/usr/local/bin/gdaldem color-relief output.nc -b 1 -alpha colorscale.txt -of PNG output.png
/usr/bin/convert output.png -transparent #FFFFFF output.png

Soy sólo se aproxima por el camino equivocado? Estoy utilizando las herramientas equivocadas?

FWIW - gdal_translate es v1.8.0

Gracias de antemano.

ACTUALIZACIÓN: "Responder" debajo de Dija - colega - es en realidad una respuesta a las observaciones sobre la pregunta en sí misma.

ACTUALIZACIÓN Accidentalmente agregó una respuesta en lugar de un comentario. Mover aquí.

He pensado que me gustaría opinar aquí. He estado trabajando en este problema por un poco de tiempo ahora. He intentado muchas iteraciones del proceso. Yo escala de estas imágenes a la 4096x4096 ya que en realidad se crea una mejor imagen de si no. Es decir, si tengo que dejar la imagen en su tamaño original, 1440x720 puedo obtener una brecha mayor en el primer meridiano (para este conjunto de datos) y en el meridiano 180º para otros conjuntos de datos. Yo escala de la imagen de salida para ser un múltiplo de 256 (para que se ajuste mejor con google) y obtener diferentes resultados.

Si yo escala a 1536x768 la brecha se muestra a la derecha del meridiano 0. Si yo escala para 1280x640 es a la izquierda. Si lo dejo solo, es bastante divide.

Me he tomado la imagen base que gdaldem escupe, cortar por la mitad y coincidir los bordes hasta...que se alinean perfectamente. Esto me lleva a pensar que hay algo raro en gdal2tiles.py

Yo escribí mi propio azulejera que me estoy preparando para activar sueltos en nuestras imágenes para ver si la brecha sigue siendo.

Voy a poner algunas capturas de pantalla para mostrar lo que estoy hablando acerca de cómo la brecha así como el informe sobre lo que mi azulejera.

Gracias por la entrada que usted pueda tener. Es apreciado.

ACTUALIZACIÓN 2 me encontré con la maptiler me escribió en contra de la trama original y funciona bien con una brecha. gdal2tiles.py parece ser el problema. Presenta un borde transparente a la última ficha que toca el meridiano 180º.

2voto

AimFirst Puntos 141

gdal_grid tiene parámetros que permiten designar a la medida de las coordenadas de su ráster de salida.

gdal_grid -txe [xcoord1] [xcoord2] -tye [ycoord1] [ycoord2]

donde la parte inferior izquierda de coordenadas sería (xcoord1,ycoord1) y la parte superior derecha sería (xcoord2,ycoord2)

GDAL va a interpolar los datos de acuerdo a su salida de las extensiones

i-Ciencias.com

I-Ciencias es una comunidad de estudiantes y amantes de la ciencia en la que puedes resolver tus problemas y dudas.
Puedes consultar las preguntas de otros usuarios, hacer tus propias preguntas o resolver las de los demás.

Powered by:

X