[GIS算法] DEM相关处理 - 窗口分析(以平均值为例)| 裁剪 - C语言实现

本文主要是介绍[GIS算法] DEM相关处理 - 窗口分析(以平均值为例)| 裁剪 - C语言实现,希望对大家解决编程问题提供一定的参考价值,需要的开发者们随着小编来一起学习吧!

#include<stdio.h>
#include<stdlib.h>typedef struct{double x0,y0; //左下角点的坐标值int dx,dy; //栅格单元大小int ycount,xcount; //行列号double **M; //矩阵
}DEM;DEM* CreateDEM(int xcount,int ycount); //创建DEM
int InitDEM(DEM *pDEM); //初始化DEM
void PrintDEM(DEM dem); //输出DEM* SubDEM(DEM dem, int row, int col, int nrow, int ncol); //裁剪
DEM* WindowAnalysis_Mean(DEM dem, int size); //窗口分析/*测试数据
5 5
1 3 5 2 0
4 3 8 3 4
8 5 0 6 3
2 1 3 5 6
9 7 2 4 1
*/
int main() {DEM *pDEM, *pMean;int xcount,ycount;printf("输入行列数:\n>>> ");scanf("%d%d", &xcount, &ycount);pDEM = CreateDEM(xcount, ycount); //创建DEMInitDEM(pDEM); //初始化PrintDEM(*pDEM); //输出printf("\n");pMean = WindowAnalysis_Mean(*pDEM, 3); //3x3窗口分析求平均值PrintDEM(*pMean);return 0;
}DEM* CreateDEM(int xcount, int ycount) {DEM *p;int i;p = (DEM *)malloc(sizeof(DEM)); if (!p) exit(0);p->xcount=xcount;p->ycount=ycount;p->M = (double **)malloc(sizeof(double *)*xcount); if (!p->M) exit(0);for (i=0; i<xcount; i++) {p->M[i] = (double *)malloc(sizeof(double)*ycount);if (!p->M[i]) exit(0);}return p;
}
int InitDEM(DEM *pDEM) {int i,j;printf("以矩阵的形式输入DEM的值:\n");for (i=0; i<pDEM->xcount; i++) {for (j=0; j<pDEM->ycount; j++) {scanf("%lf", &pDEM->M[i][j]);}}return 1;
}
DEM* SubDEM(DEM dem, int row, int col, int nrow, int ncol) {DEM *sub;int i,j;if (row+nrow>dem.xcount || col+ncol>dem.ycount || row>dem.xcount || col>dem.ycount) return NULL;sub = CreateDEM(nrow, ncol);for (i=0; i<nrow; i++) {for (j=0; j<ncol; j++) {sub->M[i][j] = dem.M[i+row][j+col];}}return sub;
}
void PrintDEM(DEM dem) {int i,j;for (i=0; i<dem.xcount; i++) {for (j=0; j<dem.ycount; j++) {printf("%lf\t", dem.M[i][j]);}printf("\n");}
}
DEM* WindowAnalysis_Mean(DEM dem, int size) {// size:窗口大小// 例如:3X3的窗口 size=3DEM *mean;int i,j;int x,y;int step;int x_step, y_step;int flag;double sum;mean = CreateDEM(dem.xcount, dem.ycount);for (i=0; i<mean->xcount; i++) {for (j=0; j<mean->ycount; j++) {sum = 0; //窗口的总和step = size/2; //位移量flag = 0; //是否越界// 计算窗口的sumfor (x_step=-step; x_step<=step && !flag; x_step++) {for (y_step=-step; y_step<=step && !flag; y_step++) {x = i+x_step;y = j+y_step;if (x<0 || y<0 || x>=mean->xcount || y>=mean->ycount) { //越界flag = 1; //越界} else { //没有越界sum += dem.M[x][y];}}}//赋值if (flag) { //越界了mean->M[i][j] = dem.M[i][j];} else { //没有越界mean->M[i][j] = sum/(size*size);}}}return mean;
}

这篇关于[GIS算法] DEM相关处理 - 窗口分析(以平均值为例)| 裁剪 - C语言实现的文章就介绍到这儿,希望我们推荐的文章对编程师们有所帮助!



http://www.chinasem.cn/article/705487

