CN108416770B - Image quality evaluation method based on visual saliency - Google Patents

Image quality evaluation method based on visual saliency Download PDF

Info

Publication number
CN108416770B
CN108416770B CN201810185626.2A CN201810185626A CN108416770B CN 108416770 B CN108416770 B CN 108416770B CN 201810185626 A CN201810185626 A CN 201810185626A CN 108416770 B CN108416770 B CN 108416770B
Authority
CN
China
Prior art keywords
image
dist
ref
gray
gradient
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.)
Active
Application number
CN201810185626.2A
Other languages
Chinese (zh)
Other versions
CN108416770A (en
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.)
Northwestern Polytechnical University
Original Assignee
Northwestern Polytechnical 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 Northwestern Polytechnical University filed Critical Northwestern Polytechnical University
Priority to CN201810185626.2A priority Critical patent/CN108416770B/en
Publication of CN108416770A publication Critical patent/CN108416770A/en
Application granted granted Critical
Publication of CN108416770B publication Critical patent/CN108416770B/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T7/00Image analysis
    • G06T7/0002Inspection of images, e.g. flaw detection
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T2207/00Indexing scheme for image analysis or image enhancement
    • G06T2207/10Image acquisition modality
    • G06T2207/10004Still image; Photographic image

Abstract

The invention provides an image quality evaluation method based on visual saliency, which relates to the technical field of image processing. The method has the advantages that the visual saliency concept is introduced into the image quality evaluation, and the visual saliency interfered in the image is described by analyzing the distribution uniformity of the interference in the image, so that the subjective influence degree of the interference in the image on human vision is reflected, and the image quality is evaluated more accurately.

Description

