OpenRadar DOA函数 Bartlett/CBF和Capon使用

2023-10-07 08:28

本文主要是介绍OpenRadar DOA函数 Bartlett/CBF和Capon使用,希望对大家解决编程问题提供一定的参考价值,需要的开发者们随着小编来一起学习吧!

Bartlett / CBF原理看这里
Capon原理看这里

这里只讲怎么调用openradar提供的aoa_bartlett和aoa_capon函数:
openradar的作者代码水平很高,但里面有很多匪夷所思的写法,比如他demo里的解析文件格式就很迷啊等等…

先获取矢量矩阵 (steering vector):

比如果你想得到[-90°,90°],分辨率为1°(就有181个角度点),8根虚拟天线的矢量矩阵,那么就会生成一个shape为[181, 8]的ndarray (也就是下面那行代码的steering_vec)。num_vec得到的是一共有多少角度(上面那个例子的话就是181),没啥用,可以不管。

num_vec, steering_vec = dsp.gen_steering_vec(ANGLE_RANGE, ANGLE_RES, VIRT_ANT)

By the way, dsp.gen_steering_vec里面用了np.complex, 高版本的numpy可能不支持,报错说和python自带的complex有歧义,把np去掉,变成complex即可。

ANGLE_RANGE 是在一侧,你要得到多少角度信息,那么一共就有 2 * ANGLE_RANGE + 1个角度信息,+1是因为0°。
ANGLE_RES 是你需要的角度分辨率
VIRT_ANT 是虚拟天线数
以下举例

ANGLE_RES = 1
ANGLE_RANGE = 90
BINS_PROCESSED = 112

意思想得到的角度为[-90°, 90°],因为有0,所以有2*ANGLE_RANGE + 1 = 181个角度格子,想要的角度分辨率为1°,处理前112个range-bin。

Bartlett/CBF:

radar_cube shape : (128, 8, 128) - (chirp, vx, adc_samples)

ra = np.zeros((ANGLE_BINS, BINS_PROCESSED), dtype=np.complex128)
# 只取前BIN_PROCESSED个range-bin处理
ra = dsp.aoa_bartlett(steering_vec, radar_cube[:, :, :BINS_PROCESSED], axis=1) 
ra = ra.sum(0) # 对所有chirp的得到的角度信息求和
plt.title('bartlet sum all chirps')
plt.imshow(np.abs(ra))

下面我尝试跟着Bartlett的原理自己写了一个函数,实测效果一致。
只取了第一个chirp,没像上面那样求和,实测这两种方法差不多,对所有chirp得到的角度信息求和(上面那种)的话效果会好一点点点点点。
radar_cube2 shape : (8, 128) - (vx, adc_samples)

ra = np.zeros((ANGLE_BINS, BINS_PROCESSED), dtype=np.complex128)
steering_vec = np.conjugate(steering_vec) # 要取共轭!
for iin range(BINS_PROCESSED): # 遍历每一个range-bin,得到每个range-bin上的角度信息tmp = steering_vec @ radar_cube2[:, i]ra[:, i] = tmp
plt.title('CBF conjugate')
# openradar bartlett源码就是把得到的每个角度的复数取模然后再平方,这样写是对的,我对比过用dsp.aoa_bartlett的range-angle图,效果一样
plt.imshow(np.abs(ra) ** 2) 

Capon

radar_cube shape : (128, 8, 256) - (chirps, vx, adc_samples)

range_azimuth = np.zeros((ANGLE_BINS, BINS_PROCESSED), dtype=np.complex128)
for i in range(BINS_PROCESSED):range_azimuth[:, i], beamWeights[:,i] = dsp.aoa_capon(radar_cube[:, :, i].T, steering_vec, magnitude=False)
plt.title('256 chirps Capon')
plt.imshow(np.abs(range_azimuth))

此外,如果你的chirp小于8,比如1个frame里只有1个chirp,那么就要把openradar的dsp.cov_matrix里的这个判断删了:

if x.shape[0] > x.shape[1]:warnings.warn("cov_matrix input should have Vrx as rows. Needs to be transposed", RuntimeWarning)x = x.T

大概意思就是他觉得你输错了,还帮你改正了,但实际上你没输错,他帮了倒忙,就是要处理一帧里面,1个chirp,8根虚拟天线的数据,anyway,如果一帧里只有小于虚拟天线个数的chirp数的话,把这段判断删了或注释掉就行。

