¿Hay alguna diferencia entre proyectar shapefiles luego rasterizar y rasterizar shapefiles luego proyectar?

8

Actualmente tengo un proyecto en el que rasterizo archivos de formas en rásteres de presencia / ausencia simples (1 y 0) y luego agrego cada uno a un total acumulado para crear un ráster de riqueza total al final. Los archivos de forma están en proyección WGS84, pero necesito que el ráster final esté en la proyección de Área igual de Behrmann.

Lo que quiero saber es si hay alguna diferencia en el ráster final si proyecto todos los archivos de forma a Behrmann y luego rasterizo y fusiono cada uno de ellos en comparación con rasterizar y fusionar, y luego proyectar el ráster final a Behrmann al final.

JPD
fuente
Gracias por la información chicos. Me he decidido a proyectar los vectores antes de la rasterización.
JPD

Respuestas:

5

La respuesta corta es: todas las cosas son iguales, rasterizar y luego transformar.

La respuesta larga depende de varios factores:

  • El tamaño de la función.
  • El error permitido
  • Si los CRS están de acuerdo en conformidad y en qué grado
  • ¿Qué elipsoide y dato se utilizan?
  • Limitaciones de tiempo de procesamiento

El problema es que las líneas rectas en un CRS no serán necesariamente rectas en otro, a menos que ambas sean conformes, por lo que la transformación de un segmento de línea posicionará sus puntos finales con precisión, pero no se agregarán nuevos puntos intermedios, por lo que la línea en su conjunto no será exacto.

En su caso, si transforma y rasteriza, y las características son de tamaño pequeño, generalmente cubren menos de un kilómetro cuadrado, no habrá un error apreciable en ninguna resolución razonable de la trama. Pero a medida que avanza hacia los polos, o cuando sus características se hacen más grandes, notará un mayor error. Por lo tanto, llegará un punto en el que la precisión de sus datos de origen es mejor que la precisión de sus datos proyectados: su error permitido.

Computacionalmente, generalmente es más rápido transformar los datos vectoriales antes de rasterizarlos, por lo que tendría que sopesar eso si tiene alguna limitación de tiempo.

Si tiene tiempo para experimentar, generaría dos rásteres usando un subconjunto de sus datos, uno usando el proceso de transformación / rasterización y otro al revés. Luego combine los dos mapas de bits con una operación exclusiva u, que resaltará las áreas donde hay un error de más de un píxel, y puede juzgar si se permite algún error. También puede comparar el tiempo que lleva procesar para informar mejor su decisión.

MerseyViking
fuente
44
Su argumento parece implicar exactamente lo contrario de su conclusión. Con una representación vectorial, los buenos SIG pueden dividir segmentos y transformarlos con precisión, pero una vez que los datos están rasterizados, no es posible realizar tales mejoras. Por lo tanto, se deduce que uno debe rasterizar el último , ¡no el primero!
whuber
¡Maldice tu sofisma! :) En mi SIG mental, que se basa en PostGIS / QGIS, los segmentos de línea nunca tienen vértices adicionales agregados, a menos que el usuario los agregue de antemano. Sin embargo, no puedo hablar por "buenos" SIG: p
MerseyViking
44
En ese caso, el orden de las operaciones debería hacer poca diferencia, excepto el remuestreo que necesariamente ocurre cuando se transforma el ráster. Si el remuestreo realiza una interpolación (a menudo es bilineal), entonces los rásteres reproyectados pueden tener algunos artefactos de interpolación a lo largo de los límites de la característica. Reproyectar las formas vectoriales y luego rasterizar no solo tenderá a ser computacionalmente más eficiente, sino que también evitará tales artefactos.
whuber