con = np.array(contours) ValueError: setting an array element with a sequence. The requested array has an inhomogeneous shape after 1 dimensions. The detected shape was (84,) + inhomogeneous part.

时间: 2023-07-24 13:00:35 浏览: 70
这个错误说明在将`contours`转换为NumPy数组时出现了问题。出现这个错误的原因是`contours`中的元素具有不同的形状,导致无法创建一个均匀形状的NumPy数组。 要解决这个问题,你可以尝试使用`np.vstack()`函数来垂直堆叠`contours`中的元素,以创建一个具有一致形状的NumPy数组。示例如下: ```python import numpy as np con = np.vstack(contours) con[:, 0] *= factor ``` 通过`np.vstack(contours)`,我们将`contours`中的元素垂直堆叠成一个NumPy数组`con`,然后可以对`con[:, 0]`进行操作。请确保`factor`是一个数值类型的变量,并且与`con[:, 0]`的形状相匹配。
相关问题

triangle = cv2.minEnclosingTriangle(max_contour) triangle = np.int0(triangle)报错:valueerror

这个错误可能是由于`cv2.minEnclosingTriangle()`函数返回值的数据类型不是整数类型,而`cv2.drawContours()`函数需要的轮廓数据类型是整数型的。 为了解决这个错误,我们需要将`cv2.minEnclosingTriangle()`函数返回的三角形数据类型转换为整数型。可以使用`np.int0()`函数将其转换为整数型。 以下是修改后的代码示例: ```python import cv2 import numpy as np # 读取图片 img = cv2.imread('example.jpg') # 转换为灰度图 gray = cv2.cvtColor(img, cv2.COLOR_BGR2GRAY) # 二值化处理 ret, thresh = cv2.threshold(gray, 127, 255, cv2.THRESH_BINARY) # 查找轮廓 contours, hierarchy = cv2.findContours(thresh, cv2.RETR_TREE, cv2.CHAIN_APPROX_SIMPLE) # 找到最大轮廓 max_contour = max(contours, key=cv2.contourArea) # 计算最小外接三角形 triangle = cv2.minEnclosingTriangle(max_contour) triangle = np.int0(triangle) # 绘制最小外接三角形 cv2.drawContours(img, [triangle], 0, (0, 0, 255), 2) # 显示结果 cv2.imshow('image', img) cv2.waitKey(0) cv2.destroyAllWindows() ``` 在此示例中,我们使用`np.int0()`函数将三角形数据类型转换为整数型。

但是运行后报错:terminate called after throwing an instance of 'cv::Exception' what(): OpenCV(4.5.0) /home/sniper/Downloads/opencv-4.5.0/modules/imgproc/src/contours.cpp:1799: error: (-215:Assertion failed) (_contours.kind() == _InputArray::STD_VECTOR_VECTOR || _contours.kind() == _InputArray::STD_VECTOR_MAT || _contours.kind() == _InputArray::STD_VECTOR_UMAT) in function 'findContours'

这个错误通常是由于传递给 `findContours()` 函数的轮廓参数类型不正确导致的。 在你的代码中,`findContours()` 函数的第一个参数是一个 `cv::Mat` 类型的矩阵,而第二个参数 `intersectionPolygon` 是一个 `std::vector<cv::Point>` 类型的向量。这里的问题是,`findContours()` 函数期望的第二个参数是一个用于存储轮廓的输出向量,而你传递了一个点的向量。因此,你需要将 `intersectionPolygon` 声明为一个 `std::vector<std::vector<cv::Point>>` 类型的向量,以存储轮廓。 更改后的代码如下所示: ``` bool isPolygonInside(const std::vector<cv::Point>& polygon1, const std::vector<cv::Point>& polygon2, double& outsideArea) { // Check if all vertices of polygon1 are inside polygon2 bool allInside = true; for (const auto& vertex : polygon1) { double distance = cv::pointPolygonTest(polygon2, vertex, true); if (distance < 0) { allInside = false; break; } } if (allInside) { return true; } // Polygon1 is partially or completely outside polygon2 std::vector<std::vector<cv::Point>> intersectionPolygon; // 修改此处 if (cv::isContourConvex(polygon1) && cv::isContourConvex(polygon2)) { cv::Mat intersectionMat; cv::intersectConvexConvex(cv::Mat(polygon1), cv::Mat(polygon2), intersectionMat); if (cv::countNonZero(intersectionMat) > 0) { cv::findContours(intersectionMat, intersectionPolygon, cv::RETR_EXTERNAL, cv::CHAIN_APPROX_SIMPLE); } } else { std::vector<cv::Point> hull1, hull2; cv::convexHull(polygon1, hull1); cv::convexHull(polygon2, hull2); std::vector<cv::Point> hullIntersection; cv::convexHull(hull1, hullIntersection, false, false); cv::fillConvexPoly(cv::Mat(hullIntersection), hull2, cv::Scalar(0), false); cv::findContours(cv::Mat(hullIntersection), intersectionPolygon, cv::RETR_EXTERNAL, cv::CHAIN_APPROX_SIMPLE); } if (intersectionPolygon.empty()) { outsideArea = 0; return false; } double intersectionArea = std::abs(cv::contourArea(intersectionPolygon[0])); // 修改此处 double polygon1Area = std::abs(cv::contourArea(polygon1)); outsideArea = polygon1Area - intersectionArea; return true; // 修改此处 } ``` 你需要将 `intersectionPolygon` 的类型从 `std::vector<cv::Point>` 改为 `std::vector<std::vector<cv::Point>>`,以及在计算交集面积时使用 `intersectionPolygon[0]` 来获取轮廓。 另外,最后一行需要返回 true,因为在计算出多边形1在多边形2外部的面积后,函数应该返回 true。

