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

グリッド全体を配列として読み取ることなく、gdalまたはpygeoprocessingを使用してポリゴン内のすべてのラスター値を取得する方法を知りたいです。 pygeoprocessingとgdalはゾーン統計を実行できますが、そのような関数から利用できるのはmin、max、mean、stdev、またはcountのみです。ゾーン統計は値にアクセスする必要があるため、同じ方法で値を抽出するのは簡単ですか？ 私はここで非常によく似た質問を見つけました：（NumPyなしのOGRポイントの下でGDALラスターのピクセル値を取得していますか？）しかし、特定の "ポイント"に対してのみです。

11 python  gdal  pygeoprocessing 

誰かがRasterioやGDALを使用して複数の.tifファイルを複数のバンドスタックにスタックする雄弁な方法を持っていますか？ 私はgdal_merge.pyのようなサブプロセスコマンドの使用を避け、それを私のpythonスクリプトの一部として持つ方法を探しています。 RasterioとGDALの両方が.tifファイルを配列として読み取ることを知っていますが、それらの配列をスタックして、結果を別々のバンドとして書き出すにはどうすればよいですか？

11 python  gdal  geotiff-tiff  multi-band  rasterio 

オープンソースのPythonツールを使用して、座標のリストからポリゴンシェープファイルを作成しようとしています。以下は、Python GDAL / OGRクックブックとこのGIS SEの回答から一緒にハッキングされた、これまでのところです。 同様の質問Pythonがあります：X、Y座標のリストからポリゴンシェープファイルを作成する方法？ただし、この質問はの使用に言及していpyshpます。ただし、gdal / ogr Pythonツールのみを使用してポリゴンシェープファイルを作成することに興味があります。 import ogr def create_polygon(coords): ring = ogr.Geometry(ogr.wkbLinearRing) for coord in coords: ring.AddPoint(coord[0], coord[1]) # Create polygon poly = ogr.Geometry(ogr.wkbPolygon) poly.AddGeometry(ring) return poly.ExportToWkt() def write_shapefile(poly, out_shp): """ /gis//a/52708/8104 """ # Now convert it to a shapefile with OGR driver = ogr.GetDriverByName('Esri Shapefile') …

10 python  gdal  polygon  ogr 

複合陰影起伏の概念は、複数の陰影起伏を異なる太陽光の方向と組み合わせて死角を回避することです。3つの陰影起伏のある画像を、それぞれ315⁰（NW日光、デフォルトレイヤー）、+355⁰レイヤー、+275⁰日光とマージする必要があります。以下に示すように、これにより詳細とエレガンスが向上します。GIS以外の完全な説明はここにあります。 任意のsrtmファイルを使用して、3つの最初の陰影起伏を次のように実行できます。 gdaldem hillshade input.tif hillshades_A.tmp.tif -s 111120 -z 5 -az 315 -alt 60 -compute_edges gdaldem hillshade input.tif hillshades_B.tmp.tif -s 111120 -z 5 -az 355 -alt 60 -compute_edges gdaldem hillshade input.tif hillshades_C.tmp.tif -s 111120 -z 5 -az 275 -alt 60 -compute_edges

10 gdal  dem  hillshade  gdal-calc 

投影座標とそれらの座標の実際のピクセル値をGeoTiffファイルから取得し、それらをnumpy配列に保存するにはどうすればよいですか？arsenci020l.tifファイルがあり、その座標はメートル単位です。以下は、私が実行したgdalinfoの要約出力です。 ~$ gdalinfo arsenci020l.tif Driver: GTiff/GeoTIFF Files: arsenci020l.tif arsenci020l.tfw Size is 10366, 7273 Coordinate System is: PROJCS["Lambert Azimuthal Equal Area projection with arbitrary plane grid; projection center 100.0 degrees W, 45.0 degrees N", 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["Lambert_Azimuthal_Equal_Area"], PARAMETER["latitude_of_center",45], PARAMETER["longitude_of_center",-100], PARAMETER["false_easting",0], PARAMETER["false_northing",0], UNIT["metre",1, AUTHORITY["EPSG","9001"]]] …

10 python  gdal  geotiff-tiff  numpy 

データを直接パイプすることは可能ogr2ogrですか？私はcurl、GeoJSONファイルにパイプしてogr2ogr、それを別の形式に変換するためにパイプします。 パイプを使用せずに、最初にファイルを書き込んでから次のコマンドを使用して変換することでこれを実現できましたが、その手順をスキップできるかどうか知りたいです。 curl "https://raw.githubusercontent.com/nvkelso/natural-earth-vector/master/geojson/ne_50m_admin_0_countries.geojson" -o "countries.geojson" && ogr2ogr -f "KML" countries.kml countries.geojson パイプする代わりにcurl出力で変数を作成しようとしましたogr2ogrが、読み取ることができませんでした。 DATA=$(curl -s "https://raw.githubusercontent.com/nvkelso/natural-earth-vector/master/geojson/ne_50m_admin_0_countries.geojson") && ogr2ogr -f "KML" countries.kml $DATA

