激光雷达水平扫描的时空误差

15195775117

4个多月过去了，今天又有了心得，跟个帖子
------------割-------------------
编写还原函数
假设前提：扫描从正北开始，顺时针旋转
那么对于扫描图上的每一个点(x,y)，我们知道激光扫到这个位置旋转的角度theta，也就知道了所花的时间t，假设的平均风场(v1,v2)×t，就是t=0时该点的位置(x0,y0)和扫到该点时的位置(x,y)之差，公式表达为：
(x0,y0)+(v1,v2)*t=(x,y)
求角度我们用向量来算：
如果x≥0，(x,y)与代表正北的(0,1)的夹角就是旋转角theta；
如果x<0，旋转角就是360-向量夹角。
代码如下：
function theta,x,y
t=acos(total([0,1]*[x,y])/sqrt(x^2+y^2))/!pi*180
return,(x ge 0)?t:(360-t)
end
根据实际扫描的图(商业机密就不放图了)，发现其中有明显的螺旋污染带，
以后有空我把数据要过来，使用这个代码，试着把图谱还原回去

15195775117

通过对水平扫描雷达图的观察，发现很多都有螺旋式的污染带，这条污染带的规律是：
越靠近扫描原点(雷达所在位置)越浓越窄，越远离原点，污染带越宽越淡