GIS剖分数据模型研究

论文作者 ： 陈东

空间数据编码化计算

---

编码真实性计算

引入$J_E (Code)$ 判断网格的真实性，1为假网格 0 为真网格
$J_E(Code) = \epsilon_{B1} \cup\epsilon_{L1}\cup\sim\epsilon_{B2}\cup\sim\epsilon_{L2}\cup\sim\epsilon_{B3}\cup\sim\epsilon_{L3}$
$\epsilon_{B1} = [(\partial C_1 \gg 1 )-(code\_B \gg 23) ]\gg 8 ;$
$\epsilon_{L1} = [\partial C_1 - (code\_L\gg23)]>>8 ;$
$\epsilon_{B2} =code\_B \ll 9 \gg 28\oplus\partial C_2 ;$
$\epsilon_{L2} =code\_L \ll 9 \gg 28\oplus\partial C_2 ;$
$\epsilon_{B3} =code\_B \ll 15 \gg 28\oplus\partial C_2 ;$
$\epsilon_{L3} =code\_L \ll 15 \gg 28\oplus\partial C_2 ;$
其中$\partial C_1 = (10110100)_2$
　　$\partial C_2 = (1111)_2$
　　$\oplus 表示异或$

逐个分析每个算子的作用：
$\epsilon_{B1}是检测纬向度级编码是否有效$
$~~~~~首先\partial C_1 = (180)_{10} 右移一位，code\_B\gg23 得到度数编码 ~~~~$，这里可以看出， 整个算法的思路是 如果纬向编码真实，即小于等于90°。则括号内的值小于$2^7=128$.继续右移八位后，结果为0,则网格为真。这里应当是利用了二进制减法中的补码原理1。 但就这一算法来讲，还存在问题，即没有考虑南北半球带来的符号位的问题。 应当增加一步消除符号位 更改过得算法如下
$\epsilon _{B1} = [(\partial C_1\gg1)-(code\_B\ll1\gg24)]\gg8;$
$同理更正\epsilon _{L1}= [\partial C_1-(code\_L\ll1\gg24)]\gg8$

以此类推 $\epsilon_{B2} 、\epsilon_{B3}分别计算了分级和秒级编码是否有效，即对应了第二次和第三次扩展$
对于分、秒 二者扩展规律相同。所以其判断规则应当类似。 这里发现，无效编码对应的为60,61,62,63 也就是二进制中的 111100 111101 111110 111111 。所以只需要筛选分、秒级编码 前四位为1111的 即为无效编码。 这里使用了异或运算符2进行判断(?)。所以在移位操作时，只需要先通过左移操作去除度级（分级）编码，再右移28位（要保留4位编码 32-4=28）。 其算法如上文所示。

针对异或运算 即保证了 如果二者相同 返回0 如果二者不同 返回1

以上就是对于$J_E(Code)$的分析与修正

2. 位移运算

暂时思路为一次一个网格进行移动，每次移动后判断网格是否存在。若不存在，则继续移动，计数不变。若存在，计数加一。
可以看出位移运算较复杂，不适应当前计算需求。需要发展新的计算方法。

3.格网距离测算

定义算符 $Dis_G$
论文中给出的计算公式：
$(\Delta C_B,\Delta C_L) = Dis_G^1(CodeA，CodeB)\\\qquad\qquad=(\sum_{code_B=min\_codeB}^{min \_codeB} \epsilon_B(code_B,\sum_{code_L=min\_codeL}^{min \_codeL} \epsilon_L(code_L)$
应当是有问题的：
需要改为：
$\qquad (\Delta C_B,\Delta C_L) = Dis_G^1(CodeA，CodeB)\\\qquad\qquad\qquad=(\sum_{code_B=min\_codeB}^{max \_codeB} \epsilon_B(code_B,\sum_{code_L=min\_codeL}^{max \_codeL} \epsilon_L(code_L)$

其中$\epsilon(code)=\epsilon_1 \cup\sim \epsilon_2 \cup\sim \epsilon_3$

其思路是先统一得出所有横纵向格子 再进行筛选。
比较复杂。

其他解决思路

• 能否通过加一个小余量进行判断？
  -能否通过重定义加减法的进制进行设计？

20160727
打卡

---

1. 补码 第一步找原码，如果是负数，就在原码的基础上全部取反后加1 ↩
2. 异或: 0 $\oplus$0 = 0, 1$\oplus$0 = 1 0$\oplus$1=1 1$\oplus$1=0 ↩