2016-10-24 8 views
0

私はマップの画像だけを使って領域の面積を推定する方法を探しています。私がこれをやっているのは、ある海面水位が上がったときに失われる領域を計算したいのです。マップ(画像形式)だけのメタデータは見つけられません。それでは、私が考えていることはEBimageパッケージを使用して、グレースケール画像にこの画像を変換することであるRを使った領域推定R

、その後数をカウントするための基準として、画素強度を使用して:ここでは、そのようなマップへのリンクです潜在的に脅威となる領域を表すピクセルの数。

私の質問は可能ですか?基準としてどのようにピクセル強度を使用できますか?そして、もしこの問題を解決する他の方法があれば?

また、そのようなマップをプロットするのに使用されるメタデータにアクセスする方法がある場合は、わからないことを教えてください。

ありがとうございます。

編集:私はrgdalパッケージを使用してグリッドデータのint型のRを読み取ることができましたhrbrmstrする

ありがとうございます。したがって領域を計算するために、rgeosパッケージを使用しようとしましたが、CRESISのデータセットにはシェイプファイルがありません。どうやってポリゴンを定義し、面積を計算できますか? この質問が愚かであれば申し訳ありません。これは私が今まで、空間データと分析を扱ったのは初めてです

+0

なぜあなたは[ESRIを使用していませんグリッドファイル](https://www.cresis.ku.edu/content/research/maps)? – hrbrmstr

+0

私はどのようにESRIファイルを扱うか分かりません。私はGISリーダーをいくつか試しますが、運はありません。 RArcinfoを試しましたが、エラーがありました。 – Mathnoob

答えて

0

データは、ESRIファイルである:

library(rgdal) 

grid_file_1m <- new("GDALReadOnlyDataset", "/full/path/to/inund1/w001001.adf") 
grid_1m <- asSGDF_GROD(grid_file_1m, output.dim=c(1000, 1000)) 

plot(grid_1m, bg="black") 

enter image description here

grid_1m_df <- as.data.frame(grid_1m) 

str(grid_1m_df) 
## 'data.frame': 2081 obs. of 3 variables: 
## $ band1: int 1 1 1 1 1 1 1 1 1 1 ... 
## $ x : num -77.2 -76.9 -76.5 -76.1 -75.8 ... 
## $ y : num 83.1 83.1 83.1 83.1 83.1 ... 
+0

ありがとうございます。これはもっと良い方法でしょう。 – Mathnoob

+0

私の提案したアプローチを使用してこの問題を解決する方法は、好奇心からだけです。 – Mathnoob

+0

これを読めば、数日間の勉強を必要とする有効な統計学的計算を行うための組み込みの方法がある場合、見た目と組み合わせた損失の多い画像操作を実行したいのですか?もしそうなら、悲しいことに、はい。 – hrbrmstr

関連する問題