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

私はgdal_proximityを使用して、米国（48州下）の最も近い主要な川までの距離を見つけています。NHD +ネットワークのフローラインをConus Albers（epsg：5070）に投影し、ストリームの順序が5を超える選択された川、ラスタ化された川を255、燃えない川を0としてラスタライズしました。これで十分です。 50km以内の場所で最も近い川へ。入力ファイルは大陸規模で30mの解像度であるため非常に大きくなりますが、変換は単純なgdal_proximityコマンドである必要があります。 gdal_proximity.bat -values 255 -distunits GEO -maxdist 50000 -nodata -999 infile.tif outfile.tif -co COMPRESS=DEFLATE -co BIGTIFF=YES -co TILED=YES これはほとんど機能しているように見えますが、出力に奇妙な幾何学模様が生成されています（画像を参照）。出力に存在するデータは正しく処理されています。なぜそんなに多くの出力が欠落しているのか、誰でも提案できますか？ 編集：これがオプションのパラメーターのいずれかによって引き起こされているかどうかをテストするために、この構成でgdal_proximityを再度実行しました。 gdal_proximity.bat H:\data\tmp\NHDplus_network_flowline_SO6plus.tif H:/data/tmp/NHDplus_network_flowline_SO6plus_proximity.tif -values 255 -maxdist 50000 -of GTiff 本質的に同じ結果が得られました： 私の唯一の考えは、それがラスターのサイズに関係しているかもしれないということです（最大100 GB非圧縮。私が知る限り、BigTiffのサイズに制限はありませんが、gdalができることには制限があるかもしれません効果的に分析しますか？

28 qgis  gdal  qgis-processing 

シェイプファイルジオメトリに従ってラスターファイルをクリップしようとしています。私は次のコードを使用しています gdalwarp -cutline INPUT.shp INPUT.tif OUTPUT.tif ただし、シェープファイルジオメトリ範囲外の黒色になります。ここで例を挙げています。クリップする最初の画像。2番目の画像は結果のラスタ​​ーですが、ジオメトリの外側がnullになる3番目のラスターが必要です。 1枚目の画像：INPUT Raster 2番目の画像：OUTPUTラスター 3枚目の画像：このような出力を取得したい

26 gdal  gdalwarp 

パターン認識のためにDEMラスターでフィルターをチェックしようとしており、常に最後の行と列が欠落しています（like..20）。PILライブラリ、イメージのロードで試しました。次に、numpyを使用します。出力は同じです。 配列の値をチェックするとき（ArcCatalogでIDを使用してピクセルを選択する）、ループに何か問題があると思ったので、ピクセル値が配列に読み込まれていないことに気付きました。 したがって、単に開いて配列に入れ、配列から画像を保存するだけです： a=numpy.array(Image.open(inraster)) #raster is .tif Float32, size 561x253 newIm=Image.new(Im.mode, Im.size) Image.fromarray(a).save(outraster) 最後の行と列を切り取ります。申し訳ありませんが、画像を投稿できません 誰もがその理由を理解するのを助けることができますか？そして、いくつかの解決策をアドバイスしますか？ 編集： だから、私はみんなの助けを借りて小さなラスターをnumpy配列にロードすることに成功しましたが、大きなイメージがあるとエラーが発生し始めます。私はそれがnumpy配列の制限についてだと思うので、配列は自動的に再形成されるか、そのようにsmthされます...だから例： Traceback (most recent call last): File "<pyshell#36>", line 1, in <module> ima=numpy.array(inDs.GetRasterBand(1).ReadAsArray()) File "C:\Python25\lib\site-packages\osgeo\gdal.py", line 835, in ReadAsArray buf_xsize, buf_ysize, buf_obj ) File "C:\Python25\lib\site-packages\osgeo\gdal_array.py", line 140, in BandReadAsArray ar = numpy.reshape(ar, [buf_ysize,buf_xsize]) File …

26 python  raster  gdal  numpy  floating-point 

Osgeoインストーラーを使用してGDALをインストールしました。ラスターレイヤーをベクターレイヤーでプログラムでクリップするにはどうすればよいですか？これに役立つGDAL APIはありますか？Pythonを使用しています。

26 python  gdal 

ESRIのファイルジオデータベースのAPIは現在、GDAL / OGRで使用することができます- http://www.gdal.org/ogr/drv_filegdb.html各個人がESRIにはログインする必要があり、個別にAPIファイルをダウンロードしているようだが。たとえば、これはここのブログ投稿およびGDALビルドノートで参照されています。 Esriが示すAPIをダウンロードする前のライセンス契約は、http：//resources.arcgis.com/node/agreement/3193にあります。 条項47の状態： ライセンシーは、Esri File Geodatabase APIを使用するソフトウェアまたはWebアプリケーションを開発して、ライセンシーのエンドユーザーに配布できます。 ライセンスに関する最良の議論は、2011年1月からのPaul Ramseyのブログ投稿です。これは、次の最終更新で終わります。 更新4：Esriから、最終ライセンスを確認して曖昧さがないことを確認し、APIが任意のアプリケーションカテゴリのアプリケーションおよび派生製品を自由に再配布可能かつロイヤリティフリーで使用できるという意図を反映することを確認するメールを受信しました。そのため、現在のライセンスに曖昧さがある限り、最終ライセンスのように赤旗と見なされるべきではありません。 Esri APIはGDALのダウンロードまたはインストールの一部として利用できるようになりましたか、いつ利用できますか？また、APIを含むアプリケーションをいつ配布できますか？