Image quality evaluation method based on visual saliency
Technical Field
The invention relates to the technical field of image processing, in particular to an image quality evaluation method.
Background
The image quality evaluation is one of basic technologies in the field of image processing, and plays an important role in aspects of algorithm analysis comparison, system performance evaluation and the like, and the full-reference image quality evaluation is a method for evaluating the quality of a distorted image by comparing the difference between the distorted image and an original image, and is widely applied to the fields of evaluating image coding compression, guiding image processing algorithms, monitoring image quality of user terminals and the like.
On the basis of the traditional image quality evaluation method, the subjective perception difference of people is described by edge features by using the perception characteristic of a visual system of human eyes which is extremely sensitive to image edges, and the method only has higher performance in the evaluation of a few compression algorithms and is still inaccurate in subjective perception description of human eyes on the whole and lacks subjective consistency.
Disclosure of Invention
In order to overcome the defects of the prior art, the invention provides the method for comprehensively evaluating the image quality by describing the significance of interference through the density degree of local structure change and combining the local structure change intensity, the change total amount, the change pixel and the contour line relationship so as to improve the evaluation accuracy.
The technical scheme adopted by the invention for solving the technical problem comprises the following specific steps:
step 1, respectively converting an image to be detected and a reference image into a gray image IdistAnd IrefFor two gray scale images IdistAnd IrefRespectively calculating phase consistency PC under frequency domaindistAnd PCref
Step 2, the gray gradient map calculation formula is as follows:
Figure BDA0001590148850000011
wherein G represents a gradient map, dxAnd dyGradient templates in both transverse and longitudinal directions:
Figure BDA0001590148850000012
wherein, a1And a2Satisfies the following conditions:
Figure BDA0001590148850000021
i is a gray image, and the gray image I can be obtained according to the formula (1)distAnd IrefGray scale gradient map GdistAnd Gref
Step 3, calculating two gray gradient images GdistAnd GrefSimilarity of gradients between SGGradient similarity map SGValue S of each pixel point p in (1)G(p) is:
Figure BDA0001590148850000022
Gdist(p) and Gref(p) respectively represent FIG. GdistAnd GrefThe gray value of the middle pixel point p; wherein, T2Is greater than 0, so as to obtain the weighted similarity S of two gray gradient imagesGPCComprises the following steps:
Figure BDA0001590148850000023
wherein, max (PC)dist,PCref) Is PCdistAnd PCrefOf the maximum value of the corresponding element, SGDetermined by the calculation of formula (2);
step 4, a gradient similarity graph S is subjectedGThe distribution uniformity S is obtainedD
First, the gradient similarity map SGDown-sampling, dividing it into n square blocks B with even side lengthiWhere i is 1,2,3 … n, each block BiThe maximum value of the gray value of the middle pixel point is used as the gray value of the pixel point at the corresponding position after down sampling, so that a down sampling image S is obtainedG(max)Of size nx×nyAnd the gray value of each pixel point in the down-sampling graph is px,yAnd x and y indicate that the pixel is at SG(max)The horizontal coordinate and the vertical coordinate in the step (2) are used for calculating the distribution uniformity SD
SD=|S1+S2|+|S3+S4|
Wherein S1,S2,S3,S4Intermediate variables are:
Figure BDA0001590148850000031
Figure BDA0001590148850000032
Figure BDA0001590148850000033
Figure BDA0001590148850000034
and 5: calculating mass fraction
Figure BDA0001590148850000035
Wherein 0 < sigma < 10, SqualityNamely the quality evaluation score of the distorted image relative to the reference image.
The method has the advantages that the visual saliency concept is introduced into the image quality evaluation, and the visual saliency interfered in the image is described by analyzing the distribution uniformity of the interference in the image, so that the subjective influence degree of the interference in the image on human vision is reflected, and the image quality is evaluated more accurately.
Drawings
FIG. 1 is a flow chart of the present invention.
Detailed Description
The invention is further illustrated with reference to the following figures and examples.
The human visual perception system is more sensitive to high-frequency information, so that the significance degree of the change of the spatial structure of image content caused by various disturbance factors is undoubtedly the core factor influencing the subjective perceptibility of image quality. At present, the existing full-reference image quality evaluation method only starts from image difference analysis, and does not pay attention to the concept of significance.
The image quality evaluation method based on visual saliency of the invention has a flow chart as shown in fig. 1, and the specific implementation method comprises the following steps:
step 1: converting the original reference image and the image to be measured obtained by distortion processing into a gray image IdistAnd IrefFourier transform is carried out on the two images and phase consistency PC is calculateddistAnd PCref
Step 2: an image gradient map is calculated. Defining transverse and longitudinal gradient templates dxAnd dyThe following were used:
Figure BDA0001590148850000041
convolution of an image with a gradient template to compute a gradient map
Figure BDA0001590148850000042
When I is taken out of I respectivelydistAnd IrefRespectively obtain GdistAnd Gref
And step 3: calculating a weighted similarity SGPC. First, the gradient similarity S is calculatedGFor each pixel point p in the gradient map, the similarity is obtained
Figure BDA0001590148850000043
Wherein T is taken out2=160。
Then, the weighted similarity S is calculatedGPC
Figure BDA0001590148850000044
Where max (PC)dist,PCref) Is PCdistAnd PCrefThe maximum value of the corresponding element in the matrix.
And 4, step 4: the distribution uniformity S is obtainedD
Firstly, similarity graph SGDown-sampling, continuously dividing the image into a plurality of blocks with the size of 8 multiplied by 8 pixels, taking the maximum value of the gray value of the pixel point in each block as the gray value of the pixel point at the corresponding position in the down-sampling image, and obtaining the down-sampling image SG(max)Of size nx×nyThe gray value of each pixel point in the graph is px,yAnd x and y indicate that the pixel is at SG(max)The horizontal coordinate and the vertical coordinate in the step (2) are used for calculating the distribution uniformity SD
SD=|S1+S2|+|S3+S4|
Wherein S is1,S2,S3,S4Intermediate variables are:
Figure BDA0001590148850000051
Figure BDA0001590148850000052
Figure BDA0001590148850000053
Figure BDA0001590148850000054
and 5: calculating a quality assessment score
Figure BDA0001590148850000055
Wherein σ is 4, SqualityNamely the quality evaluation score of the distorted image relative to the reference image.

Claims (1)

