問題タブ [qhull]

Rで非凸面をプロットする簡単な方法はありますか?

凸面に次のようなものを使用しましたが、正常に機能します。

ただし、非凸面の場合、凹面は塗りつぶされます。これは、拘束された Delaunay三角形分割や非凸オブジェクトのメッシュ生成をサポートしないQhulldelaunayn()ライブラリを使用しているため、関数の問題だと思います。

任意の提案をいただければ幸いです。

PS

ASCII ファイルとしてデータを持っていますが、3 列で 225 行の長さです。これを提供する最良の方法は何ですか？

入手可能なデータ: http://pastebin.com/R2p4Cf7d

プロットの上部は凹んでいる必要があります! persp3d()これは、サーフェスがどのように見えるかを使用して作成されたイメージです。これは、不規則なコロケーション ポイントを使用するのではなく、極座標の規則的なグリッド上でより多くのグリッド ポイントを使用して計算されています。

QHullの SciPy ラッパーを使用して、一連の点の凸包のボリュームを取得しようとしています。

QHullのドキュメントによると"FA"、総表面積と体積を取得するオプションを渡す必要があります。

ここに私が得るものがあります..何が間違っていますか?

10 次元空間で 9 点の凸包を計算しようとしています。scipy インターフェイスを介して、私は呼び出しscipy.spatial.ConvexHull(points)て取得していますQH6214 qhull input error: not enough points(9) to construct initial simplex (need 12)

次元に関係なく、凸包の定義は明確に定義されていると思います。ここで何が起こっているのですか？これを修正するために呼び出すことができる別の関数はありますか?

Windows 8.1 の VS2012 で PCL を構築しようとしているので、Kinect for Windows v2 センサーをそれに統合しようとしていますが、どうすればよいか分からない一連の CMake エラーが発生しているようです。これは部分的には、私が CMake にあまり詳しくないためです。そのため、CMake GUI を使用して依存関係を構築する方法を概説するこちらのチュートリアルに従っています。ただし、このプロセスを通じて、適切にビルドされたものは、コマンド ラインで cmake を使用してビルドしたものだけであることがわかりました。エラーは次のとおりです。

フランの場合:

VS2012 で INSTALL をビルドしようとすると (これは他のビルド エラーにはなりませんが、ここから始めましょう)、次のようになります。

エラー 21 エラー C2228: '.serialize' の左側にはクラス/構造体/共用体が必要です D:\flann-1.8.4-src\flann-1.8.4-src\src\cpp\flann\util\serialization.h 行 18 : type.serialize(ar); エラー 24 エラー LNK1181: 入力ファイル '..\lib\Release\flann.lib' を開けません D:\flann-1.8.4-src\flann-1.8.4-src\build\examples\LINK エラー 25 エラー MSB3073:コマンド "setlocal C:\CMake\bin\cmake.exe -DBUILD_TYPE=Release -P cmake_install.cmake if %errorlevel% neq 0 goto :cmEnd :cmEnd endlocal & call :cmErrorLevel %errorlevel% & goto :cmDone :cmErrorLevel exit / b %1 :cmDone if %errorlevel% neq 0 goto :VCEnd :VCEnd" コード 1 で終了しました。 C:\Program Files (x86)\MSBuild\Microsoft.Cpp\v4.0\V110\Microsoft.CppCommon.targets 134

VTK の場合:

最初に VTK_USE_QT を有効にしたにもかかわらず、VTK_USE_QTOPENGL オプションが表示されず、VS2012 のプロジェクト生成が数回クラッシュしましたが、クラッシュを乗り越えることができました (これは、実際に VTK を使用するのに良い兆候ではないのではないかと心配しています)。

クハルの場合:

FLANN と同様のエラー:

