基于matlab的运动模糊图像处理,基于matlab的运动模糊图像处理

基于matlab的运动模糊图像处理

基 于 MATLAB 的 运 动 模 糊 图 像 处理提 醒 ：我 参 考 了 文 献 里 的 书 目 和 网 上 的 一 些 代 码 而 完 成 的 ， 所 以 误 差 会 比 较 大 ， 目 前对 于 从 网 上 下 载 的 模 糊 图 片 的 处 理 效 果 很 不 好 ， 这 是 我 第 一 次 上 传 自 己 完 成 的实 验 的 文 档 ， 希 望 能 帮 到 一 些 人 吧 。研 究 目 的在 交 通 系 统 、 刑 事 取 证 中 图 像 的 关 键 信 息 至 关 重 要 ， 但 是 在 交 通 、 公 安 、 银 行 、医 学 、 工 业 监 视 、 军 事 侦 察 和 日 常 生 活 中 常 常 由 于 摄 像 设 备 的 光 学 系 统 的 失 真 、调 焦 不 准 或 相 对 运 动 等 造 成 图 像 的 模 糊 ， 使 得 信 息 的 提 取 变 得 困 难 。 但 是 相 对于 散 焦 模 糊 ， 运 动 模 糊 图 像 的 复 原 在 日 常 生 活 中 更 为 普 遍 ， 比 如 高 速 运 动 的 违规 车 辆 的 车 牌 辨 识 ， 快 速 运 动 的 人 群 中 识 别 出 嫌 疑 人 、 公 安 刑 事 影 像 资 料 中 提取 证 明 或 进 行 技 术 鉴 定 等 等 ， 这 些 日 常 生 活 中 的 重 要 应 用 都 需 要 通 过 运 动 模 糊图 像 复 原 技 术 来 尽 可 能 地 去 除 失 真 ， 恢 复 图 像 的 原 来 面 目 。 因 此 对 于 运 动 模 糊图 像 的 复 原 技 术 研 究 更 具 有 重 要 的 现 实 意 义 。图 像 复 原 原 理本 文 探 讨 了 在 无 噪 声 的 情 况 下 任 意 方 向 的 匀 速 直 线 运 动 模 糊 图 像 的 复 原 问 题 ，并 在 此 基 础 上 讨 论 了 复 原 过 程 中 对 点 扩 散 函 数 (PSF)的 参 数 估 计 从 而 依 据 自 动鉴 别 出 的 模 糊 方 向 和 长 度 构 造 出 最 为 近 似 的 点 扩 散 函 数 ， 构 造 相 应 的 复 原 模 型 ，实 现 运 动 模 糊 图 像 的 复 原 ； 在 模 糊 图 像 自 动 复 原 的 基 础 上 ， 根 据 恢 复 效 果 图 的纹 理 特 征 和 自 动 鉴 别 出 的 模 糊 长 度 和 角 度 ， 人 工 调 整 模 糊 方 向 和 长 度 参 数 ， 使得 复 原 效 果 达 到 最 佳 。实 验 过 程模 糊 方 向 的 估 计 ：对 图 1(a)所 示 的 原 始 图 像 ‘车 牌 ’图 像 做 方 向 = ， 长 度30L=20 像 素 的 匀 速 直 线 运 动 模 糊 ， 得 到 退 化 图 像 如 图 1(b)1(a) 1(b)j=imread( 车 牌 1.jpg ); figure(1),imshow(j); title( 原 图 像 ); len=20; theta=30;psf=fspecial( motion ,len,theta);j1=imfilter(j,psf, circular , conv );figure,imshow(j1);title( PSF 模 糊 图 像 ); 图 1(c)和 1(d)分 别 为 原 图 像 和 模 糊 图 像 的 二 次 傅 里 叶 变 化1(c) 1(d)J=rgb2gray(j);K=fft2(J);%傅 里 叶 变 换 M=fftshift(K);%直 流 分 量 移 到 频 谱 中 心N=abs(M);%计 算 频 谱 幅 值P=(N-min(min(N)))…/(max(max(N))-min(min(N)))*225;%归 一 化figure,imshow(P);title( 傅 里 叶 变 换 频 谱 );J1=rgb2gray(j1);K1=fft2(J1);%傅 里 叶 变 换 M1=fftshift(K1);%直 流 分 量 移 到 频 谱 中 心N1=abs(M1);%计 算 频 谱 幅 值P1=(N1-min(min(N1)))/(max(max(N1))…-min(min(N1)))*225;%归 一 化figure,imshow(P1);title( 傅 里 叶 变 换 频 谱 );利 用 图 1(d)粗 略 的 计 算 模 糊 的 方 向 ， 可 以 通 过 matlab 自 带 的 画 线 工 具 ，选 取 如 下 图 1(e)所 示 的 三 角 形 ， 计 算 a 与 c 之 间 的 夹 角Matlab 命 令 窗 口 输 ： g选 中 三 个 顶 角 后 回 车 ， 可 得 数 据ans =136.6009 146.3977137.4205 166.8895145.2074 166.8895粗 略 取 值 后 ， 经 Matlab 计 算 得 ：atan(20/9)*180/pians = 65.7723 则 运 动 方 向 为 90-ans≈ 25( 误 差 为 5) 。1(e) 模 糊 长 度 的 估 计运 动 模 糊 图 像 中 ， 在 运 动 方 向 上 大 多 数 模 糊 图 像 的 背 景 像 素 点 具 有 很 强的 相 关 性 ， 即 沿 着 运 动 模 糊 的 轨 迹 ， 背 景 像 素 点 的 灰 度 值 逐 渐 变 化 或 者不 变 。 通 过 文 献 的 学 习 ， 先 对 模 糊 图 像 进 行 一 阶 微 分 ， 然 后 进 行 自 相 关运 算 ， 可 得 一 条 鉴 别 曲 线 ， 曲 线 上 会 出 现 对 称 的 相 关 峰 ， 峰 值 为 负 ， 两相 关 峰 之 间 的 距 离 等 于 运 动 模 糊 长 度 。把 模 糊 图 像 转 换 为 灰 度 图 像 ， 采 用 Sobel 算 子 对 其 进 行 一 阶 微 分 运 算 ，Sobel 算 子10obel=2S其 自 相 关 曲 线 如 图 2(a)所 示 。 利 用 matlab 的 Data Cursor 可 以 测 得 两负 峰 之 间 的 距 离 为 20个 像 素 点 ， 此 为 运 动 模糊 的 长 度 ， 与 理 论 值 吻合 ( 比 较 精 确 ) 。 由 于长 度 的 测 量 计 较 精 确 ，可 以 根 据 长 度 对 运 动 模糊 方 向 进 行 不 同 的 取值 ， 达 到 最 好 效 果 。f1=rgb2gray(j1); f1=im2double(f1);h = fspecial( Sobel ); %Sobel 算子J = conv2(f1,h, same ); %Sobel 算子微分IP=abs(fft2(J)); %图像能量谱密度S=fftshift(real(ifft2(IP))); figure,plot(S);title( 自相关图 )