Come estrarre le breakline da un dato DEM?

Le breakline sono molto utili per descrivere un terreno, ad esempio per costruire una TIN esatta.

La mia domanda è: come estrarre queste breakline da un dato DEM?

Ho calcolato la curvatura (verticale) - raster (usando arcgis 10), dove puoi vedere le linee di discontinuità (min / max-curvatura), ma non ho capito come estrarli per ottenere le polilinee alla fine.

Qual è il flusso di lavoro per ottenere queste breakline da un DEM?

Qualcuno ha un'idea di un flusso di lavoro, di uno "strumento" o di uno script, che può essere d'aiuto?

Penso di non essere il primo / ultimo uomo che vuole svolgere questo compito :)

... nel frattempo ho trovato un documento su questo argomento che sembra essere una buona soluzione, ma la grande domanda rimane: come ridurla alla pratica?

O qualcuno ha un approccio diverso (più semplice?) Con alcuni software?

Vorrei provarlo in questo modo: cercare un valore critico nel file raster di curvatura. Al di sopra di questo valore ci si aspetta una linea di discontinuità.

Quindi userei le istruzioni condizionali nel calcolatore raster per usare "1" per i pixel sopra il valore critico e 0 sotto il valore critico. Utilizzare questo binario raster per la conversione da raster a poligono nella casella degli strumenti. troverai una risposta su come far allineare il poligono qui

dopo tutto puoi costruire un modello da questi passaggi.

Conrec lo fa. È stato originariamente scritto in fortran, quindi è stato portato in giro alcune volte.

Una libreria su cui sto lavorando in node.js, turf , utilizza questo algoritmo per creare questi contorni in base a dati vettoriali punti. Lo script per farlo è piuttosto semplice e sostanzialmente sarebbe:

var t = require('turf')
var fs = require('fs')
var z = 'elevation'
var resolution = 15
var breaks = [.1, 22, 45, 55, 65, 85,  95, 105, 120, 180]

t.load('./path/to/points.geojson', function(err, points){
  t.contour(points, z, resolution, breaks, function(err, contours){
    if(err) throw err
    fs.writeFileSync('.path/to/output.geojson', contours)
  })
})

Puoi anche generare interruzioni naturali jenks se non hai in mente interruzioni specifiche:

var t = require('turf')
var propertyName = 'elevation'
var num = 10

t.load('./path/to/points.geojson', function(err, pts){
  if(err) throw err
  t.jenks(pts, 'elevation', num, function(err, breaks){
    if(err) throw err
    console.log(breaks) // [ 11, 12, 18, 25, 29, 41, 50, 55, 76, 90, 143 ]
  })
})

Se vai sul mio sito lidar VA e cerchi un progetto con Breaklines. Ne abbiamo molti.

Ad esempio 2010/2011 Eleven County Coastal VA Lidar. Vedrai i collegamenti a tutti i rapporti e altri documenti e questo spiega come hanno fatto le discontinuità di solito in ArcGIS. Prova alcuni progetti diversi per metodi diversi, ma posso dirti che il NHD viene spesso utilizzato, così come le creste derivate da un DEM. Principalmente i report sono nei metadati. Ecco quello per sopra. CTRL F Breaklines in pdf.