• 矩阵掩膜操作


    矩阵掩膜(Mask,也被称为kernel):

    用于选定图像、图形或物体,对处理图像(全部或局部)进行遮挡,来控制图像处理区域或处理过程。用于覆盖的特定图像或者物体称为掩膜或模板。管血图像处理中可以做胶片、滤光片等

    数字图像处理中,掩膜为二维矩阵数组,有时也用多值图像。数字图像处理中,图像掩模主要用于:

    • 提取感兴趣,用预先制作的感兴趣去掩膜与待处理图像相乘,得到感兴趣区图像,感兴趣区图内图像数值保持不变,而区外图像值都为0.
    • 屏蔽作用,用掩膜对图像还是那个某些区域作屏障,使其不参与处理或不参加处理参数计算,或仅对屏蔽区域做处理或同级。
    • 结构特征提取,用相似性变量或图像匹配方法检测或提取图像中与掩膜相似的结构特征。
    • 特殊图像的制作
    • 提高对比度

    用途:掩膜是一种图像滤镜的模板,实用掩膜经常处理的是遥感图像。当提取道路或者河流或者房屋时,通过一个n*n的矩阵来对图像进行像素过滤,然后将我们需要的地物或者标志突出显示出来。这个矩阵就是一种掩膜。

       

       

    在OpenCV中,掩模操作是相对简单的。大致的意思是,通过一个掩模矩阵,重新计算图像中的每一个像素值。掩模矩阵控制了旧图像当前位置以及周围位置像素对新图像当前位置像素值的影响力度。用数学术语讲,即我们自定义一个权重表。

    代码:

     1 int cols = (src.cols - 1) * src.channels();
     2 
     3 int offsetx = src.channels();
     4 
     5 int rows = src.rows;        //图像的行数
     6 
     7 dst = Mat::zeros(src.size(), src.type());        //输出容器
     8 
     9 for(int row = 1; row < (rows - 1); row++)
    10 
    11 {
    12 
    13         const uchar* previous = src.ptr<uchar>(row - 1);
    14 
    15         const uchar* current = src.ptr<uchar>(row);
    16 
    17         const uchar* next = src.ptr<uchar>(row + 1);
    18 
    19         uchar* output = dst.ptr<uchar>(row);
    20 
    21         for (int col = offsetx; col < cols; col++)
    22 
    23         {
    24 
    25                 output[col] = saturate_cast<uchar>( 5 * current[col] - (current[col - offsetx] + current[col + offsetx] + previous[col] + next[col]));
    26 
    27         }
    28 
    29 }   

    获取图像的像素指针:

    Mat数据类型指针(Mat::ptr)

    格式:

     

     1 C++: uchar* Mat::ptr(int i=0)
     2 
     3    
     4 
     5 C++: const uchar* Mat::ptr(int i=0) const
     6 
     7    
     8 
     9 C++: template<typename _Tp> _Tp* Mat::ptr(inti=0)
    10 
    11    
    12 
    13 C++: template<typename _Tp> const _Tp*Mat::ptr(int i=0) const

    作用:表示返回指定矩阵行的指针,其中i表示基于0的索引。

    常见用法:

    例一:

     1 cv::Mat image = cv::Mat(400, 600, CV_8UC3); //宽400,长600,3通道彩色图片
     2 
     3 Vec3b* data00 = image.ptr<cv::Vec3b>(0);         //data00是指向image第一行第一个元素指针
     4 
     5 Vec3b* data10 = image.ptr<cv::Vec3b>(1);         //data10是指向image第二行第一个元素指针
     6 
     7 Vec3b* data01 = image.ptr<cv::Vec3b>(0, 1); //data01是指向image第一行第二个元素指针
     8 
     9 data01->val[0] = 0;                //修改data01即第一行第二个元素第一通道值为0
    10 
    11 data01->val[1] = 0;                //修改data01即第一行第二个元素第二通道值为0
    12 
    13 data01->val[2] = 255;         //修改data01即第一行第二个元素第三通道值为255

    结果:

    例二:

     1 Mat.ptr<uchar>(i,j)//代表第i行,第j个点的值(j的大小包含通道数),这是一个地址
     2 
     3 i = Mat.rols();
     4 
     5 j = Mat.cols()*Mat.channels();
     6 
     7 ucahr//代表这个容器存储的类型,和C++的Vector<int> test;一样的
     8 
     9    
    10 
    11 Mat.ptr<uchar>(i) //获取像素矩阵的指针,索引i表示第几行,从0开始计行数。这是一个指针
    12 
    13 const uchar* current= myImage.ptr<uchar>(row);//获得当前行指针
    14 
    15 p(row,col) = current[col]//获取当前像素点P(row, col)的像素值 这是一个值

       

       

    像素范围处理:saturate_cast<uchar>(data);

    功能:确保RGB值的范围在0~255之间。

    1 saturate_cast<uchar>(-100);        //返回值为0
    2 
    3 saturate_cast<uchar>(288);        //返回值为255        
    4 
    5 saturate_cast<uchar>(100);        //返回值为100

    矩阵掩膜函数:filter2D()函数

    格式:

    1 void filter2D( InputArray src,
    2                OutputArray dst,
    3                int ddepth,
    4                InputArray kernel,
    5                Point anchor = Point(-1,-1),
    6                double delta = 0,
    7                int borderType = BORDER_DEFAULT );               

       

    参数:

    参数

    描述

    InputArray 类型的 src

    原图像

    OutputArray 类型的 dst

    目标图像与原图尺寸和通道数相同

    Int 类型的 ddepth

    目标图像的深度

    InputArray 类型的 kernel

    卷积核(或相关核),单通道浮点矩阵;如果将不同的内核应用于不同的通道,请使用拆分将图像拆分为单独的颜色平面,然后进行处理

    Point 类型的 ancho

    内核的锚点,指示内核中过滤点的相对位置;锚应位于内核中;默认值(-1,-1)表示锚位于内核中心

    double 类型的 delta

    在将它们存储在dst中之前,将可选值添加到已过滤的像素中。类似于偏置

       

       

       

       

       

       

    int 类型的 borderType

    表示边缘类型;像素外推法(一般不管它)

    边缘处理类型:

    类型

    描述

    类型

    描述

    BORDER_CONSTANT

    填充边缘用指定像素

    BORDER_DEFAULT

    borderTye默认值,就近的元素进行映射

    BORDER_REPLICATE

    填充边缘用已知的边缘像素值(复制)

    BORDER_REFLECT_101

    对称法:以最边缘像素为轴,对称复制

    BORDER_WRAP

    用另一边的像素值来补偿填充

      

      

    代码

    1 Mat kernel = (Mat_<char>(3, 3) << 0, -1, 0,
    2 
    3                                  -1, 5, -1,
    4 
    5                                  0, -1, 0);//定义掩膜(锐化算子)
    6 
    7 filter2D(src, dst, src.depth(), kernel);

       

    获取当前CPU时间:getTickCount()函数

     1 double t = getTickCount();        //获取当前CPU时间,返回毫秒数
     2 
     3 Mat kernel = (Mat_<char>(3, 3) << 0, -1, 0,
     4 
     5                                          -1, 5, -1,
     6 
     7                                          0, -1, 0);//定义掩膜(锐化算子)
     8 
     9 filter2D(src, dst, src.depth(), kernel);
    10 
    11 double timeconsume = (getTickCount() - t) / getTickFrequency();        //时间差转换为s
    12 
    13 printf("time consume %.2f
    ", timeconsume);

       

    完整代码:

      1 #include<opencv2/opencv.hpp>
      2 
      3 #include<iostream>
      4 
      5 #include<math.h>
      6 
      7    
      8 
      9 using namespace cv;
     10 
     11 using namespace std;
     12 
     13    
     14 
     15 int main(int argc, char** argv)
     16 
     17 {
     18 
     19         Mat src, dst;
     20 
     21         src = imread("./Picture/1.jpg");
     22 
     23         if (!src.data)
     24 
     25         {
     26 
     27                 printf("cloud not load image");
     28 
     29                 return -1;
     30 
     31         }
     32 
     33         namedWindow("input", WINDOW_AUTOSIZE);        // WINDOW_AUTOSIZE 会根据图像大小显示窗口大小,无法修改窗口的小
     34 
     35         imshow("input", src);
     36 
     37    
     38 
     39         int cols = (src.cols - 1) * src.channels();        //图像的列数
     40 
     41         int offsetx = src.channels();        //返回图像src的通道数目
     42 
     43         int rows = src.rows;        //图像的行数
     44 
     45         dst = Mat::zeros(src.size(), src.type());        //创建与src大小、类型相同的零数组作为输出容器
     46 
     47         for(int row = 1; row < (rows - 1); row++)
     48 
     49         {
     50 
     51                 const uchar* previous = src.ptr<uchar>(row - 1);//previous指向src的前一行(row-1)第一个元素的指针
     52 
     53                 const uchar* current = src.ptr<uchar>(row);        //当前行
     54 
     55                 const uchar* next = src.ptr<uchar>(row + 1);//后一行
     56 
     57                 uchar* output = dst.ptr<uchar>(row);
     58 
     59                 for (int col = offsetx; col < cols; col++)
     60 
     61                 {
     62 
     63                         /*
     64 
     65                         其中 saturate_cast<uchar>(data) 的功能是确保data的值在0-255之间
     66 
     67                                 1         saturate_cast<uchar>(-100);        //返回值为0
     68 
     69                                 2         saturate_cast<uchar>(288);        //返回值为255
     70 
     71                                 3         saturate_cast<uchar>(100);        //返回值为100
     72 
     73                         */
     74 
     75                         output[col] = saturate_cast<uchar>( 5 * current[col] - (current[col - offsetx] + current[col + offsetx] + previous[col] + next[col]));
     76 
     77                 }
     78 
     79         }
     80 
     81           
     82 
     83           
     84 
     85         //double t = getTickCount();        //获取当前CPU时间
     86 
     87         //Mat kernel = (Mat_<char>(3, 3) << 0, -1, 0,
     88 
     89         //                                                                 -1, 5, -1,
     90 
     91         //                                                                 0, -1, 0);//定义掩膜(锐化算子)
     92 
     93         //filter2D(src, dst, src.depth(), kernel);
     94 
     95         //double timeconsume = (getTickCount() - t) / getTickFrequency();        //时间差
     96 
     97         //printf("time consume %.2f
    ", timeconsume);
     98 
     99    
    100 
    101    
    102 
    103         namedWindow("contrast image demo", WINDOW_AUTOSIZE);        //WINDOW_NORMAL 可以利用鼠标改变图像大小显示窗口大小
    104 
    105         imshow("contrast image demo", dst);
    106 
    107    
    108 
    109         waitKey(0);
    110 
    111         return 0;
    112 
    113 }
    View Code

       

    运行结果:

       

  • 相关阅读:
    java.math.BigDecimal
    抽象类和接口的区别:
    idea 快捷键
    java 抽象类
    mongodb3.2系统性学习——5、游标 模糊查询 findAndModify函数
    ZeroClipboard2.0 复制功能
    MySQL在远程访问时非常慢的解决skip-name-resolve 并且出现 Reading from net
    mongodb3.2系统性学习——2、write concern mongodb 写安全机制
    mysql 存储引擎 myisam innodb 区别
    MYSQL 错误 :Out of resources when opening file './datagather/mx_domain#P#p178.MYD' (Errcode: 24) 解决办法
  • 原文地址:https://www.cnblogs.com/gghy/p/12296251.html
Copyright © 2020-2023  润新知