下面的视频是我对这Bartlett/CBF和Capon和Angle-FFT在openradar的circle.bin这个数据上的效果对比:

这篇关于OpenRadar DOA函数 Bartlett/CBF和Capon使用的文章就介绍到这儿,希望我们推荐的文章对编程师们有所帮助!



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

相关文章

使用Python实现矢量路径的压缩、解压与可视化

《使用Python实现矢量路径的压缩、解压与可视化》在图形设计和Web开发中,矢量路径数据的高效存储与传输至关重要,本文将通过一个Python示例,展示如何将复杂的矢量路径命令序列压缩为JSON格式,... 目录引言核心功能概述1. 路径命令解析2. 路径数据压缩3. 路径数据解压4. 可视化代码实现详解1

Pandas透视表(Pivot Table)的具体使用

《Pandas透视表(PivotTable)的具体使用》透视表用于在数据分析和处理过程中进行数据重塑和汇总,本文就来介绍一下Pandas透视表(PivotTable)的具体使用,感兴趣的可以了解一下... 目录前言什么是透视表?使用步骤1. 引入必要的库2. 读取数据3. 创建透视表4. 查看透视表总结前言

Python 交互式可视化的利器Bokeh的使用

《Python交互式可视化的利器Bokeh的使用》Bokeh是一个专注于Web端交互式数据可视化的Python库,本文主要介绍了Python交互式可视化的利器Bokeh的使用,具有一定的参考价值,感... 目录1. Bokeh 简介1.1 为什么选择 Bokeh1.2 安装与环境配置2. Bokeh 基础2

Android使用ImageView.ScaleType实现图片的缩放与裁剪功能

《Android使用ImageView.ScaleType实现图片的缩放与裁剪功能》ImageView是最常用的控件之一,它用于展示各种类型的图片,为了能够根据需求调整图片的显示效果,Android提... 目录什么是 ImageView.ScaleType?FIT_XYFIT_STARTFIT_CENTE

Java学习手册之Filter和Listener使用方法

《Java学习手册之Filter和Listener使用方法》:本文主要介绍Java学习手册之Filter和Listener使用方法的相关资料,Filter是一种拦截器,可以在请求到达Servl... 目录一、Filter(过滤器)1. Filter 的工作原理2. Filter 的配置与使用二、Listen

Pandas使用AdaBoost进行分类的实现

《Pandas使用AdaBoost进行分类的实现》Pandas和AdaBoost分类算法,可以高效地进行数据预处理和分类任务,本文主要介绍了Pandas使用AdaBoost进行分类的实现,具有一定的参... 目录什么是 AdaBoost?使用 AdaBoost 的步骤安装必要的库步骤一:数据准备步骤二:模型

使用Pandas进行均值填充的实现

《使用Pandas进行均值填充的实现》缺失数据(NaN值)是一个常见的问题,我们可以通过多种方法来处理缺失数据,其中一种常用的方法是均值填充,本文主要介绍了使用Pandas进行均值填充的实现,感兴趣的... 目录什么是均值填充?为什么选择均值填充?均值填充的步骤实际代码示例总结在数据分析和处理过程中,缺失数

如何使用 Python 读取 Excel 数据

《如何使用Python读取Excel数据》:本文主要介绍使用Python读取Excel数据的详细教程,通过pandas和openpyxl,你可以轻松读取Excel文件,并进行各种数据处理操... 目录使用 python 读取 Excel 数据的详细教程1. 安装必要的依赖2. 读取 Excel 文件3. 读

Python的time模块一些常用功能(各种与时间相关的函数)

《Python的time模块一些常用功能(各种与时间相关的函数)》Python的time模块提供了各种与时间相关的函数,包括获取当前时间、处理时间间隔、执行时间测量等,:本文主要介绍Python的... 目录1. 获取当前时间2. 时间格式化3. 延时执行4. 时间戳运算5. 计算代码执行时间6. 转换为指

解决Maven项目idea找不到本地仓库jar包问题以及使用mvn install:install-file

《解决Maven项目idea找不到本地仓库jar包问题以及使用mvninstall:install-file》:本文主要介绍解决Maven项目idea找不到本地仓库jar包问题以及使用mvnin... 目录Maven项目idea找不到本地仓库jar包以及使用mvn install:install-file基