1.一种空间关联规则的挖掘方法，其特征在于，包括：对土地变化信息和岸线变化信息进行处理，生成多个合并网格单元和岸线变化强度；

基于所述多个合并网格单元和岸线变化强度，建立模型运算决策表；以及对所述模型运算决策表进行挖掘，以获得所述空间关联规则；

其中，在所述对土地变化信息和岸线变化信息进行处理，生成多个合并网格单元和岸线变化强度中，包括：对所述土地变化信息进行扇环形网格划分，以生成多个样本网格单元；

确定所述多个样本网格单元中每个样本网格单元的土地变化强度；

其中，在所述对所述土地变化信息进行扇环形网格划分，以生成多个样本网格单元中，包括：在平行于岸线方向上，对土地变化信息的区域岸线进行矢量化，获得矢量化数据；

基于所述矢量化数据，以预设间隔从沿海至内陆生成多个列缓冲区，其中，每个列缓冲区为所述土地变化信息对应的土地区域的一列；

根据所述多个列缓冲区，构造贴合岸线地形的扇环形，其中包括：以各土地区域的海岸的两端垂直于海岸线的边构成扇环形的两条半径，以两半径之间的岸线为扇环形的内弧、最外围列缓冲区的外边为外弧，构造所述贴合岸线地形的扇环形；

对所述扇环形进行旋转划分，获得所述多个样本网格单元，其中包括：作两半径向海洋方向的延长线并相交，并以交点为圆心，以两半径为起、止边按预设度数间隔进行旋转划分，即获得海‑陆向垂线与各缓冲区相交形成的样本网格单元。

2.根据权利要求1所述的挖掘方法，其特征在于，在对土地变化信息和岸线变化信息进行处理之前，还包括：基于土地覆被信息生成所述土地变化信息；以及

基于码头岸线信息生成所述岸线变化信息。

3.根据权利要求1所述的挖掘方法，其特征在于，在确定所述多个样本网格单元中每个样本网格单元的土地变化强度中，包括：依据所述每个样本网格单元的总面积s0和所述每个样本网格单元中的土地变化面积s1，获取所述每个样本网格单元的土地变化率f，其中包括：f＝s1/s0，f∈[0，1]；

根据所述土地变化率f，通过自然断点法获取所述每个样本网格单元的土地变化强度。

4.根据权利要求3所述的挖掘方法，其特征在于，在所述对土地变化信息和岸线变化信息进行处理，生成多个合并网格单元和岸线变化强度中，还包括：依据网格合并规则对所述多个样本网格单元进行合并，生成所述多个合并网格单元；

其中，所述网格合并规则为当所述多个样本网格单元中的一个样本网格单元与其他的至少一个样本网格单元位于同一单元行且空间位置相邻且土地变化强度相同时实现所述一个样本网格单元与其他的至少一个样本网格单元之间的合并，以生成所述多个合并网格单元中的一个合并网格单元；

所述多个合并网格单元为一处于多个在同一垂直岸线方向上的样本序列。

5.根据权利要求2所述的挖掘方法，其特征在于，在所述对土地变化信息和岸线变化信息进行处理，生成多个合并网格单元和岸线变化强度中，还包括：依据所述岸线变化信息的区域岸线总长度a0和区域岸线变化长度a1，获取所述岸线变化信息的岸线变化率d，其中包括：d＝a1/a0，d∈[0，1]；

根据所述岸线变化率d，通过自然断点法获取所述岸线变化强度。

6.根据权利要求1所述的挖掘方法，其特征在于，所述模型运算决策表为：{Lm，DnDnDn………}，其中，Lm为岸线变化强度，m为所述岸线变化强度的等级序号，Dn为土地变化强度，n为所述土地变化强度的等级序号。

7.根据权利要求1所述的挖掘方法，其特征在于，在所述对所述模型运算决策表进行挖掘以获得所述空间关联规则中，包括：通过Apriori算法对所述模型运算决策表进行挖掘，以获得所述空间关联规则。

8.根据权利要求1所述的挖掘方法，其特征在于，还包括：对所述空间关联规则进行预设最小支持度阈值和预设最小置信度阈值的判断，以确定所述空间关联规则为强关联规则。

9.一种空间关联规则的挖掘装置，以实现所述权利要求1‑8中任一项所述的方法。

10.一种电子设备，包括：

一个或多个处理器；

存储装置，用于存储一个或多个程序，

其中，当所述一个或多个程序被所述一个或多个处理器执行时，使得所述一个或多个处理器实现权利要求1‑8中任一项所述的方法。

11.一种计算机可读存储介质，其上存储有可执行指令，该指令被处理器执行时使处理器实现权利要求1‑8中任一项所述的方法。