opencv vec3f 访问的错误

Mat sample(sample_num,1,CV_32FC3);
		for(int j=0;j<sample_num;++j){
			int clusterIdx=labels.at<int>(j);
			Mat temp_cov(sample.row(j)-curr[clusterIdx].center);
			temp_cov=temp_cov.t()*temp_cov;
			temp_cov.convertTo(temp_cov,CV_64FC1);
			curr[clusterIdx].cov+=temp_cov;
			++curr[clusterIdx].mass;
		}


结果就是算出来的covariance非正定,导致马氏距离出现了负值
更改后的结果:


for(int j=0;j<sample_num;++j){
			int clusterIdx=labels.at<int>(j);
			Vec3f xxx=sample.at<Vec3f>(j,0);
			float ppp[3]={xxx[0],xxx[1],xxx[2]};
			Mat cur_pix(1,3,CV_32FC1,ppp);
			Mat temp_cov(cur_pix-curr[clusterIdx].center);
		//	cout<<temp_cov<<endl;
		//	cout<<cur_pix-curr[clusterIdx].center<<endl;
			temp_cov=temp_cov.t()*temp_cov;
			curr[clusterIdx].cov+=temp_cov;
			++curr[clusterIdx].mass;
		}


这种错误太难发现,还是通过和matlab的运算结果对比才知道哪里出了问题

评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值