相关推荐

def detect_shapes(frame): # 将图像转换为HSV颜色空间 hsv = cv2.cvtColor(frame, cv2.COLOR_BGR2HSV) # 红色范围 lower_red = np.array([0, 100, 100]) upper_red = np.array([10, 255, 255]) red_mask1 = cv2.inRange(hsv, lower_red, upper_red) lower_red = np.array([160, 100, 100]) upper_red = np.array([179, 255, 255]) red_mask2 = cv2.inRange(hsv, lower_red, upper_red) red_mask = red_mask1 + red_mask2 # 蓝色范围 lower_blue = np.array([90, 100, 100]) upper_blue = np.array([130, 255, 255]) blue_mask = cv2.inRange(hsv, lower_blue, upper_blue) # 查找轮廓 contours, _ = cv2.findContours(red_mask + blue_mask, cv2.RETR_EXTERNAL, cv2.CHAIN_APPROX_SIMPLE) for contour in contours: # 计算轮廓的近似形状 epsilon = 0.02 * cv2.arcLength(contour, True) approx = cv2.approxPolyDP(contour, epsilon, True) # 获取轮廓的外接矩形 x, y, w, h = cv2.boundingRect(approx) # 根据轮廓的顶点数和颜色进行分类 if len(approx) == 3: if np.any(red_mask[y:y+h, x:x+w]): shape_label = "Red Triangle" else: shape_label = "Blue Triangle" elif len(approx) == 4: if np.any(red_mask[y:y+h, x:x+w]): shape_label = "Red Square" else: shape_label = "Blue Square" elif len(approx) > 4: if np.any(red_mask[y:y+h, x:x+w]): shape_label = "Red Circle" else: shape_label = "Blue Circle" else: shape_label = "Unknown" # 在图像上绘制边界框和标签https://cdn-static-devbit.csdn.net/ai100/chat/imgs/icon-send-active.png cv2.rectangle(frame, (x, y), (x+w, y+h), (0, 255, 0), 2) cv2.putText(frame, shape_label, (x, y-10), cv2.FONT_HERSHEY_SIMPLEX, 0.9, (0, 255, 0),

