生成有箭头的流向Line.shp 并通过ArcMap绘制带箭头的图，根据D8算法，输入Flow Direction

最新推荐文章于 2023-07-12 09:54:34 发布
最新推荐文章于 2023-07-12 09:54:34 发布
阅读量1.3k 收藏 6
点赞数
文章标签： c++
版权声明：本文为博主原创文章，遵循 CC 4.0 BY-SA 版权协议，转载请附上原文出处链接和本声明。
本文链接：https://blog.csdn.net/qq_41644037/article/details/123904527
版权
该代码示例展示了如何使用C++结合GDAL库读取流向数据（ASC文件），并利用ArcGIS的D8算法处理，最后将处理结果保存为ESRI Shapefile（.shp）矢量图。主要涉及GDAL的环境配置、数据读取、流向计算和形状文件创建。
摘要由CSDN通过智能技术生成

简单说明：代码C++，配置GDAL环境（必须），ArcMap
附一个配置GDAL环境的教程，感谢这位大佬
VS2017编译配置GDAL库

输入流向数据的 asc文件，附网盘测试链接
链接：https://pan.baidu.com/s/1-Q35OktO8dOPbjdkI1XC8w
提取码：p2ap

代码如下：

#include <iostream>
#include <fstream>
#include <string>
#include <string.h>
#include <sstream>
#include "ogrsf_frmts.h"
#include "gdal.h"
#include "gdal_priv.h"
#include "gdal_alg.h"
using namespace std;

//ToString该函数转载自他人
template<typename T> string toString(const T& t)
{
	ostringstream oss;  //创建一个格式化输出流
	oss << t;             //把值传递如流中
	return oss.str();
}
void Draw_arrow(int fd, int row, int col, int fd_row, int fd_col, double fd_xll, double fd_yll, double fd_res, double &x1, double &x2, double &y1, double &y2);

