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

Python GDALモジュールの公式ドキュメントはありますか？ さまざまなクックブックがありますが、最も基本的な関数の各関数または引数に関する情報はほとんどありません。

9 python  gdal  documentation 

ラスターを投影しようとしています。RにはprojectRaster()これを行うための関数があります（完全に再現可能な例の下に）： # example Raster require(raster) r <- raster(xmn=-110, xmx=-90, ymn=40, ymx=60, ncols=40, nrows=40) r <- setValues(r, 1:ncell(r)) projection(r) # project to newproj <- "+init=epsg:4714" # using raster package to reproject pr1 <- projectRaster(r, crs = CRS(newproj), method = 'bilinear') それはうまくいきます。しかし、それはかなり遅いです。 速度を上げるために、gdalwarp代わりに使用します（SSD を使用すると、ディスク/ Rへの読み書きのコストはそれほど高くありません）。 ただし、をprojectRaster()使用した結果を再現できませんgdalwarp。 # using gdalwarp to reproject tf …

9 raster  gdal  r  coordinate-system  gdalwarp 

休業。この質問には、より焦点を当てる必要があります。現在、回答を受け付けていません。 この質問を改善してみませんか？質問を更新して、この投稿を編集するだけで1つの問題に焦点を当てます。 2年前休業。 OGR / GDALを使用してPythonでスクリプトを作成しています。 シェープファイルのセットとGeoTiffラスターファイルのセットがあります。 シェープファイルがラスター領域と交差しない場合、スクリプトに無視させたいのですが。 シェープファイルは長方形ではないため、layer.GetExtent（）から返されたxmin / xmax、ymin / ymaxの値を単純に比較することはできません。全体的な形状を表す実際のポリゴンが必要です。次に、そのポリゴンがラスタースクエアと交差するかどうかを判断する方法が必要です。 どういうわけか、シェープファイル内のすべてのポリゴンを1つのフィーチャにマージし、そのフィーチャのジオメトリを読み取って、その情報をラスター範囲と比較できると考えていました。ただし、具体的にどのように実行するかはわかりません。 シェープファイルからボーダーポリゴン情報を抽出する方法は？ そのポリゴンが特定の正方形の領域と交差するかどうかを判断するにはどうすればよいですか？

9 python  raster  shapefile  gdal  ogr 

