This thread has been locked.

If you have a related question, please click the "Ask a related question" button in the top right corner. The newly created question will be automatically linked to this question.

AWR1642BOOST: + DCA1000 EVM; mmwave Studio, Zero-doppler, its meaning and removal

Part Number: AWR1642BOOST

Hi Support team ,

I am using AWR1642BOOST (ES2.0) + DCA1000 EVM (FPGA version 2.8) +mmwave studio 2.1.1.0 with complex 2x ADC data format.

I am currently using advanced frame config with 2 subframes but for sake of simplicity, I will just share data for one subframe.

I have searched e2e database, but could not find my exact answer, therefore I am approaching you for guidance and help.

I request some basic but concise guidance on zero-doppler and then may be an advanced question. Please refer the following figure which shows 2d fft from mmwave studio GUI and no target (especially no moving target) in radar field of view (FoV)

Figure 1

1)    I understand that zero-doppler bin shows the target with zero radial velocity. However, I am not much clear about it though it is a very basic question. Can you please advise that what does it mean physically ?

2)   After that, can you please explain that what does zero-doppler mean mathematically ? I understand that above plot is produced by magnitude plotting of dB values of 2d fft (with fft shift anf prior flipus function to bring +ve range inupper half of plot) from the radar cube data with velocity at x-axis. But what causes those relatively high energy cells to appear in zero-doppler column  ? And more specially what and why is there a very high energy (may be single) cell at the 0 intersection of range and velocity axes ?

3)   Then a moving target passed the radar FoV. As the target speed was more than the maximum unambiguous speed (as dictated by relevant chirp configuration), therefore there was velocity ambiguity and also there was energy smear / leakage as can be see in the following figure 2.

Figure 2

4)   I am learning to use Chinese Remainder Theorem (CRT) for resolving velocity ambiguity. I understand that CRT is prone to noise and even one cell error in 2d fft Range-doppler map will produce big error but I am trying to study and implement any way. I am familiar that CRT is applied in TI mmwave SRR lab but I am not familiar with CCS and I am preferring to work on raw adc data from DCA1000 and I am working  on my own MATLAB coding so  I am working for a very basic and generic understanding. I cant understand CRT implementation in SRR therefore I request help here.

Any way, for using CRT, I need to identify, one particular cell in this 2d fft range-doppler magnitude map. I know and I repeat that we can already see the energy spread in figure and the reading will be probably not right, but that is another topic.

My question to you is that how can I obtain a single cell ( and its row (Velocity) and column (range) indices from this 2d fft range-doppler map data ? Should I use some  " maximum find " type of function to get information about cell with highest magnitude of 2d fft ? 

5)  You can see in figure 2, that it appears that, most probably, the single cell at 0 place, at intersection of range and velocity axes, has the highest energy value. So most probably, a function trying to find maximum magnitude will select this cell, but it will be wrong as I am looking for a cell with a velocity, not a  cell with zero-velocity. So what could be a solution in that case  ?

6)  Can  a solution to para 5 above could be that I just replace the center column (zero-doppler column) with zeros in the relevant data ? But as the heatmap is working on relative magnitudes, then it will result in more energy smearing in other cells.

Another option could be that I find minimum value in the 2d fft range-doppler magnitude map and then I use that minimum value to fill the full column of zero-doppler column. Can you please comment ?

Or may be you can advise something else, more better ?

Thanks in advance and regards

  • Hi,

    Let me review your questions and get back to you by end of the week

    thank you
    Cesar

  • Hi.

    Zero doppler means that the relative radial velocity is equal to zero.

    So, the high energy shows in the zero doppler bins of the heatmap for the static objects.

    Mathematically the zero doppler is the absolute value of the average values of the chirps. It is the DC component of the frequency.

    The CRT algorithm is used in the SRR demo after detection. So you would need to decide how you want to implement the detection in your matlab code.

    For more inforamation about CRT please see following

    http://staff.elka.pw.edu.pl/~jmisiure/esptr_base/lect_fmcw/kk97fm.pdf

    Thank you

    Cesar

  • Hi Cesar,

    Thank you very much for the mathematical and physical explanation of zero-doppler, it was really very helpful.

    I had already came across the research paper that you provided for CRT, i studied it again as you provided it and it was really helpful. In fact I actually carried out all the steps practically using this paper for understanding CRT.

    I am relatively comfortable with CRT now but there are some basic issues with its implementation that are, in fact, not related to CRT. Those question are, specially, about detection and / or finding the maximum cell value that would be subsequently used in CRT. 

    In view of above, I repeat some of the content from my original above post and some added material and request you to please extend guidance on it, if possible and convenient.

    1)   A moving target is made to pass the radar FoV. As the target speed is more than the maximum unambiguous speed (as dictated by relevant chirp configuration), therefore there is velocity ambiguity and also there is energy smear / leakage as can be see in the following figure 2.

    Figure 2

    2)   I am learning to use Chinese Remainder Theorem (CRT) for resolving velocity ambiguity. I understand that CRT is prone to noise and even one cell error in 2d fft Range-doppler map will produce big error but I am trying to study and implement any way. I am familiar that CRT is applied in TI mmwave SRR lab but I am not familiar with CCS and I am preferring to work on raw adc data from DCA1000 and I am working  on my own MATLAB coding so  I am working for a very basic and generic understanding. I cant understand CRT implementation in SRR therefore I request help here.

    I am familiar with the basics of CRT and I have read and fully understood the paper you have advised above.

    Any way, for using CRT, I need to identify, one particular cell in this 2d fft range-doppler magnitude map. I know and I repeat that we can already see the energy spread in figure and the reading will be probably not right, but that is another topic.

    More important issue : My question to you is that how can I obtain a single cell ( and its row (Velocity) and column (range) indices from this 2d fft range-doppler map data ? Should I use some  " maximum find " type of function to get information about cell with highest magnitude of 2d fft ? 

    May be you can reply to this question by consulting application of CRT in SRR lab , just a suggestion . . .

     

    3)  A relatively trivial issue :: You can see in figure 2, that it appears that, most probably, the single cell at 0 place, at intersection of range and velocity axes, has the highest energy value. So most probably, a function trying to find maximum magnitude will select this cell, but it will be wrong as I am looking for a cell with a velocity, not a  cell with zero-velocity. So what could be a solution in that case  ?

    4)  Can  a solution to para 3 above could be that I just replace the center column (zero-doppler column) with zeros in the relevant data ?

    But as the heatmap is working on relative magnitudes, then it will result in more energy smearing in other cells.

    Another option could be that I find minimum value in the 2d fft range-doppler magnitude map and then I use that minimum value to fill the full column of zero-doppler column. Can you please comment ?

    May be you can reply to this question by consulting application of CRT in SRR lab , just a suggestion . . .

    Or may be you can advise something else, more better ?

    Thanks in advance and regards

  • Hi,

    CFAR type algorithms are used to determine the max in the range doppler heatmap.

    https://en.wikipedia.org/wiki/Constant_false_alarm_rate

    Thank you

    Cesar