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

ラスターは、通常は画像のような形式で格納される、値の規則的なグリッドで構成されるデータ形式です。

2
arcpy.RasterToNumPyArrayを使用して空間参照を保持しますか？
ArcGIS 10.1を使用していて、2つの既存のラスターに基づいて新しいラスターを作成したい。RasterToNumPyArrayは私が適応したい良い例があります。 import arcpy import numpy myArray = arcpy.RasterToNumPyArray('C:/data/inRaster') myArraySum = myArray.sum(1) myArraySum.shape = (myArray.shape[0],1) myArrayPerc = (myArray * 1.0)/ myArraySum newRaster = arcpy.NumPyArrayToRaster(myArrayPerc) newRaster.save("C:/output/fgdb.gdb/PercentRaster") 問題は、空間参照とセルサイズを削除することです。arcpy.envを実行する必要があると考えましたが、入力ラスターに基づいてそれらを設定するにはどうすればよいですか？私はそれを把握することはできません。 ルークの答えを取り上げると、これは私の暫定的な解決策です。 Lukeのソリューションはどちらも、空間参照、範囲、セルサイズを正しく設定しました。ただし、最初の方法では配列内のデータが正しく伝達されず、出力ラスタにはどこでもnodataが入力されます。彼の2番目の方法はほとんど動作しますが、nodataの大きな領域がある場合は、ブロック状のゼロと255で埋められます。これは、nodataセルをどのように処理したかに関係している可能性があります。私が話していることの画像を含めました。 #Setting the raster properties directly import arcpy import numpy inRaster0='C:/workspace/test0.tif' inRaster1='C:/workspace/test1.tif' outRaster='C:/workspace/test2.tif' dsc=arcpy.Describe(inRaster0) sr=dsc.SpatialReference ext=dsc.Extent ll=arcpy.Point(ext.XMin,ext.YMin) # sorry that i modify calculation …
1
QGISのラスター属性テーブルにアクセスしますか？
GDALを使用してGeoTiffsのラスター属性テーブルを作成するPythonコードがいくつかあります。 ArcGISでは、これらの属性テーブルをうまく見ることができ、（分類として）好きなように動作しますが、QGISでは属性テーブルの概念全体が欠落しているようです。 QGISでRATにアクセスする方法はありますか？
2
Rを使用したラスターのランダムサンプリング？
プロセスの出力がラスタになるように、ラスタをランダムにサンプリングする簡単な方法はありますか？ r-sig-geo リストにある例を使用していsampleRandomますが、rasterパッケージ内の関数も試しました。これらのアプローチはどちらも出力を生成しますが、ラスターに変換する方法がわかりません。「SpatialPointsDataFrameラスター」のいくつかの組み合わせを検索した後、アプローチを見つけることができませんでした。 library(raster) # read in raster rasterSource <- 'landsat.TIF' r <- raster(rasterSource) # convert to spatial points data frame r.spgrd<-as(r,"SpatialPointsDataFrame") # elminate NA values r.spgrd = r.spgrd[!is.na(r.spgrd[[1]]),] # sample points selectedPoints = sample(1:length(r.spgrd[[1]]), 1000) r.sampled = r.spgrd[selectedPoints,] # try to make spgrd into a raster r.test <- raster(r.sampled) …
13 raster  r 
3
緯度/経度ポイントのPostGISラスター値
PostGIS 2.0データベースにtiffラスターがあります。指定した緯度/経度のポイントでラスターの値を取得したい。 私が持っています： postgresql-9.1 ラスターをサポートするpostgis-2.0-svn PostGISにインポートされたラスターのSRIDは3035
13 postgis  raster  point 
4
「ピクセルの長方形グリッド」をラスターまたは画像と呼ぶべきなのはいつですか？
GISの文献を読んで、著者は特定の意味を念頭に置いて2つの用語を使用していると感じることがあります。私は常に用語を交換可能に扱います。私が間違っている？ リンクの最初の行を読んでください。画像サービスがラスター（および画像）へのアクセスを提供し、2つの用語の異なる意味を示唆していることをご覧ください。
3
土地適合性のモデリングにはどのようなプラクティスが利用できますか？
現在、研究領域をセルに分割する「クラシック」なラスターベースのアプローチを使用しています。すべての入力レイヤーは、同じセル解像度でラスターに変換され、適切な評価が与えられます。 各セルの最終的な適合性評価は、各レイヤーの評価と、重要度を反映する重みを組み合わせて計算されます。 提案された土地利用に適さない水域などの領域を除外するために、最終マスクが適用されます。 このアプローチの問題は次のとおりです。 意味のある結果を提供するために大きすぎるセル解像度を選択するか、誤った精度を与える高解像度を選択します。 各入力パラメーターの重みを見つける。 土地適合性マップを作成する他の問題や代替手段はありますか？
1
Rで指定されたラスターの範囲と境界のポリゴンを作成します
2つのポリゴンを作成します。 1つは、ラスターの長方形の範囲です。r <- raster("band5.tif") e <- extent(r) plot(e)私はこれを使用してラスターの範囲を取得できることを知っ ています： 以下に示すように、ラスターの境界のポリゴンを1つ作成するにはどうすればよいですか？ これは私のラスターのようです：
13 raster  r  rgdal 
1
カーネル密度を使用してRの道路密度を計算しますか？[閉まっている]
閉じた。この質問はより集中する必要があります。現在、回答を受け付けていません。 この質問を改善したいですか？この投稿を編集するだけで1つの問題に焦点を当てるように質問を更新します。 閉じた2年前。 道路の大きな（〜70MB）シェープファイルがあり、これを各セルの道路密度を持つラスターに変換したい。理想的には、必要に応じてGDALコマンドラインツールとともにRでこれを実行したいと思います。 私の最初のアプローチは、このスレッドごとに各セルの線分セグメントの長さを直接計算することでした。これにより目的の結果が得られますが、シェープファイルが私のものよりはるかに小さい場合でも非常に遅くなります。正しいセル値が明らかな非常に単純化された例は次のとおりです。 require(sp) require(raster) require(rgeos) require(RColorBrewer) # Create some sample lines l1 <- Lines(Line(cbind(c(0,1),c(.25,0.25))), ID="a") l2 <- Lines(Line(cbind(c(0.25,0.25),c(0,1))), ID="b") sl <- SpatialLines(list(l1,l2)) # Function to calculate lengths of lines in given raster cell lengthInCell <- function(i, r, l) { r[i] <- 1 rpoly <- rasterToPolygons(r, na.rm=T) lc …
3
Rのポリゴンシェープファイルを使用してラスターからラスターを抽出する
私はRが初めてで、ラスターパッケージを使用しています。既存のラスターファイルからポリゴンを抽出する際に問題があります。私が使用する場合 extract(raster, poly_shape) ラスター上の関数は、常にデータを含むリストを作成します。本当に欲しいのは、ArcGISで再びロードできる別のラスタファイルを抽出することです。もう少し読んだ後、クロップ機能が本当に必要だと思います。しかし、この機能を使用しようとすると crop(raster, poly_shape) 私はこのエラーを受け取ります： Error in .local(x, y, ...) : extents do not overlap In addition: Warning message: In intersect(extent(x), extent(y)) : Objects do not overlap ファイルrasterとpoly_shapeは、両方の関数で同じです。ここで何が間違っているのか教えてもらえますか？トリミング機能がリストではなく別のラスターを作成するのは正しいですか？ 編集：extent（）関数は私のために動作しません。それでも同じエラーが発生します。しかし、2つのデータセットは重複していると確信しています！とともに extract(raster, poly_shape) 適切なデータを取得します。ただのリストとしてであり、私が望んでいるようなラスターとしてではありません。前にArcGISにデータセットを読み込んだところ、非常にうまく収まるため、投影を確認しませんでした。今私は試した projection(raster) # "+proj=laea +lat_0=52 +lon_0=10 +x_0=4321000 +y_0=3210000 +ellps=GRS80 +units=m +no_defs" projection(poly_shape) # "+proj=utm +zone=32 +ellps=GRS80 +units=m …
13 raster  shapefile  r  line 
2
ArcGIS Spatial AnalystのFocal Statisticsに相当するQGISを探していますか？
QGIS v2.0.1を使用して、ArcMap v10.1を使用して実行できるのと同じタスクを実行しようとしていますが、同等のツールを見つけるのに問題が発生しています。ArcMapでは、Focal Statistics（mean）Spatial Analystツールを使用していますが、QRASでGRASSまたはSAGAに同等のものがあるかどうかわかりません。 ツールが存在する場合、そのツールの名前を知っていますか？
3
隣接するラスターセルへの値の条件付き割り当て？
値ラスターがあります： m <- matrix(c(2,4,5,5,2,8,7,3,1,6, 5,7,5,7,1,6,7,2,6,3, 4,7,3,4,5,3,7,9,3,8, 9,3,6,8,3,4,7,3,7,8, 3,3,7,7,5,3,2,8,9,8, 7,6,2,6,5,2,2,7,7,7, 4,7,2,5,7,7,7,3,3,5, 7,6,7,5,9,6,5,2,3,2, 4,9,2,5,5,8,3,3,1,2, 5,2,6,5,1,5,3,7,7,2),nrow=10, ncol=10, byrow = T) r <- raster(m) extent(r) <- matrix(c(0, 0, 10, 10), nrow=2) plot(r) text(r) このラスターから、この図に従って現在のセルの8つの隣接するセルに値を割り当てる（または値を変更する）にはどうすればよいですか？このコード行から現在のセル内に赤い点を配置しました。 points(xFromCol(r, col=5), yFromRow(r, row=5),col="red",pch=16) ここで、期待される結果は次のようになります。 現在のセルの値（つまり、値ラスタの5）は0に置き換えられます。 全体として、8つの隣接セルの新しい値は次のように計算する必要があります。 新しい値=赤い長方形に含まれるセル値の平均*現在のセル（赤い点）と隣接するセル間の距離（つまり、斜めに隣接するセルの場合はsqrt（2）、そうでない場合は1） 更新 隣接セルの境界がラスターの制限を超えている場合、条件を尊重する隣接セルの新しい値を計算する必要があります。条件を尊重しない隣接セルは「NA」に等しくなります。 たとえば、[row、col]表記を使用して参照位置がc（5,5）ではなくc（1,1）である場合、右下隅の新しい値のみを計算できます。したがって、期待される結果は次のようになります。 [,1] [,2] [,3] [1,] NA NA NA [2,] NA …
12 raster  r  adjacency 
弊社のサイトを使用することにより、あなたは弊社のクッキーポリシーおよびプライバシーポリシーを読み、理解したものとみなされます。
Licensed under cc by-sa 3.0 with attribution required.