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

QGISでラスター（ピクセル値）を変更したい。 プラグインの「値ツール」は、ラスターのピクセルの値を簡単に提供できます。ただし、変更方法はわかりません。 「ラスターDEMを手動で修正する（OpenSource）」というタイトルのQ＆Aを見つけましたか？それはGRASSでそれを行うことが可能であると言います。私は、草が成功せずにどのように機能するかを理解しようとして、数時間を費やしました。 sextanteプラグインにd.rast.editツールが見つかりません。 それにアクセスする別の方法はありますか？

17 qgis  raster 

グリッド、TIF、およびIMGは、ラスタファイルの3つの形式です。 ArcGIS Desktopでは、明らかに違いはありません。 これら3つの形式の違いは何ですか？

17 raster  file-formats 

ポリゴンを含むシェープファイルがあります。また、グローバルラスターファイルがあります。シェープファイルのポリゴンをラスタグリッドにオーバーレイし、各ポリゴンの平均ラスタ値を計算します。 GDALを使用して結果をシェープファイルに書き込むにはどうすればよいですか？

17 python  raster  shapefile  gdal  spatial-statistics 

そのため、GDALは最近、S3バケットファイルのランダム読み取りを可能にする新しい機能を追加しました。ファイル全体をダウンロードせずに、画像の複数のタイルからGDAL画像をトリミングしようとしています。GDALを使用してS3バケットを設定およびアクセスする方法に関する非常にまばらなドキュメントしか見ていませんが、開始方法について少し混乱していますか？この目標を達成するために、GDALの仮想ファイルシステムをどのように設定するかについて、非常に短い例/チュートリアルを提供してくれる人がいますか？ソリューションでPythonを使用したスクリプト作成が許可されている場合、ボーナスポイント！ 明確にするために、すでにPythonでそれを行っています。Pythonの問題は、画像全体をダウンロードして操作する必要があることです。GDALの最新バージョンでは、S3バケットのマウントがサポートされているため、画像のごく一部をトリミングする必要がある場合、その小さな部分を直接操作できます。残念ながら、この機能は1月に安定版ブランチでのみリリースされたため、ドキュメントは見つかりませんでした。そのため、ソリューションでは、GDALの最新リリースのVSI3システムを使用するか、システムをスマートに使用して、ユーザーがイメージ全体をEBSドライブにダウンロードして操作する必要がないようにします。 つまり、GDALの最新バージョンにあるVSI APIを使用する回答に対して賞金が授与されるため、ファイル全体をメモリまたはディスクに読み込む必要はありません。また、使用するバケットは常に公開されているわけではないため、投稿されるHTTPトリックの多くは、多くの状況で機能しません。

17 raster  gdal  gdal-merge  amazon-ec2  amazon-s3 

これは約13 haの湖で、約10の横断線で81のサンプリングされた深度ポイントがあります。 以前は、50 haの湖と約100の深さのサンプリングポイントがあるARCgisで、このTopotoRasterツールを使用して適切な出力を得ました。ただし、QGISまたはオープンソースソフトウェアには、そのような正確な同等物は存在しないようです。 QGISでは、Raster InterpolationプラグインのTINメソッドを使用して、下の補間されたラスターを取得しました。しかし、これは幸運な結果のように思えました。繰り返してみたところ、実際にはパラメーターを変えていなかったにもかかわらず、異なる結果（不満足な結果）が得られました。 これらは湖よりも都市型の（構築された）タンクであるため、ベッドのプロファイルがかなり規則的であるため、サンプリング密度が低くなります。一部の記事では、逆距離重み付け（IDW）が最適な手法として提案されているように見えますが、それで（おそらく点密度が低いために）最悪の結果が得られるようです。 湖のサイズ、その通常のプロファイル、およびどの補間方法が適切かを判断するためのポイントサンプリング密度を考慮して、ここで使用できるヒューリスティックはありますか？（クリギング、IDW、バイリニア、キュービックコンボリューション、TIN、またはスプライン間） または、それは常に少しの試行錯誤ですか？

16 qgis  raster  interpolation  bathymetry  lakes 

のレイヤーがありcoord. ref. : NAます。 私はどのように割り当てることができますcoord. ref. : +proj=longlat +datum=WGS84 +no_defs +ellps=WGS84 +towgs84=0,0,0するにはcoord. ref. : NA？ class : RasterLayer dimensions : 127, 146, 18542 (nrow, ncol, ncell) resolution : 0.01, 0.01 (x, y) extent : -122.39, -120.93, 36.72, 37.99 (xmin, xmax, ymin, ymax) coord. ref. : NA data source : C:\2012Tif\2012ASC5min_ppt_spas1306_0001_20121222_0005_UTC.asc.tif …

16 raster  coordinate-system  r 

私はこれらの列を持つpostgresql / postgisにラスターデータベースを持っています： （ID、rast、data_of_data）。 「ラスト」は、WKT形式のラスターファイルがある列です。WGS84システム（30.424、-1.66）および2002-01-09のポイントのDN値を検索するクエリの例は次のとおりです。 SELECT st_value(rast,(st_GeomFromText('POINT(30.424 -1.66)', 4326))) as val FROM my_table WHERE date_of_data='2002-01-09' これらの種類のクエリを高速化する方法（空間インデックスなど）はありますか？

16 postgis  raster  query  performance  spatial-index 

