在python3下对数据分块(8x8大小)使用OpenCV的离散余弦变换DCT
在MATLAB中有blkproc (blockproc)对数据处理, 在python下没找到对应的Function,
这里利用numpy 的split(hsplit和vsplit) 对数据分块处理成8x8的小块, 然后在利用OpenCV的dct函数做变换, 同时利用idct 验证数据变换是否正确.
import numpy as np
import cv2
a = np.arange(256).reshape((16,16))
print("ori data: \n{}".format(a))
m,n = a.shape
print(m,n)
# Y = np.zeros(256).reshape((16,16))
hdata = np.vsplit(a,n/8) # 垂直分成高度度为8 的块
for i in range(0, n//8):
blockdata = np.hsplit(hdata[i],m/8)
#垂直分成高度为8的块后,在水平切成长度是8的块, 也就是8x8 的块
for j in range(0, m//8):
block = blockdata[j]
print("block[{},{}] data \n{}".format(i,j,blockdata[j]))
Yb = cv2.dct(block.astype(np.float))
print("dct data\n{}".format(Yb))
iblock = cv2.idct(Yb)
print("idct data\n{}".format(iblock))
以下是最后个8x8块的数据:
block[1,1] data
[[136 137 138 139 140 141 142 143]
[152 153 154 155 156 157 158 159]
[168 169 170 171 172 173 174 175]
[184 185 186 187 188 189 190 191]
[200 201 202 203 204 205 206 207]
[216 217 218 219 220 221 222 223]
[232 233 234 235 236 237 238 239]
[248 249 250 251 252 253 254 255]]
dct data
[[ 1.56400000e+03 -1.82216412e+01 0.00000000e+00 -1.90481783e+00
0.00000000e+00 -5.68239222e-01 0.00000000e+00 -1.43407825e-01]
[-2.91546259e+02 0.00000000e+00 0.00000000e+00 0.00000000e+00
0.00000000e+00 0.00000000e+00 0.00000000e+00 0.00000000e+00]
[ 0.00000000e+00 0.00000000e+00 0.00000000e+00 0.00000000e+00
0.00000000e+00 0.00000000e+00 0.00000000e+00 0.00000000e+00]
[-3.04770852e+01 0.00000000e+00 0.00000000e+00 0.00000000e+00
0.00000000e+00 0.00000000e+00 0.00000000e+00 0.00000000e+00]
[ 0.00000000e+00 0.00000000e+00 0.00000000e+00 0.00000000e+00
0.00000000e+00 0.00000000e+00 0.00000000e+00 0.00000000e+00]
[-9.09182756e+00 0.00000000e+00 0.00000000e+00 0.00000000e+00
0.00000000e+00 0.00000000e+00 0.00000000e+00 0.00000000e+00]
[ 0.00000000e+00 0.00000000e+00 0.00000000e+00 0.00000000e+00
0.00000000e+00 0.00000000e+00 0.00000000e+00 0.00000000e+00]
[-2.29452520e+00 0.00000000e+00 0.00000000e+00 0.00000000e+00
0.00000000e+00 0.00000000e+00 0.00000000e+00 0.00000000e+00]]
idct data
[[136. 137. 138. 139. 140. 141. 142. 143.]
[152. 153. 154. 155. 156. 157. 158. 159.]
[168. 169. 170. 171. 172. 173. 174. 175.]
[184. 185. 186. 187. 188. 189. 190. 191.]
[200. 201. 202. 203. 204. 205. 206. 207.]
[216. 217. 218. 219. 220. 221. 222. 223.]
[232. 233. 234. 235. 236. 237. 238. 239.]
[248. 249. 250. 251. 252. 253. 254. 255.]]
数据与原数据值大小一致.
最新文章
- linux基础命令
- poj1155 TELE (树上的背包)
- ActiveMQ 复杂类型的发布与订阅
- ASP.NET MVC3入门教程之参数(数据)传递
- linux重新设定分区大小
- spring 第一篇(1-1):让java开发变得更简单(下)转
- 如何将扩展名为.backup的文件导入postgresql中 求步骤 新手 谢谢.
- [Flex] ButtonBar系列——flex3 ButtonBar圆角菜单的运用
- C#基础性问题
- Hibernate4.x之Session--常用方法
- [AngularJS] angular-schema-form -- 1
- 高性能MySql进化论(九):查询优化器常用的优化方式
- BZOJ1653: [Usaco2006 Feb]Backward Digit Sums
- 教你如何实现微信小程序与.net core应用服务端的无状态身份验证
- PHP全路径无限分类原理
- mpvue 转小程序实践总结
- 移动端h5页面的那些坑
- 离线安装IDEA插件
- Eclipse在Debug模式下经常进入ThreadPoolExecutor解决办法
- 转error while loading shared libraries的解決方法
热门文章
- ios 一个正则表达式测试(只可输入中文、字母和数字)
- Animate a custom Dialog,自定义Dialog动画
- Hibernate更新数据报错:a different object with the same identifier value was already associated with the session: [com.elec.domain.ElecCommonMsg#297e35035c28c368015c28c3e6780001]
- java web页面跳转 总结
- 如何给MFC的菜单项添加快捷键
- 5501环路运输【(环结构)线性DP】【队列优化】
- poj1015 Jury Compromise【背包】
- SignalR循序渐进(一)简单的聊天程序
- hdu4686 简单的矩阵快速幂求前n项和
- 一个非常棒的jQuery 评分插件--好东西要分享