Parcialmente actualización de ráster (DEM)?

¿Cómo actualizar un DEM .img en ArcView 9.3?

Tengo DEM que necesita actualizarse solo si hay pocas partes. El área que se actualizará se ha producido en un DEM separado con interpolación de spline con barreras.

El área que se actualizará es en forma de S, pero recibí un DEM rectangular más grande (por cierto, esperaba recibir un valor de -9999 fuera de las barreras, pero parece que el área exterior también se interpoló).

Mi objetivo es actualizar solo el área dentro de la polilínea / barreras. ¿Como hacer eso?

Necesita usar un operador condicional . La lógica es:

Result = If {new DEM has values, use them} 
         Else {use values from the original DEM}.

Alternativamente, esto se puede expresar como

Result = If {new DEM does not have values, use original DEM}
         Else {use the new DEM}.

En ArcView, la sintaxis varía con la versión: hay una sintaxis para las versiones hasta la 3.3, otra hasta la 9.3 y una nueva sintaxis en 10.0. (Sin duda, la sintaxis continuará cambiando cada pocos años, por eso es tan importante razonar sobre SIG de forma independiente del software). La expresión alternativa es más simple en ArcView. Su sintaxis 9.3 parece

CON( IsNull([New DEM]), [DEM], [New DEM] )

En términos más generales, suponga que el nuevo DEM incluye algunos valores que no desea usar, pero tiene otra cuadrícula, que representa su polígono, que indica dónde desea que el nuevo DEM reemplace al antiguo. Llamemos a esta cuadrícula poligonal "I". Esta variante se expresaría como

CON( IsNull([I]), [DEM], [New DEM] )

Como puede ver, el procedimiento es bastante general: tiene una cuadrícula "base" (el DEM original), una nueva cuadrícula que contiene valores que desea "parchear" sobre la cuadrícula base y una cuadrícula de "zona" que indica dónde está el parcheo. debiera ser. Una expresión condicional implementa parches de cuadrícula . Este enfoque le permitirá parchear cuadrículas en cualquier SIG basado en ráster que admita álgebra de mapas.