C#实现GIS中莫顿码的计算

计算步骤：

1.将横坐标和纵坐标分别转化为二进制，需要注意的是转化后的二进制位数要相同（本文中以8为二进制码进行演示）

2.交叉排序后组合成16位2进制码

3.将16位二进制码重新转化为整形数字，此数字即为我们需要的莫顿码

C#代码如下：

public string getEightBinary(string binary)
        {
            //十进制转二进制字符串
            //   string str = Convert.ToString(row, 2);
            //输出: 1010
            string strBinary = Convert.ToString(int.Parse(binary), 2);
            //保证是8位的二进制数，位数不够用0占位
            int len = 8 - strBinary.Length;
            string prefix = "";
            for(int i = 0; i < len; i++)
            {
                prefix += "0";
            }
            strBinary = prefix + strBinary;
            return strBinary;
        }
        //获取莫顿码，输入row,col范围0-255，输出莫顿码范围0-25535
        public int getMorton(string row, string col)
        {
            if (int.Parse(row) > 255 || int.Parse(col) > 255)
            {
                return -1;
            }
            //转化为8位二进制码
            string x = getEightBinary(row);
            string y = getEightBinary(col);
            //生成二进制的M码
            string strm = "";
            //交叉算法，将两个8位二进制码合并位16位莫顿码
            for (int i = 0; i < x.Length; i++)
            {
                strm = strm + x.Substring(i, 1) + y.Substring(i, 1);
            }
            //将16为2进制码转化为整形数字，即为莫顿码
            int m = Convert.ToInt32(strm, 2);
            return m;
        }