关于widthStep造成的问题

关于widthStep造成的问题最近遇到一个很奇怪的问题,一直没有解决,就是在A图像中设置一个ROI,将其clone给B,然后对B进行二值化,输入为B,输出为C,这时二值化完后的图像C跟ROI区域的图像区域不同。通过查看发现罪魁祸首是widthStep变了。无意中解决了这个问题,做法如下:方法1:就是在A图像中设置一个ROI,将其clone给B,新建一个C,大小、位数和通道数同B,将Bclone给C,然后对C进行

大家好,又见面了,我是你们的朋友全栈君。

最近遇到一个很奇怪的问题,一直没有解决,就是在A图像中设置一个ROI,将其clone给B,然后对B进行二值化,输入

为B,输出为C,这时二值化完后的图像C跟ROI区域的图像区域不同。通过查看发现罪魁祸首是widthStep变了。

无意中解决了这个问题,做法如下:

方法1:就是在A图像中设置一个ROI,将其clone给B,新建一个C,大小、位数和通道数同B,将B clone给C,然后对C进行

InToIn的二值化操作,这个问题就解决了。

方法2:只拿二值化的阈值,在otsu函数里面不进行二值化操作,然后调用cvThreshold函数进行二值化。

现象1:异常

关于widthStep造成的问题

源代码:

// Otsu.cpp : Defines the entry point for the console application.
//

#include “stdafx.h”
#include <iostream.h>
#include <math.h>
#include <cv.h>
#include <highgui.h>
#include <stdio.h>
#include <windows.h>

int height;int width;
int otsu(IplImage* A, IplImage* B)
{

long N = height * width;
int h[256];
double p[256],u[256],w[256];
for(int i = 0; i < 256; i++)  
{   
h[i] = 0;   
p[i] = 0;   
u[i] = 0;    
w[i] = 0;  
}

for(i = 0; i < height; i++)   
for(int j = 0; j < width; j++)   
for(int k = 0; k < 256; k++)    
{      
if(((uchar*)(A->imageData + A->widthStep*i))[j] == k)    
h[k]++;  
}

for(i = 0; i < 256; i++)   
p[i] = h[i] / double(N);

int T = 0;   
double uT,thegma2fang; 
double thegma2fang_max = -10000;

for(int k = 0; k < 256; k++)  
{    
uT = 0;     
for(i = 0; i <= k; i++) 
{     
u[k] += i*p[i];   
w[k] += p[i];    
}

for(i = 0; i < 256; i++)    
uT += i*p[i];     
thegma2fang = (uT*w[k] – u[k])*(uT*w[k] – u[k]) / (w[k]*(1-w[k]));
if(thegma2fang > thegma2fang_max)  
{     
thegma2fang_max = thegma2fang;    
T = k;  
}  
}

for(i = 0; i < height; i++)  //根据最佳阈值,对原图像进行处理二值化处理    
for(int j = 0; j < width; j++)      
if(((uchar*)(A->imageData + A->widthStep*i))[j] > T)   
((uchar*)(B->imageData + B->widthStep*i))[j] = 255;   
else       
((uchar*)(B->imageData + B->widthStep*i))[j] = 0;

return T;
}

int main(int argc, char** argv)

const char* filename = argc >= 2 ?argv[1] : “lena.jpg”;  
IplImage* source = cvLoadImage( filename, 0 );
int Th;

height = source->height;   
width = source->width;
CvRect rect;
rect.x=50;
rect.y=30;
rect.width=height/2;
rect.height=width/2;
IplImage *roi = cvCreateImage(cvSize(rect.width,rect.height),IPL_DEPTH_8U,1); 
IplImage *bw = cvCreateImage(cvSize(rect.width,rect.height),IPL_DEPTH_8U,1);

cvSetImageROI(source,rect);
roi=cvCloneImage(source);
cvResetImageROI(source);

Th = otsu(roi,bw);

printf(“Th=%d\n”,Th);
cvNamedWindow( “Resource”, 1 );   
cvShowImage( “Resource”, source ); 
cvNamedWindow( “Result”, 1 );  
cvShowImage( “Result”, roi); 
cvNamedWindow( “BW”, 1 );  
cvShowImage( “BW”, bw); 
cvWaitKey(0);
return 0;
}

现象2:正常

关于widthStep造成的问题

// Otsu.cpp : Defines the entry point for the console application.
//

#include “stdafx.h”
#include <iostream.h>
#include <math.h>
#include <cv.h>
#include <highgui.h>
#include <stdio.h>
#include <windows.h>

int height;int width;
int otsu(IplImage* A, IplImage* B)
{

long N = height * width;
int h[256];
double p[256],u[256],w[256];
for(int i = 0; i < 256; i++)  
{   
h[i] = 0;   
p[i] = 0;   
u[i] = 0;    
w[i] = 0;  
}

for(i = 0; i < height; i++)   
for(int j = 0; j < width; j++)   
for(int k = 0; k < 256; k++)    
{      
if(((uchar*)(A->imageData + A->widthStep*i))[j] == k)    
h[k]++;  
}

for(i = 0; i < 256; i++)   
p[i] = h[i] / double(N);

int T = 0;   
double uT,thegma2fang; 
double thegma2fang_max = -10000;

for(int k = 0; k < 256; k++)  
{    
uT = 0;     
for(i = 0; i <= k; i++) 
{     
u[k] += i*p[i];   
w[k] += p[i];    
}

for(i = 0; i < 256; i++)    
uT += i*p[i];     
thegma2fang = (uT*w[k] – u[k])*(uT*w[k] – u[k]) / (w[k]*(1-w[k]));
if(thegma2fang > thegma2fang_max)  
{     
thegma2fang_max = thegma2fang;    
T = k;  
}  
}
return T;
}

