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

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

3
PyQGISまたはPythonを使用して多くのラスターを整列させますか？
PythonまたはPyQGISを使用して、多くのラスターをいくつかのピクセルサイズ、範囲、および投影システムで整列させたい。 最初に考えるのは、GDALを使用することです。 gdalinfo（参照ラスターから情報を検索するため） gdalwarp（投影システムとピクセルサイズを変換するため） gdal_translate（範囲を変換するため） それらは機能しますが、Pythonでこれらのツールを一緒に使用するのは簡単ではなく、この作業を完了するにはかなりの時間が必要です。また、QGISおよびAlign Rasters Tool を使用して簡単に行うことができます。 PyQGISまたはPythonを使用してこれを行うツールはありますか？（より多くのプログラミングを自動的に行いたい） 更新 私はこの質問からこのコードを見つけます： from osgeo import gdal, gdalconst inputfile = #Path to input file input = gdal.Open(inputfile, gdalconst.GA_ReadOnly) inputProj = input.GetProjection() inputTrans = input.GetGeoTransform() referencefile = #Path to reference file reference = gdal.Open(referencefile, gdalconst.GAReadOnly) referenceProj = reference.GetProjection() referenceTrans = reference.GetGeoTransform() bandreference …
8 qgis  python  raster  pyqgis  gdal 
1
Rでのリサンプルと集約および拡張を使用して、解像度と範囲が一致するラスターを作成する
解像度と範囲が異なる2つのラスターがあります。 > res(Elevation) [1] 0.002083333 0.002083333 > res(Ann_precip) [1] 0.008333333 0.008333333 > extent(Elevation) class : Extent xmin : -15.07722 xmax : -7.641806 ymin : 7.193611 ymax : 12.67694 > extent(Ann_precip) class : Extent xmin : -15.075 xmax : -7.641667 ymin : 7.191667 ymax : 12.675 私の質問は、これらの2つのラスターが一致する解像度と範囲を持つためには、 A）raster::aggregate関数を使用する > 0.008333333/0.002083333 [1] …
1
qgis2threejsプラグインのASCIIグリッドのNODATA_VALUEを処理しますか？
プロジェクトでは、地下の地層の上限と下限をモデル化しました。結果のラスタ​​ーはNODATA_VALUE、-999のASCIIファイルとして保存されます。ただし、qgis2threejsプラグインを使用してASCIIファイルを視覚化すると、ピクセルまたはセルがnodata0として表示されます（下の画像を参照）。 qgis2threejsプラグインがnodata値を持つセルを認識し、有効な値を持つラスターの一部のみを視覚化することは可能ですか？
3
GDALを使用してラスターを再分類しますか？
複数の離散値を含むラスターを再分類できるGDALのツールを探しています。これまでのところ、GDALとPythonを使用してラスタを再分類しています。ただし、ここでは単一の値しか処理できないようです。 スイートのネイティブツールのようなものはありますか？
3
英国の炭鉱のWMSレイヤーにアクセスしていますか？
UK Coal Authorityの WebサイトからWMSレイヤーにアクセスしようとしています。私は彼らのサーバーに接続して利用可能なラスターファイルを見ることができますが、それらをダウンロードすると、最初と最後のサービスからのラスターしか表示できません。 石炭鉱業報告地域-https://map.bgs.ac.uk/arcgis/services/CoalAuthority/coalauthority_coal_mining_reporting_areas/MapServer/WMSServer-石炭当局が鉱業報告書を推奨する地域。 INSPIRE準拠データ- https://map.bgs.ac.uk/arcgis/services/CoalAuthority/coalauthority_inspire/MapServer/WMSServer -石炭庁の監視ポイントとライセンス領域を備えるデータセット。 他は私のために働きません。これらは私のステップでした： 新しいWMS接続を作成し、Webサイトから名前とコピーURLの場所のみを追加します。 [ 接続]をクリックして、ラスターレイヤーを選択します。 凡例とともにレイヤーが追加されますが、マップキャンバスには何も表示されません。また、レイヤーのプロパティを確認したところ、データの行/列がないことがわかりました。 手順またはラスターに何か問題がありますか？ 疑わしいですが、私はQGIS2.16.0-Nødeboを使用しています。
2
Rでラスターを再投影すると、投影された点が有限ではないという警告が表示されますか？
1.質問 RのラスターパッケージでprojectRaster（）関数を使用して警告が発生しました。完全に再現可能な例を以下に貼り付けます。 Warning message: In rgdal::rawTransform(projto_int, projfrom, nrow(xy), xy[, 1], : 33940 projected point(s) not finite 私の質問は次のとおりです。このデータは、地上データで作業している場合に修正する必要がある問題ですか？つまり、データは「失われた」ということです。もしそうなら、これは私にとって大きな問題になるでしょう。もしそうなら、私がそれを修正できる方法があるかどうか知っていますか？ 私はこの問題の解決策をオンラインで検索したところ、ここ、ここ、およびここにいくつかの言及が見つかりましたが、この問題に対する適切な答えを提供するものはないと思います。 2.ラスターライブラリを読み込む library(raster) 3.まず、各グリッドセルに定数を使用して、世界の（長緯度）マップを作成します 各グリッドセルに定数を配置して、警告がデータに影響する可能性がある場合の問題を診断できるかどうかを確認しています。 rastertest.longlat<-raster(ncol=360, nrow=180) values(rastertest.longlat)<-c(rep(1,n=180*360)) 4.（長緯度）マップを等面積グリッドに再投影すると、警告メッセージが表示されます rastertest.eck4<-projectRaster(rastertest.longlat, res=c(100000,100000), crs="+proj=eck4",method="ngb", over=T) Warning message: In rgdal::rawTransform(projto_int, projfrom, nrow(xy), xy[, 1], : 33940 projected point(s) not finite このメッセージは基本的に、一部のグリッドセルで再投影が失敗したことを示していると思います。 5.しかし、2つのマップをプロットした場合、この警告はデータに問題を引き起こすようには見えません つまり、プロットされたデータに白いギャップはありません。私の推測では、失われた細胞は非地球上の細胞であり、世界の最上部と最果てです。何か案は？ par(mfrow=c(2,1)) plot(rastertest.longlat, col="blue") …
3
qgisラスターパレットが表示されない
メッキスタイルで新しいシングルバンドラスタのスタイルを設定しようとしていますが、「メッキ」を選択すると、色や値が表示されず、それらを追加する方法が見つかりません。 QGISエッセン、勝利7。
8 qgis  raster  style 
1
ラスター計算機が奇妙な出力を作成する
2つのラスターレイヤーを結合しようとしています（1つは夜間の衛星データ、もう1つはGPW人口密度データです）。2つを組み合わせるためにラスター計算機を使用しています（以下は除算操作ですが、加算または減算でも同じ問題が発生します）が、0またはnodataの傾向がある値のストライプを取得します。これがどのように見えるかを示すためにスクリーンショットを添付します。 なんでこんなことが起こっているの？ラスターレイヤーが同じシステム（WGS84）で地理参照されていることを確認し、ラスターの位置合わせを試みましたが、効果がありませんでした。これが新人の質問なら申し訳ありません-私はソフトウェアにかなり慣れていません。 編集：http: //beta.sedac.ciesin.columbia.edu/data/set/gpw-v4-population-density gpw is 2015 version of this data 衛星データはこちら：http : //ngdc.noaa.gov/eog/dmsp/downloadV4composites.html 特に2013年の安定したライトレイヤー
3
大きなラスターイメージを送信する効率的な方法はありますか？
かなり単純な質問：ArcSDEに格納されている2つの大きなラスターイメージがあります。これらの画像を他の誰かに送信し、ftpサイトに投稿する必要があります。これらのラスターを送信する最も効率的な方法は何ですか？ ジオデータベースを作成して画像にコピーし、データベースに貼り付けてみましたが、それには時間がかかりました。より速い方法はありますか？
1
非水平面の地理参照ラスター
（垂直）壁から生成されたラスターデータがあります。このデータをpostgisデータベースに保持し、空間参照を「最も正確な方法」でエンコードします。 現在、それらはメトリックCRSを悪用することによって保存され、壁のz座標をyとして、壁の左側からのオフセットをxとしてエンコードします。これにより、その目的のために機能するローカル参照システムが提供されますが、グローバルコンテキストは失われます。 ベクトルデータの場合、すべての頂点に3D座標を与えて（グローバル）スペースに配置するのは簡単です。これは、ラスターデータに基づいて作成する必要があります（GISユーザーインターフェイスを使用して、これらの壁の上にある関心領域をデジタル化します）。 さらに、複数の壁を隣り合わせに配置することもでき、この状況でこれらを視覚化できるはずです（同じ方位角の場合にのみ機能すれば十分です）。 これに対処する方法はいくつかあります。 実際の座標に基づいた原点を持つ垂直空間でカスタムCRSを使用します。ただし、この「元の参照」が正確にどこに格納されるかはまだ不明です。 情報をCRSに保存します（それは可能ですか？）-各参照面にいくつかの異なるCRSが必要になります。 ラインに外部キーを使用する（サンプルの赤いラインを参照）-現在の状況、冗長な情報（ラインの長さがラスターの幅に対応していない場合はどうなりますか？） 参照面として3Dポリゴンを作成-冗長情報、上記を参照 線の原点を作成します。これは、線の方位角と組み合わせて参照面にすることができます-異なる壁が同じ参照面を共有しますか？ すべてのアプローチは、どういうわけか「回避策」であるように思われ、注意が必要です。 以下の2つの画像は、状況の上面図といくつかの正面ラスター画像の構成を示しています。（それらが単一の参照面にマッピングされている場合は問題ありません） 水平空間や標高情報で地理的コンテキストを失わずに、垂直ラスターイメージをデータベースに保存する最も適切な方法は何ですか？
1
PostGISでラスターからポリゴンの平均値を計算していますか？
まず、同僚から提供された英国のNetCDF .griおよび.grdラスターファイルから始めました。私はそれをRでロンドンだけになるようにクリップし、エクスポートしてASCファイルに変換してから、Rで次のコマンドを使用してPostGISにインポートしました。 library(raster) uk_raster <- raster("AnnMean2011.grd") london_area <- extent(-720000.0,-630000.0,-50000.0,25000) london_raster <- crop(uk_raster, london_area) writeRaster(london_raster, filename="AnnMean2011.asc", format="ascii") そして、Ubuntuコマンドラインで： raster2pgsql -I -C -s 10001 -t 20x20 AnnMean2011.asc annualmean | psql -d james_traffic PostGISにラスターテーブルができました。ちなみに、10001のSRIDは次のとおりです。 INSERT INTO spatial_ref_sys(srid, auth_name, auth_srid, proj4text) VALUES (10001,'CMAQ_Urban',10001,'+proj=lcc +a=6370000 +b=6370000 +lat_1=35 +lat_2=65 +lat_0=52 +lon_0=10 +x_0=000000 +y_0=00000'); 同じデータベースに、ロンドンをカバーするポリゴンファイルSRID 27700があります。ラスターから各ポリゴン内の平均値を計算したいと思います。 私はこのようなことを試みていますが、それは正しくありません： …
2
QGISを使用してフィーチャまでの距離を持つラスターを作成しますか？
私はGIS（特にQGIS）を初めて使用します。特定の地域の海岸線までの距離を含むラスターレイヤーを作成しようとしています。 海岸線（ライン）と問題の領域（ポリゴン）の両方のベクターレイヤーがあります。 距離を測定するオブジェクトとしてラインを、マスクとしてポリゴンを使用してラスターを作成するツールはありますか？
3
Rのピクセルの塊を削除する
ラスターイメージから孤立したピクセル（または9未満のピクセルの塊）を削除したいと思います。 library(raster) # create some raster data r <- raster(ncols=12, nrows=12) set.seed(0) r[] <- round(runif(ncell(r))*0.7 ) rc <- clump(r) ピクセルの塊を削除する前<9 ピクセルの集まりを削除した後<9： Erdasには、これを行うSieveツールがありますが、Rでそれを置き換える方法は？
弊社のサイトを使用することにより、あなたは弊社のクッキーポリシーおよびプライバシーポリシーを読み、理解したものとみなされます。
Licensed under cc by-sa 3.0 with attribution required.