23 gdal  ogr  license  file-geodatabase-api 

ピラミッド/概要はGeoTIFF標準の一部ではありませんが、多くのツールがピラミッド/概要の作成をサポートしています。たとえば、vips / nip2、Orfeo Toolbox（otb）、およびossimは、これらを作成するための何らかのサポートを約束します。しかし、それらがすべて他の人がサポートする形式でファイルを作成するかどうかはわかりません。ツールのドキュメントを調べても、それについては何も言及されていません。 vipsは地理空間ツールではなく、エンドユーザーフレンドリーなドキュメントは見つかりませんでしたが、IIPImageのドキュメントでは「Tiled Pyramidal TIFF」のサポートについて言及されています：http : //iipimage.sourceforge.net/documentation/images/ otbは、「Multi Resolution Pyramids」の形式や仕様については何も言及していません：https : //www.orfeo-toolbox.org/CookBook/CookBooksu65.html ossimは、「解像度の低下したデータセット」のさまざまな出力形式をサポートしていると言いますが、実際の意味はわかりません：https : //trac.osgeo.org/ossim/wiki/img2rr GDALは、その「概要画像」に関することも実際には指定していません：http : //www.gdal.org/gdaladdo.html したがって、それらはすべてピラミッド/概要を持っていますが、それらが相互互換性があるかどうかは明確ではありません。 より一般的なページでは、次の引用を見つけました。 http://iipimage.sourceforge.net/documentation/images/は言う タイル型多重解像度（またはタイル型ピラミッド型）TIFFは、タイル型のマルチページTIFF画像であり、各解像度はTIFF内の個別のレイヤーとして保存されます。これは標準のTIFF拡張機能であり、Photoshop、GIMP、VIPS、ImageMagickなどのほとんどの画像処理アプリケーションでサポートされています。libtiffコーデックライブラリは、このような画像を完全に読み書きすることもできます。 それは誰もが使用するものですか？ 米国議会図書館はまた、いくつかの情報を持っている：http://www.digitalpreservation.gov/formats/fdd/fdd000237.shtml 彼らは注意してください。 異なるアプリケーションによって作成されたPyramid TIFFファイルは、必ずしも同じ構造ではありません。特に、JHOVEによる分析とImageMagickの識別コマンドから判断すると、AdobeのPhotoshopとImage Magickは異なる内部TIFF構造を持つファイルを生成します。どちらの場合も、TIFFを処理できるほとんどのソフトウェアは、問題なくプライマリ（フルサイズ）TIFFを認識するようです。 それで、それらのフォーマットはどこかで標準化され、指定され、文書化されていますか？他のツールと互換性のある方法でそれらを生成できるツールを見つけるにはどうすればよいですか？概要/ピラミッドに地理空間タグがありますか、それともソフトウェアを使用して画像データにタグを作成できますか？

22 gdal  geotiff-tiff  file-formats  image-pyramids  gdaladdo 

Pythonスクリプトでogr2ogrユーティリティを使用してGMLをESRIシェープファイルに変換しようとしています。 osgeoを介してGDAL / OGRパッケージを正常にインストールしましたが、現在、pythonでogr2ogrを使用するための構文の詳細を見つけて理解するのに苦労しています。 OGRをインポートした後、私が見つけたのは ogr2ogr -f "ESRI Shapefile" output.shp input.gml これをPythonでどのように使用しますか？入力ファイルと出力ファイルを割り当てるのと同じくらい簡単ですか？ 私がしようとするすべては、単に構文エラーを取得します。私を始めるための正しい方向へのポインタは素晴らしいでしょう。

22 python  shapefile  gdal  gml  ogr2ogr 

シェープファイルには1つのジオメトリタイプしか指定できませんが、他の多くの形式（kml、geojson）には複数のタイプを指定できます。 シェープファイルに変換する場合、ogr2ogrに各タイプに1つずつ複数のシェープファイルを作成するように指示することは可能ですか？ 「-where」オプションは便利かもしれませんが、マニュアルページには構文が説明されていません。 たとえば、失敗するコマンドは次のとおりです。 $ wget http://a841-tfpweb.nyc.gov/jackson-heights/wp-content/themes/tfp/kml/transit.kml $ ogr2ogr -f "ESRI Shapefile" transit.shp transit.kml Warning 6: Normalized/laundered field name: 'Description' to 'Descriptio' ERROR 1: Attempt to write non-linestring (POINT) geometry to ARC type shapefile. ERROR 1: Terminating translation prematurely after failed translation of layer Transit (use -skipfailures to skip errors)

22 shapefile  gdal  kml  geojson  ogr2ogr 

