タグ付けされた質問 「gdal」

Ubuntu 12.04でecw for QGIS 1.8を機能させる方法については、さまざまなオンラインヘルプガイドに従っています。私はそれをすべてUbuntu 10.04で動作させていましたが、12.04への移行以降、そのような運はありませんでした。 Ubuntu 12.04で機能するガイドはありますか？ 現時点では、ecwサポートがインストールに含まれているWindowsバージョンをインストールする方が簡単なようですが、Ubuntuを使い続けたいと思います。

10 qgis  gdal  ubuntu  references  ecw 

タイルベースマップを作成しようとしています。タイルを作成したい画像があり、GDALを使用しています。 タイルを作成できましたが、結果を表示すると、タイルマップが次のように見えます。 ひっくり返されます 比率が正しくありません これは地図のオーバーレイ画像です： これはタイリングの結果です： オルソフォトの画像のコーナーを比較して、座標が良いことを何度も確認しました。 もう1つは、同じ座標のGoogle Earthでは、マップ画像がGoogleマップに正確にフィットすることです。 編集済み gdalinfo map_overlay.png この出力では： ドライバー：PNG / Portable Network Graphicsファイル：map_overlay.pngサイズは3527、2494座標系は `'画像構造メタデータ： INTERLEAVE = PIXELコーナー座標：左上（0.0、0.0）左下（0.0、2494.0）右上（3527.0 、0.0）右下（3527.0、2494.0）中央（1763.5、1247.0）バンド1ブロック= 3527x1タイプ=バイト、ColorInterp =赤マスクフラグ：PER_DATASET ALPHAバンド2ブロック= 3527x1タイプ=バイト、ColorInterp =緑マスクフラグ：PER_DATASET ALPHAバンド3ブロック= 3527x1タイプ=バイト、ColorInterp =ブルー マスクフラグ：PER_DATASET ALPHAバンド4ブロック= 3527x1タイプ=バイト、ColorInterp =アルファ 次に翻訳します： gdal_translate -of VRT -a_srs EPSG：4326 -gcp 0 0 31.7431761644 35.1680410195 -gcp 3527 0 31.7493769674 …

10 gdal  tiles 

私はpython地理空間プログラミングを初めて使用します。次のスクリプトを実行して、対応するエラーメッセージを取得しました >>> import osgeo >>> import osgeo.ogr >>> shapefile = osgeo.ogr.Open("tl_2009_us_state.shp") >>> numLayers = shapefile.GetLayerCount() Traceback (most recent call last): File "<pyshell#5>", line 1, in <module> numLayers = shapefile.GetLayerCount() AttributeError: 'NoneType' object has no attribute 'GetLayerCount'

10 python  gdal 

緯度/経度座標で回転した長方形を作成しようとしています。 基本的に、私は中心の座標と中心の周りの回転角度を持っています。 形見やgdalなどを使って行うことはできますか？

10 python  gdal  geometry  shapely  affine-transformation 

QGISプリントマネージャーツールで小さな解像度のラスターファイルを表示しようとしています。各セルの値（各セルの中央付近）に小さなラベルを表示したいと思います。 どうすればできますか？このタスクを完了するのに役立つプラグインはありますか？ 私の考えは、gdal2xyzを使用し、CSVプラグインを使用してそのようなファイルをロードし、それをベクターファイルとして保存することです。残念ながら、私は多くのラスターファイルを持っているため、多くの作業が必要です。 もう1つのアイデアは、gdal2xyz.pyを使用してから、ogr2ogrを使用してベクターファイルとして保存することです。 このプロセスをバッチ処理する方法はありますか？もっと簡単にできるオープンソースプログラムはありますか？

10 qgis  raster  gdal  ogr 

多数のシェープファイルをPostGISデータベースに正常にロードしました。SQLクエリに基づいて、特定のジオメトリとレイヤーをdxfファイルとして取得できるようにしたいと思います。 WindowsのOSGeo4WスイートのPostGIS、GDAL、PostgreSQL、Pythonを使用しています。 私はogr2ogrを次のように少し使用する必要があると想像しています： $ ogr2ogr -sql "SELECT * FROM some_table" -f "dxf" outputfile.dxf -overwrite "PG:dbname=mydb user=dbusername password=dbpassword" しかし、適切な値でそれを行うと、次のような一連のエラーが発生します。 ERROR 1: DXF layer does not support arbitrary field creation, field 'objectid' not created. レイヤー内のすべての属性列でエラーが発生しsome_tableますが、dxfファイルを取得します。 .dxfファイルを既に作成していて、上記のコマンド（-overwriteオプションを含む）を再実行すると、代わりに次のエラーが発生します。 Layer sql_statement not found, and CreateLayer not supported by driver. ERROR 1: Terminating translation prematurely after …

10 postgis  gdal  ogr  dxf 

Pythonがインストールされているシステムにインストールされているgdal形式のリストを取得する方法はありますか？コマンドラインでは、次のようにして簡単に行うことができます。 gdalinfo --formats ...しかし、それはpythonで実行できますか？

10 python  raster  gdal 

私たちが持っている高解像度画像のタイルを生成するのに少し苦労しています。現在の画像は非常に大きい（+ 20GB）画像で、GeoTiffファイルとして保存されています。 TMS画像プロバイダーを使用してタイルを提供し、gdal2tilesコマンドラインユーティリティを使用してタイルを生成し、Cesiumで開いて表示します。gdalinfoを使用して、画像の詳細の一部を次に示します。 Driver: GTiff/GeoTIFF Files: image.tif Size is 52250, 56119 Coordinate System is: PROJCS["WGS 84 / UTM zone 35S", GEOGCS["WGS 84", DATUM["WGS_1984", SPHEROID["WGS 84",6378137,298.257223563, AUTHORITY["EPSG","7030"]], AUTHORITY["EPSG","6326"]], PRIMEM["Greenwich",0], UNIT["degree",0.0174532925199433], AUTHORITY["EPSG","4326"]], PROJECTION["Transverse_Mercator"], PARAMETER["latitude_of_origin",0], PARAMETER["central_meridian",27], PARAMETER["scale_factor",0.9996], PARAMETER["false_easting",500000], PARAMETER["false_northing",10000000], UNIT["meters",1], AUTHORITY["EPSG","32735"]] Origin = (606276.000000000000000,7197873.000000000000000) Pixel Size = (0.500000000000000,-0.500000000000000) Metadata: AREA_OR_POINT=Area TIFFTAG_MAXSAMPLEVALUE=13165 TIFFTAG_MINSAMPLEVALUE=1 TIFFTAG_RESOLUTIONUNIT=2 …

10 gdal  gdalwarp  cesium  gdal2tiles 

私はGDALとPythonを使用してシェープファイルを操作しています。ソースファイルを開くとき、必要なシェープファイルを使用してドライバーの「Open」メソッドを呼び出す前に、まずドライバーを「ESRIシェープファイル」に設定しました。 私はこれで問題はありませんでした。ただし、Driverインスタンスのメソッド「Open」とogr関数「Open」の違いは何なのか気になります。たとえば、次の違いは何ですか。 driver = ogr.GetDriverByName("ESRI Shapefile") source = driver.Open(“O:\myshp.shp”,0) そして source = ogr.Open(“O:\myshp.shp”,0) どちらを使用するかという特定のシナリオはありますか？

10 python  gdal 

以下に示すように、範囲/境界ボックスではなく、いくつかの単一レイヤーラスターイメージのポリゴンアウトラインを作成する必要がありますが、nodata値がない領域を作成します。。 上記の質問への回答で、Image Boundaryプラグインについて言及しましたが、Ubuntuで実行されているQGIS 1.8.0 Lisboaでは見つかりません。 ツールはまだ利用できますか？ そうでなければ、ArcMapではなく、gdal、R、QGIS、GRASSなどのオープンソースツールを使用してこれを行う方法はありますか？

10 gdal  pyqgis  qgis-plugins  footprint 

.tiffファイルから地理参照ファイルを作成したいの.pngですが。そこで.tiff、次のコマンドラインを使用して.pngを通常のファイルに変換します。 convert image.png image.tiff ここでは、すべて正常に動作.tiffします。有効に表示される有効なファイルがあります。次に、ジオリファレンスデータを追加します。そうするために、私はというツールを見つけましたtiff2geotiff、私はそれを次のように使用しています： tiff2geotiff -4 "+proj=latlong +ellps=WGS84 +datum=WGS84 +no_defs" -c none -n "-122.5575664 38.5818201 -122.5449425 38.5896175" notgeoreferenced.tiff georeferenced.tiff これにより、QGisにインポートされた黒の画像が適切な場所で正常に表示されますが、どこも黒であることを除きます（地理参照データが適切に追加されていることの証明）。 したがって、私がするとき： gdalinfo georeferenced.tiff 私は得る： Warning 1: TIFFReadDirectory:Bogus "StripByteCounts" field, ignoring and calculating from imagelength Driver: GTiff/GeoTIFF Files: araujo_tmp.tiff Size is 1585, 979 Coordinate System is: GEOGCS["WGS 84", DATUM["WGS_1984", SPHEROID["WGS …

10 qgis  gdal  php  geotiff-tiff 

30mの解像度でグローバルラスタを操作および処理しています。合計ラスターサイズは通常[1,440,000 560,000]です。スーパーコンピューターにアクセスできるので、グローバルラスタを扱いやすいチャンクに分割し、いくつかの計算を並行して実行し、それらをディスクに非常に迅速に書き込むことができるコードを作成しました。 しかし、結果を表示することに関しては、壁にぶつかっています。私は通常、地球をカバーするタイルの仮想ラスターを作成し、それをQGISに取り込みます。しかし、それは信じられないほど遅いです（もしそうならロードするのに数分）。そして、パンまたはズームしようとすると、さらに数分かかります。この問題を解決するための最初のアプローチは、gdaladdoを使用して概要を作成することでした。ただし、これらのビルドには（日単位のように）永久にかかるため、アルゴリズムの開発には役立ちません。これが私が試したことのリストと、なぜ/どのように失敗したかです。 vrtの概要を構築します。上記のように、これは8つのレベルで完了するのに2日以上かかります。それは私の目的には受け入れられません。 個々のタイルの概要を作成し、概要を含むvrtに何らかの方法でマージします。タイルのオーバービューを比較的迅速に作成できます（スーパーコンピューター）が、それらを再マージすることはできません。私は試した： 2a。概要付きのタイルのgdal_merge。ただし、概要は出力tiffに保持されません（または少なくともQGISで認識されません）。 2b。概要付きのタイルでgdalbuildvrtを使用しましたが、上記のように、概要は保持されませんでした。[これは正しくありません。編集を参照してください。] 2c。レベル1〜6のタイルの構築概要とvrt（直接オプション2b）でレベル7〜8を直接構築するハイブリッドも試してみましたが、これら2つのレベルだけではまだ時間がかかります。いくつかのテストを行ったところ、タイルの概要が実際にvrtの概要を作成するために使用されていることがわかりましたが、vrtの概要を完了するにはまだ1日ほどかかります。 だから私はここに誰かが私が次にどこに行くべきかについていくつかの提案があることを願っています。ここに私が検討しているいくつかのオプションがあります： 手動でグローバルピラミッドを手動で作成します。それらを.ovrファイルに再結合するのは注意が必要です。 マップサーバー（ジオサーバー）を使用します。私はこれについてほとんど知りませんが、プロセスに複雑さを加えている間に時間のハードルを克服できないのではないかと心配しています。 大陸またはその他の地域によってドメインを分割します。私は本当にこのオプションを避けたいです。 「なぜ地球全体を30mの解像度で表示する必要があるのですか？」1つの例：私は水のピクセルのマスクを（全体的に）取り、それをスケルトン化して川を見つけ、測定を実行します。私のスケルトン化アルゴリズムは少しの調整が必要で（枝刈り、ループの削除、一般的なクリーニングなど）、出力は必ず30mです。川や景観は世界中で多様であるため、実装した変更の影響を確認するには、周囲をパンできる必要があります。 また、QGISを調べて、巨大なラスターをより速くレンダリングするために使用できる設定がないことを確認しましたが、何も表示されませんでした。SSDドライブを購入するのに足りないので、私はそれができるだけ速く動いていると思います。（私のHDDのI / Oは〜250MB / sです）。 個々のタイルで概要を作成し、次にvrtを作成しても概要が維持されていることがわかりました。ファイルのメタデータのQGISの「Pyramid」セクションは空ですが、「Dimensions」セクションには概要の各レベルのエントリがあります（X 720000、Y 140、X 360000、Y 70など）。だから私は2bについて間違っていました。 また、すべてのタイルをQGISにプルしただけでは1分未満でレンダリングされますが、タイルを参照するvrtをプルした場合は、5分以上かかります（処理する）。 SSDを搭載したコンピューターでいくつかのテストを行ったところ、グローバルなvrts（概要なし）を適切な速度で正常に読み込み、表示、レンダリングできることがわかりました。私のコンピューターでも同じことができることを期待して、1TB PCIe SSDを注文しました。結果で更新されます。

10 qgis  gdal  performance  display  gdaladdo 

GDALで多くの問題を抱えています。ドキュメントが欠けていることもありますが、Pythonではほとんどサポートされていないようです。とにかく、これは次のように機能するバッチ再プロジェクトスクリプトです。 source_file = gdal.Open(filepath.encode("ascii")) source_wkt = source_file.GetProjectionRef() source_srs = osr.SpatialReference() source_srs.ImportFromWkt(source_wkt) reproj_file = gdal.AutoCreateWarpedVRT(source_file, source_wkt, dest_wkt) gdal.ReprojectImage(source_file, reproj_file, source_wkt, dest_wkt) reproj_attributes = reproj_file.GetGeoTransform() driver = gdal.GetDriverByName("GTiff") dest_file = driver.CreateCopy(outputpath.encode("ascii"), reproj_file, 0) dest_srsとdes_wktの部分はこのコードチャンクでは定義されていませんが、ループの外にあります（一度だけ定義する必要があるため）。それは一度動作するようです、私はそれから1つの見栄えの良いtifを得ることができ、それから私に「エラー6 WriteBlock（）はサポートされていません」とpythonがクラッシュします。それらはすべて、同じ基本データを使用して同じ方法で作成された（時間だけが異なる）すべてのGeoTIFFです。 また、（GCSからPCSへの）再投影の性質により、AutoCreateWarpedVRTは多くの空白スペースを作成する傾向がありますが、値0を与えます。これは、実際のデータ値である可能性があるため問題です。代わりにnodata値を-99に設定する方法はありますか？

10 python  gdal  coordinate-system 

緯度、経度、data_valが次の形式のASCIIファイルがあります。 35-13.643782N, 080-57.190157W, 118.6 ... GeoTiff画像ファイルを持っているので、簡単に表示できます。 ASCIIファイルにある特定の緯度/経度の位置にある画像に「ピン」（ドット/フラグ/スターまたは最も簡単なものを使用できます）を配置したいと思います。 これが私がこれまで何とかしてきたことです： 私のソース画像は次のようになります： Driver: GTiff/GeoTIFF Files: /tmp/Charlotte SEC 100.tif Size is 16867, 12358 Coordinate System is: PROJCS["Lambert Conformal Conic", GEOGCS["NAD83", DATUM["North_American_Datum_1983", SPHEROID["GRS 1980",6378137,298.2572221010042, AUTHORITY["EPSG","7019"]], AUTHORITY["EPSG","6269"]], PRIMEM["Greenwich",0], UNIT["degree",0.0174532925199433], AUTHORITY["EPSG","4269"]], PROJECTION["Lambert_Conformal_Conic_2SP"], PARAMETER["standard_parallel_1",38.66666666666666], PARAMETER["standard_parallel_2",33.33333333333334], PARAMETER["latitude_of_origin",34.11666666666667], PARAMETER["central_meridian",-78.75], PARAMETER["false_easting",0], PARAMETER["false_northing",0], UNIT["metre",1, AUTHORITY["EPSG","9001"]]] Origin = (-365041.822331817995291,240536.419747152860509) Pixel Size = (42.334586069440391,-42.334898968590878) …

9 python  gdal  latitude-longitude  geotiff-tiff  ascii 

PyPIにGDALを使用するPythonスクリプトをパッケージ化しようとしています。私は私の中に直接参照を含めることから始めましたsetup.py： install_requires=['GDAL==1.11.2'], このようにして、パッケージはテスト仮想環境にインストールできませんでした： extensions/gdal_wrap.cpp:2855:22: fatal error: cpl_port.h: No such file or directory #include "cpl_port.h" ^ compilation terminated. error: Setup script exited with error: command 'x86_64-linux-gnu-gcc' failed with exit status 1 次にpygdal、virtualenvフレンドリーバージョンとしてマークされているため、を参照して試してみました。 install_requires=['pygdal'], このようにして、インストールはエラーなしで終了します（ただし、通常のコンパイル警告のロードが伴います）。しかし、その後スクリプトを呼び出すと、次のエラーが返されます。 Traceback (most recent call last): File "/home/desouslu/.virtualenvs/test_p3/bin/hasc2gml", line 5, in <module> from pkg_resources import load_entry_point File "/home/desouslu/.virtualenvs/test_p3/lib/python3.4/site-packages/pkg_resources.py", …

9 python  gdal