1. An image quality evaluation method based on visual saliency is characterized by comprising the following steps:
step 1, respectively converting an image to be detected and a reference image into a gray image IdistAnd IrefFor two gray scale images IdistAnd IrefRespectively calculating phase consistency PC under frequency domaindistAnd PCref
Step 2, the gray gradient map calculation formula is as follows:
Figure FDA0003010674250000011
wherein G represents a gradient map, dxAnd dyGradient templates in both transverse and longitudinal directions:
Figure FDA0003010674250000012
wherein, a1And a2Satisfies the following conditions:
Figure FDA0003010674250000013
i is a gray image, and the gray image I can be obtained according to the formula (1)distAnd IrefGray scale gradient map GdistAnd Gref
Step 3, calculating two gray gradient images GdistAnd GrefSimilarity of gradients between SGGradient similarity map SGValue S of each pixel point p in (1)G(p) is:
Figure FDA0003010674250000014
Gdist(p) and Gref(p) respectively represent FIG. GdistAnd GrefThe gray value of the middle pixel point p; wherein, T2Is greater than 0, so as to obtain the weighted similarity S of two gray gradient imagesGPCComprises the following steps:
Figure FDA0003010674250000015
wherein, max (PC)dist,PCref) Is PCdistAnd PCrefOf the maximum value of the corresponding element, SGDetermined by the calculation of formula (2);
step 4, a gradient similarity graph S is subjectedGThe distribution uniformity S is obtainedD
First, the gradient similarity map SGDown-sampling, dividing it into n square blocks B with even side lengthiWhere i is 1,2,3 … n, each block BiThe maximum value of the gray value of the middle pixel point is used as the gray value of the pixel point at the corresponding position after down sampling, so that a down sampling image S is obtainedG(max)Of size nx×nyAnd the gray value of each pixel point in the down-sampling graph is px,yAnd x and y indicate that the pixel is at SG(max)The horizontal coordinate and the vertical coordinate in the step (2) are used for calculating the distribution uniformity SD
SD=|S1+S2|+|S3+S4|
Wherein S1,S2,S3,S4Intermediate variables are:
Figure FDA0003010674250000021
Figure FDA0003010674250000022
Figure FDA0003010674250000023
Figure FDA0003010674250000024
and 5: calculating mass fraction
Figure FDA0003010674250000025
Wherein 0 < sigma < 10, SqualityNamely the quality evaluation score of the distorted image relative to the reference image.
CN201810185626.2A 2018-03-07 2018-03-07 Image quality evaluation method based on visual saliency Active CN108416770B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201810185626.2A CN108416770B (en) 2018-03-07 2018-03-07 Image quality evaluation method based on visual saliency

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201810185626.2A CN108416770B (en) 2018-03-07 2018-03-07 Image quality evaluation method based on visual saliency

Publications (2)

Publication Number Publication Date
CN108416770A CN108416770A (en) 2018-08-17
CN108416770B true CN108416770B (en) 2021-06-11

Family

ID=63130303

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201810185626.2A Active CN108416770B (en) 2018-03-07 2018-03-07 Image quality evaluation method based on visual saliency

Country Status (1)

Country Link
CN (1) CN108416770B (en)

