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

ESRI ASCIIグリッドラスター形式のデジタル地形モデルファイルを実験しています。このタイプのデータに関する経験はゼロですが、以下のコードのように、Rにロードするのは非常に簡単であることがわかりました。パッケージのplot3D関数はrasterVisデフォルトで見栄えの良いプロットになりますが、モデルはかなり反映されています。それは風景であるはずなので、この光沢のある表面は本当に適切ではありません。 私はもっ​​と自然なもの（そのような単語が一連の図の3Dレンダリングについて使用できる場合）、マット、トポロジープロットを探しています。 これを回避する方法はあると思いますが、rasterいくつかのポインタを使用したり探したりした経験はありません。 プロットの反射率を低減または削除するにはどうすればよいですか？ library(raster) library(rgdal) library(rasterVis) foo <- raster(readGDAL("my.dtm.asc")) plot3D(foo)

10 raster  r  dem  visualisation  color 

ブレークラインは、たとえば正確なTINを構築するために、地形を記述するのに非常に役立ちます。 私の質問は、これらのブレークラインを特定のDEMから抽出する方法ですか？ 曲率（垂直）-ラスター（arcgis 10を使用）を計算しました。ブレークライン（最小/最大曲率）を確認できますが、それらを抽出して最終的にポリラインを取得する方法はわかりませんでした。 これらのブレークラインをDEMから取得するためのワークフローは何ですか？ 誰もが役立つワークフロー、または「ツール」、またはskriptのアイデアを持っていますか？ 私はこのタスクを実行したい最初/最後の人ではないと思います:) ...その間、私はこのトピックに関する論文を見つけましたが、それは良い解決策のようですが、大きな問題が残っています：それを実践に減らす方法は？ または、誰かがいくつかのソフトウェアで異なる（より簡単な？）アプローチを持っていますか？

10 raster  arcgis-desktop  dem  breaklines 

多数のラスターファイルを含むプロジェクトがありますが、非常に遅くなります。これをすべて処理するための最良の方法はどれでしょうか。私は単にこれらのファイルをジオリファレンスして.wldファイルを作成し（geoTiffの方が良いですか？）、ラスターとして開きます。私は世間知らずですか？

10 qgis  raster 

以下の例に示すように、地質図を利用したいと思います。これは、地図のかなり小さくて縮小された部分にすぎないため、編集/再描画が最終的に意味をなさないようです。 私の考えは、最初に岩盤と構造の断層のみを残してベクターに変換し、その後オプションで「クリーンな」ラスターファイルを使用できるようにすることです。 「方法」に関するヒントはありますか？

10 qgis  raster  vectorization 

都市化されたエリアを表すラスターデータがあります。これをトレースして、ポリゴンに変換します。 ArcMapのRaster to Polygonツールを使用して、それに近いものを取得できます。 単純化： 簡略化： ただし、コーナーを共有するポリゴンのチェーンを作成する多くの対角線（道路）があることがわかります。 これにより、ポリゴン数が大幅に増加します。これらのラスターデータセットは州レベルで処理されています。これらの線を滑らかにしてポリゴン数を減らすためのArcまたはGeoToolsツールはありますか？結果の精度がいくらか失われても問題ありません。最後に、単純化されたポリゴンで得られるものに似たものを探しています。隣接するポリゴンが結合されているだけです。

10 arcgis-desktop  arcgis-10.0  raster 

ArcGISを使用せずにラスターから標高データを取得した経験があり、情報をpython listまたはdict？ タプルのリストとしてXYデータを取得します。 xy34 =[perp_obj[j].CalcPnts(float(i.dist), orientation) for j in range (len(perp_obj))] リストをループするか、それを関数またはクラスメソッドに渡して、xyペアの対応する高度を取得します。 私はこのトピックについていくつかの調査を行い、gdal APIは有望に聞こえます。誰もが私に物事、落とし穴、サンプルコードに取り組む方法をアドバイスできますか？ 作業中のマシンでシステムパス変数を編集できないため、GDALはオプションではありません。 誰かが別のアプローチについて知っていますか？

10 python  raster  elevation 

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

10 raster  gdal  osgeo  photogrammetry  unmanned-aerial-vehicle 

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

10 python  raster  gdal  imagery 

{sp}パッケージの 'rasterize'関数を使用して、世界の水深ラスターレイヤーから関心領域の水深値を抽出しようとしています。 *編集：私が探しているもののように見える「抽出」機能を見つけました。 これは私がこれまでに行ったことです： > class(subarea0) #This is my area of interest (Eastern Canadian Arctic Sea) [1] "SpatialPolygons" attr(,"package") [1] "sp" > extent(subarea0) class : Extent xmin : -82.21997 xmax : -57.21667 ymin : 60.2 ymax : 78.16666 library(marelac) data("Bathymetry")#World bathymetric data in library (marelac) names(Bathymetry);class(Bathymetry);str(Bathymetry) [1] "x" "y" "z" …

10 raster  r  overlay 

Spatial Analystエクステンション付きのArcGIS Desktop 10を使用しています。 複数のラスターを1つに結合し、重複するセルの値から常にランダムに選択するにはどうすればよいですか？ これをよりよく説明する画像があります：

10 arcgis-desktop  arcgis-10.0  raster  spatial-analyst  map-algebra 

.ascファイルがあり、QGISにインポートしたい。 それを開く標準的な方法は、「ラスターレイヤーの追加」（「ラスターレイヤーのインポート」）からですが、別の方法はありますか？ .ascファイルを.tifまたはQGISで別のラスター形式に変換する必要があります。

