1.一种基于FPGA的球面距离定点计算方法,其特征在于,所述方法用于实现现场可编程逻辑门阵列对两定点间的球面距离的计算,所述方法包括如下内容:1)根据正弦函数的性质,在正弦函数的自变量取值区间内对正弦函数的函数值进行量化;2)基于量化后的正弦函数的函数值,构造正弦函数的基准查找表、一阶增量查找表和二阶增量查找表,并将正弦函数的基准查找表、一阶增量查找表和二阶增量查找表分别存储在现场可编程逻辑门阵列内部存储器中,以存储所有量化后的正弦函数的函数值;3)根据两点之间的球面距离计算公式,确定反正弦函数的自变量及自变量取值范围,并在反正弦函数的自变量取值区间内对反正弦函数的函数值进行非均匀量化和压缩;4)基于量化后的反正弦函数的函数值,构造反正弦函数的基准查找表、一阶增量查找表和二阶增量查找表,并将反正弦函数的基准查找表、一阶增量查找表和二阶增量查找表分别存储在现场可编程逻辑门阵列内部存储器中,以存储所有量化后的反正弦函数的函数值;5)将待计算球面距离的两定点的坐标数据和球面半径数据分别输入现场可编程逻辑门阵列,若输入的坐标数据格式为赤道坐标数据格式,现场可编程逻辑门阵列内部的数据选择器依次选取坐标数据中包含的经度数据和纬度数据,通过地址变换确定对应的正弦函数的基准查找表、一阶增量查找表和二阶增量查找表的读地址,利用确定的读地址访问现场可编程逻辑门阵列内部存储器存储的正弦函数的基准查找表、一阶增量查找表和二阶增量查找表,并通过数据变换得到经度数据和纬度数据对应的正弦函数值;若输入的坐标数据格式为直角坐标数据格式,现场可编程逻辑门阵列直接将坐标数据存储在寄存器中;现场可编程逻辑门阵列内部的数据选择器依次选取寄存器中存储的两定点的直角坐标数据中的x坐标、y坐标和z坐标计算差值,并将坐标差值送入乘法器和加法器计算坐标差值的平方和,确定对应的反正弦函数的自变量,现场可编程逻辑门阵列基于反正弦函数的自变量产生用于访问反正弦函数的基准查找表、一阶增量查找表及二阶增量查找表的读地址,利用确定的读地址访问现场可编程逻辑门阵列内部存储器存储的反正弦函数的基准查找表、一阶增量查找表和二阶增量查找表,并通过数据变换确定和输出两定点间的球面距离。
2.根据权利要求1所述的基于FPGA的球面距离定点计算方法,其特征在于,所述步骤1包括:将正弦函数的自变量取值区间[0,π/2]等分为N个长度为π/2N的子区间,且当自变量位于子区间[iπ/2N,(i+1)π/2N)(i=0,1,…,N-1)时,对应的正弦函数的函数值取为y(i)=sin((i+0.5)π/2N);其中,N的具体数值根据距离估计精度要求进行设定。
3.根据权利要求2所述的基于FPGA的球面距离定点计算方法,其特征在于,所述步骤2包括:基于正弦函数的性质并根据N的具体数值,设定三个参数B、E和M,并使N=B·E·M;分别构造长度为B·E的正弦函数的基准查找表,长度为B·(M-1)的正弦函数的一阶增量查找表,以及长度为B(E-1)·(M-1)的正弦函数的二阶增量查找表。
4.根据权利要求3所述的基于FPGA的球面距离定点计算方法,其特征在于,长度为B·E的正弦函数的基准查找表依次存储有如下B·E个数值:y(0),y(M),y(2M),…,y[(BE-1)M]。
5.根据权利要求4所述的基于FPGA的球面距离定点计算方法,其特征在于,长度为B·(M-1)的正弦函数的一阶增量查找表依次存储有B组数据,且每组数据包括有M-1个数值;其中,第k(k=0,1,…,B-1)组的M-1数据为:
6.根据权利要求5所述的基于FPGA的球面距离定点计算方法,其特征在于,长度为B(E-1)·(M-1)的正弦函数的二阶增量查找表依次存储有B组数据,且每组数据包括有(E-1)·(M-1)个数值;其中,第k(k=0,1,…,B-1)组的(E-1)·(M-1)数据为:
7.根据权利要求1所述的基于FPGA的球面距离定点计算方法,其特征在于,假设:两定点在空间直接坐标系下的坐标分别为(x 1 ,y 1 ,z 1 )和(x 2 ,y 2 ,z 2 ),球面半径为r,两定点间的球面距离表述为 所述步骤3包括:以w=(x 1 -x 2 ) 2 +(y 1 -y 2 ) 2 +(z 1 -z 2 ) 2 作为反正弦函数的自变量,对应的反正弦函数的函数值取为 并在反正弦函数的自变量取值区间[0,2r 2 ]内对反正弦函数的函数值进行非均匀量化和压缩。
8.根据权利要求7所述的基于FPGA的球面距离定点计算方法,其特征在于,在反正弦函数的自变量取值区间[0,2r 2 ]内对反正弦函数的函数值进行非均匀量化和压缩,包括:a.根据球面距离计算精度要求,设定最小量化间隔Δ,计算如下 个反正弦函数的函数值: b.对N′个函数值进行选择,包括:第1至第N 1 个函数值全部保留;第N 1 +1至N 2 个函数值进行 倍均匀抽取,对应的量化间隔为 第N 2 +1至第N 3 个函数值进行 倍均匀抽取,对应的量化间隔为 第N 3 +1至第N 4 个函数值进行 倍均匀抽取,对应的量化间隔为 按此方式不断提升量化间隔至 以涵盖全部N′个函数值,并得到压缩后的函数值数目为 其中, 以及N 1 ,N 2 ,…,N s 根据反正弦函数的性质和N′的具体数值选定。
9.根据权利要求8所述的基于FPGA的球面距离定点计算方法,其特征在于,所述步骤4包括:基于反正弦函数的性质并根据N″的具体数值,设定三个参数B″、E″和M″,并使N″=B″·E″·M″;分别构造长度为B″·E″的反正弦函数的基准查找表,长度为B″·(M″-1)的反正弦函数的一阶增量查找表,以及长度为B″(E″-1)·(M″-1)的反正弦函数的二阶增量查找表。
10.一种用于实现如权利要求1至9中任一项所述的基于FPGA的球面距离定点计算方法的球面距离定点计算装置,其特征在于,所述计算装置包括:依次电连接的坐标格式转换单元、坐标缓存单元和球面距离计算单元;所述坐标格式转换单元用于将经纬度坐标转换为空间直角坐标;所述坐标缓存单元用于接收并存储外部输入的空间直角坐标,以及接收并存储所述坐标格式转换单元输出的转换结果;所述球面距离计算单元用于接收所述坐标缓存单元输出的空间直角坐标,利用接收到的空间直角坐标计算并输出球面距离。