CN102682424A - 基于边缘方向差值的图像放大处理方法 - Google Patents
基于边缘方向差值的图像放大处理方法 Download PDFInfo
- Publication number
- CN102682424A CN102682424A CN201210003180XA CN201210003180A CN102682424A CN 102682424 A CN102682424 A CN 102682424A CN 201210003180X A CN201210003180X A CN 201210003180XA CN 201210003180 A CN201210003180 A CN 201210003180A CN 102682424 A CN102682424 A CN 102682424A
- Authority
- CN
- China
- Prior art keywords
- image
- edge direction
- processing method
- coordinate
- pixel
- 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
Links
- 230000003321 amplification Effects 0.000 title claims abstract description 9
- 238000003199 nucleic acid amplification method Methods 0.000 title claims abstract description 9
- 238000003672 processing method Methods 0.000 title abstract description 6
- 238000000034 method Methods 0.000 claims description 30
- 230000015572 biosynthetic process Effects 0.000 claims description 12
- HUTDUHSNJYTCAR-UHFFFAOYSA-N ancymidol Chemical compound C1=CC(OC)=CC=C1C(O)(C=1C=NC=NC=1)C1CC1 HUTDUHSNJYTCAR-UHFFFAOYSA-N 0.000 description 3
- 230000000694 effects Effects 0.000 description 3
- 230000000903 blocking effect Effects 0.000 description 2
- 238000004364 calculation method Methods 0.000 description 2
- 230000001965 increasing effect Effects 0.000 description 2
- 230000009466 transformation Effects 0.000 description 2
- 238000012935 Averaging Methods 0.000 description 1
- 241000220317 Rosa Species 0.000 description 1
- 230000002708 enhancing effect Effects 0.000 description 1
- 238000013178 mathematical model Methods 0.000 description 1
- 238000005070 sampling Methods 0.000 description 1
- 238000013341 scale-up Methods 0.000 description 1
- 230000011218 segmentation Effects 0.000 description 1
- 238000006467 substitution reaction Methods 0.000 description 1
Landscapes
- Image Processing (AREA)
- Editing Of Facsimile Originals (AREA)
Abstract
本发明涉及基于边缘方向插值的图像放大处理方法，将原图像放大倍数k，采用正向扫描对原图像进行扫描，确定将要构造平面所需的不在同一平面上的四顶点在二维平面上的投影坐标；判断四顶点像素值对角线的绝对值差的大小；对角线绝对值差较小的方向属于边缘方向，以对角线绝对值差较大的对角线作为分界线，分别构造平面；根据公式
计算待插像素的颜色值，并填充；本发明的方法能够解决图像放大后边缘出现锯齿或模糊现象，或者运算效率低下的问题，在图像放大时速度快，质量好，具有重大的实用价值。
Description
技术领域
本发明属于数字图像处理技术领域，涉及一种在数字图像放大处理过程中的基于边缘方向的插值方法，主要应用于智能识别中对局部小图像放大后的处理。
背景技术
图像放大是图像处理中一项很重要的基本操作，在互联网、医学图像处理、数字电视等领域以及手机、相机等诸多数码设备上有着广泛的应用。对于图像放大问题，国内外研究者根据不同的数学模型先后提出了多种类型的图像放大方法，如：最近邻插值、双线性插值、高次多项式插值、分形插值以及基于小波变换的方法等。最近邻插值最简单，处理速度最快，但是由于它只是简单地将源像素点最邻近的像素点(最邻近取样)进行复制，其结果是会产生明显的方块效应；而双线性插值法则利用了源像素点周围邻近的4个像素点的线性平均权值来计算目标像素点的值，这样虽然可以得到较好的整体效果，消除了方块效应，但图像的边缘会变得模糊，细节部分不够清晰；高次多项式插值、分形插值以及基于小波变换的方法虽然效果比最近邻插值、双线性插值要好，但运算复杂，且同样会出现边缘模糊等问题。近年来，有些学者对高次多项式、分形插值和基于小波变换的方法做出了一些改进，如夏海宏等提出了基于GPU的双二次B样条图像缩放方法，陈烨等提出了基于分形插值及小波增强的海岸带遥感图像处理方法，肖杰雄等提出了基于小波变换和PDE插值的图像超分辨率重建方法，这些方法在一定程度上都对原有方法有较好的改进，但当放大倍数增加时，运算量也急剧上升。
发明内容
本发明的目的是提供一种基于边缘方向差值的图像放大处理方法，以解决图像放大后边缘出现锯齿或模糊现象，运算效率低下的问题。
为实现上述目的，本发明的基于边缘方向插值的图像放大处理方法步骤如下：
（1）将原图像放大倍数k，采用正向扫描对原图像进行自左向右，自上而下扫描，确定将要构造平面所需的不在同一平面上的四顶点A、B、C、D在二维平面xA’y上的投影坐标；
（2）放大k倍的新图的四顶点在二维平面上的坐标分别是A’(0, 0), B’(k, 0), C’(0, k), D’(k, k)，判断四顶点像素值对角线的绝对值差的大小；
（3）对角线绝对值差较小的方向属于边缘方向，以对角线绝对值差较大的对角线作为分界线，分别构造平面；
（4）构造平面采用待定系数法，求得公式f(x,y)=ax+by+c中的a，b，c，其中，a, b, c 为待定系数，x, y对应像素点的坐标，f(x,y)为像素的颜色值；
（5）根据公式f(x,y)=ax+by+c计算待插像素的颜色值，并填充。
进一步的，所述步骤（4）中是用待定系数法求得公式f(x,y)=ax+by+c中的a，b，c。
本发明的基于边缘方向插值的图像放大处理方法能够解决图像放大后边缘出现锯齿或模糊现象，或者运算效率低下的问题，在图像放大时速度快，质量好，具有重大的实用价值。
附图说明
图1是插值空间平面示意图。
具体实施方式
如果已知空间三个点，则可以构造一个平面，本发明根据此原理进行插值。其数学表达为： f(x,y)=ax+by+c （1）
其中，a, b, c 为待定系数。x, y对应像素点的坐标，f(x,y)为像素的颜色值。
利用待插值像素对应的四个顶点,待插值像素点一定落在原图像某四邻域的网格内，也就是距待插值像素距离最近的原图像的四邻域，（在这里，我们选择计算机的设备坐标系，即水平向右为x轴正方向，垂直向下为y轴正方向，而原点则位于左上角）那么放大k倍（长宽放大倍数相同）后某待插值像素点坐标为（i，j），则它对应的原图像的四顶点（四邻域）坐标 就是（i/k,j/k）,(i/k+1,j/k）,(i/k,j/k+1),(i/k+1,j/k+1，或者这样描述，原图像四邻域坐标分别是（i，j），（i+1，j），（i，j+1），（i+1，j+1），那么放大k倍后的坐标就是（k×i，k×j），（k×（i+1），k×j），（k×i，k×（j+1）），（k×（i+1），k×（j+1）。实际像素构造插值空间平面从而进行插值运算。如图1所示的四个点，能够构造四个平面（只要四个点不在同一个平面上），其中两对相交于一条直线。由于四个像素灰度值的不同，则构造的平面的倾角不同，尤其对应于图像的边缘区域，由于灰度值变化剧烈，那么构造的两个平面夹角很大，而两个平面的交线正好对应图像边缘。对于处于边缘两侧的区域，分别根据构造的平面进行插值。
构造的空间两对平面的两条交线投影到二维图像平面上，就是对应四个实际像素的对角连线。而对角灰度值的差值就反映了图像沿±45。方向的变化率。这个45度是针对原图像的四邻域而言，即使长宽不等比例放大，那么对应的原图像仍然是四邻域，之所以选择四邻域，而不是八邻域，是由于另外四点距离中心比较远（位于斜角方向），对插值贡献不大，而且，如果考虑另外四个点的话，除增加运算量外，对于构造空间平面是没有意义，因为我们只要知道空间三个点坐标和数值，就一定能够唯一的确定一个平面。在图像边缘处，灰度变化较大的方向粗分有三种，即水平，垂直和倾斜。各种图像插值算法对水平和垂直边缘的处理效果基本相同，而对斜角方向的边缘处理则令人不尽满意，但是实际图像的边缘往往不是水平或者垂直的。
基于边缘方向插值的图像放大处理方法，步骤如下：
（1）将原图像放大倍数k，采用正向扫描对原图像进行自左向右，自上而下扫描，确定将要构造平面所需的不在同一平面上的四顶点如图1所示的A、B、C、D在二维平面xA’y上的投影坐标；由于本方法主要用于对局部小图像的放大后处理，放大倍数k为整数倍，一般为2、4。本方法所应用的背景是对图像分辨率的软件提高，采用的图像可以是灰度图像，也可以是RGB图像，如果是RGB图像，则需要对三个通道分别插值计算，原理是一样的。该方法的处理是全局性的，逐点计算，对图像边缘细节区域效果明显，有锐化作用；
（2）放大k倍的新图的四顶点在二维平面上的坐标分别是A’(0, 0), B’(k, 0), C’(0, k), D’(k, k)，判断四顶点像素值对角线的绝对值差的大小；
（3）对角线绝对值差较小的方向属于边缘方向，以对角线绝对值差较大的对角线作为分界线，分别构造平面，图1中，由于fA(0,0) – fD(k,k)的绝对值小于fB(k,0) – fC(0,k)的绝对值，这四个点是原图像四邻域灰度值，比如（x，y）为平面坐标，则f（x，y）就是灰度值，如果小于，则说明AD在空间的斜率要小于BC，也就是对应图像来说，在该处，沿BC方向的灰度变化率要比沿AD方向的变化率大，所以以AD作为构造两平面的交线。因此利用A,B,D三个点和A,C,D三个点分别构造平面，AD为两个平面的交线；
（4）构造平面采用待定系数法，求得公式 中的a，b，c，其中，a, b, c 为待定系数，x, y对应像素点的坐标，f(x,y)为像素的颜色值。因为方程中有三个未知数a，b，c，而在我们已经知道，如图1所示，ABD,ACD这两组分别点A,B,C,D的空间坐标，那么两组方程一定有解，用到的点就是对应原图像四邻域经放大后的坐标值和灰度值,灰度值经放大不变，坐标则随放大倍数而定，而放大倍数是提前已知的。因为每个待插值像素灰度值只与对应原图像的四邻域有关。如下方程组所示：
对任意一个待插值像素而言，它的落点的平面投影一定在图1对应的A’B’C’D’矩形内，我们假设它对应的原图像四邻中AD的空间斜率小于BC，则放大k倍后，一共有k×k个待插值像素点，假设其中一个为f（E），则f（E）坐标我们是知道的（u，v），根据公式f(x,y)=ax+by+c，我们只要求得a，b，c，对应于方程组（2）的第一个方程组，也就是在平面ACD内，ACD三点坐标是（0，0），（0，k），（k，k），灰度值放大后不变，则可以求解系数a1，b1，c1，然后代入公式求取待插值像素的灰度值，这里需要强调的是，由于待插值像素只跟对应的原图像四邻域有关，那么，对应原图像四邻域的坐标的选取可以选择相对坐标，这也就是为什么都是（0，0），（0，k），（k，k），（k，0）。如果选择绝对坐标，计算结果是一样，因为这是个空间线性变换，比如任意一个四邻域坐标为（x1，y1）,（x1，y1+1），（x1+1，y1+1）,（x1+1，y1）,那么放大k倍后坐标就是(k×x1,k×y1),(k×x1,k×（y1+1）),(k×（x1+1）,k×（x1+1），k×（y1+1）)，（k×（x1+1），k×y1）；
（5）根据公式f(x,y)=ax+by+c计算待插像素的颜色值，并填充。
基于边缘方向插值的图像放大处理方法对于实际应用中的数字图像，确定放大倍数后，采用正向扫描，即对原图像进行自左向右，自上而下扫描，确定将要构造平面所需点在xA’y平面上的投影坐标。判断放大k倍的新图四顶点像素值对角线的绝对值差的大小。对角线绝对值差较小的方向属于边缘方向，以对角线绝对值差较大的对角线作为分界线，分别采用待定系数法构造平面，根据公式（1）计算待插像素的颜色值。此方法能够解决了图像放大后边缘出现锯齿或模糊现象，或者运算效率低下的问题，在图像放大时速度快，质量好，具有重大的实用价值。
Claims (2)
1.一种基于边缘方向插值的图像放大处理方法，其特征在于，该方法的步骤如下：
（1）将原图像放大倍数k，采用正向扫描对原图像进行自左向右，自上而下扫描，确定将要构造平面所需的不在同一平面上的四顶点A、B、C、D在二维平面xA’y上的投影坐标；
（2）放大k倍的新图的四顶点在二维平面上的坐标分别是A’(0, 0), B’(k, 0), C’(0, k), D’(k, k)，判断四顶点像素值对角线的绝对值差的大小；
（3）对角线绝对值差较小的方向属于边缘方向，以对角线绝对值差较大的对角线作为分界线，分别构造平面；
（4）构造平面采用待定系数法，求得公式f(x,y)=ax+by+c中的a，b，c，其中，a, b, c 为待定系数，x, y对应像素点的坐标，f(x,y)为像素的颜色值；
（5）根据公式f(x,y)=ax+by+c计算待插像素的颜色值，并填充。
2.根据权利要求1所述的，其特征在于：所述步骤（4）中是用待定系数法求得公式f(x,y)=ax+by+c中的a，b，c。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201210003180.XA CN102682424B (zh) | 2012-01-07 | 2012-01-07 | 基于边缘方向差值的图像放大处理方法 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201210003180.XA CN102682424B (zh) | 2012-01-07 | 2012-01-07 | 基于边缘方向差值的图像放大处理方法 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN102682424A true CN102682424A (zh) | 2012-09-19 |
| CN102682424B CN102682424B (zh) | 2014-12-17 |
Family
ID=46814295
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201210003180.XA Expired - Fee Related CN102682424B (zh) | 2012-01-07 | 2012-01-07 | 基于边缘方向差值的图像放大处理方法 |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN102682424B (zh) |
Cited By (7)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN105160634A (zh) * | 2015-06-30 | 2015-12-16 | 上海交通大学 | 像素图像放大方法及系统 |
| CN105513011A (zh) * | 2014-09-23 | 2016-04-20 | 上海机电工程研究所 | 图像缩放方法 |
| CN105957002A (zh) * | 2016-04-20 | 2016-09-21 | 山东大学 | 基于三角形网格的图像插值放大方法和装置 |
| CN109658332A (zh) * | 2017-10-10 | 2019-04-19 | 凌云光技术集团有限责任公司 | 一种基于流式数据处理的图像放大方法及终端 |
| CN112868045A (zh) * | 2018-10-18 | 2021-05-28 | 索尼公司 | 用于基于ml的放大的帧处置 |
| CN114424268A (zh) * | 2019-09-20 | 2022-04-29 | 亚洲航测株式会社 | 超分辨率立体视觉化处理系统以及其程序 |
| CN115379180A (zh) * | 2021-05-17 | 2022-11-22 | 成都极米科技股份有限公司 | 一种投影画面调整方法、装置、设备及存储介质 |
Citations (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN101452573A (zh) * | 2007-12-04 | 2009-06-10 | 比亚迪股份有限公司 | 一种图像边缘增强方法 |
| US20110026815A1 (en) * | 2009-07-28 | 2011-02-03 | Samsung Electronics Co., Ltd. | Method and apparatus for image processing |
| CN102289799A (zh) * | 2011-08-31 | 2011-12-21 | 洛阳师范学院 | 一种平面插值和球面插值相结合的图像插值器及方法 |
-
2012
- 2012-01-07 CN CN201210003180.XA patent/CN102682424B/zh not_active Expired - Fee Related
Patent Citations (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN101452573A (zh) * | 2007-12-04 | 2009-06-10 | 比亚迪股份有限公司 | 一种图像边缘增强方法 |
| US20110026815A1 (en) * | 2009-07-28 | 2011-02-03 | Samsung Electronics Co., Ltd. | Method and apparatus for image processing |
| CN102289799A (zh) * | 2011-08-31 | 2011-12-21 | 洛阳师范学院 | 一种平面插值和球面插值相结合的图像插值器及方法 |
Non-Patent Citations (3)
| Title |
|---|
| 康牧 等: "一种Roberts自适应边缘检测方法", 《西安交通大学学报》, vol. 42, no. 10, 31 October 2008 (2008-10-31), pages 1240 - 1244 * |
| 李怀琼 等: "基于边缘保护的红外图像插值放大算法", 《兵工学报》, vol. 27, no. 4, 31 July 2006 (2006-07-31), pages 655 - 658 * |
| 龚奕刚 等: "图像放大算法的研究", 《万方学位论文数据库》, 31 March 2009 (2009-03-31), pages 1 - 42 * |
Cited By (13)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN105513011A (zh) * | 2014-09-23 | 2016-04-20 | 上海机电工程研究所 | 图像缩放方法 |
| CN105513011B (zh) * | 2014-09-23 | 2018-10-16 | 上海机电工程研究所 | 图像缩放方法 |
| CN105160634A (zh) * | 2015-06-30 | 2015-12-16 | 上海交通大学 | 像素图像放大方法及系统 |
| CN105160634B (zh) * | 2015-06-30 | 2018-09-21 | 上海交通大学 | 像素图像放大方法及系统 |
| CN105957002A (zh) * | 2016-04-20 | 2016-09-21 | 山东大学 | 基于三角形网格的图像插值放大方法和装置 |
| CN105957002B (zh) * | 2016-04-20 | 2019-03-01 | 山东大学 | 基于三角形网格的图像插值放大方法和装置 |
| CN109658332A (zh) * | 2017-10-10 | 2019-04-19 | 凌云光技术集团有限责任公司 | 一种基于流式数据处理的图像放大方法及终端 |
| CN109658332B (zh) * | 2017-10-10 | 2023-03-28 | 凌云光技术股份有限公司 | 一种基于流式数据处理的图像放大方法及终端 |
| CN112868045A (zh) * | 2018-10-18 | 2021-05-28 | 索尼公司 | 用于基于ml的放大的帧处置 |
| CN114424268A (zh) * | 2019-09-20 | 2022-04-29 | 亚洲航测株式会社 | 超分辨率立体视觉化处理系统以及其程序 |
| CN114424268B (zh) * | 2019-09-20 | 2023-09-29 | 亚洲航测株式会社 | 超分辨率立体视觉化处理系统以及存储介质 |
| CN115379180A (zh) * | 2021-05-17 | 2022-11-22 | 成都极米科技股份有限公司 | 一种投影画面调整方法、装置、设备及存储介质 |
| CN115379180B (zh) * | 2021-05-17 | 2023-09-12 | 成都极米科技股份有限公司 | 一种投影画面调整方法、装置、设备及存储介质 |
Also Published As
| Publication number | Publication date |
|---|---|
| CN102682424B (zh) | 2014-12-17 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN102682424A (zh) | 基于边缘方向差值的图像放大处理方法 | |
| CN104142157B (zh) | 一种标定方法、装置及设备 | |
| CN109272537B (zh) | 一种基于结构光的全景点云配准方法 | |
| CN103745439B (zh) | 图像放大方法以及装置 | |
| US8368714B2 (en) | Curved surface rendering system and method | |
| CN104103058B (zh) | 图像匹配系统及方法 | |
| CN111476242B (zh) | 一种激光点云语义分割方法及装置 | |
| CN103119609B (zh) | 一种确定视频前景主体图像区域的方法和装置 | |
| CN110136114B (zh) | 一种波面高度测量方法、终端设备及存储介质 | |
| CN106408596B (zh) | 基于边缘的局部立体匹配方法 | |
| US20180020205A1 (en) | Stereo matching system and method of operating thereof | |
| CN108074237B (zh) | 图像清晰度检测方法、装置、存储介质及电子设备 | |
| CN104700361B (zh) | 基于边缘检测的图像插值方法及系统 | |
| CN104851127B (zh) | 一种基于交互的建筑物点云模型纹理映射方法及装置 | |
| CN110310365A (zh) | 一种三维重建方法及装置 | |
| GB2492039A (en) | Two dimensional super resolution scaling | |
| CN102447917A (zh) | 立体图像匹配方法及其设备 | |
| CN108986024A (zh) | 一种基于网格的激光点云规则排列处理方法 | |
| CN102339478A (zh) | 一种从等高线图生成数字高程模型的方法和装置 | |
| CN105631939B (zh) | 一种基于曲率滤波的三维点云畸变校正方法及其系统 | |
| JP2009015696A (ja) | 情報処理システム | |
| CN103325109B (zh) | 适用于墙装式全景摄像机的鱼眼图像的畸变校正方法 | |
| CN108182666B (zh) | 一种视差校正方法、装置和终端 | |
| JP2012141758A (ja) | 三次元データ処理装置、方法及びプログラム | |
| CN102663805A (zh) | 一种基于投影的视锥体裁剪的方法 |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| C06 | Publication | ||
| PB01 | Publication | ||
| C10 | Entry into substantive examination | ||
| SE01 | Entry into force of request for substantive examination | ||
| C14 | Grant of patent or utility model | ||
| GR01 | Patent grant | ||
| EE01 | Entry into force of recordation of patent licensing contract | ||
| EE01 | Entry into force of recordation of patent licensing contract |
Application publication date: 20120919 Assignee: HENAN SEMI TECHNOLOGY CO.,LTD. Assignor: Henan University of Science and Technology Contract record no.: X2019980000399 Denomination of invention: Image amplification processing method based on edge direction difference Granted publication date: 20141217 License type: Exclusive License Record date: 20191104 |
|
| 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: 20141217 |