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

私はRの新人です。空間情報を持つラスターデータのNAを0に置き換えるにはどうすればよいですか？

10 raster  r  missing-data  replace 

30m解像度のフォレスト/非フォレストのバイナリラスターから240m解像度のラスターに移動したいのですが、各セルの値はフォレストされたセルの比率です。つまり、240mセルのデータ値は0から始まります（30mのラスター内のすべてのセルがゼロ/非フォレスト）から0.5（30mのセルの半分がフォレスト、半分が非フォレスト）から1（30mのラスター内のすべてのセルがフォレストされた）まで。 離散データに対して双一次補間を使用しないことが推奨されていますが、バイナリ0/1データでの結果がわかる限り、平均値（つまり、1と0の間の比率）を与えることになります。これは賢明な方法ですか、それとももっと良い方法がありますか？ Arc、QGIS、Idrisiを使用できます。

9 raster  resampling 

同じ解像度と範囲を持つ2つ以上のラスターグリッドを説明するために使用できる単語はありますか？私は「偶然」や「同じ」のような言葉で釣り回っていますが、何も正しくないと感じています。そして、「同じ範囲と解像度のラスター」よりもコンパクトなものを見つけたいと思います。

9 raster  terminology 

カラーマップが割り当てられているtiff形式のラスターがいくつかあります。 次に、画像に一致するスペックを使用してモザイクデータセットを作成し、それにすべてのTIFFを追加します（警告は無視してかまいません。これは、私が何をしたかを示すためだけです）。 カラーマップが表示されず、カラーマップをティフからインポートできません。 これらの値に従ってデータを表示できるように、カラーマップを追加/インポートする方法はありますか？

9 arcgis-desktop  raster  arcmap  mosaic-dataset  colormap 

私はlevelplot素晴らしいrasterVisパッケージの関数を使用して、値がゼロを中心に発散するマップを作成しています。赤から青へのパレットを使用してラスターファイルをプロットしていますが、マップで灰色にゼロ値を割り当てようとしています。 具体的には、この図の色を再現したいと思います。 赤から青のスケールが使用されているが、ゼロの値は灰色で着色されていることに注意してください。 現在、私の地図は次のようになっています。 そして、これはそれを再現するためのコードです（https://www.dropbox.com/s/cypfdu1eaz2fuok/r.annual.tif?dl=0で入手可能なファイル）： # Load required packages library(rasterVis) # open file r.annual <- raster("Downloads/r.annual.tif") # Set color palette myTheme=rasterTheme(region=brewer.pal('RdBu', n=11)) # Plot levelplot(annual.mask, par.settings=myTheme, margin=F) 上記のマップのすべてのゼロ値に灰色を割り当てるにはどうすればよいですか？

9 raster  r  legend  visualisation  color-ramp 

Rでベクターをラスターに変換しています。ただし、プロセスが長すぎます。より速く実行するために、スクリプトをマルチスレッドまたはGPU処理に入れる可能性はありますか？ ラスタライズされたベクトルへの私のスクリプト。 r.raster = raster() extent(r.raster) = extent(setor) #definindo o extent do raster res(r.raster) = 10 #definindo o tamanho do pixel setor.r = rasterize(setor, r.raster, 'dens_imov') r.raster クラス：RasterLayer寸法：9636、11476、110582736（nrow、ncol、ncell）解像度：10、10（x、y）範囲：505755、620515、8555532、8651792（xmin、xmax、ymin、ymax）座標。ref。：+ proj = longlat + datum = WGS84 + ellps = WGS84 + towgs84 = 0,0,0 セッター クラス：SpatialPolygonsDataFrame機能：5419範囲：505755、620515.4、8555429、8651792（xmin、xmax、ymin、ymax）座標。ref。：+ proj = utm + zone …

9 raster  r  rasterization  rgdal  parallel 

