0

同じ次元の 2 つのバイナリ ファイルがあります。1 つ目は xm と df の間の相関関係を表し、2 つ目は xm と gh データの間の相関関係も表します。これら 2 つの中から最適な相関関係を表す 1 つのマップを作成したいと考えています。例: 1- xm と df の間の相関マップの最初のピクセルと、xm と gh の間の相関マップの対応するピクセルを読み取ります。

2- 最適な相関値を取得し、それが xm と df に由来する場合は青色にします。そうでない場合は、xm と gh に由来する場合は緑色にします

3-すべてのピクセルに対して同じことを行います

4- 関連付けられたマップのようなものを取得します

以下に 2 つのファイルを示します。

  1- to read the first file correlation map:![enter image description here][1]

   conn <- file("C:\\corr.bin","rb")![enter image description here][2]
  corr<- readBin(conn, numeric(), size=4,  n=1440*720, signed=TRUE)
  y<-t(matrix((data=corr), ncol=720, nrow=1440))
    image(y)

2- 2 番目のファイルの土地被覆マップを読み取るには:

  conne <- file("C:\\cor06.bin","rb")
  over<-readBin(conne, numeric(), size=4,  n=1440*720, signed=TRUE)
  y1<-t(matrix((data=over), ncol=720, nrow=1440))
  image(y1)

3-結果を書き込む:

     to.write = file(paste("/orcomplete.bin",sep=""),"wb")
     writeBin(as.double(results), to.write, size = 4)

ここに画像の説明を入力

4

2 に答える 2

1

ベース R の例:

## example data
set.seed(1)
cor1 <- runif(100)
cor2 <- runif(100)

## find max correlation
maxCor <- pmax(cor1, cor2)
## find correct color
col <- ifelse(maxCor==cor1, "blue", "green")
于 2013-03-03T15:27:30.540 に答える
1

データの次元が同じ場合 (あなたの場合はこれが当てはまります)、次のrasterようにパッケージを使用できます。

r <-raster(t(matrix((data=corr), ncol=720, nrow=1440)))
r1 <- raster(t(matrix((data=over), ncol=720, nrow=1440)))
m <- r > r1 #Compare the two rasters
image( m , col = c("#EF8A62" , "#67A9CF" ) ) #Hexadecimal colour specification
legend( "bottomleft" , legend = c( "Y" , "Y1") , fill = c("#EF8A62" , "#67A9CF" ) , border = "#D9D9D9" , bty = "n")

ここに画像の説明を入力

于 2013-03-03T15:18:58.007 に答える