def findContours(path,isPlot=False): dcmSOPs = findSOPs(path) #path,rtFile = os.path.split(rvFileName) paths = list(map(str,path.split("\\"))) patient = paths[3] time = paths[4] rvFile = path+'\\RS.{}'.format(patient)+'.CT_{}%.dcm'.format(time) ds = pydicom.dcmread(rvFile) contours = ds.ROIContourSequence dcmFile = path+'\\CT.{}'.format(patient)+'.Image {}.dcm'.format(str(int(1))) ds = pydicom.dcmread(dcmFile) dcmOrigin = ds.ImagePositionPatient dcmSpacing = ds.PixelSpacing # GTV 为第二个轮廓 numberOfContours = len(contours[1].ContourSequence) cuts = [] # 找出包含GTV的CT minXY = 600 maxXY = -1 for k in range(0,numberOfContours): rfContent = contours[1].ContourSequence[k] # 读取该靶区所在CT切片的信息 dcmUID = rfContent.ContourImageSequence[0].ReferencedSOPInstanceUID #print(numberOfContours,len(dcmSOPs),dcmUID) #print(k,dcmSOPs.index(dcmUID)) cuts.append(dcmSOPs.index(dcmUID)) numberOfPoints = rfContent.NumberOfContourPoints # 该层靶区曲线点数 conData = np.zeros((numberOfPoints,3)) # 存储靶区曲线各点的世界坐标 pointData = np.zeros((numberOfPoints,2)) # 存储靶区曲线各点的网格体素坐标 #将靶区勾画的曲线坐标由世界坐标系转换为网格体素坐标 for jj in range(0,numberOfPoints): ii = jj*3 conData[jj,0] = rfContent.ContourData[ii+0] #轮廓世界坐标系 conData[jj,1] = rfContent.ContourData[ii+1] conData[jj,2] = rfContent.ContourData[ii+2] pointData[jj,0] = round( (conData[jj,0] - dcmOrigin[0])/dcmSpacing[0] ) #轮廓X坐标 pointData[jj,1] = round( (conData[jj,1] - dcmOrigin[1])/dcmSpacing[1] ) #轮廓Y坐标 minX = np.min(pointData[:,0]) maxX = np.max(pointData[:,0]) minY = np.min(pointData[:,1]) maxY = np.max(pointData[:,1]) if minXY>minX: minXY = minX elif minXY>minY: minXY = minY elif maxXY<maxX: maxXY = maxX elif maxXY<maxY: maxXY = maxY #print('minXY={}'.format(minXY),'maxXY={}'.format(maxXY)) cuts = np.array(cuts) writeToFile(cuts,path+'\\GTV_indexs.txt') return minXY,maxXY,cuts

def SegMentLWelder(mask_k, mask_box, img_show, contours_info, box, k_class): imgL = cv2.copyTo(mask_k, mask_box) contours_mask_k, _ = cv2.findContours( imgL, cv2.RETR_LIST, cv2.CHAIN_APPROX_NONE) (x, y), radius = cv2.minEnclosingCircle(contours_mask_k[0]) center = [int(x), int(y)] area, trg1 = cv2.minEnclosingTriangle(contours_mask_k[0]) triange_P1 = trg1[0][0] triange_P2 = trg1[1][0] triange_P3 = trg1[2][0] areaL = len(np.where(imgL == 255)[0]) triangle_cnt1 = np.array([triange_P1, center, triange_P2]) mask_triange = np.zeros_like(imgL) mask_triange = cv2.drawContours( mask_triange, [triangle_cnt1.astype(int)], 0, 255, -1) img_triange1 = cv2.copyTo(imgL, mask_triange) triangle_cnt2 = np.array([triange_P2, center, triange_P3]) mask_triange = np.zeros_like(imgL) mask_triange = cv2.drawContours( mask_triange, [triangle_cnt2.astype(int)], 0, 255, -1) img_triange2 = cv2.copyTo(imgL, mask_triange) area1 = len(np.where(img_triange1 == 255)[0]) if area1 < 0.25*areaL or area1 > 0.75*areaL: triangle_cnt1 = np.array([triange_P1, center, triange_P3]) mask_triange = np.zeros_like(imgL) mask_triange = cv2.drawContours( mask_triange, [triangle_cnt1.astype(int)], 0, 255, -1) img_triange1 = cv2.copyTo(imgL, mask_triange) else: area2 = len(np.where(img_triange2 == 255)[0]) if area2 < 0.25*areaL or area2 > 0.75*areaL: triangle_cnt2 = np.array([triange_P1, center, triange_P3]) mask_triange = np.zeros_like(imgL) mask_triange = cv2.drawContours( mask_triange, [triangle_cnt2.astype(int)], 0, 255, -1) img_triange2 = cv2.copyTo(imgL, mask_triange) FindContourSinge(k_class, img_triange1, contours_info, img_show, times=1) FindContourSinge(k_class, img_triange2, contours_info, img_show, times=1)

最新推荐

recommend-type

Dijkstra最短路径算法 - MATLAB.zip

dijkstra算法
recommend-type

文艺高逼格32.pptx

文艺风格ppt模板文艺风格ppt模板 文艺风格ppt模板 文艺风格ppt模板 文艺风格ppt模板 文艺风格ppt模板 文艺风格ppt模板 文艺风格ppt模板 文艺风格ppt模板 文艺风格ppt模板 文艺风格ppt模板 文艺风格ppt模板 文艺风格ppt模板 文艺风格ppt模板 文艺风格ppt模板
recommend-type

计算机基础知识试题与解答