私はいつものように現在のQGISマスターをコンパイルしてインストールし、postgres 9.3データベースにpostgis 2.0.2拡張子を付けていくつかのラスターデータをインポートしました。データベースマネージャーを使用してラスターデータをqgisにインポートしようとすると、次のメッセージが表示されます。 GDAL provider: Cannot get GDAL raster band: Raster layer: Provider is not valid (provider: gdal, URI: PG: dbname=gis host=127.0.0.1 user=xxx password=xxx port=5432 mode=2 schema=public table=elevation コマンドラインで、 ERROR 1: Computed PostGIS Raster dimension is invalid. You've probably specified unappropriate resolution. エラーを検索する場所がわかりません。 Postgis 2.1.2に移動したところ、同じ問題が見られます。

9 qgis  postgis  raster  gdal  raster2pgsql 

私はマイクボストックのチュートリアル「マップを作ろう」に従っています。彼はogr2ogr、機能の属性に基づいてシェープファイルをフィルタリングする方法を説明しています。QGISでシェープファイルを開き、属性テーブルを見て、それらの機能を確認できます。コマンドラインからその情報を見つける方法はありますか？

9 shapefile  gdal  ogr2ogr 

私はpythonとQGIS 2.0を使用しています。1つのポリゴンフィーチャによってフォルダー内のラスターをクリップしようとしています。私が「PyQGIS」を使用するのはこれが初めてです。以前はArcpyに慣れていました。とにかく、簡単なスクリプトが機能しないので、どんな提案でもいただければ幸いです。 import qgis.core, qgis,utils QgsApplication.setPrefixPath("C:/OSGeo4W64/apps/qgis", True) QgsApplication.initQgis() CLIP= "C:/Users/unim/Documents/Umberto/Universita/PhD/Guglielmin/Permafrost/Alta_Valtellina/Landsat_ita/study_area_foscagno.shp" INPUT_FOLDER="C:/Users/unimi/Documents/Umberto/Universita/PhD/Guglielmin/Permafrost/Alta_Valtellina/Landsat_ita/LE71930282000259EDC00" OUTPUT= "C:/Users/unim/Documents/Umberto/Universita/PhD/Guglielmin/Permafrost/Alta_Valtellina/Landsat_ita/foscagno_pyqgis/" for RASTER in INPUT_FOLDER.tif do echo "Processing $RASTER" gdalwarp -q -cutline CLIP -crop_to_cutline -of GTiff RASTER OUTPUT+ "clip_"+ RASTER done QgsApplication.exitQgis() 以下は、今までに行った改善ですが、スクリプトを機能させていませんが、近づいていると思います... import qgis.core, qgis.utils, os, fnmatch from osgeo import gdal CLIP= "C:/Users/unimi/Documents/Umberto/Universita/PhD/Guglielmin/Permafrost/Alta_Valtellina/Landsat_ita/study_area_foscagno.shp" INPUT_FOLDER= "C:/Users/unimi/Documents/Umberto/Universita/PhD/Guglielmin/Permafrost/Alta_Valtellina/Landsat_ita/LE71930282000259EDC00/DNs2Reflectance_LE71930282000259EDC00" OUTPUT= "C:/Users/unimi/Documents/Umberto/Universita/PhD/Guglielmin/Permafrost/Alta_Valtellina/Landsat_ita/Cloud_mask_AltaValtellina/clip_2_foscagno" …

9 gdal  pyqgis  gdalwarp 

大きなECW（詳細は以下）をクリップしようとしていますが、ラスターファイルが大きすぎて完全に処理できません。 以下のECWの詳細 ドライバー：ECW / ERDAS圧縮ウェーブレット（SDK 5.0） ファイルサイズ：50 GBサイズは450000、565081ピクセルサイズ：0.15 0.15 COLORSPACE = RGB COMPRESSION_RATE_TARGET = 9 VERSION = 2バンド数：4 クリップしたい領域は、元のファイルの約1/5です。 ここに私が成功せずに試した方法があります： Arcgisを使用してecwをtiff /他の形式に保存しました...（すぐにあきらめました） Qgisとそのクリッパーツールを使用しました...ファイルの作成は約40％で止まっています。 Qgis以外のオプションでOSGeo4Wのgdal_translateを使用しました。（Qgisを使用しないで一部のメモリを解放する可能性があるという考えがトリックになるだろうと考えてみました） 画像をバラバラに切り、必要なものを取得すると考えてgdal_retileを使用しました。コマンド「gdal_retile -ps 10000 10000 -of ecw -tileIndex tile.shp -targetDir input.ecwこれはさらに速くクラッシュしました」 誰かがアイデアを持っていますか？ 参考までに、16GBのRAMを搭載したi5-3470 3.2GhzでWindows 7 64ビットを実行します。

9 raster  gdal  clip  ecw  big-data 

data.gc.caはそのデータの一部をFGDB形式で配布しているため（例：Plant Hardiness Zones）、gdal / QGISで使用できるようにしたいと考えています。ただし、コンパイル/インストールは予定どおりにはいきません。 これらは私が取ったステップです： ubuntugis-unstableのソースリポジトリを有効にしました。 Linux 64ビットファイルジオデータベースAPI 1.3をダウンロードしてにインストールし/usr/local/FileGDB_API、必要なldconfig設定を追加して、テストプログラムを正常にビルドして実行しました。（ここで推奨される手順：[Qgis-user] ESRI FileGDBでgdalをコンパイルする方法） 必要なソース依存関係をインストールしました： sudo apt-get build-dep gdal ソースをプルダウン： apt-get source gdal ファイルの--with-fgdb=/usr/local/FileGDB_API \configureコマンドに行を追加しましたdebian/rules。 ソースディレクトリから発行dpkg-buildpackage -rfakeroot -uc -b（ここでアドバイスされているとおり：apt-getを使用してインストールされたパッケージのソースコードを取得するにはどうすればよいですか？-Ubuntuに問い合わせてください） ビルドはうまくいきますが、パッケージのコンパイル段階の後半で失敗し、fakerootファイルシステムにRuby 1.9ツリーが見つからないという不満があります。 … mv /home/scruss/gdal-1.10.0/debian/tmp/usr/local/lib/site_ruby/1.9/* /home/scruss/gdal-1.10.0/debian/tmp/usr/lib/ruby/1.9/. mv: cannot stat ‘/home/scruss/gdal-1.10.0/debian/tmp/usr/local/lib/site_ruby/1.9/*’: No such file or directory make: *** [install] Error 1 dpkg-buildpackage: error: fakeroot …

9 qgis  gdal  ubuntu  linux  file-geodatabase 

ほんの少し前に、私はアラスカのサイズの地域の毎日の地表面温度衛星画像を扱いました。 約1000枚の画像を処理する必要がありました。それらをダウンロードし、モザイク化し、投影し、シェープファイルの境界にクリップし、クリップされたラスターから平均を抽出し、平均を時間-温度グラフにプロットしました。時間と温度のグラフが最終的な出力です。これらはすべてArcGISで行いましたが、作業の一部を自動化するのに役立つArcPyスクリプトを作成しましたが、それでもまだ骨の折れる作業でした。 同じことを行いますが、100％自動的に実行するプログラムを作成したいと思います。関心のある領域を定義し、時間と温度のグラフである最終的な出力を生成するために、プログラムにコーナーの座標を入力したいだけです。 これは、GDALライブラリとPostGISをPythonでラスターデータストアとして使用することで可能だと思いますか

9 python  postgis  arcgis-desktop  gdal  software-recommendations 

素晴らしいPySALパッケージのマッピングコンポーネントを使用しようとしていますが、何らかの理由でogrモジュールがありません。gdalに同梱されているような印象を受けました。何か不足していますか？ import gdal import ogr --------------------------------------------------------------------------- ImportError Traceback (most recent call last) <ipython-input-29-ef9eb3418665> in <module>() ----> 1 import gdal 2 import ogr ImportError: No module named gdal 更新： 飛行機に乗るのが少し急いでいたので、いくつかの重要な情報を含めることができませんでした。ここで説明したapt-getメソッドを使用してpython-gdalを既にインストールしました。パッケージマネージャーが正しくインストールしなかった可能性はありますか？ choct155@choct155-Q550LF:~$ sudo apt-get install python-gdal Reading package lists... Done Building dependency tree Reading state information... Done python-gdal is already the newest …

9 python  gdal  ogr 

3つの個別の1バンドGeoTIFFファイルがあります。バンドの色の解釈は灰色です。3バンドRGBファイルが必要です。私はgdal_merge.pyを使用しました gdal_merge.py -separate file1.tif file2.tif file3.tif -o output_file.tif しかし、結果のoutput_file.tifの3つのバンドの色の解釈は、灰色、未定義、未定義です。それを除いて、すべてのピクセルは0です。 pct2rgb.pyを実行すると、RGBファイルが表示されます。 pct2rgb.py output_file.tif output_file_rgb.tif しかし、もちろん、ピクセルも0です。だから、私は3つの質問があります： gdal_merge.pyは、3つの1バンドファイルを1つの3バンドRGBファイルに結合するための適切なツールですか？ なぜバンドの色の解釈が未定義になるのですか？ pct2rgbは、この色の解釈を持つ3バンドファイルをRGBファイルに変換するための適切なツールですか？ 更新：ラスタにはカラーテーブルがありません。ただ色の解釈：灰色。 一方、ピクセル値は0から1023になります（これは意図的なものです） より多くのデータ：それらは回転したラスター（北を上にしない）ですが、それらはすべて同じ地理変換を持っています。 更新2：画像をワープして北に上向きにし、VRTを構築し、各バンドにColorInterpを追加しますが、出力結果で色が灰色、未定義、未定義として解釈されます。 問題は、少なくとも最初のバンドでカラーテーブルを作成する必要があることです。私が知っているそれらを作成する方法を、私は私のテーブルには、どのように多くのエントリを持っていなければならないのか分かりません。GDALラスターFAQの例に13のエントリがあるのはなぜですか？役立つ場合、すべてのピクセルの値は0〜1023です。 更新3：どうやら、TIFF形式では各バンドの色の解釈を実際に指定する方法はありません。TIFFファイルを読み取るときにGDALが色の解釈を構築する方法は、PHOTOMETRICタグとEXTRASAMPLESタグの値の組み合わせです。 これらのタグについて読む： PHOTOMETRICは、画像データの色空間を表します。値2は、ピクセル値のコンポーネントがRGBであることを意味しますが、バイトピクセルを想定しており、UInt16ピクセルを持っています（を試行し-co "PHOTOMETRIC=rgb"てエラーが発生しました）。そのため、出力ファイルにPHOTOMETRICタグを指定できません。 EXTRASAMPLESは、各ピクセルにN個の追加コンポーネントがあることを指定します。このタグを使用してマージされたファイルを作成する方法がわかりません。または、必要な場合。 それで、更新2では、ColorTableの作成をお勧めしますが、どうやって？3つの入力ファイルでは、ピクセル値は0〜1023です。それらを色と一致させる必要がありますか？1024入力のColorTableを作成する必要がありますか？どうやって？ 更新3では、マージされたファイルを作成するときにいくつかのGeoTIFFタグを使用できたようですが、実際にそれらを使用できるかどうか、またはどのように使用するかはわかりません。

9 gdal  gdal-merge 

GISのプログラミング側に不慣れな人にとっては、多くの優れたツールやリソースがアクセスできないと感じています。gdalやogrなどのライブラリでジオプロセシングを処理するためのGUIベースのウィザードに興味があります。任意のヒント？

9 gdal  ogr  gui 

私はpython gdalライブラリを使用してラスターを作成しようとしていますが、データが出力されているポイントに到達しましたが、出力データは原点のx軸で反転しています。私は何かを見落としている必要があることを知っていますが、どこが間違っているのか理解できません。何か案は？ ラスターを作成するときに、左上のx / y値を設定すると、配列には左上からインデックスが付けられ、右下に向かってインデックスが付けられているように見えます。以下のコードでは、配列に行の値を入れています。 配列を出力すると、次のようになります。 [[ 1. 1. 1. 1. 1. 1. 1. 1. 1. 1. 1. 1. 1. 1. 1. 1. 1. 1. 1. 1. 1. 1. 1. 1. 1. 1. 1. 1. 1. 1. 1. 1. 1. 1. 1. 1. 1. 1. 1. 1. 1. 1. 1. …

9 python  raster  gdal 

同僚と私は最近、Mapotecaと呼ばれるオープンソースプロジェクトをリリースしました。これは、GDAL / FWToolsインストールと統合して、空間データを探し、範囲、ジオメトリ、要素数、投影を抽出するディレクトリの完全なウォークスルーを可能にするPythonスクリプトのセットです。ユーザーがそれをcsvファイルに書き込む。関連するDBFファイルを説明する部分に取り組んでおり、各シェープファイルのサムネイル画像を生成したいと思っています。Pythonを使用してjpegまたはpngでこれらのサムネイルを生成する方法を知っていますか？

9 python  gdal  ogr  jpg 

20 GB LZW圧縮ジオティフをgdal2tilesを使用してtilepyramidに変換しようとしています。残念ながら、gdal2tilesはこれを処理できないようです。10時間実行した後、まだ最初の「ドット」に達していません。小さい（1GB）ファイルでこれを試したところ、問題なく動作しました。私の質問は、これを可能にする他のソフトウェア（商用またはオープンソース）はありますか？GlobalMapperからタイルを作成できますが、20GBのファイルをメモリにロードできません。

9 gdal  tiles  gdal2tiles