5 votos

¿Reproyección de raster con FME?

Tengo varios archivos .tiff + .tfw en EPSG:3067, que necesito reproyectar a EPSG:3878. He intentado hacerlo con FME, y lo he conseguido, más o menos. Los geotiffs se colocan bien, pero aparece una extraña zona blanca entre la cuadrícula (ver imagen adjunta) cuando añado las ortofotos a ArcMap 10.4. Utilicé CsmapReprojector, se seleccionó "ignorar alturas" y el tipo de interpolación fue Bicúbica y el Tamaño de Celda fue "Preservar Celdas". También definí manualmente las coordenadas de origen y destino a 3067 y 3878.

Entonces, ¿qué hice mal mientras 3878 reproyección se ve así:

enter image description here

Mientras que debería verse así:

enter image description here

4 votos

Mosaico en proyección original, resultado de la re proyección

6voto

DBM Puntos 51

De acuerdo con el comentario de @FelixIP, una posible solución es crear primero una trama completa utilizando RasterMosaicker transformador, luego reproyectar y volver a embalar el Raster con RasterTiler transformador. Este método puede consumir muchos recursos si tiene muchos mosaicos de entrada. Otro enfoque es aplicar los valores de los nódulos con la función RasterBandNodataSetter (Valor = 0) o RasterPaletteNodataSetter pero esto depende del modelo de color y del formato de los datos de salida.

enter image description here

El mosaico y la proyección del mosaico (y el reajuste y la aplicación del valor de los nódulos) dan como resultado

enter image description here

donde en la parte superior aplicando nodata no funcionaba correctamente, porque FME crea una banda alfa en esta zona por alguna razón.

Si se proyectan las baldosas directamente y se aplican los valores de los nódulos, se obtiene

enter image description here

donde las líneas entre los mosaicos son artefactos de visualización en el inspector de datos de FME en este caso.

0 votos

Hola, ¡muchas gracias por su respuesta! Me ha ayudado mucho. Tengo múltiples rasters, así que primero probé el "RasterBandNodataSetter". ¡Primero no funcionó, cuando se añadió después del reproyector, pero luego traté de añadir RasterPaletteNodataSetter también, justo después de RasterBandNodataSetter, y funcionó como un encanto! Ya no hay áreas blancas y reproyecta los rásteres a una velocidad adecuada.

2voto

Steven Puntos 11

No puedo añadir la imagen al comentario, así que añadiré aquí cómo pude tener éxito con la ayuda de Jochen Schwarzes.

enter image description here

He probado esto con 11 rasters, que pasaron bien. El resultado de esta reproyección en Arcmap se ve bien.

2voto

Ulf Tennfors Puntos 11

Ya que tienes ArcGIS, podrías hacer un conjunto de datos en mosaico con los rásters, y usar el conjunto de datos en mosaico como entrada a FME. Esta solución funciona bien para mí con grandes conjuntos de datos.

Ed: ArcGIS tiene la capacidad de reproyectar el conjunto de datos del mosaico sobre la marcha. Esto significa que usted tiene la opción de elegir si la reproyección debe ser hecha por ArcGIS o FME. Con un gran conjunto de datos, es probable que ArcGIS pueda reproyectar el conjunto de datos más rápido que FME, que ve el conjunto de datos como un gran raster.

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