【发布时间】:2016-12-13 22:27:22
【问题描述】:
我正在编写手工 jpeg 算法。但是,量化步骤似乎有问题: 这是我要编码的图片:Lena.png
这是量化步骤之后的结果(+解码但没有哈夫曼编码的东西):post quantization result
现在我的这一步代码:
def quantification(transformee) :
matrice=np.array([(16 , 11 , 10 , 16 , 24 , 40 , 51 , 61 ),(12 , 12 , 14 , 19 , 26 , 58 , 60 , 55 ),(14 , 13 , 16 , 24 , 40 , 57 , 69 , 56 ),(14 , 17 , 22 , 29 , 51 , 87 , 80 , 62 ),(18 , 22 , 37 , 56 , 68 , 109 , 103 , 77 ),(24 , 35 , 55 , 64 , 81 , 104 , 113 , 92 ),(49 , 64 , 78 , 87 , 103 , 121 , 120 , 101 ),( 72 , 92 , 95 , 98 , 112 , 100 , 103 , 99)],dtype=float)
quanti=[]
for i in transformee :
k=i/matrice
quanti+=[np.around(k)]
return quanti
def dequantification(quanti) :
matrice=np.array([(16 , 11 , 10 , 16 , 24 , 40 , 51 , 61 ),(12 , 12 , 14 , 19 , 26 , 58 , 60 , 55 ),(14 , 13 , 16 , 24 , 40 , 57 , 69 , 56 ),(14 , 17 , 22 , 29 , 51 , 87 , 80 , 62 ),(18 , 22 , 37 , 56 , 68 , 109 , 103 , 77 ),(24 , 35 , 55 , 64 , 81 , 104 , 113 , 92 ),(49 , 64 , 78 , 87 , 103 , 121 , 120 , 101 ),( 72 , 92 , 95 , 98 , 112 , 100 , 103 , 99)],dtype=float)
retour_transformee=[]
for k in quanti :
temp=np.zeros((8,8),dtype=float)
for i in range(8) :
for j in range(8) :
temp[i,j]=matrice[i,j]*k[i,j]
retour_transformee.append(temp)
return retour_transformee
transformee 是一个 8*8 矩阵的列表(已使用 DCT)。 量化矩阵是维基百科中给出的(JPEG)
抱歉算法中的法语部分
【问题讨论】:
标签: compression jpeg quantization