特定のラスターデータセットからタイル（jpeg）を作成する必要があるアプリに取り組んでいます。私の最初の傾向はgdal2tiles.pyを使用することでしたが、そのパフォーマンスは改善できるようです。 そこで、次の質問をすることになります。GDALを使用してラスターデータセットから「タイル化された画像」を切り取る最も速い方法は何でしょうか。また、「タイル画像」とは、単純なjpegまたはpngを意味します。 プロトタイプでは、MapServer C＃MapScriptバインディングを使用してジョブを実行します。つまり、マップオブジェクトを作成し、カットする必要があるすべてのタイル境界をループし、マップオブジェクトの範囲を設定して、結果の画像を保存します。このアプローチのパフォーマンスはgdal2tiles.pyのパフォーマンスよりも大幅に優れていますが、GDALをさらに高速にできない場合は、GDALを直接使用するのではないかと考えています。GDALメソッドで同様のワークフローを提案できる人はいますか？ 編集：本日さらに調査を行ったところ、目の前で答えが見つかりました。FWToolsをダウンロードした場合、csharp \ appsフォルダーには、C＃GDALバインディングを示すためのいくつかのクラスが含まれています。私の場合、GDALRead.csとGDALReadDirect.csは私が探していたものでした。

22 gdal  tiles 

Python + GDALを使用して、ESRIファイルジオデータベースのベクターデータセットにアクセスしようとしています。ファイルジオデータベースAPIを使用してGDALを正常にコンパイルしました。FileGDBドライバーは、 ogrinfo --formats FileGDBドライバーを表示して ogrinfo myfilegdb.gdb データベースの内容に関する正しい情報を教えてくれます。 ただし、Pythonでコンテンツ自体にアクセスする方法はわかりません。シェープファイルにアクセスするには、次のように書きます。 driver = ogr.GetDriverByName('ESRI Shapefile') ds = driver.Open('shapefile.shp', 0) FileGDBフィーチャクラスにアクセスするとき、次のコマンドを使用すると仮定します。 driver = ogr.GetDriverByName('FileGDB') ds = driver.Open('myfilegdb.gdb/feature_class', 0) しかし、これはデータセットを特定/特定できないため、機能していないようです。ESRI FileGDBから個々のフィーチャクラスを呼び出す方法を知っている人はいますか。 Ubuntu 12.04 x64でPython 2.7、GDAL 1.9.1、filegdb api 1.2を使用しています。提案をありがとう！

21 python  gdal  vector  file-formats 

GDAL pythonバインディングを使用して再プロジェクト/リサンプリングを試みていますが、コマンドラインユーティリティからの結果と比較してわずかに異なる結果が得られていますgdalwarp。 短い例については、以下の更新を参照してください このスクリプトは、Pythonのアプローチを示しています。 from osgeo import osr, gdal import numpy def reproject_point(point, srs, target_srs): ''' Reproject a pair of coordinates from one spatial reference system to another. ''' transform = osr.CoordinateTransformation(srs, target_srs) (x, y, z) = transform.TransformPoint(*point) return (x, y) def reproject_bbox(top_left, bottom_right, srs, dest_srs): x_min, y_max = top_left …

20 python  gdal  gdalwarp  numpy 

シェープファイルのフィールド名（属性）の名前を変更する必要があります。コマンドラインからこれを行う簡単な方法はありますか？

20 shapefile  gdal  convert  ogr 

PostGISサーバーからの地理データを処理するプログラムの作成を開始し、GDALを使用します。 あなたの経験によると、それを使ってプログラムするのに最適な言語はありますか？PerlとPythonでプログラミングできますが、再利用するコード、サンプル、ライブラリなどがより多く含まれている言語を知りたいです。

20 python  postgis  gdal 

Pythonで開発し、OSGEOのGDALを使用して、ラスターとシェープファイルを操作および操作しています。 ポイントフィーチャを含むシェープファイルを取得して、サーフェスラスタに補間したいです。現在、私は「RasterizeLayer」メソッドを使用しています。このメソッドは、ポイントフィーチャからラスター（すべてのnodata値で設定されている）に値を焼き付けますが、すべてのタッチされていないピクセルを「nodata」値として残します。したがって、チェッカーボードタイプのラスターが残っています。 RasterizeLayerを使用した後の内容： 最終製品に必要なもの： 私が探している関数は、arcgisscripting importの 'Spline_sa（）'として知られていると思います。 GDALには同様の機能がありますか、または希望する出力を得るための別の方法がありますか？

20 python  raster  gdal  interpolation 

Web C＃アプリケーションでシェープファイルに投影を割り当てたいため、Gdal / OGR / OSR C＃バインディングを使用し、osr_csharp.dllおよびogr_csharp.dll参照を追加します。しかし、投影をシェープファイルに割り当てるためのコードを書いているときに、いくつかの困難に直面しています。これらのライブラリのクラス/メソッドを理解するのに役立つAPIドキュメントがあります。または、誰かが投影をshapefileに割り当てるためのコードを持っているなら、それは私にとって大きな助けになるでしょう。

20 c#  shapefile  coordinate-system  gdal  ogr