"计算机基础知识试题及答案-(1).doc" 这篇文档包含了计算机基础知识的多项选择题,涵盖了计算机历史、操作系统、计算机分类、电子器件、计算机系统组成、软件类型、计算机语言、运算速度度量单位、数据存储单位、进制转换以及输入/输出设备等多个方面。 1. 世界上第一台电子数字计算机名为ENIAC(电子数字积分计算器),这是计算机发展史上的一个重要里程碑。 2. 操作系统的作用是控制和管理系统资源的使用,它负责管理计算机硬件和软件资源,提供用户界面,使用户能够高效地使用计算机。 3. 个人计算机(PC)属于微型计算机类别,适合个人使用,具有较高的性价比和灵活性。 4. 当前制造计算机普遍采用的电子器件是超大规模集成电路(VLSI),这使得计算机的处理能力和集成度大大提高。 5. 完整的计算机系统由硬件系统和软件系统两部分组成,硬件包括计算机硬件设备,软件则包括系统软件和应用软件。 6. 计算机软件不仅指计算机程序,还包括相关的文档、数据和程序设计语言。 7. 软件系统通常分为系统软件和应用软件,系统软件如操作系统,应用软件则是用户用于特定任务的软件。 8. 机器语言是计算机可以直接执行的语言,不需要编译,因为它直接对应于硬件指令集。 9. 微机的性能主要由CPU决定,CPU的性能指标包括时钟频率、架构、核心数量等。 10. 运算器是计算机中的一个重要组成部分,主要负责进行算术和逻辑运算。 11. MIPS(Millions of Instructions Per Second)是衡量计算机每秒执行指令数的单位,用于描述计算机的运算速度。 12. 计算机存储数据的最小单位是位(比特,bit),是二进制的基本单位。 13. 一个字节由8个二进制位组成,是计算机中表示基本信息的最小单位。 14. 1MB(兆字节)等于1,048,576字节,这是常见的内存和存储容量单位。 15. 八进制数的范围是0-7,因此317是一个可能的八进制数。 16. 与十进制36.875等值的二进制数是100100.111,其中整数部分36转换为二进制为100100,小数部分0.875转换为二进制为0.111。 17. 逻辑运算中,0+1应该等于1,但选项C错误地给出了0+1=0。 18. 磁盘是一种外存储设备,用于长期存储大量数据,既可读也可写。 这些题目旨在帮助学习者巩固和检验计算机基础知识的理解,涵盖的领域广泛,对于初学者或需要复习基础知识的人来说很有价值。
recommend-type

管理建模和仿真的文件

管理Boualem Benatallah引用此版本:布阿利姆·贝纳塔拉。管理建模和仿真。约瑟夫-傅立叶大学-格勒诺布尔第一大学,1996年。法语。NNT:电话:00345357HAL ID:电话:00345357https://theses.hal.science/tel-003453572008年12月9日提交HAL是一个多学科的开放存取档案馆,用于存放和传播科学研究论文,无论它们是否被公开。论文可以来自法国或国外的教学和研究机构,也可以来自公共或私人研究中心。L’archive ouverte pluridisciplinaire
recommend-type

【进阶】音频处理基础:使用Librosa

