問題タブ [gdal]

gdal / ogrを使用してベクターマップ/シェープファイルをロードし、Qtを使用して設計されたGUIに表示するアプリケーションを作成しています。私はベクターファイルを扱うのが初めてで、GUIでそれらをレンダリングする方法を決めることができません。私の頭に浮かぶいくつかのアプローチは次のとおりです。-ベクターファイルをラスターに変換してラスターとして表示する。-ベクターファイルの各形状を1つずつ描画する（これを行う方法がわからない、それができるかどうかさえわからない） GDAL / OGRを使用してこれを行うことも可能です。私の場合、ベクターファイルは非常に複雑になる可能性があるため、これはどれほど複雑になるでしょう）

ベクターファイルをレンダリングする方法の決定は私にとって非常に重要です。ベクターファイルを表示した後、投影の変更、インタラクティブな選択、地理参照などのいくつかの高度な機能をサポートする必要があるためです。

ベクトル/地理情報システムに取り組んだことがある人なら誰でも、これが通常どのように行われるか、そしてそのようなデータを表示するためにどのQtウィジェットが使用されるかを教えてください。

Windows 用の OSGeo4W の最新のコンパイルをインストールしました。利用可能なさまざまな GDAL ツールを使用したいと考えています。OSGeo4Wコマンド シェルを開くと、 Python と入力して次のように記述します。

次のエラーが表示されます。

私のPATHには、次の OSGeo4W エントリが含まれています。

Web で読んだことによると、このエラーは互換性のないバージョンに関するものである可能性があります。OSGeo4W パッケージはそのまま動作するようにコンパイルされていると想定しているため、以前にコンピューターにインストールされた他のバージョンまたは干渉する残留物である可能性があります。しかし、干渉しているものとその情報がどこに保存されているかを判断する方法がわかりません。

この問題を解決するにはどうすればよいですか?

GDAL を使用してグリッド (ASCII グリッド) でラスター値を操作しようとしています。しかし、これに進む前に、新しい値をファイルに書き込むのに問題があります。これらのエラー メッセージは、slopeband.WriteArray(s) が呼び出されたときに表示されます。

エラー 6: 勾配.asc、バンド 1: このデータセットでは WriteBlock() はサポートされていません。

エラー 1: スロープ.asc、バンド 1: ダーティ ブロックの書き込み中にエラーが発生しました

これが非常に基本的なものである場合は申し訳ありませんが、特に python と GDAL にはまだ慣れていません。Mac OS X 10.6.8 および Python 2.7 で GDAL 1.9.0 を使用しています。ありがとうございました！

TIF ファイルから PostGIS データベースにサブデータセットをロードしたいと考えています。

gdalinfo は、TIF ファイルに関する次の情報を提供してくれます。

2 番目のサブデータセット (392x514352) を読み込みたいと思います。

現在、raster2pgsql とファイル名を入力として使用すると、最初のサブデータセットのみが読み込まれ、2 番目のサブデータセットは無視されます。バンドは 1 つしかないため、-b スイッチを使用しても機能しません。何かアドバイス？

どうもありがとう

トーマス

地理空間データを取得していくつかのステップを実行するアルゴリズムを作成しました。入力データは、大規模なラスター調査エリア (~1 億 5000 万ピクセル) のポリゴンと共変量ラスターのシェープファイルです。手順は次のとおりです。

1. シェープファイルのポリゴン内からのサンプル ポイント
2. サンプリング ポイントごとに、共変量ラスターから値を抽出します。
3. サンプリング ポイントで予測モデルを構築する
4. ターゲット グリッド ポイントの共変量を抽出する
5. 予測モデルをターゲット グリッドに適用する
6. 一連の出力グリッドに予測を書き込む

プロセス全体を数回 (たとえば 100 回) 反復する必要がありますが、現在、連続して処理すると、各反復に 1 時間以上かかります。各反復で最も時間がかかるのはステップ 4 と 5 です。ターゲット グリッドが非常に大きいため、一度に 1 ブロック (たとえば 1000 行) ずつ処理しています。

私は 32 Gb の RAM を備えた 6 コアの CPU を使用しているため、各反復内でオブジェクトを含む Python のmultiprocessingモジュールを使用しPoolて、多数のブロックを同時に処理し (ステップ 4 と 5)、出力 (予測) を書き込んでみました。グローバル出力書き込み関数を呼び出すコールバック関数を使用して、出力グリッドの共通セットに。これは機能しているように見えますが、各ブロックを連続して処理するよりも高速ではありません (実際には、おそらく低速です)。

だから私の質問は、それを行うためのより効率的な方法はありますか? multiprocessing モジュールのQueueクラスに興味がありますが、どのように機能するのかよくわかりません。たとえば、ステップ 4 と 5 を実行するキューを作成してから、ステップ 6 を実行する別のキューに結果を渡す方が効率的かどうか疑問に思っています。

任意のポインタをいただければ幸いです。

TIFF形式の地理データのない画像（マップ）があります。
画像からGeoTIFFファイルを取得する必要があります。地図の各隅に緯度と経度があります。地理データをGoogle空間参照の画像に追加して、geotiffを取得するにはどうすればよいですか？私はGDALがそれを手伝ってくれることを知っています。誰かが私がコマンドを作成するのを手伝ってくれますか

Rails 3 プロジェクトで GDAL を使用しています。システムルビーをデフォルトに設定せずに正しく動作させる方法を理解できないようです。Rails は、「gdal/org」がインストールされているにもかかわらず見つからないと文句を言います。RVM でこのリンクを作成し、Ruby gemset を引き続き使用するにはどうすればよいですか?

過去 2 日間、python スクリプトを取得して、属性を持つ 289995 ポイントのシェープファイルの作成を完了しようとしました。ポイントは作成できますが、スクリプトは属性を完成しません。コードは、ループまで正しく実行されます: for j, p in enumerate(wCoords):(2 番目のループの下のコードを参照) しばらくすると、セグメンテーション違反が発生します。特定の間隔でプロセスを停止するifステートメントを追加して、ループサイクルが発生したときにその場所を見つけることができるかどうかを確認しようとしました。ループは 1000 サイクルまで問題なく循環しますが、10,000 サイクルまではフィードバックなしで停止し、エンドレス ループのように見えます。このプログラムは、木のポイントを作成し、そのポイントに木の高さ属性をアタッチします。コードは以下のとおりです。

gdal/ogr については、これ以上の情報を提供できるほど詳しくありません。助けてください。ジャック

RgoogleMaps パッケージの GetMap() 関数を使用して「静的マップ」を派生させました。それ（MyMap）をPNGとしてハードドライブに保存できます。ただし、空間参照が失われます。

そのような PNG から (GDAL で読み取り可能なデータ形式という意味で) 空間オブジェクトを作成することに成功した人はいますか?

pip install gdal仮想環境(Ubuntu)内にpipからgdalをインストールしようとしています。見つからないので失敗cpl_port.h

ただし、GDAL は正しくインストールされており、ヘッダー ファイルは にあります/usr/include/gdal/cpl_port.h。pipがヘッダーファイルを見つけるために設定する必要があるGDALの環境変数はありますか?