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

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

15 python  raster  gdal 

Landsat TM画像から派生した土地被覆分類があります-ENVIを使用してこの分類を作成しました。航空写真と比較して出力の精度を向上させるために、分類のいくつかの領域でセルの個々のグループとグループを編集したいと思います。 ArcGIS Desktopを使用して、セルの大きなグループを選択し、それらを新規または既存のクラスに割り当てることにより、手動で編集する方法はありますか？ クラス統計に基づいてメソッドをマージすることは知っていますが、分類している領域は薄い雲に覆われているためスペクトルの混乱が非常に大きくなっています。

15 arcgis-desktop  raster  arcgis-10.0  editing 

PostGIS2.0を使用して、ラスター/ポリゴンの交差を行います。どの操作を使用する必要があるのか​​、これを実行する最も速い方法は何かを理解するのが困難です。私の問題は次のとおりです。 ポリゴンとラスターがあります ポリゴン内にあるすべてのピクセルを見つけて、ピクセル値の合計を取得したい （更新された問題）：クエリを実行すると、元のラスターに存在しないいくつかのピクセルの大きな値を取得しています ST_Intersects()またはを使用する必要があるかどうかがわかりませんST_Intersection()。また、ピクセルを合計するための最良のアプローチが何であるかわかりません。これが私が試した最初のアプローチです（＃1）： SELECT r.rast FROM raster as r, polygon as p WHERE ST_Intersects(r.rast, p.geom) これはrast値のリストを返しますが、これをどうすればいいのかわかりません。を使用して要約統計を計算しようとしましたST_SummaryStats()が、これがポリゴン内にあるすべてのピクセルの加重和であるかどうかはわかりません。 SELECT (result).count, (result).sum FROM ( SELECT ST_SummaryStats(r.rast) As result FROM raster As r, polygon As p WHERE ST_Intersects(r.rast, p.geom) ) As tmp 私が試した他のアプローチ（＃2）はST_Intersection()次を使用します： SELECT (gv).geom, (gv).val FROM ( SELECT ST_Intersection(r.rast, p.geom) AS …

15 raster  polygon  sql  spatial-database  postgis-2.0 

R内でシェープファイルをラスターに変換しようとしています。私のアプローチは、ラスターを次のように読み込むことです。 library(rgdal) # Loads SP package by default demo <- readOGR('F:/data/', 'shapefile') # Creates a SpatialPolygonsDataFrame class (sp) これはうまく機能し、プロットできます。しかし、それは大きなシェープファイルであり、ラスターに変換したいです。私は次を試しました： r <- raster(ncol=180, nrow=180) Demo_ras = rasterize(r, demo, 'pop') # pop is an integer here エラーメッセージ： Error in function (classes, fdef, mtable) : unable to find an inherited method for function …

15 raster  shapefile  r 

これは私が使用しているラスターのきれいな写真です-エフェクトは適用されていません-カラーランプとヒストグラムを調整して、すべてのレリーフを表示します-非常に2Dですが！ 陰影起伏をデフォルトとして表示したいラスターレイヤーがあります。これはコードで行うことができ、実際には、ArcMapで陰影起伏オプションを選択するのと同じです（以下を参照）。 ただし、これを行うと、期待はずれに見えます。 特に、ラスターにズームインすると、あらゆる種類のグリッドアーティファクトが表示されます。 ただし、Spatial Analystのラスターを使用して適切な陰影起伏レイヤーを作成し、85％の透明度でオーバーレイすると、かなり見栄えがよく、非常に異なって見えます！： ArcMapでこの効果を実現し、Hillshadeオプションのデフォルトの機能を無効にするにはどうすればよいでしょうか（これは不可能か、ラスターが正しく表示されないために見落とした設定があると思います）。正しい設定を上書きまたは選択できる場合-ArcObjectsで実行できますか。 オプション1が不可能な場合、Hillshadeレイヤーを作成して（Spatial Analystを使用せずに）エフェクトを作成し、それを元のラスターに効果的に融合できます。つまり、ArcMapのコンテンツウィンドウで2つのレイヤーを表示したくない。1つのレイヤーだけが必要ですが、選択/構成の既定の陰影起伏レンダリングが必要です。これもプログラムで行う必要があります。つまり、ArcObjectsです。 私が学んだことの1つは、正しい空間参照および測定ユニットにDEMを配置する必要があること、または結果の陰影が黒すぎて誇張に見えることです-詳細については、このビデオをチェックしてください-非常に良い説明です。 だから、高品質の陰影起伏を生成するために何をする必要があるかがわかったと思いますが、ArcObjectsコードを微調整する必要があると思います。自分のグラフィックレンダリングを考え出し、これをArcMapにプラグインする必要があると考えています（グリッドのビットマップをクールな方法でレンダリングするdllがあります）。これを示す良い場所は、IRasterRendererのようなものでDraw（）イベントをオーバーライドし、独自の画像レンダリングを挿入することだと思います。もちろん、再投影などで少し複雑になるかもしれません。 ArcMapでこのような低レベルのグラフィックレンダリングを行った人はいますか？もしそうなら、ガイダンスや知っておくべきことを提供できますか？

15 arcobjects  arcmap  raster  hillshade 

最新のArcGISバージョン（10）には、ファジーロジック用の2つのarctoolboxツール（Fuzzy Membership、およびFuzzy Overlay）が含まれています。誰がどのように適用できるかを示す例を見つけることができますか？ （Desktop 10ヘルプから。）

15 arcgis-10.0  raster  spatial-analyst 

このソースから複数のDEMラスターを1つに結合する必要があります。 http //srtm.csi.cgiar.org/SELECTION/inputCoord.aspなので、Mosaic to New Rasterツールを使用します。 ツールをセットアップしました： ラスターを入力します（すべて同じソース、同じサイズ、投影なし）。 ピクセルタイプを設定します（元のラスターのように32ビット浮動小数点）。 セルサイズを元のラスターとして設定します。 元のラスターとしてのバンド数= 1。モザイクオペレーター、私はBLENDとMEANでそれを行いました（同じ結果が得られます）。 私が抱えている問題は、結果のモザイクが、個々のラスターの最大値と最小値とは異なる最大値と最小値の範囲を示すことです。たとえば、ラスター1（-5123.8、23.25）、ラスター2（-5974.6、40.09）、ラスター3 （-57770.2、38）、ラスター4（-2534.3、23.55）、および最終モザイクラスター（-5975.8、81.1）。 この解決策は正しくないと思います。少なくともそれを期待していませんでした。これが大丈夫かどうか、そしてそうでない場合、誰でもそれを解決し、正しい最大値と最小値を持つ適切なモザイクラスターを取得する方法を知っていますか？ ArcGIS 10.2.2 for Desktopを使用しています。

14 arcgis-desktop  raster  arcgis-10.2  dem  mosaic 

私は地理データの絶対的な初心者なので、質問が適切でない場合はご容赦ください。 NCDC NARRからデータをダウンロードし、rasterパッケージを使用してRにロードすることに成功しました。緯度、経度、値のリストを取得したいのですが。私はそれrasterToPoints()が私が望むことを正確に行うべきであることを理解していますが、私の緯度と経度の値は奇妙に見えます： r <- raster(myfile) data_matrix <- rasterToPoints(r) head(data_matrix) x y value [1,] -5405401 4347242 70 [2,] -5372938 4347242 88 [3,] -5340475 4347242 76 [4,] -5308012 4347242 85 [5,] -5275549 4347242 87 [6,] -5243086 4347242 88 私は、現在Lambert Conformal Conic（LCC）である投影法で何かをすべきだと思います。ラスターに関する詳細情報は次のとおりです。 > r class : RasterLayer dimensions : 277, 349, 96673 …

14 raster  r 

特定の日の旅行時間の等時線を表すポリゴンのテーブルがあります。各原点には、5つのアイソクロネジオメトリがあります（別々の行に格納されます）。各原点について、5つのアイソクロナス（バイナリNULLまたは1）をラスタライズし、それらを1つのラスタレイヤーに結合します。このラスターレイヤーには単純なマップ代数sum / 5が必要です。そのため、各原点は、最終的に[NULL、0.2、0.4、0.6、0.8、1]の値を持つ単一のラスターレイヤーに関連付けられます。構成レイヤーはオーバーラップします。それは確率曲面です。 私のデータはすべてPostgres 9.3（PostGISを使用）に保存されます。私の問題は、PostGISラスタの使用方法を学びたいが、非常に急な学習曲線を持っているように思え、すべての例が単一のラスタレイヤに対処できることです。例では、このレイヤーはポリゴンオーバーレイの一部として使用され、各ポリゴンのラスターの値を平均化しています。以下を組み合わせた複製可能な例は見つかりませんでした。a）ベクトル->ラスターb）マップ代数。およびc）最初の段落ごとのGROUP BY属性。 このタスクを実行するために必要な場合、GDALまたはGRASSを使用しても問題ありませんが、これはPostGISが処理できるもののようです。入力データがすでにPostGISジオメトリである場合、そうするのが便利です。そして、私は本当に PostGISラスターに同意したいと思っています。 サンプルデータ構造： areaid time date isogeom (polygon) 1000 07:15:00 2014-05-05 xxx 1000 07:15:00 2014-05-06 xxy ... 1006 07:15:00 2014-05-05 zzz ラスタライズし、areaidでグループ化し、マップ代数演算を実行して次のことを行います。 areaid isorast (raster) 1000 aaa 1006 bbb 私はこれをPostGISに含めることに成功していません。私のアプローチは、ベクターをラスターに変換し、ラスターを配列にダンプし、psycopg2を介してnumpy配列との組み合わせを実行してから、GeoTIFFに書き込む（PostGISに戻す）ことです。理想的ではありませんが、実行可能です。

14 postgis  raster  map-algebra  postgis-raster 

道路ネットワークを表すラインシェープファイルがあります。このデータをラスター化して、ラスターの結果値がラスターセル内に収まるラインの全長を表示するようにします。 データはBritish National Gridプロジェクションにあるため、単位はメートルになります。 理想的には、を使用してこの操作を実行したいと思いますが、パッケージRのrasterize機能がrasterこれを達成するのに役立つと推測していますが、適用される機能がどうあるべきかわかりません。

14 raster  r  line  rasterization 

私は困惑しており、かなりの数の検索が役に立たなかった。 最終的には3つのラスターレイヤーと3つのマーカーレイヤー（ラスターとその関連データを切り替えるコントロール、データのポップアップ/ツールチップを含む）を含むマップの最初のステップにいますが、できません。 Tilemillタイルを取得して適切にレンダリングします。（ラスターをQGISで作成し、GDALを使用してそれらを色付けして再投影し、TileMillに配置しました。） 要するに、 それはChromeにあります。これはSafariです。 そして、これは私にはまったく無害に見えるコードです（この時点でタイルを除くすべてをコメントアウトしたので、これですべてが含まれました）。 <html> <head> <title> - Leaflet demo</title> <script type='text/javascript' src="http://cdn.leafletjs.com/leaflet-0.7.3/leaflet.js"></script> <style type='text/css'> html, body{ width:100%; height: 100%, margin:0; padding:0; } #map{ width:700px; height:600px } img {margin:0px;padding:0px}; </style> <script type='text/javascript'>//<![CDATA[ window.onload=function(){ //one base layer on map var map = new L.Map( "map" ).setView([38.173,-94.384],5); L.tileLayer('http://a.tiles.mapbox.com/v3/mtgmediagroup.beef2013/{z}/{x}/{y}.png', { attribution: …

14 raster  javascript  leaflet  tilemill  browser 

整数値を持つラスターイメージ（GeoTIFF）があります。カラーマップなどを使用してQGISでラスター値を視覚化する方法は知っていますが、代わりに各カラー値にシンボルを割り当てたいと思います-カラーマップと同様に、シンボルを使用します。 例： 質問： これは可能ですか？

14 qgis  raster  symbology  visualisation 

私は1x1メートルのDEMと.lasの元のLiDAR点群を持っています。ベクトルフィーチャ（ポイント、ポリライン）に川の堤防（堤防の最高点）を抽出する必要があります。 アルゴリズムまたは既存のツールのアイデアはありますか？ 最初の画像では堤防が明るい灰色で表示され、その下にある堤防のある領域の点群画像が下にあります。青の川の中心線。 DEMから線形フィーチャを抽出するためのツールはありますか？ ENVIの空間特徴抽出モジュール（7ページ）のようなものが必要だと思いますが、ENVIの部分はありません。

14 raster  dem  lidar  algorithm  feature-extraction 

私はGISデータの使用に非常に慣れており、Rをあまり経験していません。spatial-analyst.netPDFブックを使用して空間データを分析する方法について読んでいるので、完全に失われたわけではありませんが、私の問題と人々はアイデアを提案するかもしれません。 さまざまな緯度/経度座標で約2000の測定値を含むデータセットがありますが、データが3年間にわたって収集され、条件が時間とともに変化したため、おそらくこのデータセットを細分化します。測定対象の変数を「IP」と呼びましょう。 サンプルデータに対してKrigingまたは他の補間方法を使用して、問題のある領域全体にIPのマップを作成したい。次に、さまざまなIPバケットの土地の量を測定するヒストグラムを作成します。また、各バケット内のサンプル数を示すヒストグラムを作成する必要があります（サンプルは、クリギングがその土地に対して予測するものよりも実際のIPが大きい場合も小さい場合もあります）。 データをSpatialPointsDataFrameに読み込み、クリギング分析を実行する方法に従いますが、その問題を抱えているのは、ヒストグラム分析を行うためにそのデータをグリッドデータフレームに変換する方法です。 ポイントをグリッドに変換するための提案はありますか？

14 raster  r  kriging 

ベクターレイヤーのポリゴンごとにラスターポイント（マルチバンド）の合計、平均などを計算するにはどうすればよいですか。これは「ゾーン統計」と呼ばれると言われました。最初にQGISで試しました。 それを行う方法はありますが、それは非常に遅すぎます（ラスターをベクターに変換し、2番目のベクターレイヤーと交差し、ジオメトリを計算し、数値をエクスポートし、スプレッドシートまたは他のプログラムで統計を計算し、結果を再インポートし、永遠に時間がかかります350.000ラスターポイント）。 私はまた、サガギを使用するヒントを与えられました。これには「ゾーン統計」がありますが、これらはベクターレイヤーのポリゴンではなく、ラスターバンドのカテゴリに基づいています。したがって、これを使用するには、ベクターレイヤーをラスターに変換してから統計情報を計算する必要があります。 これはこれを解決する間違った方法のようです。2つ以上のポリゴンに属しているラスタポイントは、ポリゴン境界と交差しているため、それらを考慮する方法はありません。ポリゴンベースの統計でこれを処理できるはずだと思うので、正しいモジュールがまだ見つからなかったと思います。 Saga-gisには本当に多くのモジュールがあります。このアプリケーションに適したものを教えてください。

14 qgis  raster  statistics  saga