10 gdal  ogr2ogr 

pythonおよびを使用してラスターを地理参照したいGDAL。私の現在のアプローチは、地上管制点の醜いリストを呼び出しgdal_translateてgdalwarp使用os.systemすることです。私はこれを内でネイティブに行う方法を本当に望んでいますpython。 これは私が使用している現在のプロセスです： import os os.system('gdal_translate -of GTiff -gcp 1251.92 414.538 -7.9164e+06 5.21094e+06 -gcp 865.827 107.699 -7.91651e+06 5.21104e+06 "inraster.tif" "outraster1.tif"') os.system('gdalwarp -r bilinear -tps -co COMPRESS=NONE "outraster2.tif" "outraster3.tif"') インポート後にアクセスできると述べている2012年からの以前の質問と回答gdal_translateがありますgdal。が廃止されているのか、それが間違っているのかfrom osgeo import gdalはわかりませんが、実行しgdal.gdal_translateてもオプションとして表示されません。 それが存在するかどうかはわかりませんが、ラスタをpythonicの方法で変換および再投影できたら嬉しいです。例えば： # translate gcp_points = [(1251.92, 414.538), (-7.9164e+06, 5.21094e+06)] gdal.gdal_translate(in_raster, gcp_points, out_raster1) # warp gdal.gdalwarp(out_raster1, out_raster2, 'bilinear', args*) そのようなアプローチは可能ですか？

10 python  raster  gdal  coordinate-system  gdalwarp 

休業。この質問には、より焦点を当てる必要があります。現在、回答を受け付けていません。 この質問を改善してみませんか？質問を更新して、この投稿を編集するだけで1つの問題に焦点を当てます。 2年前休業。 私はこのように14のgeotiffをマージしようとします： 各geotiffは約50Mbです。出力にgeotiffが必要です 私のワークフロー： gdalbuildvrt -input_file_list list.txt test.vrt （私のリストにはtifの名前が含まれています） 次に： gdal_translate -of Gtiff test.vrt test.tif Input file size is 79841, 59955 動作しますが、結果は13,3 Gbのgeotiffです！14個のファイル（それぞれ50 Mb）について、13 Gbではなく700 Mbのgeotiffを試みました。 gdalはデフォルトでは圧縮しないことを知っているので、次のコマンドを試しました。 gdal_translate -of Gtiff -co COMPRESS=JPEG test.vrt test_compressed.tif しかし、ファイルの「マージ」は、JPEG圧縮には大きすぎます。 Input file size is 79841, 59955 0ERROR 1: JPEGPreEncode:Strip/tile too large for JPEG …

10 raster  gdal  geotiff-tiff  gdal-translate  gdalbuildvrt 

不要な列が含まれているGMLデータがあります。ファイルサイズを小さくするために、これらの列を削除します。（明確にするために：ファイルには多くのポイントがあり、これらのポイントのそれぞれにメタデータが添付されていることを意味します。メタデータのすべてが必要なわけではありません。） 多くのファイルがあるので、コマンドラインからもこれを実行します。 GMLファイルをを使用してさまざまな形式に変換できるogr2ogrので、私の質問は次のとおりです。コマンドラインから地理データファイルから列を削除する最も簡単な方法は何ですか。 私はPythonプログラマーなので、本能はファイルをGeoJSONに変換し、PythonのJSON解析メソッドを使用してJSONファイルの読み取りと再書き込みを行うことです。 しかし、GDALツールなどを使用するより直接的な方法があるかもしれません。

10 gdal  data  geojson  gml  file-size 

英国の雨の日数のデータを含むNetCDFファイルがあります。それをGeoTIFFに変換したいのですが、GDALを正しくジオリファレンスできません。以下のgdalinfoからのダンプを参照してください。 「raindays_tmean_abs」変数を示すシングルバンドTIFFが必要です。「lat」および「lon」変数には、英国にデータを配置するための適切な最小値/最大値があるようですが、gdal_translateで変換すると、アフリカの西海岸に到達します。 最初のディメンションが「時間または垂直」ディメンションではないという警告があります。私はgdalがデータが特定の順序であることを期待していることを読んだので、おそらくこれが問題ですが、それを変更する方法がわかりません。 このデータを変換/再投影するにはどうすればよいですか？ gdalinfo NETCDF：RainDays1_1961-1990_LTA_25km.nc Warning 1: dimension #0 (meaning_period) is not a Time or Vertical dimension. Driver: netCDF/Network Common Data Format Files: RainDays1_1961-1990_LTA_25km.nc RainDays1_1961-1990_LTA_25km.nc.aux.xml Size is 39, 52 Coordinate System is `' Origin = (-17.650008352179277,11.329999312758446) Pixel Size = (0.219999614514803,-0.219999998807907) Metadata: axis_0#units= meaning_period#comment=Meaning periods are all months (12) beginning …

