R言語の画像復元プログラムの解読を行っているのですが、文中の(1),(2)式のように、要素が[i,j]から[i,j,1,1]へと変化した理由を教えていただきたいです。行パターンと列パターンの2系統が必要になったからだと考えているのですが、どうでしょうか。
R
1 2library(survey) 3library(imager) 4 5n <- 256 #64 6n2 <- n*n 7img = load.image("House256.jpg") #( "insects.jpg" ) 8grayimg = grayscale(img) 9subi<-imsub(grayimg,x>=90) 10subi<-imsub(subi,x<=256) 11subi<-imsub(subi,y>=100) 12subi<-imsub(subi,y<=256) 13#subi<-imsub(grayimg,x>= 101) 14#subi<-imsub(subi,x<=64) 15#subi<-imsub(subi,y>=141) 16#subi<-imsub(subi,y<=64) 17 18subih<-subi 19subiv<-subi 20rsubih<-subi 21rsubiv<-subi 22rsubi<-subi 23#plot(subi) 24H <- hadamard(n) 25T <- matrix(0, nrow=n, ncol=n) 26suim <- matrix(0, nrow=n, ncol=n) 27suiv <- matrix(0, nrow=n, ncol=n) 28suih <- matrix(0, nrow=n, ncol=n) 29rsuiv <- matrix(0, nrow=n, ncol=n) 30rsuih <- matrix(0, nrow=n, ncol=n) 31rsui <- matrix(0, nrow=n, ncol=n) 32 33for (i in 1:n){ 34 for(j in 1:n){ 35 if(H[i,j]==0) {T[i,j]=-1} else {T[i,j]=1} 36 (1)suim[i,j]=subi[i,j,1,1] 37 } 38} 39 40suiv = T %*% suim / n 41 42suih = t( T %*% t(suim) ) / n 43 44 45# Fukugen 46 47rsuiv = T %*% suiv 48 49rsuih = t( T %*% t(suih) ) 50 51rsui = (rsuiv+rsuih)/2 52 53for (i in 1:n){ 54 for(j in 1:n){ 55 (2)rsubi[i,j,1,1]=rsui[i,j] 56 } 57} 58 59#rsubi=(rsubiv+rsubih)/2 60
もととなる論文を見て考えているのですが、確信が得られませんでした。
同じ内容の質問を何度も繰り返し投稿しないように。
あなたの回答
tips
プレビュー