![【进阶】音频处理基础:使用Librosa](https://picx.zhimg.com/80/v2-a39e5c9bff1d920097341591ca8a2dfe_1440w.webp?source=1def8aca) # 2.1 Librosa库的安装和导入 Librosa库是一个用于音频处理的Python库。要安装Librosa库,请在命令行中输入以下命令: ``` pip install librosa ``` 安装完成后,可以通过以下方式导入Librosa库: ```python import librosa ``` 导入Librosa库后,就可以使用其提供的各种函数
recommend-type

设置ansible 开机自启

Ansible是一个强大的自动化运维工具,它可以用来配置和管理服务器。如果你想要在服务器启动时自动运行Ansible任务,通常会涉及到配置服务或守护进程。以下是使用Ansible设置开机自启的基本步骤: 1. **在主机上安装必要的软件**: 首先确保目标服务器上已经安装了Ansible和SSH(因为Ansible通常是通过SSH执行操作的)。如果需要,可以通过包管理器如apt、yum或zypper安装它们。 2. **编写Ansible playbook**: 创建一个YAML格式的playbook,其中包含`service`模块来管理服务。例如,你可以创建一个名为`setu
recommend-type

计算机基础知识试题与解析

"计算机基础知识试题及答案(二).doc" 这篇文档包含了计算机基础知识的多项选择题,涵盖了操作系统、硬件、数据表示、存储器、程序、病毒、计算机分类、语言等多个方面的知识。 1. 计算机系统由硬件系统和软件系统两部分组成,选项C正确。硬件包括计算机及其外部设备,而软件包括系统软件和应用软件。 2. 十六进制1000转换为十进制是4096,因此选项A正确。十六进制的1000相当于1*16^3 = 4096。 3. ENTER键是回车换行键,用于确认输入或换行,选项B正确。 4. DRAM(Dynamic Random Access Memory)是动态随机存取存储器,选项B正确,它需要周期性刷新来保持数据。 5. Bit是二进制位的简称,是计算机中数据的最小单位,选项A正确。 6. 汉字国标码GB2312-80规定每个汉字用两个字节表示,选项B正确。 7. 微机系统的开机顺序通常是先打开外部设备(如显示器、打印机等),再开启主机,选项D正确。 8. 使用高级语言编写的程序称为源程序,需要经过编译或解释才能执行,选项A正确。 9. 微机病毒是指人为设计的、具有破坏性的小程序,通常通过网络传播,选项D正确。 10. 运算器、控制器及内存的总称是CPU(Central Processing Unit),选项A正确。 11. U盘作为外存储器,断电后存储的信息不会丢失,选项A正确。 12. 财务管理软件属于应用软件,是为特定应用而开发的,选项D正确。 13. 计算机网络的最大好处是实现资源共享,选项C正确。 14. 个人计算机属于微机,选项D正确。 15. 微机唯一能直接识别和处理的语言是机器语言,它是计算机硬件可以直接执行的指令集,选项D正确。 16. 断电会丢失原存信息的存储器是半导体RAM(Random Access Memory),选项A正确。 17. 硬盘连同驱动器是一种外存储器,用于长期存储大量数据,选项B正确。 18. 在内存中,每个基本单位的唯一序号称为地址,选项B正确。 以上是对文档部分内容的详细解释,这些知识对于理解和操作计算机系统至关重要。
recommend-type

"互动学习:行动中的多样性与论文攻读经历"

多样性她- 事实上SCI NCES你的时间表ECOLEDO C Tora SC和NCESPOUR l’Ingén学习互动,互动学习以行动为中心的强化学习学会互动,互动学习,以行动为中心的强化学习计算机科学博士论文于2021年9月28日在Villeneuve d'Asq公开支持马修·瑟林评审团主席法布里斯·勒菲弗尔阿维尼翁大学教授论文指导奥利维尔·皮耶昆谷歌研究教授:智囊团论文联合主任菲利普·普雷教授,大学。里尔/CRISTAL/因里亚报告员奥利维耶·西格德索邦大学报告员卢多维奇·德诺耶教授,Facebook /索邦大学审查员越南圣迈IMT Atlantic高级讲师邀请弗洛里安·斯特鲁布博士,Deepmind对于那些及时看到自己错误的人...3谢谢你首先,我要感谢我的两位博士生导师Olivier和Philippe。奥利维尔,"站在巨人的肩膀上"这句话对你来说完全有意义了。从科学上讲,你知道在这篇论文的(许多)错误中,你是我可以依
recommend-type

【基础】网络编程入门:使用HTTP协议

![【基础】网络编程入门:使用HTTP协议](https://img-blog.csdnimg.cn/direct/4fbc6b5a6d744a519429654f56ea988e.png) # 1. 网络编程基础** **1.1 网络基础知识** 网络是连接计算机和设备的系统,允许它们交换数据和资源。它由节点(计算机、服务器、路由器等)和连接它们的数据链路组成。网络可以是局域网(LAN)、广域网(WAN)或互联网。 **1.2 TCP/IP协议栈** TCP/IP协议栈是一组通信协议,用于在网络上传输数据。它分为四层: * **链路层:**处理物理连接和数据传输。 * **网络层:
recommend-type

时间序列大模型的研究进展

时间序列大模型是近年来自然语言处理领域的一个热门研究方向,它们专注于理解和生成基于时间顺序的数据,比如文本序列、音频或视频信号。这类模型通常结合了Transformer架构(如BERT、GPT等)与循环神经网络(RNNs, LSTM)的长短期记忆能力,以捕捉序列数据中的时间依赖性。 近期的研究进展包括: 1. 长序列建模:研究人员一直在努力提高模型能够处理长序列的能力,例如M6和Turing-NLG,这些模型扩展了序列长度限制,增强了对长期依赖的理解。 2. 结合外部知识:一些模型开始融合外部知识库,如ProphetNet和D-PTM,以提升对复杂时间序列的预测精度。 3. 强化学习和