void main()
{
	int i, j, k;
	int fd_row = 359;
	int fd_col = 719;
	double fd_xll = 112.6;
	double fd_yll = 32.500277777778;
	double fd_res = 0.00027777777799999;
	string flow_path_file = "E:\\arcgis\\fdr.asc";
	string line_shp_path = "E:\\arcgis\\fdr.shp";
	string flow_file_path;
	int **fd_value = new int*[fd_row*fd_col];
	for (i = 0; i < fd_row*fd_col; i++)
	{
		fd_value[i] = new int[8];
	}
	int **order = new int*[fd_row*fd_col];
	for (i = 0; i < fd_row; i++)
	{
		order[i] = new int[fd_col];
	}
	for (i = 0; i < fd_row*fd_col; i++)
		for (j = 0; j < 8; j++)
			fd_value[i][j] = 0;

	ifstream FDR(flow_path_file.c_str());
	if (!FDR)
	{
		cout << "no in FDR" << endl;
		system("pause");
	}
	string str;
	for (i = 0; i < 6; i++)
		getline(FDR, str);
	//cout << "1" << endl;
	//system("pause");
	int tmp;
	for (i = 0; i < fd_row; i++)
	{
		for (j = 0; j < fd_col; j++)
		{
			FDR >> tmp;
			if (tmp == -9999 || tmp == 0)
			{
				continue;
			}
			// 这里根据ArcGIS的D8算法流向修改了一下，不过问题不大，个人比较喜欢这么计算
			if (tmp == 32)
			{
				fd_value[i*fd_col + j][0] = 1;
			}
			else if (tmp == 64)
			{
				fd_value[i*fd_col + j][1] = 1;
			}
			else if (tmp == 128)
			{
				fd_value[i*fd_col + j][2] = 1;
			}
			else if (tmp == 16)
			{
				fd_value[i*fd_col + j][3] = 1;
			}
			else if (tmp == 1)
			{
				fd_value[i*fd_col + j][4] = 1;
			}
			else if (tmp == 8)
			{
				fd_value[i*fd_col + j][5] = 1;
			}
			else if (tmp == 4)
			{
				fd_value[i*fd_col + j][6] = 1;
			}
			else if (tmp == 2)
			{
				fd_value[i*fd_col + j][7] = 1;
			}
		}
	}
	GDALAllRegister();
	const char * pFileName = line_shp_path.c_str();
	const char * pszDriverName = "ESRI Shapefile";
	GDALDriver * poDriver = GetGDALDriverManager()->GetDriverByName(pszDriverName);
	if (poDriver == NULL)
		return;
	GDALDataset * poDs = poDriver->Create(pFileName, 0, 0, 0, GDT_Unknown, NULL);
	if (poDs == NULL)
		return;
	//wkbPolygon表示绘制多边形；wkbLinearRing表示绘制线环
	OGRLayer * poLayer = poDs->CreateLayer("ring", NULL, wkbLineString, NULL);//不删
	if (poLayer == NULL)
	{
		return;
	}
	int line_num;
	int count_line;
	double x1, x2, y1, y2;
	//	system("pause");
	for (i = 0; i < fd_row; i++)
	{
		for (j = 0; j < fd_col; j++)
		{
			count_line = 0;
			for (k = 0; k < 8; k++)
			{
				if (fd_value[i*fd_col + j][k] == 1)
					count_line++;

			}
			if (count_line == 0)
				continue;
			if (count_line == 1)
			{
				for (k = 0; k < 8; k++)
				{
					if (fd_value[i*fd_col + j][k] == 1)
					{
						Draw_arrow(k + 1, i, j, fd_row, fd_col, fd_xll, fd_yll, fd_res, x1, x2, y1, y2);
						OGRFeature * poFeature = OGRFeature::CreateFeature(poLayer->GetLayerDefn());
						OGRLineString line;

						line.setNumPoints(2);
						line.setPoint(0, x1, y1, 0);//第一个输入变量是序号，第二、三、四分别是线XYZ坐标
						line.setPoint(1, x2, y2, 0);
						poFeature->SetGeometry(&line);
						if (poLayer->CreateFeature(poFeature) != OGRERR_NONE)
						{
							return;
						}
						OGRFeature::DestroyFeature(poFeature);
					}
				}
			}
			else if (count_line > 1)
			{
				for (k = 0; k < 8; k++)
				{
					if (fd_value[i*fd_col + j][k] == 1)
					{	// k+1
						Draw_arrow(k + 1 + 8, i, j, fd_row, fd_col, fd_xll, fd_yll, fd_res, x1, x2, y1, y2);
						OGRFeature * poFeature = OGRFeature::CreateFeature(poLayer->GetLayerDefn());
						OGRLineString line;
						line.setNumPoints(2);
						line.setPoint(0, x1, y1, 0);//第一个输入变量是序号，第二、三、四分别是线XYZ坐标
						line.setPoint(1, x2, y2, 0);
						poFeature->SetGeometry(&line);
						if (poLayer->CreateFeature(poFeature) != OGRERR_NONE)
						{
							return;
						}
						OGRFeature::DestroyFeature(poFeature);
					}
				}
			}
		}
	}
	GDALClose(poDs);
	OGRCleanupAll();
	for (i = 0; i < fd_row; i++)
		delete[]fd_value[i];
	delete[]fd_value;
	cout << "***** Finished *****" << endl;
	//system("pause");
}

