R言語の画像処理におけるエッジ検出プログラムの画素値について

Question

R言語を使った画像処理においてエッジ検出プログラムの解読を行っているのですが、文中の最後の部分で、なぜ0.03を境に画素値を0と1に変換しているかという問題なのですが、この画像の平均輝度が0.03だからという理由で良いのでしょうか？

library(imager)

n <- 256    #64
n2 <- n*n

m <- n/4
img = load.image("House.jpg")    #( "insects.jpg" )
grayimg = grayscale(img)
subi<-imsub(grayimg,x>=90)
subi<-imsub(subi,x<=256)
subi<-imsub(subi,y>=100)
subi<-imsub(subi,y<=256)

#subi<-imsub(grayimg,x>= 1)
#subi<-imsub(subi,x<=256)
#subi<-imsub(subi,y>=55)
#subi<-imsub(subi,y<=256)

#subi<-imsub(grayimg,x>= 385)
#subi<-imsub(subi,y>=100)
#subi<-imsub(subi,y<=256)

#subi<-imsub(grayimg,x>= 101)
#subi<-imsub(subi,x<=64)
#subi<-imsub(subi,y>=141)
#subi<-imsub(subi,y<=64)

subih<-subi
subiv<-subi

# Vertical edge
for (j in 1:n){
   for (i in 1:m){
      tmp=subiv[(i-1)*4+1,j,1,1]+subiv[(i-1)*4+2,j,1,1]-subiv[(i-1)*4+3,j,1,1]-subiv[(i-1)*4+4,j,1,1]
      subiv[(i-1)*4+1,j,1,1]=abs(tmp/4)
      subiv[(i-1)*4+2,j,1,1]=abs(tmp/4)
      subiv[(i-1)*4+3,j,1,1]=abs(tmp/4)
      subiv[(i-1)*4+4,j,1,1]=abs(tmp/4)
   }
}

gsubiv= subiv + 0.9
gsubiv[ gsubiv > 1 ] = 1

# Horizontal edge
for (i in 1:n){
   for (j in 1:m){
      tmp=subih[i,(j-1)*4+1,1,1]+subih[i,(j-1)*4+2,1,1]-subih[i,(j-1)*4+3,1,1]-subih[i,(j-1)*4+4,1,1]
      subih[i,(j-1)*4+1,1,1]=abs(tmp/4)
      subih[i,(j-1)*4+2,1,1]=abs(tmp/4)
      subih[i,(j-1)*4+3,1,1]=abs(tmp/4)
      subih[i,(j-1)*4+4,1,1]=abs(tmp/4)
   }
}

#gsubih= subih + 0.9
#gsubih[ gsubih > 1 ] = 1

gsubi= (subiv+subih)/2

########################

#gsubi= gsubi + 0.9
#gsubi[ gsubi > 1 ] = 1

########################

gsubi[ gsubi > .03] = 1        #.07
gsubi[ gsubi <= .03 ] = 0
plot( gsubi )

Answer 1

R言語は全く知りませんので，
この処理の結果として「エッジか否かの２値画像」が得られている，という仮定の上での話になりますが…

その0.03というのは単なる閾値でしょう．
つまり，エッジ検出処理の目的に合わせて適当に調整すべき処理パラメータなので，
「今現在の値が何故0.03なのか？」に特別な理由を求めなくても良いと思います．

試しにこの値を増やしたり減らしたりして結果を比較してみてはいかがでしょうか．

Comments

[masa_0705]

回答ありがとうございます。値を増減した結果、理由がはっきりと分かりました。