10 gdal  coordinate-system  netcdf 

次の質問は、ArcGISでの航空写真の処理に関するものでした。 表示のみの航空写真を管理するための最も効果的なフォーマット 航空写真の保存、再サンプリング、再投影などの2つの主要なオプションがあるようです。 JP2000 / JP2 / JPEG 2000（最近、GDAL処理用の5つのコード） ECW（ERDAS圧縮ウェーブレット（.ecw）） 他に見逃したことはありますか？ 両方のQGISバージョンに応じて私が理解したことは、通常、いくつかの追加ライブラリをインストールする必要があります。ECWにはいくつかの制限があります-圧縮のためにライセンスを購入する必要がありますか？ 大きなファイルには使用できないjpegをテストしました（最大ディメンションの制限）。また、大きいディメンションでも遅いです。 答えには以下が含まれている必要があります。 QGIS 2.0.1デスクトップやOSGEOでデフォルトで利用できるものは何ですか？ 大きなファイルでどのように機能するか-ズームイン/ズームアウト（ピラミッド）？ ある作成オプションは - 決議 JP2ピラミッドのために？

10 qgis  gdal  ecw  aerial-photography  jpeg-2000 

arcpyとgdalモジュールの両方を使用するArcGIS 10.1のPythonウィンドウを使用して、いくつかのPythonコードを実行しようとしています。ただし、gdalモジュールをインポートしようとすると、エラーが発生します。 ImportError: No module named osgeo 明らかにモジュールを見つけることができないので、メインのpython site-packagesディレクトリをsys.pathリストに追加しました： sys.path.append(r"C:\Python27\lib\site-packages") しかし、私は試してみて、実行したときimport osgeoやfrom osgeo import gdal、私は次のエラーを取得します： Runtime error Traceback (most recent call last): File "<string>", line 1, in <module> File "C:\Python27\lib\site-packages\osgeo\__init__.py", line 21, in <module> _gdal = swig_import_helper() File "C:\Python27\lib\site-packages\osgeo\__init__.py", line 17, in swig_import_helper _mod = imp.load_module('_gdal', fp, pathname, description) ImportError: …

10 arcgis-desktop  arcgis-10.1  gdal  importerror  python-window 

クエリが失敗した後にOGRでPostgreSQL出力エラーメッセージをキャッチすることは可能ですか？ 例えば conn.ExecuteSQL('SELECT * FROM non_existing_table;') Postgresではこれはエラーメッセージを返します： "Error: relation 'non_existing_table' does not exist." ogrでは標準を取得します。 "Error: current transaction is aborted, commands ignored until end of transaction block" Postgresのエラーメッセージをキャッチ/解析できれば、おまけです。私はSWIG / CPLについて読んだことを考えると、これは少し長いショットだと思いますが、誰かがアイデアを持っているかどうかを確認するように頼む価値があります。理想的にはこれはPythonで行われますが、C / C ++コメントも使用します。 ありがとう。

10 python  postgis  gdal  postgresql  ogr 

どうすれば画像を取得して、オープンソースソフトウェアで自動的にオルソ補正できますか？各画像に関連する情報は次のとおりです。 カメラスペック カメラの緯度、経度、高度 カメラの向き：UAVを下向きにして、ロールピッチとヨーを指定 geotiffファイルを出力し、うまくいけば画像をつなぎ合わせて地図を作成したいと思います。

10 raster  gdal  osgeo  photogrammetry  unmanned-aerial-vehicle 

いくつかのテストを行うために、Ceroからマルチスペクトル画像を作成したいと考えています。塩とコショウのノイズが乗った5つの完全に均一なバンド、または中央に異なる値の2乗のような本当にシンプルなもの。明らかに、これは行列のスタック、つまり多次元配列であり、生成は非常に簡単です。私はpythonとgdalを使用してこれを実現したいのですが、gdalはかなり気密性が高く、コツがつかめません。geotiffファイルを作成するのが理想的です。誰かがこれを手伝ってくれませんか？いくつかのポインタまたは非常に穏やかなgdalチュートリアル？皆さん、ありがとうございました。

10 python  raster  gdal  imagery