QGIS 2.4.0での作業Chugiakラスター画像のクリッピングについて質問があります。シェイプファイルの範囲を使用してクリップしたい航空写真があります。これまでのところ、私はラスター/抽出/クリッパーツールを使用してこれを行いました。問題は、出力ラスターのファイルサイズが元のファイル（つまり、それぞれ725.283および249.693 kb）のほぼ3倍であることです。 ファイルが大きいのはなぜですか？そして、出力ファイルのサイズを「設定」する方法はありますか？ クリッパーのヘルプ機能は次のページを参照しています：http: //www.gdal.org/gdal_translate.html私はそれから多くを作ることはできませんが。 私はQGISに比較的慣れていません。

9 qgis  raster  clip  file-size 

ラスターを投影しようとしています。RにはprojectRaster()これを行うための関数があります（完全に再現可能な例の下に）： # example Raster require(raster) r <- raster(xmn=-110, xmx=-90, ymn=40, ymx=60, ncols=40, nrows=40) r <- setValues(r, 1:ncell(r)) projection(r) # project to newproj <- "+init=epsg:4714" # using raster package to reproject pr1 <- projectRaster(r, crs = CRS(newproj), method = 'bilinear') それはうまくいきます。しかし、それはかなり遅いです。 速度を上げるために、gdalwarp代わりに使用します（SSD を使用すると、ディスク/ Rへの読み書きのコストはそれほど高くありません）。 ただし、をprojectRaster()使用した結果を再現できませんgdalwarp。 # using gdalwarp to reproject tf …

9 raster  gdal  r  coordinate-system  gdalwarp 

.adfファイルをRにロードしたい。データは次のページからのものである：http ://www.fao.org/geonetwork/srv/en/metadata.show?id=14057 インターネットでの調査の結果、次のコードを見つけました。問題は、RasterLayerクラスで、そこにあるべきではないネガティブな値を取得することです。なぜこれが起こるのかわかりませんので、誰かが私を助けてくれるといいのですが！？ コード： library(rgdal) library(RColorBrewer) dpath<- path... x <- new("GDALReadOnlyDataset", dpath) getDriver(x) getDriverLongName(getDriver(x)) xx<-asSGDF_GROD(x) r <- raster(xx) 「r」の出力は次のとおりです。 rクラス：RasterLayer寸法：2160、4320、9331200（nrow、ncol、ncell）解像度：0.08333333、0.08333333（x、y）範囲：-180、180、-90、90（xmin、xmax、ymin、ymax）座標。ref。：+ proj = longlat + ellps = WGS84 + towgs84 = 0,0,0,0,0,0,0 + no_defsデータソース：メモリ名内：band1値：-997、16（最小、最大） 値の「16」は、成長期間の長さの16クラスを指します。しかし、それらの「-997」はどこから来たのでしょうか。たぶん、コーディネートに問題があります。ref？ ここにも「xx」のデータの概要があります。 データの概要：最小。第1四半期 中央値第3四半期 マックス。NAの-997 3 5 -9 8 16 7123158 そして、xxのデータをさらに詳しく見てみると、 テーブル（xx $ band1） -997 1 …

9 raster  data  r  esri-grid-format 

指定されたフォルダーからすべての複数のモザイクラスタ（各モザイクは7つのレイヤで構成されています）をR にインポートします。次に、それらを単一のマルチレイヤラスタとしてアクセスします。 だから私はしました： # read all mosaics named "mos....img" in R raster_data<-list.files(path=getwd(), pattern="mos.*.img$") # read files as rasters s <- stack(raster_data) # check my imported rasters p.ex. raster n°8 from "s" raster stack s[[8]] また、ラスターs [[8]]には1つのレイヤーしか含まれていないため、モザイク全体がインポートされていません。 nlayers(s[[8]]) [[1]] 各モザイクを個別に読んだ場合、それは機能します： # read 1 mosaic (composed by 7 bands) mosaic1<-brick("mosaic1.img") # extract …

9 raster  r  remote-sensing  landsat 

