关于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)


相关推荐

  • html+css网页开发 之 头部导航条(logo、导航栏、搜索框)

    html+css网页开发 之 头部导航条(logo、导航栏、搜索框)页面布局整体思路:确定页面的版心(可视区),测量可知。 分析页面中的行模块,以及每个行模块中的列模块。 一行中列模块常用浮动布局,先确定每个列的大小,之后确定列的位置。 制作HTML结构。遵循先有结构,后有样式的原则。头部制作1号是版心盒子header1200*42的盒子水平居中对齐 版心盒子内包含2号盒子logo 版心盒子内包含3号盒子nav导航栏 版心盒子内包含4号盒子search搜索框 版心盒子内包含5号盒子user个人信息 注意4个盒子都必须是浮动style.c..

  • response contentType值的问题

    response contentType值的问题response,contentType,UTF-8,ISO-8859-1

  • Adaptive thresholding using the integral image

    Adaptive thresholding using the integral image#include #include #include#includeusingnamespacestd; voidvvAdaptiveThreshold(IplImage*inImg,IplImage*outImg){intS=inImg->width>>5;intT=10;char*input,*bin;

  • lseek函数的用法

    lseek函数的用法所有打开的文件都有一个当前文件偏移量(currentfileoffset),以下简称为cfo。cfo通常是一个非负整数,用于表明文件开始处到文件当前位置的字节数。读写操作通常开始于cfo,并且使cfo增大,增量为读写的字节数。文件被打开时,cfo会被初始化为0,除非使用了O_APPEND 。   使用lseek函数可以改变文件的cfo。       #in…

  • QIIME 2教程. 01简介和安装 Introduction & Install(2020.11)

    QIIME 2教程. 01简介和安装 Introduction & Install(2020.11)QIIME2https://qiime2.org/简介QIIME2是微生物组分析软件QIIME(截止17.7.13被引7771次)的全新版(不是升级版),全部python3全新编写,并于明年全面接替QIIME,是代表末来的分析方法标准(大牛们制定方法标准,我们跟着用就好了)。优点更易于安装:QIIME1的安装让无数生信人竞折腰,现在官方发布了docker,下载即可运行;使用方法多样:支持命令

  • day72Django

    day72Django

发表回复

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

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