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

今日、私はISSからヨーロッパを見下ろす写真を偶然見つけました。それは本当に美しいです。 だから、私は写真にさらに追加できると考え始めました。世界の国境、都市、道路網、すべてのもの。楽しみのために。 問題は、proj4衛星投影を使用したセンサーの方位角または高さについてほとんど何も知らないことです その画像の投影を割り当てるのを手伝ってもらえますか？ ここにソースがあります（大きい！）

11 gdal  coordinate-system  georeferencing  satellite 

HDFからGeoTiffへの変換に問題があります。 このWebページに記載されている手順に従っています。問題は、gdaltranslateコマンドを実行するためのNDVIサブデータセット名（XXXX）が見つからないことです。 gdal_translate -of GTiff HDF4_EOS:EOS_GRID:"MOD13C2.A2001001.005.2007078152825.hdf" :MOD_Grid_monthly_CMG_VI:XXXX modis_ndvi01.tif これは、ファイル（MOD13C2.A2001001.005.2007078152825.hdf）の1つのgdalinfo出力です。 E:\GDAL>gdalinfo MOD13C2.A2001001.005.2007078152825.hdf Driver: HDF4/Hierarchical Data Format Release 4 Files: MOD13C2.A2001001.005.2007078152825.hdf Size is 512, 512 Coordinate System is `' Metadata: HDFEOSVersion=HDFEOS_V2.9 LOCALGRANULEID=MOD13C2.A2001001.005.2007078152825.hdf PRODUCTIONDATETIME=2007-03-19T19:28:25.000Z DAYNIGHTFLAG=Both REPROCESSINGACTUAL=reprocessed LOCALVERSIONID=5.2.1 REPROCESSINGPLANNED=further update is anticipated SCIENCEQUALITYFLAG=Not Investigated AUTOMATICQUALITYFLAGEXPLANATION=No automatic quality assessment is performed in the PGE AUTOMATICQUALITYFLAG=Passed …

11 gdal  modis 

QGISの「マスター」（ナイトリービルド）バージョンを使用して、コンパイルせずにfgdbサポートを取得できることを望んでいました（私は開発者ではなく、率直に何も運のないものをコンパイルしようとする生産的な時間を無駄にしていました） 。 gdal 1.9はfgdbをサポートしているように見えるので、開発者バージョンをダウンロードできると思った。ただし、最新バージョンのQGIS（191a229）では、gdal 1.9でもgdbファイルを追加できないようです。 マスタービルドに含まれているgdalバージョンは、fgdbサポートでコンパイルされていないだけですか？新しいバージョンをコンパイルせずにこれを行う方法があるかどうか誰もが知っていますか？

11 qgis  gdal  file-geodatabase 

2つのラスターを追加し、1つのラスター出力を生成するPythonメソッドを書いています。私が制御できない理由により、入力ラスタの範囲は異なりますが、重複しています。 出力ラスタ範囲が2つの入力の交差領域のみになるように出力を生成するために、2つの重複領域で重複する入力ラスタピクセルのみを操作することは可能ですか？

11 python  raster  gdal 

GeoTIFFドライバーを使用して、GDALでいくつかの浮動小数点ラスターを作成しています。結果の画像をQGISまたはArcに読み込むと、デフォルトのシンボルでは最小値が-3.40282e38、最大値が3.40282e38になるため、ラスターは正確に灰色に見えます。実際の範囲をGeoTIFFに直接書き込む方法はありますか？GISプログラムにロードすると、ヒストグラムが自動的にうまくスケーリングされますか？私はこのようなデフォルトのヒストグラムを作成しようとしました： rasterMin, rasterMax = raster.GetRasterBand(1).ComputeRasterMinMax() raster.GetRasterBand(1).SetDefaultHistogram(rasterMin, rasterMax, 255) しかし、バインドの問題のようなエラーが表示されます。 File "/usr/lib/python2.6/dist-packages/osgeo/gdal.py", line 846, in SetDefaultHistogram return _gdal.Band_SetDefaultHistogram(self, *args) TypeError: not a sequence 私が何ができるか、または間違えた場合の提案はありますか？

11 python  raster  gdal  typeerror 

[編集：私の元の投稿では、タイプミスに基づくエラーを受け取っていたので、質問のその部分を廃棄しましたが、まだ当てはまる質問を残しました。] でogr2ogrドキュメントは言います： -nlt type: Define the geometry type for the created layer. One of NONE, GEOMETRY, POINT, LINESTRING, POLYGON, GEOMETRYCOLLECTION, MULTIPOINT, MULTIPOLYGON or MULTILINESTRING. Add "25D" to the name to get 2.5D versions. このパート：Add "25D" to the name to get 2.5D versionspostgresへのデータの読み込みと、postgresから.shpファイルへのデータのエクスポートの両方に適用されますか？ さらに、25Dは、xy座標の各ペアに対応するz値を持つことができると想定していますが（PolygonZシェープファイルタイプの場合）、これらの座標は重複しない場合があります。これは正しいです？この場合、2.5Dと3Dの違いは何ですか？ ありがとう

11 postgis  gdal  3d 

私は質問がかなり曖昧であることを知っていますが、w / meを負担してください。私は、Google / Bingマップタイルを作成するために使用したさまざまな方法論について、人々がどのような製品パフォーマンス（特にタイミング）を見てきたかを把握しようとしています。これを実行する方法は多数あります（gdal2tiles、FME、maptilerなど）。まともなLinuxサーバー上で、単に大きなPNGを取得し、imagemagickを使用してタイルを作成するという最初の試みでは、かなり長い処理時間が生じました。新しいタイルは少なくとも毎日生成する必要があるため、これに要する時間は非常に重要です。 唯一の実際の要件は、Linuxサーバーで実行できることです。明らかに、無料の方が良いのですが、私はそれだけに制限したくありません。入力には、生のグリッド/ラスターデータまたは大きな画像を使用できます。出力は、GoogleまたはBingマップでそのまま使用できる画像タイルである必要があります。 比較のためだけに、タイミングはGoogleマップのズームレベル7に合わせるべきだと言います。 皆様のご協力に感謝します。この質問がいかに曖昧であるかをおaび申し上げます。 更新：入力に関しては、現在、NetCDF、GRIB、GRIB2などのさまざまな形式の複数の（生の）データソースがあります。生データ自体に加えて、そのデータの非常に大きな画像を生成し、その画像をスライス/タイル化することもできます。 理想的には、画像を切り刻むだけですが、最速の結果が得られるものであれば何でも試してみたいと思います。

11 google-maps  gdal  tiles  bing-maps  optimization 

ジオサーバーストアImageMosaic、WorldImage、およびImagePyramidの違いは何ですか？ImagePyramidが最新のGeoserverに含まれていないことに気付きました。これらの1つは他のものよりも効率的ですか？

11 raster  geoserver  gdal  tiles 

4コアのdebianマシンでgdal2tilesを実行します。 sudo /usr/bin/gdal2tiles.py -r cubic -a 0,0,0 -z 10-15 /home/adm/topo/ekb.vrt /var/www/tiles ではhtopIを参照してください。 gdalが作業でコアを1つだけ使用するのはなぜですか？Gdalバージョンは2.1.2です。

11 gdal 

特定のシェープファイルを使用して、NetCDFからデータをサブセット化する必要があります。データは、1/4度の解像度での海面温度と海の色です。米国を表す4つのポリゴンがあります。北東大陸棚の大規模な海洋生態系とそのデータを抽出するために使用する必要があるサブコンポーネント。私は1982-2014の月次複合ファイルで作業しているので、このデータ抽出ルーチンは自動化する必要があります。ファイルはすでに[35、45、-80、-60]のおおよその作業領域グリッドにサブセット化されています。 以前は、RでHDF5データファイルをラスターに変換してこの方法で処理していましたが、この方法は非常に効率が悪く、現在のNetCDFファイルを使用するpythonでより良い解決策があると確信しています。 これまでは、GDALとフィオナを使用して形状ファイルを読み込み、NetCDF4を使用してデータファイルをロードしてきました。データをサブセット化する方法がわかりません。私はこれを見つけました： GDAL for Python：NetCDFファイルからサブドメインを抽出しますか？ しかし、私はこれらのポリゴンが確かにそうではない単純なバウンディングボックス以外の何かを使用してnetcdfファイルをサブセット化する方法についての最も曖昧な考えを持っていません。 ポリゴンルーチンのポイントはおそらく機能するのに永遠を要しますが、これらの形状に合うように回転された小さな境界ボックスを使用してデータをサブセット化し、最初の開始点としてポリイン検索を行うことができます。 lon / lat境界ボックスを使用して、曲線のnetCDFファイル（ROMSモデル出力）をサブセット化します。 何か案は？ 編集1： OpenClimateGISパッケージに出くわしました。これは法案に完全に適合しているようです... http://ncpp.github.io/ocgis/examples。 html＃advanced-subsetting

11 python  shapefile  gdal  polygon  netcdf 

私はNumPyで処理するためにラスターレイヤーを配列に変換するツールを開発しています。理想的には、すべてをエクスポートせずに.gdbにパッケージ化されたラスターで作業できるようにしたいと思います（特に、ArcGISまたはArcPyを使用する必要がある場合） ）。 これはOpenFileGDBドライバーで可能ですか？ 私が言うことができることから、このドライバーはラスターレイヤーをベクターレイヤーと同じように扱うようです。これにより、レイヤーに関するいくつかのデータにアクセスできますが、ReadAsArray機能は利用できません。

11 python  gdal  file-geodatabase 

1バンドのGeoTIFFファイルがあります。バンドの色の解釈は灰色です。私は私の変換するために、カラーテーブルを必要Col.tiffとしますCol.jpegかCol.png。 私は.vrtを介してこのtiffにカラーテーブルを入れてみました： gdal_translate -of VRT Col.tiff Col.vrt 次に、.vrtを追加します。 <ColorInterp>Palette</ColorInterp> <ColorTable> <Entry c1="0" c2="0" c3="255" c4="240"/> <Entry c1="0" c2="0" c3="254" c4="240"/> <Entry c1="0" c2="0" c3="253" c4="240"/> <Entry c1="0" c2="0" c3="252" c4="240"/> </ColorTable> そしてもう一度翻訳してください： gdal_translate Col.vrt Col.tiff しかし、それは答えました： Warning 1: Unable to export color table to GeoTIFF file. Color tables can only …

11 gdal  gdal-translate  colortable 

私はgdalwarpを使用して、静止からランバート正角に再マッピングしようとしています。入力データはnetcdfにあり、地理座標（度）にあります。再マップしたデータをnetcdfに出力したいと思います。入力netcdfデータに対応するvrtファイルを作成しました。Gdalwarpはnetcdfファイルを出力しますが、出力データはすべてゼロであり、次のエラーを受け取ります。 Creating output file that is 5120P x 5120L. Processing input file netcdf.vrt. ERROR 1: Too many points (441 out of 441) failed to transform, unable to compute output bounds. Warning 1: Unable to compute source region for output window 0,0,5120,5120, skipping. 0...10...20...30...40...50...60...70...80...90...100 - done. 私は次のコマンドを試みました： /usr/bin/gdalwarp -s_srs "+proj=geos +h=35785831 …

11 gdal  coordinate-system  proj  gdalwarp  netcdf 

QGIS 3をインストールし、「ふるい」「マージ」「ポリゴン化」などのツールを実行してみました。それらすべては私にエラーメッセージを与えました： GDAL command output: /bin/sh: (scriptname).py: command not found との両方/Username/miniconda3/binでスクリプトが見つかりました/Username/miniconda3/lib/python3.6/site-packages/GDAL-2.2.2-py3.6-macosx-10.9-x86_64.egg-info/scripts。 ターミナルでenvを確認miniconda3/binすると、パスにあることがわかります。 PATH=/usr/local/bin:/usr/bin:/bin:/usr/sbin:/sbin://Users/Username/miniconda3/bin これらのスクリプトの実行が失敗する原因は何ですか？

11 qgis  gdal  osx 

Rの同じGeoPackageに複数のレイヤーを書き込もうとしていますが、エラーが発生しますCreation of output file failed。RGDALを使用した.gpkgファイルの読み取りと書き込みのドキュメントを検索してみました。たとえば、writeOGR()実際に複数のレイヤーをサポートしているかどうかを調べたところ、ほとんど成功しませんでした。これは可能であれば、どのように行うのですか？最小限の作業例： library(sp) library(maptools) library(rgdal) data(wrld_simpl) norway <- wrld_simpl[wrld_simpl$NAME == "Norway", ] sweden <- wrld_simpl[wrld_simpl$NAME == "Sweden", ] file <- tempfile("scandinavia", fileext = c(".gpkg")) writeOGR(norway, dsn = file, layer = "norway", driver = "GPKG") writeOGR(sweden, dsn = file, layer = "sweden", driver = "GPKG") ogrListLayers(file) トリックを実行するogr2ogr シェルコマンド（ヒントmdsumner）があるようですが、これをR関数でラップできます。ただし、sfパッケージのwriteOGR（）やst_write（）にこれが組み込まれているとlayer_options便利です。GDALに依存していると思いますが、GDALの …

11 gdal  r  rgdal  geopackage  sf