特定のDEMからブレークラインを抽出する方法は？

ブレークラインは、たとえば正確なTINを構築するために、地形を記述するのに非常に役立ちます。

私の質問は、これらのブレークラインを特定のDEMから抽出する方法ですか？

曲率（垂直）-ラスター（arcgis 10を使用）を計算しました。ブレークライン（最小/最大曲率）を確認できますが、それらを抽出して最終的にポリラインを取得する方法はわかりませんでした。

これらのブレークラインをDEMから取得するためのワークフローは何ですか？

誰もが役立つワークフロー、または「ツール」、またはskriptのアイデアを持っていますか？

私はこのタスクを実行したい最初/最後の人ではないと思います:)

...その間、私はこのトピックに関する論文を見つけましたが、それは良い解決策のようですが、大きな問題が残っています：それを実践に減らす方法は？

または、誰かがいくつかのソフトウェアで異なる（より簡単な？）アプローチを持っていますか？

私はこの方法で試してみます。曲率ラスターファイルで重要な値を検索します。この値を超えると、ブレークラインが予想されます。

次に、ラスター計算機の条件ステートメントを使用して、臨界値を超えるピクセルと臨界値を下回る0のピクセルに「1」を使用します。ツールボックスのラスターからポリゴンへの変換には、このバイナリラスターを使用します。ここにポリゴンを線で結ぶ方法の答えがあります

結局のところ、これらの手順からモデルを構築できます。

Conrecはこれを行います。もともとはFortranで書かれていたため、何度か移植されています。

私がnode.jsで作業しているライブラリturfは、このアルゴリズムを使用して、ポイントベクトルデータに基づいてこれらの輪郭を作成します。これを行うスクリプトは非常に単純で、基本的に次のようになります。

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)
  })
})

特定の休憩がない場合は、ジャンク自然休憩を生成することもできます。

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 ]
  })
})

私のVA LIDAR サイトにアクセスして、ブレークラインのあるプロジェクトを探したとします。たくさんあります。

たとえば、2010/2011 Eleven County Coastal VA Lidarです。すべてのレポートおよびその他のドキュメントへのリンクが表示され、通常はArcGISでブレークラインが作成された方法が説明されます。いくつかの異なるプロジェクトをさまざまな方法で試してみてください。ただし、DEMから派生したリッジのように、NHDがよく使用されていることがわかります。ほとんどのレポートはメタデータに含まれています。ここで1以上のため。PDFのCTRL Fブレークライン。