エラー 1 エラー MSB3073: コマンド "setlocal C:\CMake\bin\cmake.exe -DBUILD_TYPE=Debug -P cmake_install.cmake if %errorlevel% neq 0 goto :cmEnd :cmEnd endlocal & call :cmErrorLevel %errorlevel% & goto : cmDone :cmErrorLevel exit /b %1 :cmDone if %errorlevel% neq 0 goto :VCEnd :VCEnd" コード 1 で終了しました。C:\Program Files (x86)\MSBuild\Microsoft.Cpp\v4.0\V110\Microsoft. CppCommon.targets 134

しかし、残りはうまく構築されます。どんなアドバイスも大歓迎です。

衛星画像からデータを補間しています。私の初期データは通常のグリッドではありません。データは、たとえば、および でch1指定された位置を参照します。最近傍補間を試みると、得られる結果は良好です。longitudelatitude

ただし、双一次補間を要求すると、Qhull エラーが発生します。

私が得るエラーは次のとおりです。

griddata ランタイム エラー -- Python/SciPy (Interpolation)の投稿を読みましたが、私の場合、一方の方法でデータを補間しますが、もう一方の方法では補間しません。

私は何を間違っていますか？

下の画像では、最終点の位置 (赤) をプロットし、初期位置 (青) を重ねています。

私のコードでは、一連の点の Convexhull を計算する必要があり、何らかの理由でqhullライブラリを使用する必要があります。このライブラリには、必要なことを正確に実行するメソッドqconvexがあります。ターミナルでこのコマンドを実行して、必要なものを取得できます。たとえば、次のような入力があるとしますpoints.txt。

これらのコマンドを端末で実行して、結果を取得できます:qconvex Fx < points.txtまたはcat points.txt | qconvex -Fx、出力は次のとおりです。

ここで私の質問は、C++ コードで入力に対して反復的にこのコマンドを呼び出す方法forですfloat **rs_tmp;。毎回これらの 10 ポイントの qconvex。コードで実行し、入力としてqconvexパイプする方法を教えてください。rs_tmpコードを超高速にする必要があるため、一時ファイルに書き込むことを避けrs_tmp、そこから読み取ることを好みます。

ライブラリ Scipy と ConvexHull を使用して凸包を作成しようとしています。私の知る限り、それは QHull を呼び出します。

追加したい点に「全次元」がない場合に問題が発生します。例：

出力用:

ただし、凸包が完全な次元になるように、ポイントを追加すると、次のようになります。

その後、すべてが機能します。1 つの例と他の例の違い (他の多くの例を行ったので、確信があります) は、最初のケースの凸包が 2 次元空間で 1 次元であるのに対し、2 番目のケースでは 2- 2 次元空間の次元 (つまり、全次元)。

何か案は？回答で言及されているように、ドキュメントが示しているので、いくつかのqhull_optionsを渡すと思いました：

QHullError 解決するオプションが有効になっていないときの幾何学的縮退など、Qhull がエラー状態に遭遇したときに発生します。

ただし、QHull の多くのオプションを読みましたが、この問題に対処しているようには見えません。私はそれらのいくつかを無作為に試しましたが、ほとんど成功しませんでした.

どんな助けでも役に立ちます。私は何百ものこれらの船体を作成するプログラムに取り組んでおり、そのうちのいくつかは完全な次元ではありません.

Visual Studio と C++ の基本は知っていますが、上級ユーザーでもプログラマーでもありません。私がやりたいことは、ドロネー三角形分割を実行し、qhull ライブラリを使用してボロノイ グリッドを生成することです (Windows で Visual Studio C++ を使用)。計算に必要なデータの準備は問題ありません。私の質問は次のとおりです。

1. qhull からインクルードするのに最低限必要なファイルは何ですか? (公開されているものから、qhull は非常に強力で、膨大な量のファイルが含まれています)

2. ドローネ三角形分割を簡単に実行し、ボロノイ グリッドを生成するにはどうすればよいですか? (Matlab の delaunayn.m ファイルのように、" t = qhullmx(x', 'd ', opt); " を呼び出してから、delaunay 三角形分割を取得します。)

あなたの経験を共有してくれてありがとう！