python在5000*4000二维矩阵中求20*20二维移动窗口内极大值的坐标

时间: 2023-04-02 11:04:54 浏览: 50
可以使用numpy库中的函数实现,具体代码如下: import numpy as np # 生成一个500*400的随机矩阵 matrix = np.random.rand(500, 400) # 定义移动窗口的大小 window_size = 20 # 使用maximum_filter函数求出每个窗口内的最大值 max_matrix = np.maximum_filter(matrix, size=(window_size, window_size)) # 使用argwhere函数找到最大值的坐标 max_coords = np.argwhere(matrix == max_matrix) # 输出第一个最大值的坐标 print(max_coords[])
相关问题

python在5000*4000二维矩阵中求每个20*20二维移动窗口内极大值的坐标

可以使用numpy库中的stride_tricks函数来实现这个功能。具体步骤如下: 1. 将二维矩阵转换为一维数组,方便后续操作。 2. 使用stride_tricks函数将一维数组转换为二维数组,每个子数组的大小为20*20。 3. 对每个子数组求最大值的坐标。 4. 将每个子数组的最大值坐标转换为在原始二维矩阵中的坐标。 代码示例: import numpy as np # 生成一个500*400的随机二维矩阵 matrix = np.random.rand(500, 400) # 定义移动窗口大小 window_size = 20 # 将二维矩阵转换为一维数组 flat_matrix = matrix.ravel() # 使用stride_tricks函数将一维数组转换为二维数组 shape = (matrix.shape[]-window_size+1, matrix.shape[1]-window_size+1, window_size, window_size) strides = (matrix.strides[], matrix.strides[1], matrix.strides[], matrix.strides[1]) sub_matrices = np.lib.stride_tricks.as_strided(flat_matrix, shape=shape, strides=strides) # 对每个子数组求最大值的坐标 max_indices = np.unravel_index(np.argmax(sub_matrices, axis=(2,3)), shape[:2]) # 将每个子数组的最大值坐标转换为在原始二维矩阵中的坐标 max_indices = (max_indices[]+window_size//2, max_indices[1]+window_size//2) print(max_indices)

python二维矩阵求二维移动窗口内极大值的坐标

可以使用numpy库中的函数numpy.lib.stride_tricks.as_strided来实现二维移动窗口的操作,具体实现可以参考以下代码: import numpy as np def max_in_window(matrix, window_size): # 计算窗口的步长 stride = matrix.strides stride = (stride[], stride[1], stride[], stride[1]) # 利用as_strided函数生成移动窗口的视图 window_view = np.lib.stride_tricks.as_strided(matrix, shape=( matrix.shape[] - window_size + 1, matrix.shape[1] - window_size + 1, window_size, window_size), strides=stride) # 计算每个窗口内的最大值的坐标 max_index = np.unravel_index(np.argmax(window_view, axis=(2, 3)), (window_view.shape[], window_view.shape[1])) # 将坐标转换为在原矩阵中的位置 max_index = np.array(max_index).T max_index[:, ] += window_size // 2 max_index[:, 1] += window_size // 2 return max_index # 测试代码 matrix = np.random.randint(, 10, size=(5, 5)) print(matrix) print(max_in_window(matrix, 3))

相关推荐

最新推荐

recommend-type

Python二维数组实现求出3*3矩阵对角线元素的和示例

今天小编就为大家分享一篇Python二维数组实现求出3*3矩阵对角线元素的和示例,具有很好的参考价值,希望对大家有所帮助。一起跟随小编过来看看吧
recommend-type

Python:二维列表下标互换方式(矩阵转置)

今天小编就为大家分享一篇Python:二维列表下标互换方式(矩阵转置),具有很好的参考价值,希望对大家有所帮助。一起跟随小编过来看看吧
recommend-type

python中字符串变二维数组的实例讲解

下面小编就为大家分享一篇python中字符串变二维数组的实例讲解,具有很好的参考价值,希望对大家有所帮助。一起跟随小编过来看看吧
recommend-type

python保存二维数组到txt文件中的方法

今天小编就为大家分享一篇python保存二维数组到txt文件中的方法,具有很好的参考价值,希望对大家有所帮助。一起跟随小编过来看看吧
recommend-type

python NumPy ndarray二维数组 按照行列求平均实例

今天小编就为大家分享一篇python NumPy ndarray二维数组 按照行列求平均实例,具有很好的参考价值,希望对大家有所帮助。一起跟随小编过来看看吧
recommend-type

zigbee-cluster-library-specification

最新的zigbee-cluster-library-specification说明文档。
recommend-type

管理建模和仿真的文件

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

【实战演练】MATLAB用遗传算法改进粒子群GA-PSO算法

![MATLAB智能算法合集](https://static.fuxi.netease.com/fuxi-official/web/20221101/83f465753fd49c41536a5640367d4340.jpg) # 2.1 遗传算法的原理和实现 遗传算法(GA)是一种受生物进化过程启发的优化算法。它通过模拟自然选择和遗传机制来搜索最优解。 **2.1.1 遗传算法的编码和解码** 编码是将问题空间中的解表示为二进制字符串或其他数据结构的过程。解码是将编码的解转换为问题空间中的实际解的过程。常见的编码方法包括二进制编码、实数编码和树形编码。 **2.1.2 遗传算法的交叉和
recommend-type

openstack的20种接口有哪些

以下是OpenStack的20种API接口: 1. Identity (Keystone) API 2. Compute (Nova) API 3. Networking (Neutron) API 4. Block Storage (Cinder) API 5. Object Storage (Swift) API 6. Image (Glance) API 7. Telemetry (Ceilometer) API 8. Orchestration (Heat) API 9. Database (Trove) API 10. Bare Metal (Ironic) API 11. DNS
recommend-type

JSBSim Reference Manual

JSBSim参考手册,其中包含JSBSim简介,JSBSim配置文件xml的编写语法,编程手册以及一些应用实例等。其中有部分内容还没有写完,估计有生之年很难看到完整版了,但是内容还是很有参考价值的。