相关文章

Android kotlin中 Channel 和 Flow 的区别和选择使用场景分析

《Androidkotlin中Channel和Flow的区别和选择使用场景分析》Kotlin协程中,Flow是冷数据流,按需触发,适合响应式数据处理;Channel是热数据流,持续发送,支持... 目录一、基本概念界定FlowChannel二、核心特性对比数据生产触发条件生产与消费的关系背压处理机制生命周期

c++ 类成员变量默认初始值的实现

《c++类成员变量默认初始值的实现》本文主要介绍了c++类成员变量默认初始值,文中通过示例代码介绍的非常详细,对大家的学习或者工作具有一定的参考学习价值,需要的朋友们下面随着小编来一起学习学习吧... 目录C++类成员变量初始化c++类的变量的初始化在C++中,如果使用类成员变量时未给定其初始值,那么它将被

Qt使用QSqlDatabase连接MySQL实现增删改查功能

《Qt使用QSqlDatabase连接MySQL实现增删改查功能》这篇文章主要为大家详细介绍了Qt如何使用QSqlDatabase连接MySQL实现增删改查功能,文中的示例代码讲解详细,感兴趣的小伙伴... 目录一、创建数据表二、连接mysql数据库三、封装成一个完整的轻量级 ORM 风格类3.1 表结构

基于Python实现一个图片拆分工具

《基于Python实现一个图片拆分工具》这篇文章主要为大家详细介绍了如何基于Python实现一个图片拆分工具,可以根据需要的行数和列数进行拆分,感兴趣的小伙伴可以跟随小编一起学习一下... 简单介绍先自己选择输入的图片,默认是输出到项目文件夹中,可以自己选择其他的文件夹,选择需要拆分的行数和列数,可以通过

Python中将嵌套列表扁平化的多种实现方法

《Python中将嵌套列表扁平化的多种实现方法》在Python编程中,我们常常会遇到需要将嵌套列表(即列表中包含列表)转换为一个一维的扁平列表的需求,本文将给大家介绍了多种实现这一目标的方法,需要的朋... 目录python中将嵌套列表扁平化的方法技术背景实现步骤1. 使用嵌套列表推导式2. 使用itert

Python使用vllm处理多模态数据的预处理技巧

《Python使用vllm处理多模态数据的预处理技巧》本文深入探讨了在Python环境下使用vLLM处理多模态数据的预处理技巧,我们将从基础概念出发,详细讲解文本、图像、音频等多模态数据的预处理方法,... 目录1. 背景介绍1.1 目的和范围1.2 预期读者1.3 文档结构概述1.4 术语表1.4.1 核

Python使用pip工具实现包自动更新的多种方法

《Python使用pip工具实现包自动更新的多种方法》本文深入探讨了使用Python的pip工具实现包自动更新的各种方法和技术,我们将从基础概念开始,逐步介绍手动更新方法、自动化脚本编写、结合CI/C... 目录1. 背景介绍1.1 目的和范围1.2 预期读者1.3 文档结构概述1.4 术语表1.4.1 核

在Linux中改变echo输出颜色的实现方法

《在Linux中改变echo输出颜色的实现方法》在Linux系统的命令行环境下,为了使输出信息更加清晰、突出,便于用户快速识别和区分不同类型的信息,常常需要改变echo命令的输出颜色,所以本文给大家介... 目python录在linux中改变echo输出颜色的方法技术背景实现步骤使用ANSI转义码使用tpu

Python使用python-can实现合并BLF文件

《Python使用python-can实现合并BLF文件》python-can库是Python生态中专注于CAN总线通信与数据处理的强大工具,本文将使用python-can为BLF文件合并提供高效灵活... 目录一、python-can 库:CAN 数据处理的利器二、BLF 文件合并核心代码解析1. 基础合

Python使用OpenCV实现获取视频时长的小工具

《Python使用OpenCV实现获取视频时长的小工具》在处理视频数据时,获取视频的时长是一项常见且基础的需求,本文将详细介绍如何使用Python和OpenCV获取视频时长,并对每一行代码进行深入解析... 目录一、代码实现二、代码解析1. 导入 OpenCV 库2. 定义获取视频时长的函数3. 打开视频文