- iOS/Objective-C 元类和类别
- objective-c - -1001 错误,当 NSURLSession 通过 httpproxy 和/etc/hosts
- java - 使用网络类获取 url 地址
- ios - 推送通知中不播放声音
我正在将一些代码从 matlab 转换为 opencv。我尝试在 opencv 中使用 Sobel,但 opencv 和 matlab 的输出完全不同,这可能是什么原因。如何使 opencv 的输出与 matlab 相同?我的 MATLAB 代码是:
[sobel_edges,T,V,H] = edge(rgb2gray(im),'sobel',0.03);
sobel_angles = atan2(V,H);
sobel_weights = (V.*V+H.*H).^0.5;
其中 0.03 是阈值。在 opencv 中,当我使用预建的 Sobel 滤波器时,输出与 matlab 完全不同,甚至在 openc 中计算的 engle 和 magnitude 也不同。opencv代码为:
Mat gray_img=Mat::zeros(img.size(),CV_8U);
Mat gradientX=Mat::zeros(gray_img.size(),CV_64F);
Mat gradientY=Mat::zeros(gray_img.size(),CV_64F);
Mat sobel_edge=Mat::zeros(gray_img.size(),CV_64F);
cvtColor(img, gray_img, CV_BGR2GRAY);
Sobel(gray_img, gradientX, gradientX.type(), 1, 0, 3);
Sobel(gray_img, gradientY, gradientY.type(), 0, 1, 3);
Sobel(gray_img,sobel_edge,sobel_edge.type(),1,1,3);
sobel_edge.convertTo(sobel_edge,CV_8U);
sobel_edge.convertTo(sobel_edge,CV_64F);
sobel_edge=sobel_edge/255.0; //I divided this my 255 becuz in MATLAB the output is between 0 to 1
imshow("Sobel",sobel_edge);
Mat magnitude(gray_img.size(), CV_64F, cv::Scalar(0.0));
Mat angles=Mat::zeros(gradientX.size(),CV_64F);
bool anglesInDegrees = true;
cartToPolar(gradientX, gradientY, magnitude, angles, anglesInDegrees);
sobel 边缘本身不同,大小和角度也不同,我试图通过查看 matlab 中的边缘函数在 opencv 中手动转换 sobel 但输出仍然不同,因为它证明了 opencv 和 imfilter 的 filter2D在 matlab 中返回不同的输出。如何在matlab和opencv中获得相同的sobel输出???手动将matlab的sobel转换为opencv的代码是:
Mat gray_img=Mat::zeros(img.size(),CV_32FC1);
cvtColor(img, gray_img, CV_RGB2GRAY);
double minVal,maxVal;
cv::Mat gray = cv::Mat(gray_img.size(),CV_32FC1);
gray_img.convertTo(gray_img,CV_32FC1);
gray=gray_img/255.0;
cout<<gray<<endl<<"End";
double data[]={1,2,1,0,0,0,-1,2,-1};
Mat op=Mat(3,3,CV_64F,data).clone();
op=op/8;
Mat x_mask;
transpose(op,x_mask);
cout<<x_mask<<endl;
Mat y_mask=op.clone();
int scale=4;
int offset[]={0,0,0,0};
double sobel_thresh=0.03;
Mat bx,by,bx_mul,by_mul,b;
Point anchor(0,0);
float delta = 0.0;
cv::filter2D(gray, bx, CV_32FC1, x_mask, anchor, delta, BORDER_REPLICATE);
bx=abs(bx);
imshow("f1",bx);
cv::filter2D(gray, by, CV_32FC1, y_mask, anchor, delta, BORDER_REPLICATE);
by=abs(by);
imshow("by",by);
pow(bx,2,bx_mul);
imshow("f2",bx_mul);
pow(by,2,by_mul);
b= bx_mul+by_mul;
imshow("f3",b);
double cut_off;
cut_off=pow(sobel_thresh,2);
Mat sobel_edge(gray.size(),CV_32FC1);
for(int i=0;i<b.rows;i++)
{
for(int j=0;j<b.cols;j++)
{
if((b.at<float>(i,j))>cut_off)
{
sobel_edge.at<float>(i,j)=1;
}
else
{
sobel_edge.at<float>(i,j)=0;
}
}
}
imshow("Sobel_edge",sobel_edge);
最佳答案
这段代码给出了与 MATLAB 代码相同的结果:
int main(int argc, char* argv[])
{
namedWindow("result");
Mat img=imread("D:\\ImagesForTest\\1.tiff",0);
img.convertTo(img,CV_32FC1,1.0/255.0);
Mat h,v,g;
cv::Sobel(img,h,-1,1,0,3,1.0/8.0);
cv::Sobel(img,v,-1,0,1,3,1.0/8.0);
cv::magnitude(h,v,g);
// Check extremums
double m,M;
cv::minMaxLoc(g,&m,&M);
cout << m << ":" << M << endl;
cv::minMaxLoc(h,&m,&M);
cout << m << ":" << M << endl;
cv::minMaxLoc(v,&m,&M);
cout << m << ":" << M << endl;
imshow("result",g);
cv::waitKey(0);
}
OpenCV 从不缩放卷积结果,所以要小心。
关于c++ - opencv 和 Matlab 的 Sobel 滤波器输出不同,我们在Stack Overflow上找到一个类似的问题: https://stackoverflow.com/questions/18613377/
我正在使用 OUTFILE 命令,但由于权限问题和安全风险,我想将 shell 的输出转储到文件中,但出现了一些错误。我试过的 #This is a simple shell to connect t
我刚刚开始学习 Java,我想克服在尝试为这个“问题”创建 Java 程序时出现的障碍。这是我必须创建一个程序来解决的问题: Tandy 喜欢分发糖果,但只有 n 颗糖果。对于她给第 i 个糖果的人,
你好,我想知道我是否可以得到一些帮助来解决我在 C++ 中打印出 vector 内容的问题 我试图以特定顺序在一个或两个函数调用中输出一个类的所有变量。但是我在遍历 vector 时收到一个奇怪的错误
我正在将 intellij (2019.1.1) 用于 java gradle (5.4.1) 项目,并使用 lombok (1.18.6) 来自动生成代码。 Intellij 将生成的源放在 out
编辑:在与 guest271314 交流后,我意识到问题的措辞(在我的问题正文中)可能具有误导性。我保留了旧版本并更好地改写了新版本 背景: 从远程服务器获取 JSON 时,响应 header 包含一
我的问题可能有点令人困惑。我遇到的问题是我正在使用来自 Java 的 StoredProcedureCall 调用过程,例如: StoredProcedureCall call = new Store
在我使用的一些IDL中,我注意到在方法中标记返回值有2个约定-[in, out]和[out, retval]。 当存在多个返回值时,似乎使用了[in, out],例如: HRESULT MyMetho
当我查看 gar -h 的帮助输出时,它告诉我: [...] gar: supported targets: elf64-x86-64 elf32-i386 a.out-i386-linux [...
我想循环遍历一个列表,并以 HTML 格式打印其中的一部分,以代码格式打印其中的一部分。所以更准确地说:我想产生与这相同的输出 1 is a great number 2 is a great
我有下面的tekton管道,并尝试在Google Cloud上运行。集群角色绑定。集群角色。该服务帐户具有以下权限。。例外。不确定需要为服务帐户设置什么权限。
当尝试从 make 过滤非常长的输出以获取特定警告或错误消息时,第一个想法是这样的: $ make | grep -i 'warning: someone set up us the bomb' 然而
我正在创建一个抽象工具类,该类对另一组外部类(不受我控制)进行操作。外部类在某些接口(interface)点概念上相似,但访问它们相似属性的语法不同。它们还具有不同的语法来应用工具操作的结果。我创建了
这个问题已经有答案了: What do numbers starting with 0 mean in python? (9 个回答) 已关闭 7 年前。 在我的代码中使用按位与运算符 (&) 时,我
我写了这段代码来解析输入文件中的行输入格式:电影 ID 可以有多个条目,所以我们应该计算平均值输出:**没有重复(这是问题所在) import re f = open("ratings2.txt",
我需要处理超过 1000 万个光谱数据集。数据结构如下:大约有 1000 个 .fits(.fits 是某种数据存储格式)文件,每个文件包含大约 600-1000 个光谱,其中每个光谱中有大约 450
我编写了一个简单的 C 程序,它读取一个文件并生成一个包含每个单词及其出现频率的表格。 该程序有效,我已经能够在 Linux 上运行的终端中获得显示的输出,但是,我不确定如何获得生成的显示以生成包含词
很难说出这里要问什么。这个问题模棱两可、含糊不清、不完整、过于宽泛或夸夸其谈,无法以目前的形式得到合理的回答。如需帮助澄清此问题以便重新打开,visit the help center . 关闭 1
1.普通的输出: print(str)#str是任意一个字符串,数字··· 2.格式化输出: ?
我无法让 logstash 正常工作。 Basic logstash Example作品。但后来我与 Advanced Pipeline Example 作斗争.也许这也可能是 Elasticsear
这是我想要做的: 我想让用户给我的程序一些声音数据(通过麦克风输入),然后保持 250 毫秒,然后通过扬声器输出。 我已经使用 Java Sound API 做到了这一点。问题是它有点慢。从发出声音到
我是一名优秀的程序员,十分优秀!