2つのレイヤーがあります。多くのタイルを含むポリゴンシェイプレイヤーと、カラーマップの多くのカテゴリを含むCORINE 2006土地被覆を含むラスターレイヤー。シェープレイヤーの各ポリゴンについて、ラスターレイヤーの各土地被覆カテゴリの合計を取得します。 たとえば、IDが '2'のポリゴンがあり、このポリゴンに対して次のような属性が必要です（パーセントまたは平方メートル）。 耕地：15％ フォレスト：11％ 通り：2％（...など） 私は草、qgis（機能なし）、saga（すべてを合計値に合計する）r（合計）でそれをしようとしましたが、まだ解決策が見つかりませんでした。ほとんどのプラグイン（qgisのゾーン統計）は、0-1ラスターレイヤーのみをサポートします。v.rast.statsも役に立たなかった。Iamは、優れたスマートソリューションにオープンです。たぶん私は間違ったアプローチを使用したか、間違いを犯した。 Arcgisでは、このタスクは非常に簡単です。覚えていれば正しいですが、私はあなたの日常のLinuxユーザーにとって良い解決策がまだありません。 私はDebian Linuxシステムを実行しているため、このOS用のプログラムしか使用できません。 編集：この質問にはまだ多くのビューと訪問者がいるので、QGISプラグインを作成しました。これは、ラスターレイヤーの土地被覆を計算することもできます。ポリゴンオーバーレイはまだコーディングしていませんが、間違いなく計画されています。ここでプラグインを見つけて、最初にScipyライブラリをインストールします。

16 qgis  raster  polygon  r  statistics 

作業中の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 

DEM（デジタル標高モデル）とライン（x1、y1）->（x2、y2）（ペアの座標（x、y））が与えられた場合、そのラインに投影された標高プロファイルまたは断面をどのように計算できますか？ 私は自分のプロジェクトで使用できるオープンソースコード、またはガイドとして使用できる擬似コードを探しています。

16 raster  dem  open-source-gis  elevation  profile 

QGISで最も近い点までの距離を計算する方法はありますか？ArcGISには、ユークリッド距離と呼ばれる関数があり、入力ベクトルまたはラスターを取得し、距離のラスターを生成できます。 QGISに同様の機能はありますか？ドキュメント内で見つけることができませんでした。RasterQGIS の関数とMMGISHub Distanceからの計算（非ラスターソリューションの場合）を調査しましたが、どちらを使用してこれを行うかはわかりませんでした。ベクトルアプローチが唯一の解決策である場合、それでも問題はありません。

16 qgis  raster  distance  point 

私は本当にPythonが初めてで、ArcGIS 10のPythonを使用して、ピクセルごとのラスターのセル値と座標（各ピクセルの中心のXY座標をマップ）を取得する簡単な方法があるかどうかを知りたいですか？ これをさらに説明するには、最初のピクセルのマップX、マップY、およびセル値を取得し、それらの3つの値を3つの変数に割り当て、他のピクセルの残りに対してこの手順を繰り返す必要があります（ラスター全体をループします）。 私の質問をもっと説明する必要があると思います。問題は、最初のラスタのピクセルのXY位置を取得し、そのXY位置に対応するいくつかの他のラスタのセル値を取得する必要があることです。このプロセスは、中間ポイントシェープファイルを作成せずに、最初のラスターのすべてのピクセルをループ処理する必要があります。80億ピクセル近くのラスターを処理する必要があるため、本当に時間がかかります。また、ArcGIS 10でPythonを使用してこれを行う必要があります。 @JamesS：ご提案ありがとうございます。はい、これは1つのラスターで機能しますが、他のいくつかのラスターのセル値も収集する必要があります。問題は、最初のラスタの最初のピクセルのX座標とY座標を取得した後、最初のラスタのX、Y位置、3番目のラスタなどに対応する2番目のラスタのセル値を取得する必要があることです。だから、最初のラスターをループするとき、ピクセルのXとYの位置を取得し、その位置に対応する他のラスターのセル値を取得することは同時に行われるべきだと思いますが、わかりません。これは、最初のラスターをポイントシェープファイルに変換し、ArcGIS 10でポイント値への複数値の抽出を実行することで実行できますが、 @hmfly：ありがとう、はい、このメソッド（RastertoNumpyarray）は、配列の既知の行と列の値の座標を取得できれば機能します。 @whuber：計算を実行したくないので、テキストファイルにXY座標とセル値を書き込むだけで十分です。

16 arcgis-10.0  python  raster 

QGISのラスター内のセル値のカウントの概要を取得するにはどうすればよいですか？レイヤーのプロパティに要約統計とヒストグラムを提供するために、QGISがこれを行っていることを知っています。要約統計もヒストグラム画像も必要ではなく、個々の値のカウントが必要です。 ある時点で、このための特別な機能を追加する試みがあったようです（機能＃4321、http：//hub.qgis.org/issues/4321）。 私は本当に明白な何かを見逃していますか？

16 qgis  raster 

どのCRSあなたは歴史的な世界をジオリファレンスする推薦でのもののようにマッピングしピープルズアトラス- - 1920年ロンドン国土地理院、例えば マップには、コントロールポイントを選択するために使用する緯度/経度グリッドが表示されます。私はWGS84 EPSG：4326を使用して、いくつかの合理的な結果を得ることができましたが、より良いオプションがないのではないかと思います。

15 raster  georeferencing 

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

15 raster  gdal  band