問題タブ [gdal]

簡単なスクリプトを実行してジオメトリを手動で作成し、そのジオメトリに対してバッファを実行しようとしています。ただし、バッファ操作を実行しようとすると、python が完全にクラッシュします。

Windows 7 64ビットを実行しています。

GDAL と Python を使用して画像分類の最小距離アルゴリズムを実装しようとしています。サンプル領域の平均ピクセル値を計算し、それらを配列のリスト ("sample_array") に格納した後、画像を "values" という配列に読み込みました。次のコードでは、この配列をループします。

大きな画像の場合、この手順には非常に時間がかかります。誰かがより速い方法を知っているかどうか尋ねたいのはそのためです。Pythonのさまざまな変数のアクセス速度についてはあまり知りません。または、誰かが私が使用できるライブラリを知っているかもしれません。前もってありがとう、マリオ

Windows 7 ワークステーションに GDAL と Python をインストールしようとしています。

• Python バージョン: 2.6.6
• グダル 1.8
• (Visual C++ Studio 2010 Express)

http://ucsb.pbworks.com/w/page/1936549/Using-GDAL-with-Pythonの指示に従い、システム環境変数に PYTHONPATH と GDAL_DATA を追加しました。さらに、GDAL バイナリ フォルダーを追加して PATH 変数を変更しました。

ダウンロードした GDAL Python バインディング: http://vbkto.dyndns.org/sdk/PackageList.aspx?file=release-1600-gdal-1-8-mapserver-5-6.zip

" from osgeo import gdal " により、次のエラー メッセージが表示されます。

前もってありがとう、Martwig

Delphi と GDAL18 を使用してアイソラインを作成しようとしています。そのために、次のコードを使用しています。

GDALContourGenerate コマンドは、「サポートされていないジオメトリ タイプです」というエラーを返します。

次の方法で gdal18.dll を含めました。

他のジオメトリ タイプ (wkbLineString25D など) も試しましたが、これは役に立ちませんでした。何か提案があれば幸いです。どうもありがとう、マリオ

[編集] 「レイヤー」(GDALContourGenerate 内) を「nil」で再生すると、同じエラーが発生することがわかりました。したがって、問題は別の場所にある可能性があります。[/編集]

ディレクトリ内のいくつかのテキストファイルを繰り返し処理し、データ操作を行い、RのrgdalパッケージのwriteOGRを使用してESRIシェープファイルに書き込むスクリプトでエラーが発生しました。特定の数のファイルの後でエラーが発生し続けます書き込まれます（各シェープファイルが多数のファイルを作成するため、約43〜44個の個別のファイル、または約14個の完全なシェープファイルのように見えます）。

エラーは通常次のようなものです

さらに、スクリプトがクラッシュした後、次のようなエラーが発生します

別のファイルを読んだり、何かをしようとしたとき。

この時点では、開いているファイルは表示されません。

問題を再現するコードを次に示します。

OSXでR2.12.2を使用しています。エラーは64ビットバージョンと32ビットバージョンの両方で発生します。

GDAL ツールで ogrinfo [-spat xmin ymin xmax ymax] オプションを使用した人はいますか? 形状ファイルで - sql クエリを実行して回答/形状を取得することはできますが、スパットで緯度/経度の値を使用すると、結果が得られません (ただし、エラーは発生しません)。例がまったく見つかりませんでした。

python、gdal、numpyを使用して標高/高さフィールドラスターを作成したいと思います。私はnumpy（そしておそらくpythonとgdal）で立ち往生しています。

numpyでは、私は次のことを試みてきました：

osgeoimportgdalからgdalconstimportから*

シンプルなものが欠けていると思いますので、アドバイスをお待ちしております。

ありがとう、

クリス

（後で続く）

• terragendataset.cpp、v 1.2 *
  • プロジェクト：Terragen（tm）TERドライバー
  • 目的：TerragenTERドキュメントのリーダー
  • 著者：Ray Gardener、Daylon Graphics Ltd. *
  • GDALドライバーから派生したこのモジュールの一部
  • フランク・ウォーマーダム、http： //www.gdal.orgを参照

RayGardenerとFrankWarmerdamに事前にお詫び申し上げます。

Terragenフォーマットノート：

書き込みの場合：SCAL =グリッドポストの距離（メートル単位）hv_px = hv_m / SCAL span_px = span_m/SCALオフセット=TerragenDataset:: write_header（）を参照scale = TerragenDataset :: write_header（）を参照物理hv =（hv_px-offset）* 65536.0 / scale

発信者に次のように伝えます。

これは、上記のWriteArray（somearray）の前に、GeoTransformとSetProjectionおよびSetUnitTypeの両方を（潜在的に）機能させるために設定する必要があることを示しています。

GDAL APIチュートリアルから：Python import osr import numpy

過度に長い投稿と告白を作成してしまったことをお詫びします。これが正しければ、すべてのメモを1か所（長い投稿）にまとめておくと便利です。

告白：

以前に写真（jpeg）を撮り、それをgeotiffに変換し、タイルとしてPostGISデータベースにインポートしました。現在、画像をドレープするための標高ラスターを作成しようとしています。これはおそらくばかげているように思えますが、私が尊敬したいアーティストがいますが、同時にこれらの優れたツールの作成と改善に熱心に取り組んできた人々を怒らせません。

アーティストはベルギー人なので、メーターが適切です。彼女はニューヨークのロウアーマンハッタン、UTM18で働いています。上記の写真はw=3649 / h = 2736です、これについて少し考えなければなりません。

