莫顿码Morton code的c语言实现

十进制 Morton码可以使用栅格单元的行列号计算

###原理部分
遵循C语言规范，矩阵的第一行为“0”行、第一列为“0”列，先将十进制的行列号转换成二进制数，进行“位”运算操作，即行号和列号的二进制数两两交叉，得到以二进制数表示的MD码，再将其转换为十进制数。
例如图中第二行和第三列对应的栅格单元，其二进制的行列号分别为:I＝0010，J＝0011；得到的MD码为:MD＝（0001101）2＝（13）10；用类似的方法，也可以由MD码反求栅格单元的行列号。

###基本算法

// 输入 x 和 y, 交叉后为 xyxyxy。
unsigned short x;
// x在奇数位
unsigned short y;
// y在偶数位
unsigned int md = 0; // 存储生成的Md码
for (int i = 0; i < sizeof(x) * CHAR_BIT; i++) // 核心交叉算法，基本思想为移位和位运算..
{
md |= (x & 1U << i) << (i + 1) | (y & 1U << i) << i;
}

###完整代码
#####代码一：

#include<stdio.h>
#include<limits.h>	//定义CHAR_BIT
int main(void)
{
unsigned short i,j;	//行号、列号
unsigned short n;	//栅格大小
unsigned int md[n][n];	//每个码位的值
printf("请输入栅格的行列数（0<=n<=65535)：");
if(scanf("%hd",&n) == 0)
return 0;
for(i = 0;i < n;i++)	//行循环
{
for(j = 0;j < n;j++)	//列循环
{
md[i][j] = 0;
for(int x = 0;x < sizeof(unsigned short) * CHAR_BIT;x++)	//行列号计算
md[i][j] |= (i & 1U << x) << (x + 1) | (j & 1U << x) << x;
printf("%5u ",md[i][j]);	//打印输出
}
printf("n");
}
return 1;
}

#####代码二：

#include<stdio.h>
#include<limits.h>
void to_binary1(unsigned short n);
void to_binary2(unsigned int n);
int main(void)
{
unsigned short x,y;
unsigned int md;
printf("请输入行号和列号（q to quit)：");
while(scanf("%hd %hd",&x,&y) != 0)
{
printf("行号I=(%hd)10=(",x);
to_binary1(x);
printf(")2n");
printf("列号J=(%hd)10=(",y);
to_binary1(y);
printf(")2n");
md = 0;
for(int i = 0;i < sizeof(x) * CHAR_BIT;i++)
md |= (x & 1U << i) << (i + 1) | (y & 1U << i) << i;
printf("莫顿码MD=(%u)10=(",md);
to_binary2(md);
printf(")2n");
printf("n请输入行号和列号（q to quit)：");
}
printf("Done.n");
return 0;
}
void to_binary1(unsigned short n)
{
int r;
r = n % 2;
if (n >= 2)
to_binary1(n / 2);
putchar(r == 0 ? '0' : '1');
return;
}
void to_binary2(unsigned int n)
{
int r;
r = n % 2;
if (n >= 2)
to_binary2(n / 2);
putchar(r == 0 ? '0' : '1');
return;
}

最后

以上就是安详棒棒糖为你收集整理的莫顿码Morton code的c语言实现的全部内容，希望文章能够帮你解决莫顿码Morton code的c语言实现所遇到的程序开发问题。

如果觉得靠谱客网站的内容还不错，欢迎将靠谱客网站推荐给程序员好友。