void Draw_arrow(int fd, int row, int col, int fd_row, int fd_col, double fd_xll, double fd_yll, double fd_res, double &x1, double &x2, double &y1, double &y2)
{
	int i = row;
	int j = col;
	if (fd == 1)
	{
		x1 = fd_xll + (j + 1)* fd_res;
		x2 = fd_xll + (j)*fd_res;
		y1 = fd_yll + (fd_row - i - 1)*fd_res;
		y2 = fd_yll + (fd_row - i) * fd_res;
	}
	else if (fd == 2)
	{
		x1 = fd_xll + (2 * j + 1)* fd_res / 2.0;
		x2 = fd_xll + (2 * j + 1)*fd_res / 2.0;
		y1 = fd_yll + (fd_row - i - 1)*fd_res;
		y2 = fd_yll + (fd_row - i) * fd_res;
	}
	else if (fd == 3)
	{
		x1 = fd_xll + (j)* fd_res;
		x2 = fd_xll + (j + 1)*fd_res;
		y1 = fd_yll + (fd_row - i - 1)*fd_res;
		y2 = fd_yll + (fd_row - i) * fd_res;
	}
	else if (fd == 4)
	{
		x1 = fd_xll + (j)*fd_res;
		x2 = fd_xll + (j + 1)* fd_res;
		y1 = fd_yll + (fd_row - (2 * i + 1) / 2.0)*fd_res;
		y2 = fd_yll + (fd_row - (2 * i + 1) / 2.0) * fd_res;
	}
	else if (fd == 8)
	{
		x1 = fd_xll + (j + 1)* fd_res;
		x2 = fd_xll + (j)*fd_res;
		y1 = fd_yll + (fd_row - (2 * i + 1) / 2.0)*fd_res;
		y2 = fd_yll + (fd_row - (2 * i + 1) / 2.0) * fd_res;
	}

	else if (fd == 5)
	{
		x1 = fd_xll + (j)*fd_res;
		x2 = fd_xll + (j + 1)* fd_res;
		y1 = fd_yll + (fd_row - i) * fd_res;
		y2 = fd_yll + (fd_row - i - 1)*fd_res;
	}

	else if (fd == 6)
	{
		x1 = fd_xll + (2 * j + 1)*fd_res / 2.0;
		x2 = fd_xll + (2 * j + 1)* fd_res / 2.0;
		y1 = fd_yll + (fd_row - i) * fd_res;
		y2 = fd_yll + (fd_row - i - 1)*fd_res;
	}

	else if (fd == 7)
	{
		x1 = fd_xll + (j + 1)*fd_res;
		x2 = fd_xll + (j)* fd_res;
		y1 = fd_yll + (fd_row - i) * fd_res;
		y2 = fd_yll + (fd_row - i - 1)*fd_res;
	}

	//以下是别的，可忽略；
	else if (fd == 9)
	{
		x1 = fd_xll + (j + 1)* fd_res - 0.5*fd_res;
		x2 = fd_xll + (j)*fd_res;
		y1 = fd_yll + (fd_row - (i + 1))*fd_res + 0.5*fd_res;
		y2 = fd_yll + (fd_row - i) * fd_res;
	}
	else if (fd == 10)
	{
		x1 = fd_xll + (2 * j + 1)* fd_res / 2.0;
		x2 = fd_xll + (2 * j + 1)*fd_res / 2.0;
		y1 = fd_yll + (fd_row - (i + 1))*fd_res + 0.5*fd_res;
		y2 = fd_yll + (fd_row - i) * fd_res;
	}
	else if (fd == 11)
	{
		x1 = fd_xll + (j)* fd_res + 0.5*fd_res;
		x2 = fd_xll + (j + 1)*fd_res;
		y1 = fd_yll + (fd_row - i - 1)*fd_res + 0.5*fd_res;
		y2 = fd_yll + (fd_row - i) * fd_res;
	}
	else if (fd == 12)
	{
		x1 = fd_xll + (j)*fd_res + 0.5*fd_res;
		x2 = fd_xll + (j + 1)* fd_res;
		y1 = fd_yll + (fd_row - (2 * i + 1) / 2.0)*fd_res;
		y2 = fd_yll + (fd_row - (2 * i + 1) / 2.0) * fd_res;
	}
	else if (fd == 16)
	{
		x1 = fd_xll + (j + 1)* fd_res - 0.5*fd_res;
		x2 = fd_xll + (j)*fd_res;
		y1 = fd_yll + (fd_row - (2 * i + 1) / 2.0)*fd_res;
		y2 = fd_yll + (fd_row - (2 * i + 1) / 2.0) * fd_res;
	}

	else if (fd == 13)
	{
		x1 = fd_xll + (j)*fd_res + 0.5*fd_res;
		x2 = fd_xll + (j + 1)* fd_res;
		y1 = fd_yll + (fd_row - i) * fd_res - 0.5*fd_res;
		y2 = fd_yll + (fd_row - i - 1)*fd_res;
	}

	else if (fd == 14)
	{
		x1 = fd_xll + (2 * j + 1)*fd_res / 2.0;
		x2 = fd_xll + (2 * j + 1)* fd_res / 2.0;
		y1 = fd_yll + (fd_row - i) * fd_res - 0.5*fd_res;
		y2 = fd_yll + (fd_row - i - 1)*fd_res;
	}

	else if (fd == 15)
	{
		x1 = fd_xll + (j + 1)*fd_res - 0.5*fd_res;
		x2 = fd_xll + (j)* fd_res;
		y1 = fd_yll + (fd_row - i) * fd_res - 0.5*fd_res;
		y2 = fd_yll + (fd_row - i - 1)*fd_res;
	}
}

将代码生成的fdr.shp加载入ArcMap
在这里插入图片描述

对其形状标志选择向右的箭头
在这里插入图片描述
放大看就完成了
在这里插入图片描述

宇天y
关注
  • 0
    点赞
  • 踩
  • 6
    收藏
    觉得还不错? 一键收藏
  • 打赏
    打赏
  • 0
    评论