別の試み：

明らかに近づいていますが、SetUTM（18,1）がピックアップされたかどうかは不明です。これはハドソン川の4x4ですか、それともLocal_CS（座標系）ですか？サイレントフェイルとは何ですか？

もっと読む

4x4メートルはかなり小さな論理スパンです。

だから、おそらくこれはそれが得るのと同じくらい良いです。SetGeoTransformは単位を正しく取得し、スケールを設定して、Terragenワールドスペースを作成します。

最終的な考え：私はプログラムしませんが、ある程度は従うことができます。そうは言っても、私は最初に私の小さなTerragen World Spaceでデータがどのように見えるか疑問に思いました（http://www.gis.usu.edu/~chrisg/python/2009/第4週に感謝します）：

ですから、これは満足のいくものです。上記で使用したnumpycとこの結果の違いは、この非常に小さな論理スパン全体にFloat16を適用するアクションにあると思います。

そして「hf2」へ：

私はラコンコルディアペルーにいるように見えますが、ほぼ完全に満足しています。だから私は、ニューヨークの北のように、もっと北のように言う方法を理解する必要があります。SetUTMは、「北または南の距離」を示唆する3番目の要素を取りますか。昨日、赤道でC、ニューヨークエリアでTまたはSのような文字ラベルゾーンがあるUTMチャートに出くわしたようです。

私は実際、SetGeoTransformが本質的に左上の北と東を確立していると思っていました。これは、SetUTMの「南北の距離」の部分に影響を与えていました。gdal-devに移動します。

そして後でまだ：

パディントンベアはチケットを持っていたのでペルーに行きました。そこに行きたいと言ったのでそこに着きました。Terragenは、そのように機能して、ピクセルスペースを与えてくれました。その後のsrsの呼び出しは、hf2（SetUTM）に作用しましたが、東と北はTerragenの下で確立されたため、UTM 18が設定されましたが、赤道のバウンディングボックスにありました。十分です。

gdal_translateは私をニューヨークに連れて行ってくれました。私はウィンドウズにいるので、コマンドラインです。そしてその結果。

だから、私たちはニューヨークに戻ってきました。これらすべてにアプローチするためのより良い方法がおそらくあります。numpyからのデータセットも仮定/即興していたので、Createを受け入れるターゲットが必要でした。作成を許可する他の形式を調べる必要があります。GeoTiffでの標高は可能性です（私は思います）。

すべてのコメント、提案、および適切な読書に向けた穏やかな微調整に感謝します。Pythonでマップを作成するのは楽しいです！

クリス

USDA Crop Data Layer (データはこちらから入手できます。大きなファイルに注意してください: http://www.nass.usda.gov/research/Cropland/SARS1a.htm ) をダウンロードし、分析のために R にインポートしたいと考えています。他の大きな GTIFF ファイル (LANDSAT データなど) をインポートしましたが、関数を呼び出すと次のエラーが発生します。

ここで、「ファイル」は USDA が提供する .tif ファイルのいずれかです (中西部地域のコンパイルを使用しています)。

問題は、GTiff が真の序数/スペクトル データではなく、カテゴリ データ (土地利用カテゴリ) であることだと思いますが、画像エディター (プレビューなど) で正常にレンダリングされ、R / GDAL でボークする理由がわかりません。 .

助けてくれてありがとう。

編集: これが GIS フォーラムにある場合は、私に知らせてください。一般的に、空間問題に対するプログラミングのアプローチはここにあるようです。

注: もう少し読みやすくするために、各コマンドから他の 5 つのファイル名を削除しました。私はバッチファイルを作成するPythonツールに取り組んでいますが、この.batを実行しようとすると、最初のコマンドが正常に実行されてからコンソールに戻るだけです。FWTools シェルで実行しています。16 進エディタでバッチ ファイルを開き、余分な改行文字と改行文字を検索しましたが、通常とは異なるものは何も表示されないため、Python は正しく記述しているようです。

簡単に：

.prj WKTファイルをパラメーター化して、7ポイントの変換を実行できるようにする方法（wiki）。パラメータがどのように機能するかは知ってfalse_easting and false_northingいますが、スケールを調整するにはどうすればよいですか？scale_factor'を意味するのではありません

それが問題の説明です：

非GIS環境（輸送モデリングソフトウェア）に輸送ネットワーク（ベクトル層）を保存しています。ネットワークは、ノード（ポイント）とポリライン（道路リンク）で構成されます。投影や座標などに関係なく、ほとんどの場合、ランダムな背景から実行されます。

ネットワークに適切な予測を設定する必要があります。

.prjファイルにアクセスできます（たとえば、WGS84プロジェクションを使用している場合は、他のプロジェクションに切り替えることができます）

それが私が試していることです
。7ポイントのヘルメルト変換（http://proj.maptools.org/gen_parms.html）を試します。.prjファイルのWKTパラメーターとしてtowgs84変換を使用します。ここでは、回転行列がゼロであると想定し（そうすることはできますか？）、delta_x、delta_y、およびscaleparamのみを計算します。

ただし、動作しません。これは私の.prjです。TOWGS84のパラメータは変換に影響しません。

だから私はパラメータを使用しようとしましたfalse_norting、false_eastingそしてそれらはうまく機能し、そして私のネットワークを適切に変換します、しかし：

それは私のネットワークの規模を変えることはなく、位置だけを変えます。では、.prjファイルを使用してネットワークを再スケーリングするにはどうすればよいですか？

ヒントをありがとう