休業。この質問には、より焦点を当てる必要があります。現在、回答を受け付けていません。 この質問を改善してみませんか？質問を更新して、この投稿を編集するだけで1つの問題に焦点を当てます。 2年前休業。 OGR / GDALを使用してPythonでスクリプトを作成しています。 シェープファイルのセットとGeoTiffラスターファイルのセットがあります。 シェープファイルがラスター領域と交差しない場合、スクリプトに無視させたいのですが。 シェープファイルは長方形ではないため、layer.GetExtent（）から返されたxmin / xmax、ymin / ymaxの値を単純に比較することはできません。全体的な形状を表す実際のポリゴンが必要です。次に、そのポリゴンがラスタースクエアと交差するかどうかを判断する方法が必要です。 どういうわけか、シェープファイル内のすべてのポリゴンを1つのフィーチャにマージし、そのフィーチャのジオメトリを読み取って、その情報をラスター範囲と比較できると考えていました。ただし、具体的にどのように実行するかはわかりません。 シェープファイルからボーダーポリゴン情報を抽出する方法は？ そのポリゴンが特定の正方形の領域と交差するかどうかを判断するにはどうすればよいですか？

9 python  raster  shapefile  gdal  ogr 

スタックされたラスターから最大セル値を見つけるにはどうすればよいですか。 Rmax <- maxValue(RAD1998.all[[1]]) 正常に動作しますが Rmax <- maxValue(RAD1998.all[[2]]) NAを与えます。 もちろん、スタックされたラスターではありません。 これが私のコードです： RAD1998 <- raster(paste(getwd(), "/1998bil/1998ASC5min_ppt_spas1214_0001_19980202_0810_UTC.asc.bil", sep = "")) list.ras <- mixedsort(list.files(paste(getwd(), "/1998bil/", sep = ""), full.names = T, pattern = ".asc.bil")) RAD1998.all <- stack(list.ras)

9 raster  r 

カーネル密度出力の各レベルの面積を測定できるArcMapのツールまたはメソッドはありますか？ これは単なるカーネル密度プロットであることに注意してください（シェープファイルやポリゴンではありません）。カーネル密度分析（濃い緑（10％）から赤（90％））（1） いくつか試しましたが、正確である必要があります。

9 arcgis-desktop  raster  arcmap  kernel-density 

私はいつものように現在のQGISマスターをコンパイルしてインストールし、postgres 9.3データベースにpostgis 2.0.2拡張子を付けていくつかのラスターデータをインポートしました。データベースマネージャーを使用してラスターデータをqgisにインポートしようとすると、次のメッセージが表示されます。 GDAL provider: Cannot get GDAL raster band: Raster layer: Provider is not valid (provider: gdal, URI: PG: dbname=gis host=127.0.0.1 user=xxx password=xxx port=5432 mode=2 schema=public table=elevation コマンドラインで、 ERROR 1: Computed PostGIS Raster dimension is invalid. You've probably specified unappropriate resolution. エラーを検索する場所がわかりません。 Postgis 2.1.2に移動したところ、同じ問題が見られます。

9 qgis  postgis  raster  gdal  raster2pgsql 

同じ領域の2つの.ecwレイヤーが重なり合っていますが、1つは中央がクリッピングされ、白いエッジが付いています。背景の値を255、255、255に設定し、それを中空として定義すると、特定のRGB値のみが中空になります。実際には、レイヤーに次のような他の「白い」値があります。 254,254,254; 253,253,253 .. 249,249,249まで。 ディスプレイの背景の値を定義して、複数の一意のRGBセットの値を含める方法はありますか？ これを回避できることはわかっていますが（クリッピング、ラスター計算機など）、重複を生成せずに元の.ecwファイルを1つだけ保持したいので、Arcmapシンボルまたはイメージレンダリングツールのみを使用します。 私が必要とするものはCustom transparency options、QGISのとまったく同じですが、ArcMap用です。

9 raster  arcmap  arcgis-10.2  ecw  display