[发明专利]一种基于局部网格的动目标活动范围记录方法

权利要求书

1.一种基于局部网格的动目标活动范围记录的方法，其特征在于，包括以下步骤：

步骤1、选择地理坐标系统，剖分地理空间形成初始网格；

步骤2、旋转初始网格；

步骤3、设置网格分界线；

步骤4、计算生成网格编码；

步骤5、基于剖分后的网格区域，实时记录动目标经过的网格单元编码以及时间信息；

步骤1包括：选择地理坐标系统，确定需剖分的地理空间范围：动目标活动区域的起始位置G₁的坐标为(Lng_G1,Lat_G1)，终止位置G₂的坐标为(Lng_G2,Lat_G2),G₁和G₂两个点在投影坐标系下y轴方向距离L₁千米，x轴方向距离L₂千米；网格起始点记为A，其经纬度坐标为(Lng_A,Lat_A)，其中A点处于G1在方位角90度距离L₂/2的位置，以网格起始点A经度作为网格中心线经度，以固定Dis千米为基本网格单元向中心线两侧分别扩展m列网格线，从网格中心点纬度递增方向扩展n行网格线，形成(2m)*(n-1)个相等地理距离单元组成的初始网格，其中动目标的平均速度为V千米/小时，动目标在网格单元的停留时间为u秒，Dis＝V*u/3600，L₁1，L₂1，m≥1，n1，m＝(L₂/2/Dis+1)，n＝L₁/Dis+1；

步骤2包括：初始网格中心线是一条垂直于赤道的地理垂直直线，不与网格中心点所在经度线重合，通过计算旋转角度，再根据选择角度旋转网格使网格中心线与经度线重合；

步骤2中，旋转角度计算过程如下：选择网格起始点A同一经度线上的另一个点B，点B的经纬度坐标(Lng_B,Lat_B)，其中Lng_A＝Lng_B,Lat_BLat_A，获取A点与B点投影坐标系下对应的平面直角坐标分别为(x_A,y_A)和(x_B,y_B)，取另外一点C的坐标为(x_A,y_B)，旋转角度θ计算公式如下：

步骤3包括：从网格区域外轮廓的上边界开始，沿着网格线取点，直到外轮廓的下边界为止，形成网格分界线，网格分界线将整个网格区域一分为二，其中网格分界线的起始点、终止点分别在网格外轮廓的上方和下方；

步骤3包括：从网格区域外轮廓的下边界开始，沿着网格线取点，直到外轮廓的上边界为止，形成网格分界线，网格分界线将整个网格区域一分为二，其中网格分界线的起始点、终止点分别在网格外轮廓的下方和上方；

步骤4包括：一个网格单元编码由编码头以及网格编号组成，其中网格编号根据所述网格单元所在行列号以及分界线所在位置确定，网格编码组成结构为：编码头+网格单元行号+网格单元列号；设定网格分界线左侧网格区域的编码头为W，右侧网格区域的编码头为Y，则一个网格单元Z编码的计算过程如下：

步骤4-1，网格单元Z的4个角点的经纬度坐标已知，通过4个角点经纬度坐标计算得到网格单元Z的中心点坐标(Lng_center,Lat_center)，对应的直角坐标为(x_center,y_center)，

则旋转前的网格中心点直角坐标(x_c,y_c)计算公式如下所示：

步骤4-2，以网格中心线作为0列，网格外轮廓下边界线作为0行，则网格单元Z所在的行列号计算公式如下：

其中Dis为网格间距；row_z,col_z分别为网格单元Z的行号和列号；offset_row为行偏移量，offset_col为列偏移量，计算网格单元行列号时offset_row＝1，offset_col＝1；

根据公式(3)与公式(4)计算网格分界线每一个点坐标对应的网格行列号集合{(row₁,col₁)(row₂,col₂)......(row_k,col_k)}，此时offset_row＝0.5，offset_col＝0.5，其中row_k和col_k分别表示第k个点坐标对应的网格行号和列号；

从网格行列号集合中得到分界线最大与最小的列号分别为col_large和col_small；

步骤4-3，分界线在网格区域每一行的列号集合为：{(1,col₁)(2,col₂)......(n-1,col_n-1)}，则网格单元Z编码头Head_z确认公式如下：

其中，i为网格行号，col_i为第i行在分界线的网格列号；

步骤4-4，当网格单元Z编码头为R时，相对于网格分界线所在的网格列号col_real计算公式如下：

当网格单元Z编码头为L时，相对于网格分界线所在网格列号col_real计算公式如下：

最终网格单元Z的编码组成形式为Head_z+row_z+col_real。

2.根据权利要求1所述的方法，其特征在于，步骤5包括：将动目标活动的日期以及编号作为唯一标识该动目标的字段，动目标经过的网格单元编码集合作为位置信息，将动目标进入与离开网格单元的时间作为时间信息，{动目标标识，网格编码，时间信息}组成了记录动目标活动范围的结构体。