Arcgis中怎样设置调查路线线型（带箭头的虚线），附带1:1万地形图符号库
weixin_40406267的博客
11-09 831
Arcgis中怎样设置调查路线线型（带箭头的虚线），附带1:1万地形图符号库
Arcgis-在Arcgis中添加箭头、线段、图框的方法
热门推荐
Along6的博客
05-12 2万+
问题描述：一般在研究区图中，我们经常为了能清晰表示自己的研究区所在地，会用到全国-省市-研究区所在地，从大到小的不同范围的研究区图，这就会用到箭头指示，线段链接，矩形图框圈定等，具体如何使用下面介绍两种不同的方法：
ArcGIS教程：“流向”的工作原理
07-21 1777
获取表面的水文特征的关键之一是能够确定从栅格中的每个像元流出的方向。这可通过流向工具来完成。 　　该工具把表面作为输入，然后输出一个显示从每个像元流出方向的栅格。如果选择了输出下降率栅格数据选项，则会创建一个以百分比的形式表示的输出栅格，显示从沿流向的每个像元到像元中心之间的路径长度的高程的最大变化率。如果选择了强制所有边缘像元向外流动选项，则表面栅格边缘处的所有像元都将从表面栅格向外流动。
arcgis箭头绘制
08-28
arcgis 箭头绘制、箭头绘制、箭头绘制、箭头绘制、箭头绘制、箭头绘制、箭头绘制、箭头绘制、箭头绘制、箭头绘制
ArcGIS基于DEM计算水流方向的方法（D8算法）
weixin_34302561的博客
09-16 6815
ArcGIS采用D8算法计算水流方向（9.3.1后新增），输入数据应首先完成了洼地填充处理： One of the keys to deriving hydrologic characteristics about a surface is the ability to determine the direction of flow from every cell in the raster. ...
Arcgis实现路线等线方向箭头表示
12-21
Arcgis实现路线等线方向箭头表示
国界线.zip_arcmap 国界 shp_shp 国界线_中国国界gis图_国界数据_国界线数据网
07-13
中国国界线GIS数据，shp格式，完整版
China.shp.rar_bou2_4p.shp_shp 中国地图_shp中国_中国shp_地图shp
07-13
中国GIS地图，shp格式，ARCview可以使用
C#不使用ArcEngine生成.shp文件
08-09
然而，如果你不想依赖像ArcEngine这样的大型库，而是希望通过C#编程语言直接生成`.shp`文件，那么这里将介绍一种不使用ArcEngine的方法。 首先，我们需要了解`.shp`文件的结构。`.shp`文件是基于ESRI Shapefile格式...
提取.shp文件中信息并绘图保存
04-21
这可以通过遍历`shapes`数组并根据特定条件（比如地理位置或形状大小）来实现。例如，你可以根据经纬度坐标将数据划分为多个区域： ```matlab minLon = -180; maxLon = 180; minLat = -90; maxLat = 90; deltaLon =...
长江经济带地级市shp图.zip
07-09
【标题】"长江经济带地级市shp图.zip"是一个包含有关长江经济带地级市地理信息的数据压缩包。这个文件集合提供了对中国长江流域内地级行政区域的详细地图数据，采用的是Esri的Shapefile（shp）格式，这是一种广泛...
arcgis api for js绘制箭头图
qq_42430948的博客
09-17 2461
//人才需求箭头图 rcxqToRcgy(data){ const _this = this; let png = this.graduatPng if (_this.rcxqToRcgyLayer != "") { // console.log(123456); _this.rcxqToRcgyLayer.clear(); } //如果传的data为[] 之前又有人才需求箭头流向图层 先清除图层 在ret.
arcgis for android: 绘制带流向箭头的线
qq_37532493的博客
07-10 3765
arcgis for android: 绘制带流向箭头的线 需求描述 项目中需要绘制如下图所示的线,但arcgis for android中linesymbol只提供了单箭头/颜色/实虚线等简单配置项,无法通过配置更改线型. 思路 既然没法直接配置线型,那我们就把箭头画上去! 利用PictureMarkerSymbol把箭头画到线上,我们需要这些东西: 1.箭头图标 2.每条线上各自需要多少箭头,以及相应的坐标位置 3.由于线的倾斜程度各不相同,我们需要算出每条线上箭头的旋转角度 得到这些东西 我们就可以
arcgis符号方向_ArcGIS制图表达-河流渐变与符号旋转
weixin_34237125的博客
01-14 2459
在制图表达原理的章节中已经介绍过制图表达的两个字段，rule_ID和override ，之后的章节都是围绕这规则来讲解(也就是rule_ID记录的规则)，那么本章节就开始介绍override 字段所显示的内容。原则中，override 记录的是规则覆盖的内容。先来看看河流线渐变和测站符号垂直河流旋转的效果。河流渐变效果符号自动旋转效果河流线渐变首先来看河流线渐变效果。图上的效果是一种夸大的效果，主...
【Matlab】智能优化算法_流向算法FDA
最新发布
敲代码两年半的练习生博客
07-12 1198
在流域中，过量或有效降雨量是指流过地表且未渗入土壤的降雨量。事实上，降水和损失（如截留、蒸发和渗透）后留在地表的水量被称为直接径流。到目前为止，已经提出了几种确定直接径流的方法，其中之一是φ指数法。单位为厘米/小时的指数φ代表降雨过程中的平均失水量，超过此限度的降雨量将直接转化为径流。换言之，指数φ的高位和等于直接径流的高度。从降雨量中减去每个时间间隔的指数φ，即可得到直接径流量。图1说明了φ指数的概念。关系式（1）表示直接径流计算方法。
手把手教会你使用Arcgis画流场图
MarineWitch的博客
07-22 1万+
手把手教会你使用Arcgis画流场图 手把手教会你使用Arcgis画流场图 我以前一直都是用MATLAB绘制流场图的，但是工作原因让我不得不使用Arcgis画流场图。在网上找到了一些教程，但是感觉讲得不是很清晰，让Arcgis小白的我走了很多弯路。好不容易琢磨出来了，于是给大家分享一下如何用Arcgis画流场图。方法过于不方便，希望大家永远不需要用Arcgis画流场。 1.准备数据 （1）这里我使用的是HYCOM流场数据，大家手中有什么流场数据都可以使用。因为我的lon、lat数组均是1维：lon(301X
基于ArcMap的地图符号制作(详细教程)
jingss_3的专栏
01-06 1万+
本帖来源于http://bbs.esrichina-bj.cn/ESRI/thread-40793-1-1.html1 地图符号的具体制作过程　　在我们制作以上几种符号之前，首先要单独创建一个自己的样式库即符号库，具体的做法是：启动ArcMap，点击Tools/styles/ style manager菜单项。接着在弹出的style manager 对话框中点击style 的下拉按钮，在弹
ArcGIS教程：制作风或水流速流向图
11-11 1万+
风速风向图或流速流向图相信大家都已经见过不少，但不知道有多少人会制作这样炫的专题图，下面这边文章向我们展示了当基本数据U和V矢量被存储时，怎样计算风或水流的速度和方向和对其进行符号化。 　　 　　 　　为了获取风或水流的速度和方向，风速计或多普勒流速剖面仪在两个相互垂直的方向U和V量测了风或水的速度。U是指向东的方向，V是指向北的方向。有时候，U指的是纬向速度，V指的是径向速
GIS开发：胖箭头的画法
地理信息技术杂谈 博客，欢迎技术交流
04-23 1250
之前分享过在两类地图上的箭头画法，效果是单线，顶部加了箭头的效果，这里分享一下在地图上胖箭头的画法。 胖箭头是在箭头的首尾两点上，横向进行了扩充，最后以面状图层的形式画出来。 胖箭头根据首尾两点生成面的方式： –将首尾两点连接成直线，计算直线在地图中的北方向夹角，可以使用turf.js进行计算； –计算和首尾直线垂直的6个点，分别是箭头尾部的两个点和箭头头部的4个样式点，turf.js中有方法计算这6个点，可以参见文章开始的简单箭头中的方法； –将以上的6个点按照顺序连接起来，形成面状的效果，就能在地图
评论
添加红包

请填写红包祝福语或标题

个

红包个数最小为10个

元

红包金额最低5元

当前余额3.43元 前往充值 >
需支付：10.00元
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包

打赏作者

宇天y

你的鼓励将是我创作的最大动力

¥1 ¥2 ¥4 ¥6 ¥10 ¥20
扫码支付：¥1
获取中
扫码支付

您的余额不足，请更换扫码支付或充值

打赏作者

实付元
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明：

1.余额是钱包充值的虚拟货币，按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载，可以购买VIP、付费专栏及课程。

余额充值