int main(int argc, char** argv)

const char* filename = argc >= 2 ?argv[1] : “lena.jpg”;  
IplImage* source = cvLoadImage( filename, 0 );
int Th;

height = source->height;   
width = source->width;
CvRect rect;
rect.x=0;
rect.y=0;
rect.width=height/2;
rect.height=width/2;
IplImage *roi = cvCreateImage(cvSize(rect.width,rect.height),IPL_DEPTH_8U,1); 
IplImage *bw = cvCreateImage(cvSize(rect.width,rect.height),IPL_DEPTH_8U,1);

cvSetImageROI(source,rect);
roi=cvCloneImage(source);
cvResetImageROI(source);

bw = cvCloneImage(roi);
Th = otsu(bw,bw);
cvThreshold(roi,bw,Th,255,CV_THRESH_BINARY);

printf(“roi->widthStep=%d,bw->widthStep=%d\n”,roi->widthStep,bw->widthStep);
printf(“Th=%d\n”,Th);
cvNamedWindow( “Resource”, 1 );   
cvShowImage( “Resource”, source );
cvNamedWindow( “Result”, 1 );  
cvShowImage( “Result”, roi);
cvNamedWindow( “BW”, 1 );  
cvShowImage( “BW”, bw); 
cvWaitKey(0);
return 0;
}

版权声明:本文内容由互联网用户自发贡献,该文观点仅代表作者本人。本站仅提供信息存储空间服务,不拥有所有权,不承担相关法律责任。如发现本站有涉嫌侵权/违法违规的内容, 请发送邮件至 举报,一经查实,本站将立刻删除。

发布者:全栈程序员-用户IM,转载请注明出处:https://javaforall.cn/132961.html原文链接:https://javaforall.cn

【正版授权,激活自己账号】: Jetbrains全家桶Ide使用,1年售后保障,每天仅需1毛

【官方授权 正版激活】: 官方授权 正版激活 支持Jetbrains家族下所有IDE 使用个人JB账号...

(0)
blank

相关推荐

  • 分布式事务TCC方案Hmily——springcloud + feign + mybatis

    分布式事务TCC方案Hmily——springcloud + feign + mybatisTCC理论:分布式事务基础理论——TCCHmily介绍:分布式事务TCC方案——Hmily金融级柔性分布式事务解决方案介绍本文demo代码:GitHub依赖<dependency><groupId>org.dromara</groupId><artifactId>hmily-springcloud</artifactId><vers

  • pycharm如何编译程序_xc8编译器

    pycharm如何编译程序_xc8编译器现在轮到PyCharmPyCharm安装教程PyCharm是一款功能强大的Python编辑器,具有跨平台性,鉴于目前最新版PyCharm使用教程较少,为了节约时间,来介绍一下PyCharm在Windows下是如何安装的。第一步:进入PyCharm官网这是PyCharm的下载地址:http://www.jetbrains.com/pycharm/download/#section=windows进入该网站后,我们会看到如下界面:professional表示专业版,co

  • Ext4文件系统fsck后损坏修复方法-linux数据恢复案例[通俗易懂]

    Ext4文件系统fsck后损坏修复方法-linux数据恢复案例[通俗易懂]在数据恢复案例开始之前有几个概念需要了解块组:Ext4文件系统的全部空间被划分为若干个块组,每个块组内的结构都是大致相同的。块组描述符表:每个块组都对应一个块组描述符,这些块组描述符统一放在文件系统的前部,称为块组描述符表。每个块组描述符大小为32字节,其主要描述块位图、i-节点位图及i-节点表的地址等信息。超级块(Superblock):用于存储文件系统的配置参数(如块大小、总块数、i-节点数…

  • Kali更新源签名无效问题[通俗易懂]

    Kali更新源签名无效问题[通俗易懂]一:Kali更新源签名无效问题1.1:问题展示·更换中科大源报错root@kali:~#apt-getupdate获取:1https://repo.huaweicloud.com/kalikali-rollingInRelease[30.5kB]错误:1https://repo.huaweicloud.com/kalikali-rollingInRelease下列签名无效:EXPKEYSIGED444FF07D8D0BF6KaliLinuxReposito

  • 计算机 修改 虚拟ip,怎么样在电脑中设置虚拟IP地址?

    计算机 修改 虚拟ip,怎么样在电脑中设置虚拟IP地址?满意答案wtc69812020.03.01采纳率:56%等级:9已帮助:114人更改IP地址广域IP:1、如果是PPOE上网只需断开连接再重新连上就好了,服务器会从IP地址池中随机分配一个IP地址给你。2、固定IP上网那你要找运营商更改了,这样改是快不了的。局域IP:网络邻居右键属性,我的连接右键属性,打开tcp/ip更改。改网卡物理地址:1、直接改:可以在桌面上的“网上邻居”图标上单击右…

    2022年10月20日
  • jquery下载教程[通俗易懂]

    jquery下载教程[通俗易懂]jquery下载教程建议使用google浏览器下载,不建议使用IE浏览器(没有办法生成.js文件,只能复制粘贴)第一种方法第二种方法

发表回复

您的电子邮箱地址不会被公开。

关注全栈程序员社区公众号