发明名称 一种脉冲雷达成像中的坐标转换方法
摘要 本发明公开了一种脉冲雷达成像中的坐标转换方法,建立了新型查找表,通过附加信息标示死点情况,在坐标输出时按附加信息对死点进行补点处理,提高了雷达成像质量;建立的查找表只记录了0°~45°极坐标到直角坐标的转换,圆周内其他极坐标到直角坐标的转换通过映射公式、查表法和转换公式处理完成,缩减了查找表需要的储存空间;通过使用片外存储器,为处理器节约了片内存储空间;本发明可应用于基于FPGA硬件平台的脉冲雷达图像显示、基于PC的脉冲雷达图像显示,还可以应用于各种极坐标软硬件显示系统。
申请公布号 CN102313886B 申请公布日期 2013.04.03
申请号 CN201110241144.2 申请日期 2011.08.22
申请人 宁波成电泰克电子信息技术发展有限公司 发明人 李宏波;李浩;田丹;罗长阳;鄢林;敬洁;冉元进;陈闹;陶吉怀
分类号 G01S13/89(2006.01)I;G01S7/02(2006.01)I 主分类号 G01S13/89(2006.01)I
代理机构 宁波奥圣专利代理事务所(普通合伙) 33226 代理人 程晓明;周珏
主权项 1.一种脉冲雷达成像中的坐标转换方法,其特征在于包括以下处理步骤:s1.雷达扫描一周的扫描线总数为N,每条雷达扫描线上有扫描点M个,n为第n条扫描线,在雷达扫描圆周的0°~45°范围内,将所有雷达扫描线的扫描点的极坐标(r,θ)通过公式<maths num="0001"><![CDATA[<math><mfenced open='{' close=''><mtable><mtr><mtd><msub><mi>x</mi><mi>i</mi></msub><mo>=</mo><mi>r</mi><mi>cos</mi><mrow><mo>(</mo><mi>&theta;</mi><mo>)</mo></mrow></mtd></mtr><mtr><mtd><msub><mi>y</mi><mi>i</mi></msub><mo>=</mo><mi>r</mi><mi>sin</mi><mrow><mo>(</mo><mi>&theta;</mi><mo>)</mo></mrow></mtd></mtr></mtable></mfenced></math>]]></maths>计算得到对应的直角坐标(x<sub>i</sub>,y<sub>i</sub>),组成数组{(x<sub>i</sub>,y<sub>i</sub>)},其中,N,M,n,i,r都为整数,且n∈[0,N/8],i∈[0,MN/8),r∈[0,M),<img file="FDA0000085056140000012.GIF" wi="195" he="84" />s2.将位于直角坐标系第一象限0°~45°范围内且满足条件<maths num="0002"><![CDATA[<math><mfenced open='{' close=''><mtable><mtr><mtd><mi>y</mi><mo>&le;</mo><mi>x</mi></mtd></mtr><mtr><mtd><msup><mi>x</mi><mn>2</mn></msup><mo>+</mo><msup><mi>y</mi><mn>2</mn></msup><mo>&le;</mo><msup><mi>M</mi><mn>2</mn></msup></mtd></mtr></mtable></mfenced></math>]]></maths>的所有直角坐标点(x,y)找出,组成数组{(x,y)},再将数组{(x,y)}中的每一个直角坐标点分别与数组{(x<sub>i</sub>,y<sub>i</sub>)}中的每一个直角坐标点比较,若满足条件<maths num="0003"><![CDATA[<math><mrow><mfenced open='{' close=''><mtable><mtr><mtd><mi>y</mi><mo>=</mo><msub><mi>y</mi><mi>i</mi></msub></mtd></mtr><mtr><mtd><mi>x</mi><mo>=</mo><msub><mi>x</mi><mi>i</mi></msub></mtd></mtr></mtable></mfenced><mo>,</mo></mrow></math>]]></maths>则将此坐标点(x,y)从数组{(x,y)}中删除,最后将筛选过的数组{(x,y)}重新定义为数组{(x′,y′)};s3.选定数组{(x<sub>i</sub>,y<sub>i</sub>)}中的第i个直角坐标点,i=0,1,2,...,MN/8-1,再遍历数组{(x′,y′)},找出满足条件<maths num="0004"><![CDATA[<math><mfenced open='{' close=''><mtable><mtr><mtd><msup><mi>y</mi><mo>&prime;</mo></msup><mo>&le;</mo><msub><mi>y</mi><mi>i</mi></msub></mtd></mtr><mtr><mtd><msup><mi>x</mi><mo>&prime;</mo></msup><mo>&le;</mo><msub><mi>x</mi><mi>i</mi></msub></mtd></mtr></mtable></mfenced></math>]]></maths>的所有直角坐标组成新的数组{(x<sub>ij</sub>′,y<sub>ij</sub>′)}并从数组{(x′,y′)}中删除,其中(x<sub>ij</sub>′,y<sub>ij</sub>′)表示数组{(x<sub>i</sub>,y<sub>j</sub>)}中第i个坐标点周围满足条件<maths num="0005"><![CDATA[<math><mfenced open='{' close=''><mtable><mtr><mtd><msup><mi>y</mi><mo>&prime;</mo></msup><mo>&le;</mo><msub><mi>y</mi><mi>i</mi></msub></mtd></mtr><mtr><mtd><msup><mi>x</mi><mo>&prime;</mo></msup><mo>&le;</mo><msub><mi>x</mi><mi>i</mi></msub></mtd></mtr></mtable></mfenced></math>]]></maths>的第j个坐标点,j为整数且j∈[0,J),J表示数组{(x<sub>ij</sub>′,y<sub>ij</sub>′)}中所含元素个数;用直角坐标(x<sub>i</sub>,y<sub>i</sub>)的横坐标x<sub>i</sub>分别减去数组{(x<sub>ij</sub>′,y<sub>ij</sub>′)}中的每一个直角坐标的横坐标x<sub>ij</sub>′,直角坐标(x<sub>i</sub>,y<sub>i</sub>)的纵坐标y<sub>i</sub>分别减去数组{(x<sub>ij</sub>′,y<sub>ij</sub>′)}中的每一个直角坐标的纵坐标y<sub>ij</sub>′得到(x<sub>i</sub>-x<sub>i1</sub>′,y<sub>i</sub>-y<sub>i1</sub>′)、(x<sub>i</sub>-x<sub>i2</sub>′,y<sub>i</sub>-y<sub>i2</sub>′)...(x<sub>i</sub>-x<sub>ij</sub>′,y<sub>i</sub>-y<sub>ij</sub>′),再用一个宽度为2J*k bit的二进制数add来表示(x<sub>i</sub>-x<sub>i1</sub>′,y<sub>i</sub>-y<sub>i1</sub>′)、(x<sub>i</sub>-x<sub>i2</sub>′,y<sub>i</sub>-y<sub>i2</sub>′)...(x<sub>i</sub>-x<sub>ij</sub>′,y<sub>i</sub>-y<sub>ij</sub>′),表示方法为:二进制数add的0~(k-1)bit表示数值x<sub>i</sub>-x<sub>i1</sub>′、二进制数add的k~(2k-1)bit表示数值y<sub>i</sub>-y<sub>i1</sub>′,2k~(3k-1)bit表示数值x<sub>i</sub>-x<sub>i2</sub>′、二进制数add的3k~(4k-1)bit表示数值y<sub>i</sub>-y<sub>i2</sub>′,...,二进制数add的(2(j-1)*k)~((2j-1)*k-1)bit表示数值x<sub>i</sub>-x<sub>ij</sub>′、二进制数add的((2j-1)*k)~((2j*k)-1)bit表示数值y<sub>i</sub>-y<sub>ij</sub>′,其中k为所有x<sub>i</sub>-x<sub>ij</sub>′和y<sub>i</sub>-y<sub>ij</sub>′数值中的最大值化为二进制数后的位宽;将直角坐标(x<sub>i</sub>,y<sub>i</sub>)及其相应的add值保存为(x,y,add)<sub>i</sub>,(x,y,add)<sub>i</sub>代表了直角坐标(x<sub>i</sub>,y<sub>i</sub>)和数组{(x<sub>ij</sub>′,y<sub>ij</sub>′)};s4.重复执行s3,其中i依次取值0,1,2,...,MN/8-1,再将所有的(x,y,add)<sub>i</sub>保存到查找表;s5.处理器在接收雷达扫描圆周范围内第n<sub>1</sub>条扫描线及其相应灰度信息后,将第n<sub>1</sub>条扫描线通过映射公式转换到0°~45°范围内相应扫描线n;其中,n<sub>1</sub>∈[0,N),映射公式为<maths num="0006"><![CDATA[<math><mrow><mfenced open='{' close=''><mtable><mtr><mtd><mi>n</mi><mo>=</mo><msub><mi>n</mi><mn>1</mn></msub><mo>,</mo></mtd><mtd><mn>0</mn><mo>&le;</mo><msub><mi>n</mi><mn>1</mn></msub><mo>&le;</mo><mi>N</mi><mo>/</mo><mn>8</mn></mtd></mtr><mtr><mtd><mi>n</mi><mo>=</mo><mi>N</mi><mo>/</mo><mn>4</mn><mo>-</mo><msub><mi>n</mi><mn>1</mn></msub><mo>,</mo></mtd><mtd><mi>N</mi><mo>/</mo><mn>8</mn><mo>&lt;</mo><msub><mi>n</mi><mn>1</mn></msub><mo>&le;</mo><mi>N</mi><mo>/</mo><mn>4</mn></mtd></mtr><mtr><mtd><mi>n</mi><mo>=</mo><msub><mi>n</mi><mn>1</mn></msub><mo>-</mo><mi>N</mi><mo>/</mo><mn>4</mn><mo>,</mo></mtd><mtd><mi>N</mi><mo>/</mo><mn>4</mn><mo>&lt;</mo><msub><mi>n</mi><mn>1</mn></msub><mo>&le;</mo><mn>3</mn><mi>N</mi><mo>/</mo><mn>8</mn></mtd></mtr><mtr><mtd><mi>n</mi><mo>=</mo><mi>N</mi><mo>/</mo><mn>2</mn><mo>-</mo><msub><mi>n</mi><mn>1</mn></msub><mo>,</mo></mtd><mtd><mn>3</mn><mi>N</mi><mo>/</mo><mn>8</mn><mo>&lt;</mo><msub><mi>n</mi><mn>1</mn></msub><mo>&le;</mo><mi>N</mi><mo>/</mo><mn>2</mn></mtd></mtr><mtr><mtd><mi>n</mi><mo>=</mo><msub><mi>n</mi><mn>1</mn></msub><mo>-</mo><mi>N</mi><mo>/</mo><mn>2</mn><mo>,</mo></mtd><mtd><mi>N</mi><mo>/</mo><mn>2</mn><mo>&lt;</mo><msub><mi>n</mi><mn>1</mn></msub><mo>&le;</mo><mn>5</mn><mi>N</mi><mo>/</mo><mn>8</mn></mtd></mtr><mtr><mtd><mi>n</mi><mo>=</mo><mn>3</mn><mi>N</mi><mo>/</mo><mn>4</mn><mo>-</mo><msub><mi>n</mi><mn>1</mn></msub><mo>,</mo></mtd><mtd><mn>5</mn><mi>N</mi><mo>/</mo><mn>8</mn><mo>&lt;</mo><msub><mi>n</mi><mn>1</mn></msub><mo>&le;</mo><mn>3</mn><mi>N</mi><mo>/</mo><mn>4</mn></mtd></mtr><mtr><mtd><mi>n</mi><mo>=</mo><msub><mi>n</mi><mn>1</mn></msub><mo>-</mo><mn>3</mn><mi>N</mi><mo>/</mo><mn>4</mn><mo>,</mo></mtd><mtd><mn>3</mn><mi>N</mi><mo>/</mo><mn>4</mn><mo>&lt;</mo><msub><mi>n</mi><mn>1</mn></msub><mo>&le;</mo><mn>7</mn><mi>N</mi><mo>/</mo><mn>8</mn></mtd></mtr><mtr><mtd><mi>n</mi><mo>=</mo><mi>N</mi><mo>-</mo><msub><mi>n</mi><mn>1</mn></msub><mo>,</mo></mtd><mtd><mn>7</mn><mi>N</mi><mo>/</mo><mn>8</mn><mo>&lt;</mo><msub><mi>n</mi><mn>1</mn></msub><mo>&lt;</mo><mi>N</mi></mtd></mtr></mtable></mfenced><mo>;</mo></mrow></math>]]></maths>通过寻址查找表得到扫描线n对应的M个扫描点数据,将这M个扫描点数据组成数组<img file="FDA0000085056140000022.GIF" wi="254" he="51" />n<sub>1</sub>∈[0,N);s6.将s5得到的数组<img file="FDA0000085056140000023.GIF" wi="241" he="51" />保存到处理器的第一存储器,将第n<sub>1</sub>条扫描线上每个扫描点对应的灰度值保存到处理器的第二存储器;s7.从第一存储器中取出第n<sub>1</sub>条扫描线上的第r个扫描点数据<img file="FDA0000085056140000024.GIF" wi="255" he="51" />r∈[0,M),并从第二存储器读取与<img file="FDA0000085056140000025.GIF" wi="220" he="51" />对应的灰度值;根据<img file="FDA0000085056140000026.GIF" wi="220" he="51" />中的直角坐标<img file="FDA0000085056140000027.GIF" wi="133" he="51" />及add值得到数组<img file="FDA0000085056140000028.GIF" wi="260" he="50" />将直角坐标<img file="FDA0000085056140000029.GIF" wi="134" he="51" />及数组<img file="FDA00000850561400000210.GIF" wi="230" he="50" />中的每一个坐标元素代入转换公式得到直角坐标<img file="FDA00000850561400000211.GIF" wi="166" he="51" />和<img file="FDA00000850561400000212.GIF" wi="284" he="50" />其中,转换公式为<maths num="0007"><![CDATA[<math><mrow><mfenced open='{' close=''><mtable><mtr><mtd><msub><mi>x</mi><mn>1</mn></msub><mo>=</mo><mi>x</mi><mo>,</mo><msub><mi>y</mi><mn>1</mn></msub><mo>=</mo><mi>y</mi><mo>,</mo></mtd><mtd><mn>0</mn><mo>&le;</mo><msub><mi>n</mi><mn>1</mn></msub><mo>&le;</mo><mi>N</mi><mo>/</mo><mn>8</mn></mtd></mtr><mtr><mtd><msub><mi>x</mi><mn>1</mn></msub><mo>=</mo><mi>y</mi><mo>,</mo><msub><mi>y</mi><mn>1</mn></msub><mo>=</mo><mi>x</mi><mo>,</mo></mtd><mtd><mi>N</mi><mo>/</mo><mn>8</mn><mo>&lt;</mo><msub><mi>n</mi><mn>1</mn></msub><mo>&le;</mo><mi>N</mi><mo>/</mo><mn>4</mn></mtd></mtr><mtr><mtd><mrow><msub><mi>x</mi><mn>1</mn></msub><mo>=</mo><mo>-</mo><mi>y</mi><mo>,</mo><msub><mi>y</mi><mn>1</mn></msub><mo>=</mo><mi>x</mi><mo>,</mo></mrow></mtd><mtd><mi>N</mi><mo>/</mo><mn>4</mn><mo>&lt;</mo><msub><mi>n</mi><mn>1</mn></msub><mo>&le;</mo><mn>3</mn><mi>N</mi><mo>/</mo><mn>8</mn></mtd></mtr><mtr><mtd><msub><mi>x</mi><mn>1</mn></msub><mo>=</mo><mo>-</mo><mi>x</mi><mo>,</mo><msub><mi>y</mi><mn>1</mn></msub><mo>=</mo><mi>y</mi><mo>,</mo></mtd><mtd><mn>3</mn><mi>N</mi><mo>/</mo><mn>8</mn><mo>&lt;</mo><msub><mi>n</mi><mn>1</mn></msub><mo>&le;</mo><mi>N</mi><mo>/</mo><mn>2</mn></mtd></mtr><mtr><mtd><msub><mi>x</mi><mn>1</mn></msub><mo>=</mo><mo>-</mo><mi>x</mi><mo>,</mo><msub><mi>y</mi><mn>1</mn></msub><mo>=</mo><mo>-</mo><mi>y</mi><mo>,</mo></mtd><mtd><mi>N</mi><mo>/</mo><mn>2</mn><mo>&lt;</mo><msub><mi>n</mi><mn>1</mn></msub><mo>&le;</mo><mn>5</mn><mi>N</mi><mo>/</mo><mn>8</mn></mtd></mtr><mtr><mtd><msub><mi>x</mi><mn>1</mn></msub><mo>=</mo><mo>-</mo><mi>y</mi><mo>,</mo><msub><mi>y</mi><mn>1</mn></msub><mo>=</mo><mo>-</mo><mi>x</mi><mo>,</mo></mtd><mtd><mn>5</mn><mi>N</mi><mo>/</mo><mn>8</mn><mo>&lt;</mo><msub><mi>n</mi><mn>1</mn></msub><mo>&le;</mo><mn>3</mn><mi>N</mi><mo>/</mo><mn>4</mn></mtd></mtr><mtr><mtd><mrow><msub><mi>x</mi><mn>1</mn></msub><mo>=</mo><mi>y</mi><mo>,</mo><msub><mi>y</mi><mn>1</mn></msub><mo>=</mo><mo>-</mo><mi>x</mi><mo>,</mo></mrow></mtd><mtd><mn>3</mn><mi>N</mi><mo>/</mo><mn>4</mn><mo>&lt;</mo><msub><mi>n</mi><mn>1</mn></msub><mo>&le;</mo><mn>7</mn><mi>N</mi><mo>/</mo><mn>8</mn></mtd></mtr><mtr><mtd><msub><mi>x</mi><mn>1</mn></msub><mo>=</mo><mi>x</mi><mo>,</mo><msub><mi>y</mi><mn>1</mn></msub><mo>=</mo><mo>-</mo><mi>y</mi><mo>,</mo></mtd><mtd><mn>7</mn><mi>N</mi><mo>/</mo><mn>8</mn><mo>&lt;</mo><msub><mi>n</mi><mn>1</mn></msub><mo>&lt;</mo><mi>N</mi></mtd></mtr></mtable></mfenced><mo>;</mo></mrow></math>]]></maths>将直角坐标<img file="FDA0000085056140000032.GIF" wi="166" he="51" />与<img file="FDA0000085056140000033.GIF" wi="220" he="51" />对应的灰度值一起输出到显存;将数组<img file="FDA0000085056140000034.GIF" wi="254" he="51" />中的每一个直角坐标分别与<img file="FDA0000085056140000035.GIF" wi="220" he="51" />对应的灰度值一起输出到显存。
地址 315040 浙江省宁波市高新区扬帆路999弄5号7-4室