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

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

相关文章

利用Python调试串口的示例代码

《利用Python调试串口的示例代码》在嵌入式开发、物联网设备调试过程中,串口通信是最基础的调试手段本文将带你用Python+ttkbootstrap打造一款高颜值、多功能的串口调试助手,需要的可以了... 目录概述:为什么需要专业的串口调试工具项目架构设计1.1 技术栈选型1.2 关键类说明1.3 线程模

Python Transformers库(NLP处理库)案例代码讲解

《PythonTransformers库(NLP处理库)案例代码讲解》本文介绍transformers库的全面讲解,包含基础知识、高级用法、案例代码及学习路径,内容经过组织,适合不同阶段的学习者,对... 目录一、基础知识1. Transformers 库简介2. 安装与环境配置3. 快速上手示例二、核心模

Java的栈与队列实现代码解析

《Java的栈与队列实现代码解析》栈是常见的线性数据结构,栈的特点是以先进后出的形式,后进先出,先进后出,分为栈底和栈顶,栈应用于内存的分配,表达式求值,存储临时的数据和方法的调用等,本文给大家介绍J... 目录栈的概念(Stack)栈的实现代码队列(Queue)模拟实现队列(双链表实现)循环队列(循环数组

使用Python从PPT文档中提取图片和图片信息(如坐标、宽度和高度等)

《使用Python从PPT文档中提取图片和图片信息(如坐标、宽度和高度等)》PPT是一种高效的信息展示工具,广泛应用于教育、商务和设计等多个领域,PPT文档中常常包含丰富的图片内容,这些图片不仅提升了... 目录一、引言二、环境与工具三、python 提取PPT背景图片3.1 提取幻灯片背景图片3.2 提取

Python实现word文档内容智能提取以及合成

《Python实现word文档内容智能提取以及合成》这篇文章主要为大家详细介绍了如何使用Python实现从10个左右的docx文档中抽取内容,再调整语言风格后生成新的文档,感兴趣的小伙伴可以了解一下... 目录核心思路技术路径实现步骤阶段一:准备工作阶段二:内容提取 (python 脚本)阶段三:语言风格调

使用Java将DOCX文档解析为Markdown文档的代码实现

《使用Java将DOCX文档解析为Markdown文档的代码实现》在现代文档处理中,Markdown(MD)因其简洁的语法和良好的可读性,逐渐成为开发者、技术写作者和内容创作者的首选格式,然而,许多文... 目录引言1. 工具和库介绍2. 安装依赖库3. 使用Apache POI解析DOCX文档4. 将解析

一文详解如何在Python中从字符串中提取部分内容

《一文详解如何在Python中从字符串中提取部分内容》:本文主要介绍如何在Python中从字符串中提取部分内容的相关资料,包括使用正则表达式、Pyparsing库、AST(抽象语法树)、字符串操作... 目录前言解决方案方法一:使用正则表达式方法二:使用 Pyparsing方法三:使用 AST方法四:使用字

C++使用printf语句实现进制转换的示例代码

《C++使用printf语句实现进制转换的示例代码》在C语言中,printf函数可以直接实现部分进制转换功能,通过格式说明符(formatspecifier)快速输出不同进制的数值,下面给大家分享C+... 目录一、printf 原生支持的进制转换1. 十进制、八进制、十六进制转换2. 显示进制前缀3. 指

使用Python实现全能手机虚拟键盘的示例代码

《使用Python实现全能手机虚拟键盘的示例代码》在数字化办公时代,你是否遇到过这样的场景:会议室投影电脑突然键盘失灵、躺在沙发上想远程控制书房电脑、或者需要给长辈远程协助操作?今天我要分享的Pyth... 目录一、项目概述:不止于键盘的远程控制方案1.1 创新价值1.2 技术栈全景二、需求实现步骤一、需求

Java中Date、LocalDate、LocalDateTime、LocalTime、时间戳之间的相互转换代码

《Java中Date、LocalDate、LocalDateTime、LocalTime、时间戳之间的相互转换代码》:本文主要介绍Java中日期时间转换的多种方法,包括将Date转换为LocalD... 目录一、Date转LocalDateTime二、Date转LocalDate三、LocalDateTim