Families Citing this family (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN109447980B (en) * 2018-11-12 2021-12-10 公安部第三研究所 Method for implementing image quality evaluation control, computer readable storage medium and processor
CN112991362B (en) * 2021-03-17 2022-11-01 合肥高晶光电科技有限公司 Color sorter adhesion material image segmentation method based on Gaussian mixture model

Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN104091343A (en) * 2014-07-22 2014-10-08 西北工业大学 Image quality evaluation method based on sparse structure
CN105825503A (en) * 2016-03-10 2016-08-03 天津大学 Visual-saliency-based image quality evaluation method
CN106504230A (en) * 2016-10-11 2017-03-15 华侨大学 Complete based on phase equalization refers to color/graphics image quality measure method
CN107578395A (en) * 2017-08-31 2018-01-12 中国地质大学(武汉) The image quality evaluating method that a kind of view-based access control model perceives
CN106023214B (en) * 2016-05-24 2018-11-23 武汉大学 Image quality evaluating method and system based on central fovea view gradient-structure similitude

Patent Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN104091343A (en) * 2014-07-22 2014-10-08 西北工业大学 Image quality evaluation method based on sparse structure
CN105825503A (en) * 2016-03-10 2016-08-03 天津大学 Visual-saliency-based image quality evaluation method
CN106023214B (en) * 2016-05-24 2018-11-23 武汉大学 Image quality evaluating method and system based on central fovea view gradient-structure similitude
CN106504230A (en) * 2016-10-11 2017-03-15 华侨大学 Complete based on phase equalization refers to color/graphics image quality measure method
CN107578395A (en) * 2017-08-31 2018-01-12 中国地质大学(武汉) The image quality evaluating method that a kind of view-based access control model perceives

Non-Patent Citations (4)

* Cited by examiner, † Cited by third party
Title
An Image Quality Assessment Index Based on Visual Saliency and Gradient Amplitude for Telemedicine Application;Ke Pang等;《2017 4th International Conference on Information Science and Control Engineering (ICISCE)》;20171116;第172-176页 *
Full-reference quality assessment of stereoscopic images using disparity-gradient-phase similarity;Feifei Li等;《2015 IEEE China Summit and International Conference on Signal and Information Processing (ChinaSIP)》;20150903;第658-662页 *
基于视觉注意的全参考彩色图像质量评价方法;温阳等;《计算机测量与控制》;20170630;第25卷(第6期);第279-281页 *
基于视觉特征的图像质量评价算法;陈玉坤;《中国优秀硕士学位论文全文数据库信息科技辑》;20131215;第2013年卷(第S2期);I138-1172 *

Also Published As

Publication number Publication date
CN108416770A (en) 2018-08-17

Similar Documents

Publication Publication Date Title
Feichtenhofer et al. A perceptual image sharpness metric based on local edge gradient analysis
Li et al. No-reference image blur assessment based on discrete orthogonal moments
CN108109147B (en) No-reference quality evaluation method for blurred image
CN103353985B (en) A kind of Measurement Method of image Gaussian Blur
CN109523506B (en) Full-reference stereo image quality objective evaluation method based on visual salient image feature enhancement
JP2005527028A (en) Method and system for estimating a sharpness index based on a statistical distribution of local edges
CN116758077B (en) Online detection method and system for surface flatness of surfboard
CN103454282A (en) Gray-gradient-based measurement method of mottle of printed product
CN102036098B (en) Full-reference type image quality evaluation method based on visual information amount difference
US10332244B2 (en) Methods and apparatuses for estimating an ambiguity of an image
CN104036485A (en) Method about image resampling tampering detection
CN108416770B (en) Image quality evaluation method based on visual saliency
Nevriyanto et al. Image enhancement using the image sharpening, contrast enhancement, and Standard Median Filter (Noise Removal) with pixel-based and human visual system-based measurements
CN107895366B (en) Imaging method and system for color evaluation and computer readable storage device
CN104318577A (en) Multi-scale and full-reference image quality evaluation method
CN107944497A (en) Image block method for measuring similarity based on principal component analysis
CN111445435B (en) Multi-block wavelet transform-based reference-free image quality evaluation method
CN112801141B (en) Heterogeneous image matching method based on template matching and twin neural network optimization
Samani et al. Transform domain measure of enhancement—TDME—for security imaging applications
CN116883412B (en) Graphene far infrared electric heating equipment fault detection method
CN116612112A (en) Visual inspection method for surface defects of bucket
CN106611405B (en) Image interpolation method and device
CN111325651B (en) Quantization watermarking method and device based on perceptual JND model
CN111325720B (en) Multi-visual feature integrated reference-free super-resolution image quality evaluation method
Saad et al. Natural DCT statistics approach to no-reference image quality assessment

Legal Events

Date Code Title Description
PB01 Publication
PB01 Publication
SE01 Entry into force of request for substantive examination
SE01 Entry into force of request for substantive examination
GR01 Patent grant
GR01 Patent grant