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

現在、次のスクリプトを使用して、テーブルからいくつかの属性データを多数の個々のシェープファイルに追加しています。 library(rgdal) specieslist <- read.csv("SpeciesList1.txt", header=F) attdata <- read.table("TestAtt.csv", sep = ",", header=T) for (n in 1:dim(specieslist)[1]) { speciesname <- specieslist[n,1] shp <- readOGR("Mesoamerica_modified_polygons", speciesname) shp$ENGL_NAME<-attdata[n,2] writeOGR(shp, "PolygonsV2", speciesname, driver="ESRI Shapefile") } 最後に次の警告が表示されます。 1: In writeOGR(shp, "PolygonsV2", speciesname, driver = "ESRI Shapefile") : Field names abbreviated for ESRI Shapefile driver …

16 shapefile  gdal  r  ogr  attribute-table 

同じマップに表示したいgdal2tilesを使用して生成された複数のタイルソースがあります。そのため、タイルを提示するときは、どのソースからタイルを提供するかを決定し、境界をチェックする必要があります。 （gdal2tilesによって生成されたファイル構造から）ズーム、x、およびyのみに基づいて単一のタイルの緯度経度境界を計算する方法を知っていますか？ところで：計算を支援するためにAPIからいくつかの機能を呼び出す必要がある場合に備えて、私はGoogle Maps API v3を使用しています。 私が純粋にズーム、x、およびyに制限されている理由は、タイルがマップ上のオーバーレイによって呼び出されるだけでなく、マップ外に印刷できるカスタム印刷機能もあるためです。

16 gdal  google-maps  gdal2tiles 

