中位中心栅格遍历法

最新推荐文章于 2023-03-03 17:48:15 发布

Yjiokm

最新推荐文章于 2023-03-03 17:48:15 发布

阅读量1.4k

点赞数 1

分类专栏： GIS 简单算法文章标签： GIS 欧氏距离变换栅格化

本文链接：https://blog.csdn.net/qq_34581491/article/details/51489429

版权

本文详细介绍了中位中心栅格遍历法，这是一种在GIS领域中常用的技术，用于处理和分析地理空间数据。通过这种方法，可以有效地进行欧氏距离变换和栅格化的操作，提高数据处理效率。

摘要由CSDN通过智能技术生成

中位中心栅格遍历法

#include <iostream>
#include <cmath>

namespace geo
{
	struct POINT
	{
		double x;
		double y;
	};

	//栅格法寻找中位中心
	POINT GetMedianCentre(const POINT * pts, int n, double d/*设置栅格化格网宽度*/)
	{
		double min_x = pts[0].x;
		double min_y = pts[0].y;
		double max_x = pts[0].x;
		double max_y = pts[0].y;

		for(int i = 0; i < n; i++)
		{
			if(pts[i].x < min_x)
			{
				min_x = pts[i].x;
			}
			else if(pts[i].x > max_x)
			{
				max_x = pts[i].x;
			}
			if(pts[i].y < min_y)
			{
				min_y = pts[i].y;
			}
			else if(pts[i].y > max_y)
			{
				max_y = pts[i].y;
			}
		}

		int raster_row = (max_y - min_y) / d + 1;//注意 + 1
		int raster_column = (max_x - min_x) / d + 1;

		struct INTPT {int row; int col;};//行列数对
		INTPT * pts_o = new INTPT [n];//原有点集栅格化后的行列数对
		for(int i = 0; i < n; i++)
		{