Todo depende de dónde trazar la línea. Independientemente, este problema parece que puede ser fácilmente abordado con la morfológicas de las funciones disponibles en arcgis Spatial Analyst, especialmente umbralización (realizado con "<" y ">" operaciones locales) y "RegionGroup" para identificar y extraer los componentes.
A pesar de que no tienen acceso a la DEM para ilustrar, la primera imagen es casi lo suficientemente bueno para hacer el trabajo de todos modos. Por ejemplo, aquí es una secuencia de regiongrouped versiones de la banda roja (con valores de escala de 0.0=negro a 1.0=blanco), a partir de un umbral de 0 y se mueve de izquierda a derecha, de arriba a abajo, en incrementos de 0.02:
![Dikes]()
(Estas imágenes se encogió de reproducción aquí: todo el análisis fue hecho en la resolución de la imagen original.)
Los diques surgen alrededor de la frontera de la región negra por el inicio de la segunda fila (umbral = 0.08). En el inicio de la tercera fila (umbral = 0.16) los diques de formar sus propios componentes (en azul oscuro) y en ese punto puede ser fácilmente extraído por separado como rejillas o polígonos (y sus límites pueden ser extraídos como polilíneas después de una detección de los bordes de paso). Sólo la más amplia diques, siendo el último de la fila (umbral = 0,24 y mayor). Usted tendrá que seleccionar un umbral adecuado con el fin de obtener exactamente lo que usted desea.
En el original de la DEM, elevaciones de jugar el papel de intensidad, por lo que estos procedimientos deben ser igualmente eficaz con los DEM sí mismo. Si el DEM tiene un gran alcance (incluyendo características muy lejos del río o diques), las características pertinentes puede ser seleccionado como aquellas adyacentes a la componente en el cual el río se cuentan mentiras.
Delgado no deseados regiones que se encuentran por RegionGroup puede ser eliminada mediante la aplicación de una pequeña erosión (búfer negativo), seguido por un equivalente de la dilatación positiva (buffer). Los pequeños no deseados regiones pueden ser excluidos en un criterio de tamaño (recuento celular total o área). Los puntos más altos de cada dique (si eso es lo que realmente es necesario) se puede encontrar mediante la comparación de la DEM a zonal máximo cuadrículas (el uso de los diques como zonas).
Por cierto, esta secuencia de imágenes fue producido en Mathematica 8. Aquí están los comandos para aquellos que deseen realizar esta opción.
original = Import["http://i.stack.imgur.com/gV7Du.jpg"];
{r, g, b} = ColorSeparate[original];
frames = ParallelTable[Colorize[MorphologicalComponents[Binarize[r, t/100]],
ColorFunction -> "ThermometerColors"], {t, 0, 30, 2}];
Rasterize @ TableForm[Partition[frames, 4]]
Total de tiempo de reloj (después de importar la imagen) fue de 0.94 segundos, de los cuales la mitad era necesario para reducir y exportación de todos los 16 de imágenes: morfológicas de las operaciones tienden a ser rápido (lo cual es bueno, porque LIDAR DEMs puede ser enorme).