作業中のDEMの大きなセットを受け取ったばかりで、それらから輪郭を生成したいと思います。DEMの解像度は1m、サイズは1kmx1kmです。 gdalinfoからの出力： Driver: AAIGrid/Arc/Info ASCII Grid Files: 380000_6888000_1k_1m_DEM_ESRI.asc Size is 1000, 1000 Coordinate System is `' Origin = (380000.000000000000000,6888000.000000000000000) Pixel Size = (1.000000000000000,-1.000000000000000) Corner Coordinates: Upper Left ( 380000.000, 6888000.000) Lower Left ( 380000.000, 6887000.000) Upper Right ( 381000.000, 6888000.000) Lower Right ( 381000.000, 6887000.000) Center ( 380500.000, 6887500.000) Band …

16 raster  gdal  dem  contour 

例として、次のタイルhttp://a.tile.openstreetmap.org/3/4/2.pngを取得し、「4_2.png」として保存します。 このタイルのWGS84座標を計算するか、読み取ることができますが、対応するタイルをクリックして： 0 66.51326044311185 45 40.97989806962013 (West North East South) タイルを正しくジオリファレンスする方法（gdalを使用してジオティフまたは別のジオリファレンス形式を生成する）： ビットマップを引き伸ばす必要はありません（=ジオティフのピクセルは元のビットマップのピクセルとまったく同じです） 結果の画像は、GISビューア/エディタの適切な場所で開かれますか（TatukGIS Free Viewerなど）？ （2011年9月19日に編集して、質問を明確にし、結論を含めました） 私の結論： 私は最初に、3番目のアイデア（以下を参照）が正しいものであると考えました。GISビューアーでジオティフを開き、表示された座標を予想と比較しました。2番目のアイデアのジオティフは、2ピクセル北にシフトしているようです。だからこそ、私はアイデア3（または4）を正しいものと考えました。 しかし、はるかに高いズームレベルでタイルを試してみると、アイデア3からのジオティフは決定的に南にシフトします。ズームレベル3のタイル上の座標を比較するのはばかげていました。そのようなズームレベルでの国境は、比較が良い結果をもたらさないように単純化されています。 Dan S.は正しかった、タイル画像はすでにEPSG：3857にあります。2番目のアイデアは正しいアイデアです（高ズームレベルでも良い結果が得られます） 最初のアイデア：EPSG：4326 WGS84座標のEPSGコードはEPSG：4326です。そのため、単にgdal_translateを使用してWGS84座標を使用して、タイルをgeotifとしてジオリファレンスします。 gdal_translate -of Gtiff -co tfw=yes -a_ullr 0 66.51326044311185 45 40.97989806962013 -a_srs EPSG:4326 4_2.png t4326.tif 結果のマップは適切な場所に表示されますが、投影が正しくなく、タイルの中央にシフトがある可能性があります。gdalwarpを使用してマップを再投影して確認するために長い時間を試した後、Global Mapperのデモバージョンをダウンロードしました。EPSG：4326座標を使用できるように、画像をストレッチする必要があります。 2番目のアイデア：EPSG：3857 このタイルは、「web mercator」プロジェクション（エイリアスgoogleマッププロジェクション）を使用します。EPSG コードはEPSG：3857（エイリアスEPSG：900913）になりました。私は単にgdaltransformを使用して座標を変換します： gdaltransform -s_srs EPSG:4326 -t_srs EPSG:3857 0 66.51326044311185 …

16 web-mapping  gdal  epsg  tiles 

次の質問がやや馬鹿げている場合はおologiesび申し上げますが、私はこのGISのことをまったく知りません。 Pythonでgdalを使用して、投影されたgeoTiff画像をWGS84に変換しようとしています。次のようなものを使用して、投影されたGeoTiffs内のポイントを変換するプロセスの概要を示す投稿を見つけました。 from osgeo import osr, gdal # get the existing coordinate system ds = gdal.Open('path/to/file') old_cs= osr.SpatialReference() old_cs.ImportFromWkt(ds.GetProjectionRef()) # create the new coordinate system wgs84_wkt = """ GEOGCS["WGS 84", DATUM["WGS_1984", SPHEROID["WGS 84",6378137,298.257223563, AUTHORITY["EPSG","7030"]], AUTHORITY["EPSG","6326"]], PRIMEM["Greenwich",0, AUTHORITY["EPSG","8901"]], UNIT["degree",0.01745329251994328, AUTHORITY["EPSG","9122"]], AUTHORITY["EPSG","4326"]]""" new_cs = osr.SpatialReference() new_cs .ImportFromWkt(wgs84_wkt) # create a transform object …

16 python  gdal  coordinate-system  wgs84 

これは簡単なことのように思えますが、私の人生ではそれを理解することはできません。 単一のバンドを含む既存のGeoTiffがある場合、そのGeoTiffに別のバンドを追加するにはどうすればよいですか？GDALCreate（）を使用して追加のバンドを持つ新しい出力ファイルを作成し、既存のバンドを元のファイルから新しいファイルにコピーする唯一の方法はありますか？これが私がこれまでに考案した唯一のソリューションですが、既存のデータセットにバンドを追加する簡単な方法があるはずです。

15 raster  gdal  band 

OK VRT形式について読んでいます。どうか、これについてのあなたの意見をお願いします。また、私が提案することが理にかなっている場合は： VRTは基本的に、ピクセルサイズ、ジオロケーションなど、実際のラスターファイルのさまざまなプロパティを記述するメタデータXMLファイルです。 ステップ1：すべてのGeoTIFFファイルをVRTに変換しました（実際にはVRTを生成しただけではありません）。 これで、VRTをマージできることがわかったので、すべてのメタデータを保持する1つの大きなVRTのVRTメタデータファイルのモザイクを作成できると思います。素晴らしい、つまり、メタデータはおそらくレンダリングの最適化に使用されます。QGISは、どのタイルが現在のキャンバス範囲に収まり、どのタイルが収まらないかを（仮定して）決定するためです。 ステップ2：VRTモザイクを作成する さて、次の利点は、QGISの目次に1つのエントリを追加するだけで、必要に応じてスケール依存のレンダリングやnullデータ値の透過性など、さまざまなものを設定できることです。議論の余地があります。GeoTIFFをQGISにロードしてグループに入れることはできますが、共通のSRSを設定する以外はグループを制御できません。 オルソフォトを小さなスケール（大きなズームアウト）で表示したいが、レンダリングを高速化したい場合は、元のラスターのリサンプリングセットのみでステップを繰り返し、いくつかの大きなVRTデータセットを作成します。 これは理にかなっていますか？:)

15 qgis  gdal  vrt 

短い： ラスターの解像度を変更し、以下の画像に示すようにグレーのランプを滑らかにします。GDAL、PIL、またはNumpyを使用することをお勧めします。 説明： 私は、高性能地球統計ライブラリーを使用して、20メートルの出力解像度でラスターにポイントをクリギングしています。補間時間が指数関数的に増加するため、出力解像度を変更したくありません。 この解像度では、出力画像は見苦しくなります（ピクセル化およびエイリアス化）。概念的に正しいかどうかはわかりませんが、次の例のように画像をより滑らかにしたいです。これは、画像をより良い解像度に「再補間」するようなものです。私はpythonを使用しているので、私の好みはGDAL、Python Imaging LibraryまたはNumpyです。答えは、アルゴリズム名やこの種の操作の概念を指摘するような理論的なものです。 ソース： 先： gdalwarpキュービックスプラインを使用した結果の編集：

15 python  raster  gdal 

シェープファイルをインポートするためにPythonでGDALを使用します。ファイルのフィールド名を知りたいのですが、現在の方法は次のとおりです。 fields = [] for i in range(1, layer.GetFeature(0).GetFieldCount()): field = layer.GetFeature(0).GetDefnRef().GetFieldDefn(i).GetName() fields.append(field) しかし、この方法で、最初のレイヤーの機能を取得しています。レイヤーごとに機能が異なる可能性があるということですか？ そうでない場合、これほど深くに入る代わりに、フィールド名を一度に取得することは可能ですか？はいの場合、フィールド名を取得する簡単な方法はありますか？

15 python  shapefile  gdal 

Windows XPでPython 2.7.3のgdalをインポートしようとしています： >>> import gdal Traceback (most recent call last): File "<stdin>", line 1, in <module> File "C:\Python27\lib\site-packages\gdal.py", line 2, in <module> from osgeo.gdal import deprecation_warn File "C:\Python27\lib\site-packages\osgeo\__init__.py", line 21, in <mo _gdal = swig_import_helper() File "C:\Python27\lib\site-packages\osgeo\__init__.py", line 17, in swit_helper _mod = imp.load_module('_gdal', fp, pathname, description) ImportError: DLL …

15 python  gdal  windows-xp 

Pythonで機能ごとにシェープファイルを分割することは可能ですか？（最良の解決策は、結果のベクトルオブジェクトを一時的にディスクではなくメモリに保存できるソリューションです）。 理由：シェープファイルのいくつかの異なるサブセットでgdal rasterizeLayer関数を使用したい。この関数には、osgeo.ogr.Layerオブジェクトが必要です。 まあ、私は少し試してみましたが、次のように機能します。次のように、フィーチャごとにgdalレイヤーオブジェクトのジオメトリを取得できます。 # Load shape into gdal shapefile=str(vectorPath) layer_source = ogr.Open(shapefile) lyr = layer_source.GetLayer(0) for i in range(0,lyr.GetFeatureCount()): feat = lyr.GetFeature(i) ge = feat.geometry() このジオメトリに基づいてosgeo.ogr.layerオブジェクトを作成する方法を知る必要があります。 明確にするために。プレーンogr / gdalコードの関数が必要です！これは他の人にも興味があるようであり、セカンダリモジュールなしのソリューションが引き続き必要です（ただし、ここから来るソリューションは無料のqgis​​プラグインで使用されます）。

15 python  gdal  rasterization  osgeo  split-by-attribute 

QGISによって動的に生成された画像を表示するWebポータルを開発したいのですが、QGISが提供するTIFF画像のサイズが大きすぎるため、画像ビューアーを使用してWebサイトに表示することはできません。そのため、これらのTIFF画像をjpg / jpegに変換してサイズを縮小し、ブラウザーフレンドリーにすることで、これらの画像をWebサイトに簡単に表示できるようにしました。しかし、gdalコマンドを使用してTIFFイメージをjpg / jpegに変換する方法がわかりません。 この点で私を助けてください。

15 qgis  gdal  geotiff-tiff 

ソースからビルドせずにCentOS7にGDAL 2.Xをインストールしようとしています。別のサーバーでGDAL 2.Xを構成する必要があり、ソースコードによるインストールは非常に時間と手間がかかる作業です。 opengeoやELGISなどのリポジトリをいくつか見つけましたが、CentOS7用です。EPELリポジトリでは、GDAL 2.xが利用できません。 [root@localhost ~]# yum repo-pkgs epel list gdal* (...) Available Packages gdal.x86_64 1.11.4-1.el7 @epel gdal-libs.x86_64 1.11.4-1.el7 @epel gdal-python.x86_64 1.11.4-1.el7 @epel gdal-devel.x86_64 1.11.4-1.el7 @epel gdal-doc.noarch 1.11.4-1.el7 @epel gdal-java.x86_64 1.11.4-1.el7 @epel gdal-javadoc.noarch 1.11.4-1.el7 @epel gdal-perl.x86_64 1.11.4-1.el7 @epel バージョン1.11は、1.xメジャーバージョンの最新版ではありません。 任意のリポジトリからインストールすることはできますか、それともソースからのみインストールできますか？そうでない場合、RPMとしてビルドして複数のサーバーに簡単にインストールできますか？

15 gdal  centos  repository 

それぞれ約100 MBの約20個のdem tiffをマージしたかったので、gdalwarpを使用しました（仮想ファイルオプションについて知っています）。 gdalwarp $(list_of_tiffs) merged.tiff ただし、このコマンドには非常に長い時間がかかり、約15分後にセットから5番目の画像を処理していました。終了します。 次に、gdal_merge.pyスクリプトがモザイクにも使用できることを発見し、試してみました。 gdal_merge.py $(list_of_tiffs) 3分以内に完了しました。 両方のコマンドで同じ結果が得られると思っていたので、結果が同じであるとgdalwarpに時間がかかるのはなぜですか？

14 gdal  gdalwarp  gdal-merge 

pypyを使用してpythonコードを最適化しようとしています。ほとんどのライブラリをインポートして実行できますが、GDALライブラリをインポートすることはできません。Python 2.7でインポートした場合、同じライブラリは問題なく実行されます。 私はGDALの場所をPYTHONPATHに追加しようとしましたが、GDAL PythonバインディングをPyPyディレクトリにダウンロードしようとしても、pypyはGDALを認識したくないため、次のエラーが発生します。 Traceback (most recent call last): File "app_main.py", line 51, in run_toplevel File "ImportENC-P.py", line 50, in <module> from osgeo import gdal, ogr, osr File "B:\PYPY\site-packages\osgeo\__init__.py", line 21, in <module> _gdal = swig_import_helper() File "B:\PYPY\site-packages\osgeo\__init__.py", line 13, in swig_import_helper import _gdal ImportError: No module named _gdal

14 python  gdal