佛像控制点质心提取(附代码)

2024-03-22 10:50

本文主要是介绍佛像控制点质心提取(附代码),希望对大家解决编程问题提供一定的参考价值,需要的开发者们随着小编来一起学习吧!

本文主要包含一些矩阵运算的代码(矩阵转置、求逆、相乘),如果你只对矩阵运算这部分感兴趣的话请直接移步到下面Matrix_function.cpp和Matrix_function.h处,需要注意的是矩阵运算代码是与你的矩阵形式息息相关的,比如3*4的矩阵,有人用二维方式存储,有人用一维方式存储,我这儿用的是二维的。
计算机视觉老师让写空间后方交会求解相机外部参数的代码。具体就是老师会给一张相片,相片上有黑色控制点,这些控制点带有编号,老师会给出控制点对应的世界坐标以及相机内方位元素。然后我们需要自己用PS或者用程序提取出该控制点的像素坐标,在这里我会给出用OpenCV代码提取自动提取控制点像素坐标的方法,但是事实上程序提取的像素坐标不是很准,因为是提取的质心,采用这个像素坐标解算出来的方位元素与老师提供的标准答案是一点点差距,所以事实上我本人是用PS工具(或者选择用VisualStudio 打开图片)刺探出来的,人工刺探的像素坐标很准,解算结果很不错。因为没有啥底子,所以用了摄影测量里空间后方交会的方法,但是本质都是一样的,无伤大雅。

1.数据:
(1)消除畸变的影像2张(整个过程不考虑畸变);
(2)控制点编号方式和控制点坐标;
(3)相机内部参数,相机外部参数的初始值。
实习内容:
(1)在影像上提取控制点并编号:
(2)分别计算每张影像的外部参数(后方交会);
2.平差模型:
本次实习采用共线条件方程:

OpenCV提取质心的方法,是依据这篇文章改的,https://blog.csdn.net/yuyangyg/article/details/70888944
[](https://blog.csdn.net/yuyangyg/article/details/70888944)
,感谢大佬提供的参考!
我改之后的代码如下:```cpp
#include <opencv2/opencv.hpp>
#include <opencv2/imgproc/imgproc.hpp>
#include<iostream>
using namespace cv;
using namespace std;int g_nMedianBlurValue = 2;
Mat element = getStructuringElement(MORPH_RECT, Size(3, 3));
vector<Point> point_seq;Mat selectChannel(Mat src, int channel);
bool objectDetection(Mat  src, int threshold_vlaue, int areasize, int channel)
{int i;Mat displayImage = src.clone();Mat imageSC = selectChannel(src, channel);medianBlur(imageSC, imageSC, g_nMedianBlurValue * 2 + 1);//中值滤波threshold(imageSC, imageSC, threshold_vlaue, 255, CV_THRESH_BINARY);if (1)threshold(imageSC, imageSC, threshold_vlaue, 255, CV_THRESH_BINARY_INV);    //cvNot(imageSC,imageSC);//把元素的每一位取反dilate(imageSC, imageSC, element);//膨胀Scalar color = Scalar(155, 155, 155);vector<vector<Point>> Contours;vector<Vec4i> Hierarchy;findContours(imageSC, Contours, Hierarchy, CV_RETR_TREE, CV_CHAIN_APPROX_SIMPLE, Point(0, 0));vector<Moments> mu(Contours.size());vector<Point2f> mc(Contours.size());Mat drawing = Mat::zeros(src.size(), CV_8UC3);for (int i = 0; i < Contours.size(); i++) {mu[i] = moments(Contours[i], false);}for (int i = 0; i < Contours.size(); i++) {mc[i] = Point2d(mu[i].m10 / mu[i].m00, mu[i].m01 / mu[i].m00);}for (int i = 0; i < Contours.size(); i++) {double tmparea = fabs(contourArea(Contours[i]));if (tmparea > areasize) {drawContours(displayImage, Contours, i, color, 2, 8, Hierarchy, 0, Point());//you can change 1 to CV_FILLEDif (1) {char tam[10000];sprintf(tam, "(%0.0f,%0.0f)", mc[i].x, mc[i].y);cout << "质心x " << mc[i].x << " 质心y  " << mc[i].y  << endl;circle(src, mc[i], 100, cvScalar(255, 255, 255), 2);  //画圆,第三个参数是半径,最后一个参数是线宽putText(src, tam, Point(mc[i].x, mc[i].y), FONT_HERSHEY_COMPLEX, 1.3, Scalar(0, 0, 255), 2, 1); //画出重心的坐标,最后一个参数是线宽}}}namedWindow("控制点坐标", WINDOW_NORMAL);imshow("控制点坐标", src);waitKey(0);return false;
}Mat selectChannel(Mat src, int channel)
{///select image channelMat image, gray, hsv;image = src.clone(); //not directly operate on source imagecvtColor(image, gray, CV_BGR2GRAY);cvtColor(image, hsv, CV_BGR2HSV);vector<Mat> imageRGBORHSV;Mat imageSC;switch (channel){case 1://cvSplit(image,imageSC,0,0,0);split(image, imageRGBORHSV);imageSC = imageRGBORHSV[0];break;case 2://cvSplit(image,0,imageSC,0,0);split(image, imageRGBORHSV);imageSC = imageRGBORHSV[1];break;case 3://cvSplit(image,0,0,imageSC,0);split(image, imageRGBORHSV);imageSC = imageRGBORHSV[2];break;case 4://cvSplit(hsv,imageSC,0,0,0);split(hsv, imageRGBORHSV);imageSC = imageRGBORHSV[0];break;case 5://cvSplit(hsv,0,imageSC,0,0);split(hsv, imageRGBORHSV);imageSC = imageRGBORHSV[1];break;case 6://cvSplit(hsv,0,0,imageSC,0);split(hsv, imageRGBORHSV);imageSC = imageRGBORHSV[2];break;default://cvCopy( gray, imageSC, 0 );imageSC = gray;}return imageSC;
}int main() {//无畸变影像 - 左.bmpMat src = imread("去掉佛像.jpg", 1);//Mat imageSC=selectChannel(src, 1);objectDetection(src, 100, 100, 1);}

如果图像有的部分有干扰,不妨用Photoshop套索去掉

这篇关于佛像控制点质心提取(附代码)的文章就介绍到这儿,希望我们推荐的文章对编程师们有所帮助!



http://www.chinasem.cn/article/835149

相关文章

C#使用SQLite进行大数据量高效处理的代码示例

《C#使用SQLite进行大数据量高效处理的代码示例》在软件开发中,高效处理大数据量是一个常见且具有挑战性的任务,SQLite因其零配置、嵌入式、跨平台的特性,成为许多开发者的首选数据库,本文将深入探... 目录前言准备工作数据实体核心技术批量插入:从乌龟到猎豹的蜕变分页查询:加载百万数据异步处理:拒绝界面

用js控制视频播放进度基本示例代码

《用js控制视频播放进度基本示例代码》写前端的时候,很多的时候是需要支持要网页视频播放的功能,下面这篇文章主要给大家介绍了关于用js控制视频播放进度的相关资料,文中通过代码介绍的非常详细,需要的朋友可... 目录前言html部分:JavaScript部分:注意:总结前言在javascript中控制视频播放

Spring Boot 3.4.3 基于 Spring WebFlux 实现 SSE 功能(代码示例)

《SpringBoot3.4.3基于SpringWebFlux实现SSE功能(代码示例)》SpringBoot3.4.3结合SpringWebFlux实现SSE功能,为实时数据推送提供... 目录1. SSE 简介1.1 什么是 SSE?1.2 SSE 的优点1.3 适用场景2. Spring WebFlu

java之Objects.nonNull用法代码解读

《java之Objects.nonNull用法代码解读》:本文主要介绍java之Objects.nonNull用法代码,具有很好的参考价值,希望对大家有所帮助,如有错误或未考虑完全的地方,望不吝赐... 目录Java之Objects.nonwww.chinasem.cnNull用法代码Objects.nonN

SpringBoot实现MD5加盐算法的示例代码

《SpringBoot实现MD5加盐算法的示例代码》加盐算法是一种用于增强密码安全性的技术,本文主要介绍了SpringBoot实现MD5加盐算法的示例代码,文中通过示例代码介绍的非常详细,对大家的学习... 目录一、什么是加盐算法二、如何实现加盐算法2.1 加盐算法代码实现2.2 注册页面中进行密码加盐2.

python+opencv处理颜色之将目标颜色转换实例代码

《python+opencv处理颜色之将目标颜色转换实例代码》OpenCV是一个的跨平台计算机视觉库,可以运行在Linux、Windows和MacOS操作系统上,:本文主要介绍python+ope... 目录下面是代码+ 效果 + 解释转HSV: 关于颜色总是要转HSV的掩膜再标注总结 目标:将红色的部分滤

在C#中调用Python代码的两种实现方式

《在C#中调用Python代码的两种实现方式》:本文主要介绍在C#中调用Python代码的两种实现方式,具有很好的参考价值,希望对大家有所帮助,如有错误或未考虑完全的地方,望不吝赐教... 目录C#调用python代码的方式1. 使用 Python.NET2. 使用外部进程调用 Python 脚本总结C#调

Java时间轮调度算法的代码实现

《Java时间轮调度算法的代码实现》时间轮是一种高效的定时调度算法,主要用于管理延时任务或周期性任务,它通过一个环形数组(时间轮)和指针来实现,将大量定时任务分摊到固定的时间槽中,极大地降低了时间复杂... 目录1、简述2、时间轮的原理3. 时间轮的实现步骤3.1 定义时间槽3.2 定义时间轮3.3 使用时

Java中&和&&以及|和||的区别、应用场景和代码示例

《Java中&和&&以及|和||的区别、应用场景和代码示例》:本文主要介绍Java中的逻辑运算符&、&&、|和||的区别,包括它们在布尔和整数类型上的应用,文中通过代码介绍的非常详细,需要的朋友可... 目录前言1. & 和 &&代码示例2. | 和 ||代码示例3. 为什么要使用 & 和 | 而不是总是使

Java强制转化示例代码详解

《Java强制转化示例代码详解》:本文主要介绍Java编程语言中的类型转换,包括基本类型之间的强制类型转换和引用类型的强制类型转换,文中通过代码介绍的非常详细,需要的朋友可以参考下... 目录引入基本类型强制转换1.数字之间2.数字字符之间引入引用类型的强制转换总结引入在Java编程语言中,类型转换(无论