> restart: > Digits:=14: > st:= time(): > with(LinearAlgebra): > with(PolynomialTools): > with(DiscreteTransforms): > with(Groebner): > with(ImageTools): > read"../Src/One_RGB.mpl": > infolevel[S_stls]:=1: > interface(rtablesize=190): > printlevel:=0: > m:=4.2:k:=6: > Ph:=Read("../Result/Figure.4/ORGB_original.jpg"): > Vx:=Vector(175,i->x^(i-1)):Vy:=Vector(133,i->y^(i-1)): > P1:=RandomMatrix(175,133,density=1,generator=-10..10): > pert_f:=expand(LinearAlgebra:-Transpose(Vx).P1.Vy):#noise with m=4 > mean *10^(-4) > P1:=(x+x^6+y^5+x^6*y^6+y+1)/k:# the blurring function > Pho1:=RGB_Blurred(Ph,P1,10^(-m)*pert_f): > View([Ph,Pho1]); Initializing Java runtime environment. > lcoeff(P1); 1/6 > Write("../Result/Figure.4/ORGB_blur.jpg",Pho1): > RGB_SNR(Ph,P1,10^(-m)*pert_f);# the SNR 63.375752892174 > st:= time()-st; st := 6.734 > st:= time(): > e:=5: > Pho:=RGB_SSD_ImageDeconvolution(Pho1,e):# for one RGB blurred version RGB_SSD_ImageDeconvolution: Before we compute gcd using computing time is st1 = 0.438 evlu1 = 3/10 evlu2 = 1/2 evlu3 = 1/5 evlu1 = 9/10 evlu2 = 3/5 evlu3 = 4/5 Get_AB: the degree of GCD is degx = 6 degy = 6 st = 0.140 Get_AB: evaluation time 1D f and g is st1 = 0. deg = 6 Ap1D_gcd: the time for gcd from cofactor t = 0.016 Get_AB: the GCD's computing time 1D f and g is st1 = 0.266 deg = 6 Get_AB: evaluation time 1D f and g is st1 = 0.031 deg = 6 Ap1D_gcd: the time for gcd from cofactor t = 0.016 Get_AB: the GCD's computing time 1D f and g is st1 = 0.610 deg = 6 Get_AB: evaluation time 1D f and g is st1 = 0.031 deg = 6 Ap1D_gcd: the time for gcd from cofactor t = 0. Get_AB: the GCD's computing time 1D f and g is st1 = 0.391 deg = 6 Get_AB: evaluation time 1D f and g is st1 = 0.031 deg = 6 Ap1D_gcd: the time for gcd from cofactor t = 0.016 Get_AB: the GCD's computing time 1D f and g is st1 = 0.344 deg = 6 Get_AB: evaluation time 1D f and g is st1 = 0.125 deg = 6 Ap1D_gcd: the time for gcd from cofactor t = 0.016 Get_AB: the GCD's computing time 1D f and g is st1 = 0.359 deg = 6 Get_AB: evaluation time 1D f and g is st1 = 0.031 deg = 6 Ap1D_gcd: the time for gcd from cofactor t = 0.015 Get_AB: the GCD's computing time 1D f and g is st1 = 0.406 deg = 6 Get_AB: evaluation time 1D f and g is st1 = 0.032 deg = 6 Ap1D_gcd: the time for gcd from cofactor t = 0.109 Get_AB: the GCD's computing time 1D f and g is st1 = 0.593 deg = 6 Get_AB: the first time GCD's computing time is st1 = 3.250 Get_AB: evaluation time 1D f and g is st1 = 0. deg = 6 Ap1D_gcd: the time for gcd from cofactor t = 0.016 Get_AB: the GCD's computing time 1D f and g is st1 = 0.454 deg = 6 Get_AB: evaluation time 1D f and g is st1 = 0.015 deg = 6 Ap1D_gcd: the time for gcd from cofactor t = 0.015 Get_AB: the GCD's computing time 1D f and g is st1 = 0.703 deg = 6 Get_AB: evaluation time 1D f and g is st1 = 0.016 deg = 6 Ap1D_gcd: the time for gcd from cofactor t = 0. Get_AB: the GCD's computing time 1D f and g is st1 = 1.000 deg = 6 Get_AB: evaluation time 1D f and g is st1 = 0.031 deg = 6 Ap1D_gcd: the time for gcd from cofactor t = 0.016 Get_AB: the GCD's computing time 1D f and g is st1 = 0.719 deg = 6 Get_AB: evaluation time 1D f and g is st1 = 0.031 deg = 6 Ap1D_gcd: the time for gcd from cofactor t = 0.015 Get_AB: the GCD's computing time 1D f and g is st1 = 0.875 deg = 6 Get_AB: evaluation time 1D f and g is st1 = 0.016 deg = 6 Ap1D_gcd: the time for gcd from cofactor t = 0.110 Get_AB: the GCD's computing time 1D f and g is st1 = 0.797 deg = 6 Get_AB: evaluation time 1D f and g is st1 = 0.031 deg = 6 Ap1D_gcd: the time for gcd from cofactor t = 0. Get_AB: the GCD's computing time 1D f and g is st1 = 0.828 deg = 6 Get_AB: the second time GCD's computing time is st1 = 5.516 SD_Dividen: we get Polynomials' coefficient matrix using computing time is st1 = 0. SD_Dividen: The first FFT time is st1 = 0.078 SD_Dividen: The second FFT time is st1 = 0.110 SD_Dividen: The divide time is st1 = 0.562 SD_Dividen: The inverse FFT time is st1 = 0.141 SD_Dividen: we get Polynomials' coefficient matrix using computing time is st1 = 0. SD_Dividen: The first FFT time is st1 = 0.078 SD_Dividen: The second FFT time is st1 = 0.109 SD_Dividen: The divide time is st1 = 0.609 SD_Dividen: The inverse FFT time is st1 = 0.125 SD_Dividen: we get Polynomials' coefficient matrix using computing time is st1 = 0. SD_Dividen: The first FFT time is st1 = 0.234 SD_Dividen: The second FFT time is st1 = 0.110 SD_Dividen: The divide time is st1 = 0.640 SD_Dividen: The inverse FFT time is st1 = 0.125 RGB_SSD_ImageDeconvolution: AFTER we compute gcd using computing time is st1 = 3.969 > View([Ph,Pho1,Pho]);#the original , blurred and deblurring images > Write("../Result/Figure.4/ORGB_deblur.jpg",Pho): > st:= time()-st;#the time is st := 13.328 > >