10 qgis  raster  geotiff-tiff  ascii 

Rでラスターデータを操作したいのですが、ArcMapからラスターデータセットをエクスポートできません。「データのエクスポート」/ GRIDを使用して、*。aux拡張子の付いたファイルを取得します。必要なのは* .grdです。 誰かが私を手伝ってくれる？

10 arcmap  raster  r 

たとえば、ArcMap 10ドキュメント（関連付けられたtfwを含む.tiff）に読み込まれた長方形の地理参照ラスターファイルがある場合、その中心点を簡単に見つけてその点をポイントベクターレイヤーに保存するにはどうすればよいですか？ また、ArcMapドキュメントにそのようなラスターが複数ある場合、それらすべてにプロセスを適用するにはどうすればよいですか？ 残念ながら、Pythonの経験はありません。したがって、プログラムによるソリューションは問題ありませんが、既存のスクリプトをArcGIS 10にロードして、問題のラスターで実行する方法についての具体的な指示が必要になります。（ところで、ラスターはすべて個別のレイヤーにあります）。

10 arcgis-desktop  raster  arcgis-10.0  point  geotiff-tiff 

カテゴリカルデータをグリッド化する方法を探しています。海図とフィールドシートから、海底の性質を特定するポイントのコレクションを抽出しました。これらのデータは、数値ではなくカテゴリであり、定期的またはランダムにサンプリングされることもありません。航海図は、ナビゲーションとアンカーを支援するために作成されます。それらは生息地をマッピングするために作成されていません。そのため、比較的浅い水深が航行に危険をもたらす可能性があり、船が停泊する傾向がある海岸近くで、より多くの探査が行われます。岸からさらに遠く、航行に十分な深さがあり、アンカーが実用的でない場合、探査はそれほど頻繁に行われません。 他の誰かが海図から格子状の基質マップを作成しようとしましたか？ 私はティーセン（ボルノイ）ポリゴンを調べましたが、海岸に沿った音響の集中は、海岸沿いの細かい「ハニカム」、オフショア、オフショアに伸びる長いパイ形ポリゴンの間にある大きなポリゴンにつながります。最近傍を使用したグリッドは、ほぼ同じ結果になります。 岸に近い浅いポイントの影響を制限する方法が必要です。長いパイ形のポリゴンを制限する方法です。深い海では、海底の性質が岸に近い底部の続きになるとは思わない。私は2つの線に沿って考え始めました-両方とも深さを使用しています。1つは、グリッドセルと隣接する点の間の深さの違いを使用して、「最も近い」隣接の選択に重みを付けることです。もう1つは、特定の許容差よりも深さが異なる隣接ポイントの選択を解除することです。または、事前に指定された許容値ではなく、深度範囲をビン化して、隣接するポイントの選択を同じ深度範囲またはビン内のポイントに制限することもできます。 これらの2つのオプションのいずれかを実装する方法についての考えはありますか？ 他のフォーラムで同僚と話し合って以来、私は他にもいくつかのアプローチを検討してきました。1つ目は、岸壁のデータの影響を制限するために、バリア（深度100mのコンター）を使用することです。このアプローチの課題は、バリアを使用できるすべてのESRI補間ルーチンが、不連続データではなく連続データを処理するように設計されていることです。ティーセンポリゴンを作成する前に、バリアを使用してポイントを岸の近くの浅いポイントとより深いポイントに分割できます。ただし、ArcGISは複雑な領域ではなく長方形の領域に対してティエセンポリゴンを作成するため、エッジ効果が蔓延すると予想しています。 2番目のアプローチ-複数の同僚によって提案された-はクリギングでした。私は、継続的なデータのためにそれを検討したことがないので、最初は手に負えないクリギングを却下しました。クリギングの課題は、クリギングもカテゴリカルデータ用に設計されていないことです。今、私は表面の深さと性質を使ったcokrigingを見ていますが、どのようなタイプのクリギングでも、表面の性質に整数の数値コードを使用する必要があります。その後、結果の浮動小数点数値コードは、元の整数コーディングに戻す必要があります。きれいじゃない。 他の方針に従うことを提案できますか？（おそらく、地形分析を使用できます。たとえば、安息角より急な斜面は堆積物にはなりません。もっと単純なものを探しています。とにかく、十分な空間解像度のデータがありません。） よろしく、

10 raster  interpolation  kriging  voronoi-thiessen 

Rの「ラスター」パッケージを使用して、Windowsで.DEM形式のラスターファイルを読み取ろうとしています。 Windows 7のRにデータをロードするときにNA値に関する問題が発生しますが、OSX Lionを搭載したMacでは問題がありません。Windowsでは、NA値が正しく読み取られないようです。問題は、なぜこれが起こるのかということです。 使用されたラスターファイルは、次のRコードを使用してUSGSからダウンロードされました。 download.file('http://edcftp.cr.usgs.gov/pub/data/gtopo30/global/e020n90.tar.gz', 'e020n90.tar.gz') untar('e020n90.tar.gz') 次に、「ラスター」パッケージを使用してラスターをRに読み込みます。OSX LionおよびR64バージョン2.13.1では、NA値が認識されます。 > onMac <- raster('E020N90.DEM') > onMac class : RasterLayer dimensions : 6000, 4800, 28800000 (nrow, ncol, ncell) resolution : 0.008333333, 0.008333333 (x, y) extent : 20, 60, 40, 90 (xmin, xmax, ymin, ymax) coord. ref. : +proj=longlat +ellps=WGS84 +towgs84=0,0,0,0,0,0,0 +no_defs …

10 raster  r  windows  osx