CN102722895A - 基于分数傅立叶变换的形状描述方法 - Google Patents

基于分数傅立叶变换的形状描述方法 Download PDF

Info

Publication number
CN102722895A
CN102722895A CN2012101500265A CN201210150026A CN102722895A CN 102722895 A CN102722895 A CN 102722895A CN 2012101500265 A CN2012101500265 A CN 2012101500265A CN 201210150026 A CN201210150026 A CN 201210150026A CN 102722895 A CN102722895 A CN 102722895A
Authority
CN
China
Prior art keywords
image
fourier transform
fractional fourier
matrix
order
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Granted
Application number
CN2012101500265A
Other languages
English (en)
Other versions
CN102722895B (zh
Inventor
郑丽颖
石大明
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Nanhai Innovation And Development Base Of Sanya Harbin Engineering University
Original Assignee
Harbin Engineering University
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by Harbin Engineering University filed Critical Harbin Engineering University
Priority to CN201210150026.5A priority Critical patent/CN102722895B/zh
Publication of CN102722895A publication Critical patent/CN102722895A/zh
Application granted granted Critical
Publication of CN102722895B publication Critical patent/CN102722895B/zh
Expired - Fee Related legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Landscapes

  • Image Analysis (AREA)

Abstract

本发明提供的是一种基于分数傅立叶变换的形状描述方法。其步骤是：第一步，估计图像中目标物体的质心和主轴方向；第二步，将图像旋转和平移；第三步，去除旋转和平移之后的图像中的目标外接矩形之外的部分，得到目标图像，并将目标图像的大小进行归一化；第四步，计算归一化图像的不同阶数的分数傅立叶变换；第五步，舍弃分数傅里叶变换的幅值，仅保留相位；第六步，将归一化图像的各阶分数傅立叶变换的相位直方图作为物体的形状特征。本发明所提出的形状描述方法具有旋转和平移不变性以及利用分数傅里叶变换相位的特点，实现了图像在不同尺度上的滤波，保证了使得所描述的特征既和物体的边界轮廓有关，又与物体的内部区域有关。

Description

基于分数傅立叶变换的形状描述方法
技术领域
本发明涉及一种形状描述方法，特别是涉及一种数字图像中的形状描述方法。
背景技术
目标的形状特征被广泛用于目标识别中，对于物体形状的描述是计算机视觉的重要任务之一。现有的目标形状描述方法可大体上分为两大类：基于边界的形状描述和基于区域的形状描述。其中，基于边界的方法描述目标的形状仅考虑目标边界上的像素点。与之相反，基于区域的方法则从整个目标图像的区域来抽取形状的特征。目前比较常用的基于边界的目标形状描述技术主要有傅里叶描述子、小波描述子、小波-傅里叶描述子等。傅里叶描述子因其计算简单、定义清晰、能量向低频集中的特点，被广泛应用于形状描述领域。傅里叶变换描述子的主要思想是将目标轮廓转化为一维信号，并对此信号进行傅里叶变换，将其从空间域变换到频率域,用图像的频谱特征来描述目标的形状特征。但是传统的傅里叶变换采用全局性的基函数决定其只能用来处理确定的平稳信号，对于时变的非平稳信号则无能为力。
分数傅立叶变换是傅立叶变换的推广形式，它将信号表示为chirp信号的组合。分数傅立叶变换克服了傅立叶变换将时域信息完全摒弃的缺点，通过一种简单的方式实现了信号的从纯粹时间域到纯频率域的全过程的综合描述，能够展示出信号从纯时间域到纯频率域的所有变化特征。目前也有学者尝试用分数傅立叶变换的方法来描述目标的形状特征。张婉婷、王斌提出一种应用分数阶傅里叶变换的形状描述方法，利用目标一维轮廓的分数傅里叶变换的幅度信息来描述目标形状（应用分数阶傅里叶变换的形状描述方法研究.计算机工程与应用.2011.07.20）。
已有的研究工作表明，图像的相位谱中含有比幅值谱更多的信息。随着变换阶数取不同的值，在图像的分数阶相位谱中会包含不同程度的低频和高频信息，对于比较小的变换阶数甚至可以完全从相位谱恢复的图像中看到原图像的大部分信息。
发明内容
本发明的目的在于提供一种综合利用目标轮廓信息和区域信息的基于分数傅立叶变换的形状描述方法。
本发明的目的是这样实现的：
第一步，估计图像中目标物体的质心和主轴方向；第二步，将图像旋转和平移；第三步，去除旋转和平移之后的图像中的目标外接矩形之外的部分，得到目标图像，并将目标图像的大小进行归一化；第四步，计算归一化图像的不同阶数的分数傅立叶变换；第五步，舍弃分数傅里叶变换的幅值，仅保留相位；第六步，将归一化图像的各阶分数傅立叶变换的相位直方图作为物体的形状特征。
本发明的有益效果是：1）所提出的形状描述方法具有旋转和平移不变性；2）利用分数傅里叶变换相位的特点，实现了图像在不同尺度上的滤波，形状特征既和物体的边界轮廓有关，又与物体的内部区域有关。
附图说明
图1是目标的主轴方向。
具体实施方式
下面结合附图举例对本发明做详细的描述：
本发明所提出的形状描述方法具体实现步骤如下：
设I为的输入图像矩阵，I中只包含一个目标，且目标为黑色，背景为白色。
第一步，估计图像中物体的质心和主轴方向，具体为：根据（1）计算目标的质心：
Oc＝[M10/M00,M01/M00]T    (1)
其中，
Figure BDA00001640930100021
变量k和l的取值为0或者1，Mkl表示图像的(k+l)阶矩，Oc为图像中目标的质心，W和H分别为图像的宽度和高度；I(h,w)表示输入图像矩阵I的第h行第w列的元素；
结合图1，根据目标的最小外接矩形估计物体主轴方向；1是目标，2是最小外接矩形，3是主轴；
第二步，根据主轴方向和质心建立新的坐标系，其中坐标原点为质心，主轴方向为x轴，按照新的坐标系将输入图像矩阵I平移和旋转，得到经平移和旋转操作之后的图像矩阵I′；
第三步，将图像的大小进行归一化，具体为：去除图像I′中外接矩形之外的部分，将剩余部分用矩阵L表示，将矩阵L归一化为H′×W′，归一化之后的矩阵用N表示；
第四步，计算归一化图像的不同阶数的分数傅立叶变换，具体为：令变换阶数a=0.2,0.4,…,1.8，利用快速离散分数傅立叶变换算法计算矩阵N的(a,a)阶二维分数傅里叶变换，得到变换结果其中N(a，a)(u,v)为N的(a,a)阶二维份数傅里叶变换，u=1,2,…H′,v=1,2,…,W′，为虚数单位，exp(.)表示指数函数，|N(a，a)(u,v)|为幅值，表示相位；
第五步，舍弃分数傅里叶变换的幅值|N(a，a)(u,v)|，仅保留相位
Figure BDA00001640930100025
第六步，将归一化图像的各阶分数傅立叶变换的相位直方图作为物体的形状特征，具体为：对于每一个变换阶数a，将量化为20个等级，并计算
Figure BDA00001640930100032
的直方图，设
Figure BDA00001640930100033
的直方图为ga(i),其中i=1,2,…,20表示量化等级,a=0.2,0.4,…,1.8表示变换阶数，将直方图{ga(i)|i＝1,..,20;a＝0.2,0.4...,1.8}作为输入图像矩阵I中的目标的形状特征。

Claims (1)

1.一种基于分数傅立叶变换的形状描述方法，其特征在于包括以下步骤：
设I为的输入图像矩阵，I中只包含一个目标，且目标为黑色，背景为白色，
第一步，估计图像中物体的质心和主轴方向，具体为：根据（1）计算目标的质心：
Oc＝[M10/M00,M01/M00]T    (1)
其中，
Figure FDA00001640930000011
变量k和l的取值为0或者1，Mkl表示图像的(k+l)阶矩，Oc为图像中目标的质心，W和H分别为图像的宽度和高度；I(h,w)表示输入图像矩阵I的第h行第w列的元素；
根据目标的最小外接矩形估计物体主轴方向；
第二步，根据主轴方向和质心建立新的坐标系，其中坐标原点为质心，主轴方向为x轴，按照新的坐标系将输入图像矩阵I平移和旋转，得到经平移和旋转操作之后的图像矩阵I′；
第三步，将图像的大小进行归一化，具体为：去除图像I′中外接矩形之外的部分，将剩余部分用矩阵L表示，将矩阵L归一化为H′×W′，归一化之后的矩阵用N表示；
第四步，计算归一化图像的不同阶数的分数傅立叶变换，具体为：令变换阶数a=0.2,0.4,…,1.8，利用快速离散分数傅立叶变换算法计算矩阵N的(a,a)阶二维分数傅里叶变换，得到变换结果
Figure FDA00001640930000012
其中N(a，a)(u,v)为N的(a,a)阶二维分数傅里叶变换，u=1,2,…H′,v=1,2,…,W′，为虚数单位，exp(.)表示指数函数，|N(a，a)(u,v)|为幅值，表示相位；
第五步，舍弃分数傅里叶变换的幅值|N(a，a)(u,v)|，仅保留相位
第六步，将归一化图像的各阶分数傅立叶变换的相位直方图作为物体的形状特征，具体为：对于每一个变换阶数a，将
Figure FDA00001640930000016
量化为20个等级，并计算
Figure FDA00001640930000017
的直方图，设的直方图为ga(i),其中i＝1,2,…,20表示量化等级,a=0.2,0.4,…,1.8表示变换阶数，将直方图{ga(i)|i＝1,..,20;a＝0.2,0.4...,1.8}作为输入图像矩阵I中的目标的形状特征。
CN201210150026.5A 2012-05-15 2012-05-15 基于分数傅立叶变换的形状描述方法 Expired - Fee Related CN102722895B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201210150026.5A CN102722895B (zh) 2012-05-15 2012-05-15 基于分数傅立叶变换的形状描述方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201210150026.5A CN102722895B (zh) 2012-05-15 2012-05-15 基于分数傅立叶变换的形状描述方法

Publications (2)

Publication Number Publication Date
CN102722895A true CN102722895A (zh) 2012-10-10
CN102722895B CN102722895B (zh) 2014-03-26

Family

ID=46948641

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201210150026.5A Expired - Fee Related CN102722895B (zh) 2012-05-15 2012-05-15 基于分数傅立叶变换的形状描述方法

Country Status (1)

Country Link
CN (1) CN102722895B (zh)

Cited By (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN107851088A (zh) * 2015-07-20 2018-03-27 学校法人冲绳科学技术大学院大学学园 用于实时应用的具有同时边缘伪像去除的二维离散傅里叶变换
CN113379777A (zh) * 2021-05-26 2021-09-10 西安电子科技大学 一种基于最小外接矩形垂直内距离比例的形状描述与检索方法

Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101944235A (zh) * 2009-09-18 2011-01-12 哈尔滨工程大学 基于分数傅立叶变换的图像压缩方法
CN102054272A (zh) * 2010-11-03 2011-05-11 浙江大学 一种基于联合分数傅里叶变换的位移矢量探测方法及装置

Patent Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101944235A (zh) * 2009-09-18 2011-01-12 哈尔滨工程大学 基于分数傅立叶变换的图像压缩方法
CN102054272A (zh) * 2010-11-03 2011-05-11 浙江大学 一种基于联合分数傅里叶变换的位移矢量探测方法及装置

Non-Patent Citations (2)

* Cited by examiner, † Cited by third party
Title
LIYING ZHENG,ET AL.: "Comments on the complete generalized fractional Fourier transform", 《OPTICS COMMUNICATIONS》 *
张婉婷等: "应用分数阶傅里叶变换的形状描述方法研究", 《计算机工程与应用》 *

Cited By (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN107851088A (zh) * 2015-07-20 2018-03-27 学校法人冲绳科学技术大学院大学学园 用于实时应用的具有同时边缘伪像去除的二维离散傅里叶变换
US10121233B2 (en) 2015-07-20 2018-11-06 Okinawa Institute Of Science And Technology School Corporation 2D discrete fourier transform with simultaneous edge artifact removal for real-time applications
CN107851088B (zh) * 2015-07-20 2019-03-01 学校法人冲绳科学技术大学院大学学园 对对象图像数据执行离散傅立叶变换的方法和系统
CN113379777A (zh) * 2021-05-26 2021-09-10 西安电子科技大学 一种基于最小外接矩形垂直内距离比例的形状描述与检索方法

Also Published As

Publication number Publication date
CN102722895B (zh) 2014-03-26

Similar Documents

Publication Publication Date Title
CN108122213B (zh) 一种基于YCrCb的低对比度图像增强方法
CN104200437A (zh) 一种图像去雾处理方法
CN104574293A (zh) 基于有界运算的多尺度Retinex图像清晰化算法
CN110472521B (zh) 一种瞳孔定位校准方法及系统
CN104182941A (zh) 高光谱图像条带噪声去噪方法
CN112308872B (zh) 基于多尺度Gabor一阶导数的图像边缘检测方法
CN107784639B (zh) 一种无人机遥感图像改进多边滤波去噪方法
CN103559697A (zh) 基于fft的碎纸片纵切拼接复原算法
CN116468958B (zh) 通信铁塔安全检测方法及系统
CN104537624B (zh) 基于ssim校正聚类稀疏表示的sar图像降斑方法
Gongwen et al. On medical image segmentation based on wavelet transform
CN111091107A (zh) 一种人脸区域边缘检测方法、装置及存储介质
CN108932492A (zh) 一种基于非采样剪切波变换的图像指纹提取方法
CN102722895A (zh) 基于分数傅立叶变换的形状描述方法
CN104978720A (zh) 一种视频图像雨滴去除方法及装置
Mustaghfirin et al. The comparison of iris detection using histogram equalization and adaptive histogram equalization methods
CN104751470A (zh) 一种图像快速匹配方法
Thriveni Satellite image enhancement using discrete wavelet transform and threshold decomposition driven morphological filter
Luo et al. Monogenic signal theory based feature similarity index for image quality assessment
Huang et al. An improved retinex algorithm for underwater image enhancement based on HSV model
CN111652806B (zh) 一种图像去阴影的方法和系统
CN112927169A (zh) 一种基于小波变换和改进的加权核范数最小化的遥感影像去噪方法
Lei et al. Cement pavement surface crack detection based on image processing
CN112101213A (zh) 一种获取指纹方向信息的方法
Mahajan et al. Artefact removal and contrast enhancement for dermoscopic images using image processing techniques

Legal Events

Date Code Title Description
C06 Publication
PB01 Publication
C10 Entry into substantive examination
SE01 Entry into force of request for substantive examination
GR01 Patent grant
GR01 Patent grant
TR01 Transfer of patent right
TR01 Transfer of patent right

Effective date of registration: 20201202

Address after: Area A129, 4th floor, building 4, Baitai Industrial Park, Yazhou Bay science and Technology City, Yazhou District, Sanya City, Hainan Province, 572024

Patentee after: Nanhai innovation and development base of Sanya Harbin Engineering University

Address before: 150001 Heilongjiang, Nangang District, Nantong street,, Harbin Engineering University, Department of Intellectual Property Office

Patentee before: HARBIN ENGINEERING University

CF01 Termination of patent right due to non-payment of annual fee
CF01 Termination of patent right due to non-payment of annual fee

Granted publication date: 20140326

Termination date: 20210515