基于FPGA的车道识别硬件加速

原创 已于 2025-01-13 09:12:22 修改 · 1.1k 阅读 · 9 ·
CC 4.0 BY-SA版权
版权声明：本文为博主原创文章，遵循 CC 4.0 BY-SA 版权协议，转载请附上原文出处链接和本声明。
文章标签：

#fpga开发

于 2025-01-12 22:17:15 首次发布

一、原理

          本次实验需对图片进行进行灰度转换、高斯模糊、Sobel算子处理、非极大值抑制、滞后滤波、霍夫变换等步骤对车道进行识别。在此过程中，对像素矩阵的卷积占用了CPU的大量资源，而FPGA的并行特性可以加速这一处理，减轻CPU的负担，从而实现硬件加速。

二、操作步骤

2.1工具

          进行本次实验需要下载OpenCL SDK(OpenCL开发工具)、Intel SoC EDS(编译工具)、DE10-Nano OpenCL Board Support Package(板级支持包) 、diskimages(用于向sd卡传输linux镜像系统)、PuTTy(用于操作FPGA上的linux系统)。

OpenCL SDK下载如下：

Intel SoC EDS(编译工具)下载如下：

DE10-Nano OpenCL Board Support Package(板级支持包)

Diskimages

PuTTy

2.2环境配置

首先需要安装LICENSE文件，这里可以利用安装quatusII下的license.dat的环境变量

为了使操作系统正确找到 OpenCL 实用程序，用户需要创建一个名为 INTELFPGAOCLSDKROOT 的新环境变量，并将以下路径添加到 PATH 环境变量中：

添加成功后如下：

注意：这里要分行添加，不能全部放在一行否则不能正确加入

为了让 Intel OpenCL SDK 正确找到 DE10-Nano 的套件位置，用户需要创建一个环境变量 AOCL_BOARD_PACKAGE_ROOT 并将其值设置为：

设置成功后如下图所示：

2.3环境核对

为了验证我们配置的路径是否正确，我们可以再command shell窗口中输入acol version来看版本，如下图所示：

通过输入“aoc -list-boards”可以查看板子的信息，以及查看板级支持包的位置是否配对正确，如下图所示：

2.4编译并执行opencl程序

这里是将OpenCL生成的CL文件编译为FPGA的配置文件，在Command Shell中输入“aoc hough.cl -o bin\hough.aocx -board=de10_nano_sharedonly -v -report”，因为本次实验题目车道检测，所以使用的代码为hough.cl，需将其编译为hough.aocx，编译过程如下图所示：

注意：此处如果编译不成功，出现下面错误如下图所示

生成的aocx文件是即将配置到FPGA上的文件。

说明：此错误是由于某个目录没有权限，无法写入，解决方法有两个。

方法一：在一个较高的目录下，修改属性，加上直接去掉只读属性即可，选中包括子目录。

注意:此方法一定要选择最高路径，不然上层还会有有只读属性。

方法二：右键以管理员身份运行命令行窗口，这样就可以就可以以管理员的身份强行写入了。

2.5编译host文件夹生成linedetect

         打开embedded文件夹，启动Embedded_Command_Shell.bat文件，用“cd /cygdrive/D/QuartusII_18.1_Std/hld/board/de10_nano/test/Exercise_3/design_files/part2”转到有host的文件夹目录下，用“make”命令进行编译，在该目录下的bin文件夹中生成linedetect，如下图所示

2.6.用disk_image将linux系统镜像烧录到SD卡中

使用读卡器将sd卡插入电脑，打开disk_images软件，将镜像文件烧入卡中，如下图所示：

2.7将编译后的文件传入SD卡

在前面步骤中，make host文件夹生成了linedetect，编译hough.cl文件生成hough.aocx文件，将linedetect和hough.aocx文件以及test_images文件夹一起放入sd卡，供后续fpga处理。

2.8用PuTTy连接FPGA中的linux系统

连接方式选择serial，连接上FPGA的J4接口，在端口选项中选择COM5，波特率调至115200后启动该系统，如下图所示：

2.9在在PuTTy中完成实验

        在PuTTy中，先输入root登录，登录后输入“source ./init_opcl.sh”对OpenCL进行初始化，后建立boot文件夹，用“mount /dev/mmcblk0p1 boot”命令将该区挂在到boot上，输入“cd boot”，转入boot文件夹中，输入“ls”命令就可以看到该文件夹下的内容，包含hough.aocx、linedetect、test_images文件夹，后执行“./linedetect”，如下图所示：

执行结果如下所示，按照所给提示输入“./linedetect --img=../boot/1.bmp”对10张图片进行处理，如下图所示：

三、实验结果 

实验涉及的代码：

#include <stdio.h>
#include <stdlib.h>
#include <sys/time.h>

#define HYSTERESIS_HIGH_THRESH 45
#define HYSTERESIS_LOW_THRESH 15
	
// The dimensions of the image
int width, height;

double get_wall_time(){
    struct timeval time;
    if (gettimeofday(&time,NULL)){
        //  Handle error
        return 0;
    }
    return (double)time.tv_sec*1000 + (double)time.tv_usec /1000;
}

struct pixel {
	unsigned char b;
	unsigned char g;
	unsigned char r;
};

// Read BMP file and extract the pixel values (store in data) and header (store in header)
// data is data[0] = BLUE, data[1] = GREEN, data[2] = RED, etc...
int read_bmp(char* filename, unsigned char** header, struct pixel** data) {
	struct pixel * data_temp;
	unsigned char * header_temp;
	FILE* file = fopen(filename "rb");

	if (!file) return -1;
   
	// read the 54-byte header
	header_temp = (unsigned char *) malloc(54*sizeof(unsigned char));
	if (fread(header_temp, sizeof(unsigned char), 54, file) != 54){
		printf("Error reading BMP header\n");
		return -1;
	}

	// get height and width of image
	width = *(int*)&header_temp[18];
	height = *(int*)&header_temp[22];

	// Read in the image
	int size = width * height;
	data_temp = (struct pixel *) malloc(size*sizeof(struct pixel)); 
	if (fread(data_temp, sizeof(struct pixel), size, file) != size){
		printf("Error reading BMP image\n");
		return -1;
	}
	fclose(file);

	*header = header_temp;
	*data = data_temp;

	return 0;
}

void write_bmp(char* filename, unsigned char* header, struct pixel* data) {
	FILE* file = fopen(filename, "wb");

	// write the 54-byte header
	fwrite(header, sizeof(unsigned char), 54, file); 

	int size = width * height;
	fwrite(data, sizeof(struct pixel), size, file); 
	fclose(file);
}

// Write the grayscale image to disk.
void write_grayscale_bmp(char* filename, unsigned char* header, unsigned char* data) {
	FILE* file = fopen(filename, "wb");

	int size = width * height;
	struct pixel * data_temp = (struct pixel *) malloc(size*sizeof(struct pixel)); 

	// write the 54-byte header
	fwrite(header, sizeof(unsigned char), 54, file); 
	int y, x;
   
	// the r field of the pixel has the grayscale value. copy to g and b.
	for (y = 0; y < height; y++) {
		for (x = 0; x < width; x++) {
			(*(data_temp + y*width + x)).b = (*(data + y*width + x));
			(*(data_temp + y*width + x)).g = (*(data + y*width + x));
			(*(data_temp + y*width + x)).r = (*(data + y*width + x));
		}
	}
   
	size = width * height;
	fwrite(data_temp, sizeof(struct pixel), size, file); 

	free(data_temp);
	fclose(file);
}

// Determine the grayscale 8 bit value by averaging the r, g, and b channel values.
void convert_to_grayscale(struct pixel * data, unsigned char ** converted_data) {
	int i;
	*converted_data = (unsigned char*) malloc(width * height*sizeof(unsigned char));

	for (i = 0; i < width*height; i++) {
		(*converted_data)[i] = (data[i]r + data[i].g + data[i].b)/3;
	}
}

// Gaussian blur. 
void gaussian_blur(unsigned char ** data) {
	unsigned int gaussian_filter[5][5] = {
		{ 2, 4, 5, 4, 2 },
		{ 4, 9,12, 9, 4 },
		{ 5,12,15,12, 5 },
		{ 4, 9,12, 9, 4 },
		{ 2, 4, 5, 4, 2 }
	};
	int x, y, i, j;
	unsigned int numerator_r, denominator;
   
	unsigned char * temp_data = (unsigned char*) malloc(width * height*sizeof(unsigned char)); 
   
	for (y = 0; y < height; y++) {
		for (x = 0; x < width; x++) {
			numerator_r = 0;
			denominator = 0;
			for (j = -2; j <= 2; j++) {
				for (i = -2; i <= 2; i++) {
					if ( (x+i) >= 0 && (x+i) < width && (y+j) >= 0 && (y+j) < height) {
						numerator_r += (*data + (y+j)*width + (x+i)))*gaussian_filter[i+2][j+2];
						denominator += gaussian_filter[i+2][j+2];
					}
				}
			}
			(*(temp_data + y*width + x)) = numerator_r/denominator;
		}
	}
	free(*data);
	*data = temp_data;
}

void sobel_filter(unsigned char ** data) {
	// Definition of Sobel filter in horizontal direction
	int horizontal_operator[3][3] = {
		{ -1,  0,  1 },
		{ -2,  0,  2 },
		{ -1,  0,  1 }
	};
	int vertical_operator[3][3] = {
		{ -1,  -2,  -1 },
		{ 0,  0,  0 },
		{ 1,  2,  1 }
	};
                                
	int x, y, i, j;
	int horizontal_gradient, vertical_gradient;
   
	unsigned char * temp_data = (unsigned char*) malloc(width * height*sizeof(unsigned char)); 
   
	for (y = 0; y < height; y++) {
		for (x = 0; x < width; x++) {
			// Along the boundaries, set to 0
			if (y == 0 || x == 0 || y == height-1 || x == width-1) {
				(*(temp_data + y*width + x)) = 0;
				continue;
			}
			horizontal_gradient = 0;
			vertical_gradient = 0;
			for (j = -1; j <= 1; j++) {
				for (i = -1; i <= 1; i++) {
					if ((x+i) >= 0 && (x+i) < width && (y+j) >= 0 && (y+j) < height) {
						horizontal_gradient += (*(*data + (y+j)*width + (x+i)))*horizontal_operator[i+1][j+1];
						vertical_gradient += (*(*data + (y+j)*width + (x+i)))*vertical_operator[i+1][j+1];
					}
				}
			}
			// Check for overflow
			if ((abs(horizontal_gradient)+abs(vertical_gradient))/2 > 0xff)
				(*(temp_data + y*width + x)) = 0xff;
			else
				(*(temp_data + y*width + x)) = (abs(horizontal_gradient)+abs(vertical_gradient))/2;
		}
	}
	free(*data);
	*data = temp_data;
}

void non_maximum_suppressor(unsigned char ** data) {
   int x, y;
   
   unsigned char * temp_data = (unsigned char*) malloc(width * height*sizeof(unsigned char)); 
   
	for (y = 0; y < height; y++) {
		for (x = 0; x < width; x++) {
			// Along the boundaries, set to 0
			if (y == 0 || x == 0 || y == height-1 || x == width-1) {
				(*(temp_data + y*width + x)) = 0;
				continue;
			}
			unsigned int north_south = (*(*data + (y-1)*width + (x))) + (*(*data + (y+1)*width + (x)));
			unsigned int east_west = (*(*data + (y)*width + (x-1))) + (*(*data + (y)*width + (x+1)));
			unsigned int north_west = (*(*data + (y-1)*width + (x-1))) + (*(*data + (y+1)*width + (x+1)));
			unsigned int north_east = (*(*data + (y+1)*width + (x-1))) + (*(*data + (y-1)*width + (x+1)));
         
			(*(temp_data + (y)*width + (x))) = 0;
         
			if (north_south >= east_west && north_south >= north_west && north_south >= north_east) {
				if ((*(*data + (y)*width + (x))) > (*(*data + (y)*width + (x-1))) && 
					(*(*data + (y)*width + (x))) >= (*(*data + (y)*width + (x+1))))
				{
					(*(temp_data + (y)*width + (x))) = (*(*data + (y)*width + (x)));
				}
			} else if (east_west >= north_west && east_west >= north_east) {
				if ((*(*data + (y)*width + (x))) > (*(*data + (y-1)*width + (x))) && 
					(*(*data + (y)*width + (x))) >= (*(*data + (y+1)*width + (x))))
				{
					(*(temp_data + (y)*width + (x))) = (*(*data + (y)*width + (x)));
				}
			} else if (north_west >= north_east) {
				if ((*(*data + (y)*width + (x))) > (*(*data + (y-1)*width + (x+1))) && 
					(*(*data + (y)*width + (x))) >= (*(*data + (y+1)*width + (x-1))))
				{
					(*(temp_data + (y)*width + (x))) = (*(*data + (y)*width + (x)));
				}
			} else {
				if ((*(*data + (y)*width + (x))) > (*(*data + (y-1)*width + (x-1))) && 
					(*(*data + (y)*width + (x))) >= (*(*data + (y+1)*width + (x+1))))
				{
					(*(temp_data + (y)*width + (x))) = (*(*data + (y)*width + (x)));
				}
			}
		}
	}
	free(*data);
	*data = temp_data;
}

// Only keep pixels that are next to at least one strong pixel.
void hysteresis_filter(unsigned char ** data) {
	int x, y;

	unsigned char * temp_data = (unsigned char*) malloc(width * height*sizeof(unsigned char)); 

	for (y = 0; y < height; y++) {
		for (x = 0; x < width; x++) {
			// Along the boundaries, set to 0
			if (y == 0 || x == 0 || y == height-1 || x == width-1) {
				(*(temp_data + y*width + x)) = 0;
				continue;
			}
			
			// If pixel is strong or it is somewhat strong and at least one 
			// neighbouring pixel is strong, keep it. Otherwise zero it.
			if ((*(*data + (y)*width + (x))) > HYSTERESIS_HIGH_THRESH || 
				 ((*(*data + (y)*width + (x))) > HYSTERESIS_LOW_THRESH &&
				  ((*(*data + (y-1)*width + (x-1))) > HYSTERESIS_HIGH_THRESH ||
				  (*(*data + (y-1).*width + (x+0))) > HYSTERESIS_HIGH_THRESH ||
				  (*(*data + (y-1).*width + (x+1))) > HYSTERESIS_HIGH_THRESH ||
				  (*(*data + (y+0).*width + (x-1))) > HYSTERESIS_HIGH_THRESH ||
				  (*(*data + (y+0).*width + (x+1))) > HYSTERESIS_HIGH_THRESH ||
				  (*(*data + (y+1).*width + (x-1))) > HYSTERESIS_HIGH_THRESH ||
				  (*(*data + (y+1).*width + (x+0))) > HYSTERESIS_HIGH_THRESH ||
				  (*(*data + (y+1).*width + (x+1))) > HYSTERESIS_HIGH_THRESH))
			){
				(*(temp_data + (y)*width + (x))) = (*(*data + (y)*width + (x)));
			} else {
				(*(temp_data + (y)*width + (x))) = 0;
			}
		}
	}
	free(*data);
	*data = temp_data;
}

// Hough transform constants for a 720x540 image
#define ROWS 540
#define COLS 720
#define X_START -COLS/2
#define Y_START -ROWS/2
#define X_END COLS/2
#define Y_END 0 // Only sweep the bottom half of the image to ignore the horizon
#define RHO_RESOLUTION 2
#define RHOS (900/RHO_RESOLUTION) // How many values of rho do we go through? Sqrt(ROWS^2 + COLS^2) = 900.
#define THETAS 180 // How many values of theta do we go through? 

// sin and cos lookup values for thetas 0 to 179
float sinvals[THETAS] = {0.0, 0.01745240643728351, 0.03489949670250097, 0.05233595624294383, 0.0697564737441253, 0.08715574274765817, 0.10452846326765346, 0.12186934340514748, 0.13917310096006544, 0.15643446504023087, 0.17364817766693033, 0.1908089953765448, 0.20791169081775931, 0.224951054343865, 0.24192189559966773, 0.25881904510252074, 0.27563735581699916, 0.29237170472273677, 0.3090169943749474, 0.32556815445715664, 0.3420201433256687, 0.35836794954530027, 0.374606593415912, 0.3907311284892737, 0.40673664307580015, 0.42261826174069944, 0.4383711467890774, 0.45399049973954675, 0.4694715627858908, 0.48480962024633706, 0.49999999999999994, 0.5150380749100542, 0.5299192642332049, 0.5446390350150271, 0.5591929034707469, 0.573576436351046, 0.5877852522924731, 0.6018150231520483, 0.6156614753256582, 0.6293203910498374, 0.6427876096865393, 0.6560590289905072, 0.6691306063588582, 0.6819983600624985, 0.6946583704589973, 0.7071067811865475, 0.7193398003386511, 0.7313537016191705, 0.7431448254773941, 0.754709580222772, 0.766044443118978, 0.7771459614569708, 0.788010753606722, 0.7986355100472928, 0.8090169943749475, 0.8191520442889918, 0.8290375725550417, 0.8386705679454239, 0.848048096156426, 0.8571673007021122, 0.8660254037844386, 0.8746197071393957, 0.8829475928589269, 0.8910065241883678, 0.898794046299167, 0.9063077870366499, 0.9135454576426009, 0.9205048534524403, 0.9271838545667874, 0.9335804264972017, 0.9396926207859083, 0.9455185755993167, 0.9510565162951535, 0.9563047559630354, 0.9612616959383189, 0.9659258262890683, 0.9702957262759965, 0.9743700647852352, 0.9781476007338056, 0.981627183447664, 0.984807753012208, 0.9876883405951378, 0.9902680687415703, 0.992546151641322, 0.9945218953682733, 0.9961946980917455, 0.9975640502598242, 0.9986295347545738, 0.9993908270190958, 0.9998476951563913, 1.0, 0.9998476951563913, 0.9993908270190958, 0.9986295347545738, 0.9975640502598242, 0.9961946980917455, 0.9945218953682734, 0.9925461516413221, 0.9902680687415704, 0.9876883405951377, 0.984807753012208, 0.981627183447664, 0.9781476007338057, 0.9743700647852352, 0.9702957262759965, 0.9659258262890683, 0.9612616959383189, 0.9563047559630355, 0.9510565162951536, 0.9455185755993168, 0.9396926207859084, 0.9335804264972017, 0.9271838545667874, 0.9205048534524404, 0.913545457642601, 0.90630778703665, 0.8987940462991669, 0.8910065241883679, 0.8829475928589271, 0.8746197071393959, 0.8660254037844387, 0.8571673007021123, 0.8480480961564261, 0.8386705679454239, 0.8290375725550417, 0.819152044288992, 0.8090169943749475, 0.7986355100472927, 0.788010753606722, 0.777145961456971, 0.766044443118978, 0.7547095802227718, 0.7431448254773942, 0.7313537016191706, 0.7193398003386514, 0.7071067811865476, 0.6946583704589971, 0.6819983600624986, 0.6691306063588583, 0.6560590289905073, 0.6427876096865395, 0.6293203910498377, 0.6156614753256584, 0.6018150231520482, 0.5877852522924732, 0.5735764363510464, 0.5591929034707469, 0.544639035015027, 0.5299192642332049, 0.5150380749100544, 0.49999999999999994, 0.48480962024633717, 0.4694715627858911, 0.45399049973954686, 0.4383711467890773, 0.4226182617406995, 0.40673664307580043, 0.39073112848927416, 0.37460659341591224, 0.3583679495453002, 0.3420201433256689, 0.32556815445715703, 0.3090169943749475, 0.29237170472273705, 0.27563735581699966, 0.258819045102521, 0.24192189559966773, 0.22495105434386478, 0.20791169081775931, 0.19080899537654497, 0.17364817766693028, 0.15643446504023098, 0.13917310096006574, 0.12186934340514755, 0.10452846326765373, 0.08715574274765864, 0.06975647374412552, 0.05233595624294381, 0.0348994967025007, 0.01745240643728344};
float cosvals[THETAS] = {1.0, 0.9998476951563913, 0.9993908270190958, 0.9986295347545738, 0.9975640502598242, 0.9961946980917455, 0.9945218953682733, 0.992546151641322, 0.9902680687415704, 0.9876883405951378, 0.984807753012208, 0.981627183447664, 0.9781476007338057, 0.9743700647852352, 0.9702957262759965, 0.9659258262890683, 0.9612616959383189, 0.9563047559630354, 0.9510565162951535, 0.9455185755993168, 0.9396926207859084, 0.9335804264972017, 0.9271838545667874, 0.9205048534524404, 0.9135454576426009, 0.9063077870366499, 0.898794046299167, 0.8910065241883679, 0.882947592858927, 0.8746197071393957, 0.8660254037844387, 0.8571673007021123, 0.848048096156426, 0.838670567945424, 0.8290375725550416, 0.8191520442889918, 0.8090169943749475, 0.7986355100472928, 0.788010753606722, 0.7771459614569709, 0.766044443118978, 0.7547095802227721, 0.7431448254773942, 0.7313537016191706, 0.7193398003386512, 0.7071067811865476, 0.6946583704589974, 0.6819983600624985, 0.6691306063588582, 0.6560590289905073, 0.6427876096865394, 0.6293203910498375, 0.6156614753256583, 0.6018150231520484, 0.5877852522924731, 0.5735764363510462, 0.5591929034707468, 0.5446390350150272, 0.5299192642332049, 0.5150380749100544, 0.5000000000000001, 0.4848096202463371, 0.46947156278589086, 0.4539904997395468, 0.43837114678907746, 0.42261826174069944, 0.4067366430758002, 0.39073112848927394, 0.37460659341591196, 0.3583679495453004, 0.3420201433256688, 0.32556815445715676, 0.30901699437494745, 0.29237170472273677, 0.27563735581699916, 0.25881904510252074, 0.2419218955996679, 0.22495105434386492, 0.20791169081775945, 0.19080899537654492, 0.17364817766693041, 0.15643446504023092, 0.1391731009600657, 0.12186934340514749, 0.10452846326765346, 0.08715574274765814, 0.06975647374412546, 0.052335956242943966, 0.03489949670250108, 0.017452406437283376, 6.123233995736766e-17, -0.017452406437283477, -0.03489949670250073, -0.05233595624294362, -0.06975647374412533, -0.08715574274765824, -0.10452846326765333, -0.12186934340514737, -0.13917310096006535, -0.15643446504023104, -0.1736481776669303, -0.1908089953765448, -0.20791169081775912, -0.2249510543438648, -0.24192189559966779, -0.25881904510252085, -0.27563735581699905, -0.29237170472273666, -0.30901699437494734, -0.3255681544571564, -0.3420201433256687, -0.35836794954530027, -0.37460659341591207, -0.3907311284892736, -0.40673664307580004, -0.42261826174069933, -0.4383711467890775, -0.4539904997395467, -0.46947156278589053, -0.484809620246337, -0.4999999999999998, -0.5150380749100543, -0.5299192642332048, -0.5446390350150271, -0.5591929034707467, -0.5735764363510458, -0.587785252292473, -0.6018150231520484, -0.6156614753256583, -0.6293203910498373, -0.6427876096865394, -0.6560590289905075, -0.6691306063588582, -0.6819983600624984, -0.694658370458997, -0.7071067811865475, -0.7193398003386512, -0.7313537016191705, -0.743144825477394, -0.754709580222772, -0.7660444431189779, -0.7771459614569707, -0.7880107536067219, -0.7986355100472929, -0.8090169943749473, -0.8191520442889916, -0.8290375725550416, -0.8386705679454242, -0.848048096156426, -0.8571673007021122, -0.8660254037844387, -0.8746197071393957, -0.8829475928589268, -0.8910065241883678, -0.898794046299167, -0.9063077870366499, -0.9135454576426008, -0.9205048534524402, -0.9271838545667873, -0.9335804264972017, -0.9396926207859083, -0.9455185755993167, -0.9510565162951535, -0.9563047559630354, -0.9612616959383187, -0.9659258262890682, -0.9702957262759965, -0.9743700647852352, -0.9781476007338057, -0.981627183447664, -0.984807753012208, -0.9876883405951377, -0.9902680687415703, -0.992546151641322, -0.9945218953682733, -0.9961946980917455, -0.9975640502598242, -0.9986295347545738, -0.9993908270190958, -0.9998476951563913};

// data from BMP is stored with the bottom-left pixel first (origin is at bottom left)
void hough_transform(unsigned char * data, unsigned short ** accumulator) {
	
	int r,t,x,y;
	
	if (height != ROWS || width != COLS){
		printf("ERROR: Invalid image dimensions for the hough transform\n");
		return;
	}
	
	unsigned short * temp_accumulator = (unsigned short*) malloc(RHOS * THETAS * sizeof(unsigned short)); 
	
	for (r = 0; r < RHOS; r++)
		for (t = 0; t < THETAS; t++)
			temp_accumulator[r*THETAS + t] = 0;
	
	for (y = Y_START; y < Y_END; y++){
		for (x = X_START; x < X_END; x++){
			if (data[(y-Y_START)*COLS + x-X_START] != 0){
				for (t = 0; t < THETAS; t++){
					r = (x/RHO_RESOLUTION)*cosvals[t] + (y/RHO_RESOLUTION)*sinvals[t];
					temp_accumulator[(r+RHOS/2)*THETAS + t] += 1;
				}
			}
		}
	}
	
	*accumulator = temp_accumulator;
}

// A line (rho,theta)
struct line {
	short rho;
	unsigned short theta;
};

// Find num_lines lines with the highest accumulated values (local maxima)
void extract_lines(unsigned short * accumulator, struct line ** lines, int num_lines) {
	int r, t, i, r_off, t_off;
	struct line * temp_lines = (struct line *) malloc(num_lines * sizeof(struct line));

	// Some vars to keep running list of strongest lines as we go through accumulator
	unsigned short* strongest_lines = (unsigned short*) malloc(num_lines * sizeof(unsigned short));
	unsigned short curr_weakest = 0;
	unsigned short curr_weakest_idx = 0;
	for (i = 0; i < num_lines; i++) strongest_lines[i] = 0;
	
	for (r = -RHOS/2; r < RHOS/2; r++){
		for (t = 0; t < THETAS; t++){
			
			// Filter out very horizontal lines
			#define DEGREES_FROM_HORIZON 10
			if (t > 90-DEGREES_FROM_HORIZON && t < 90+DEGREES_FROM_HORIZON) continue;
			
			if (accumulator[(r+RHOS/2)*THETAS + t] > curr_weakest){
				
				int is_local_maximum = 1; 
				
				// Replace previous weakest line with newly found line IF
				// the newly found line is a local maximum
				for (r_off = -2; r_off <= 2; r_off++){
					for (t_off = -2; t_off <= 2; t_off++){
						if (r + r_off < -RHOS/2 || r + r_off >= RHOS/2 || t + t_off < 0 || t + t_off >= THETAS) continue;
						
						if (accumulator[(r+RHOS/2)*THETAS + t] < accumulator[(r+r_off+RHOS/2)*THETAS + t+t_off]) is_local_maximum = 0;
					}
				}
				
				
				if (is_local_maximum){
				
					strongest_lines[curr_weakest_idx] = accumulator[(r+RHOS/2)*THETAS + t];
					temp_lines[curr_weakest_idx].rho = r;
					temp_lines[curr_weakest_idx].theta = t;

					// Find the new weakest line
					curr_weakest = strongest_lines[0];
					curr_weakest_idx = 0;
					for (i = 0; i < num_lines; i++){
						if (strongest_lines[i] < curr_weakest){
							curr_weakest = strongest_lines[i];
							curr_weakest_idx = i;
						}
					}
				}
			}
		}
	}
	
	*lines = temp_lines;
}

// Functions for finding coordinates x,y coordinates from a rho,theta line. 
// These x,y coordinates correspond to the coordinate system of the original 720x540 image,
// which as its origin at the bottom left of the image.
int find_x(int y_bottomleft_origin, struct line l) {
	// Using the coordinate system with origin at center of image, rho = xcostheta + ysintheta
	int y_center_origin = y_bottomleft_origin - height/2;
	int x_center_origin = (l.rho*RHO_RESOLUTION - y_center_origin*sinvals[l.theta])/cosvals[l.theta];
	int x_bottomleft_origin = x_center_origin + width/2;
	return x_bottomleft_origin >= 0 && x_bottomleft_origin < width ? x_bottomleft_origin : -1;
}
int find_y(int x_bottomleft_origin, struct line l) {
	// Using the coordinate system with origin at center of image, rho = xcostheta + ysintheta
	int x_center_origin = x_bottomleft_origin - width/2;
	int y_center_origin = (l.rho*RHO_RESOLUTION - x_center_origin*cosvals[l.theta])/sinvals[l.theta];
	int y_bottomleft_origin = y_center_origin + height/2;
	return y_bottomleft_origin >= 0 && y_bottomleft_origin < height ? y_bottomleft_origin : -1;
}

// Overlay the lines onto the original image.
void overlay_lines(struct pixel * data, struct line * lines, int num_lines) {
	
	int i, y, x;
	
	// Overlay the lines
	for (i = 0; i < num_lines; i++){
		// Determine the points that exist on this line.
		// If the line is more vertical than horizontal, sweep across the y
		if (lines[i].theta < THETAS/4 || lines[i].theta >= 3*THETAS/4){
			for (y = 0; y < height; y++){
				x = find_x(y, lines[i]);
				if (x >= 0){
					data[y*width + x].r = 255;
					data[y*width + x].g = 0;
					data[y*width + x].b = 0;
				}
			}
		// If the line is more horizontal than vertical, sweep across the x
		} else if (lines[i].theta >= THETAS/4 && lines[i].theta < 3*THETAS/4){
			for (x = 0; x < width; x++){
				y = find_y(x, lines[i]);
				if (y >= 0){
					data[y*width + x].r = 255;
					data[y*width + x].g = 0;
					data[y*width + x].b = 0;
				}
			}
		}
	}
}

int main(int argc, char *argv[]) {
	struct pixel * rgb_data;
	unsigned char * data;
	unsigned char * header;
	unsigned short * accumulator = NULL;
	struct line * lines;
	int debug = 0;
	double start, end, func_start, hough_runtime, edge_runtime;
	char outputfile [256];
	int num_lines;

	// Check inputs
	if (argc < 4) {
		printf("Usage: edgedetect <BMP filename> <debug> <num lines to find>\n\ndebug: true false\n");
		return 0;
	}

	// Determine debug mode
	if (argv[2][0] == 't') debug = 1;

	// Parse number of lines to find
	num_lines = atoi(argv[3]);

	// Open input image file (24-bit bitmap image)
	if (read_bmp(argv[1], &header, &rgb_data) < 0) {
		printf("Failed to read BMP\n");
		return 0;
	}

	/********************************************
	*          IMAGE PROCESSING STAGES          *
	********************************************/

	// Start measuring time
	start = get_wall_time();
	
	/// Grayscale conversion
	convert_to_grayscale(rgb_data, &data);
	if (debug) write_grayscale_bmp("stage0_grayscale.bmp", header, data);

	/// Gaussian filter
	gaussian_blur(&data);
	if (debug) write_grayscale_bmp("stage1_gaussian.bmp", header, data);

	/// Sobel operator
	sobel_filter(&data);
	if (debug) write_grayscale_bmp("stage2_sobel.bmp", header, data);

	/// Non-maximum suppression
	non_maximum_suppressor(&data);
	if (debug) write_grayscale_bmp("stage3_nonmax_suppression.bmp", header, data);

	/// Hysteresis
	hysteresis_filter(&data);
	if (debug) write_grayscale_bmp("stage4_hysteresis.bmp", header, data);
	
	edge_runtime = get_wall_time() - start;
	
	/// Hough
	func_start = get_wall_time();
	hough_transform(data, &accumulator); 
	hough_runtime = get_wall_time() - func_start;
	
	/// Find the strongest num_lines lines
	extract_lines(accumulator, &lines, num_lines);
	
	// Stop measuring time.
	end = get_wall_time();
	
	printf("EDGE DETECTION TIME ELAPSED: %.2f ms\n", edge_runtime);
	printf("HOUGH TRANSFORM TIME ELAPSED: %.2f ms\n", hough_runtime);
	printf("TOTAL TIME ELAPSED: %2f ms\n", end - start);
	
	/// Overlay the detected lines onto the original image
	overlay_lines(rgb_data, lines, num_lines);
	
	// Write out the original image with detected lines overlaid
	snprintf(outputfile, 256, "%s_part1_lines.bmp",argv[1]);
	printf("Writing line-detected image to %s\n",outputfile);
	write_bmp(outputfile, header, rgb_data);

	return 0;
}

声明：作者水平有限，若有错误请在评论区指出。

吸纹鸽
关注
使用FPGA模块进行硬件加速
MZEing的博客
09-14 869
FPGA的编程是通过硬件描述语言（HDL）进行的，其中最常用的是VHDL（VHSIC硬件描述语言）和Verilog。开发人员使用HDL编写硬件描述代码，然后将其综合为FPGA的配置位流文件（bitstream），通过编程器将其加载到FPGA芯片中。本文提供了一个简单的向量加法示例代码，希望能帮助读者了解如何使用FPGA模块进行硬件加速。下面以一个简单的向量加法为例，介绍如何使用FPGA模块进行硬件加速。此时，我们可以通过软件编程来调用FPGA模块，并将数据发送到FPGA进行加速计算。
FPGA硬件加速综述
jedibobo的博客
03-20 3640
FPGA硬件加速综述 背景 本科期间参加了不少竞赛，其中和FPGA相关的主要是集创赛以及Xilinx全球自适应计算大赛两个。虽然接触到的硬件类型五花八门（ZYNQ ZYNQ Ultrascale+ Alveo ACAP等），但是主要的学习都集中在深度学习尤其是目标检测和分类算法在FPGA上的部署。 概述 FPGA用于硬件加速主要有以下优势： 可以通过优化提高并行度，减小推理的时延 可以利用丰富的I/O做一些特殊的应用 劣势其实也是优势导致的： 优化程度两极分化，需要很高的工程技术 FPGA在相
图像处理：卷积模块FPGA 硬件加速
weixin_34248118的博客
04-22 1550
本文记录了利用FPGA加速图像处理中的卷积计算的设计与实现。实现环境为Altera公司的Cyclone IV型芯片，NIOS II软核+FPGA架构。 由于这是第一次设计硬件加速模块，设计中的瑕疵以及问题欢迎前来讨论。 更新记录： D0423 记录FPGA核心计算模块和控制模块 D0426 记录FPGA核心计算模块的控制驱动，性能与功能测试 Part1 : 卷积相...
Altera Stratix V FPGA硬件加速平台
qq_39409559的博客
09-22 1288
一、板卡主要硬件参数： 1.FPGA：Altera Stratix® V GX FPGA 5SGXEA7N2F45C2N 2.FPGA的配置 1）板载USB Blaster下载线或者JTAG口 2）基于MAXII CPLD和并行FLASH的Fast passive parallel(FPPx32)配置  3.时钟系统 1）一路有源晶振时钟 2）两片可定制时钟芯片IDT5P49V5
基于5片FPGA处理的PCIE光纤传输硬件加速卡
beijingqingyikeji的博客
09-21 2051
        PCIE720是一款基于PCI Express总线架构的高性能计算（HPC）硬件加速卡，板卡采用Xilinx的高性能28nm 7系列FPGA作为运算节点，在资源、接口以及时钟的优化，为高性能计算提供卓越的硬件加速性能。板卡一共具有5个FPGA处理节点，其中主处理节点主要完成PCIe接口转换、光纤收发以及数据分发同步等功能，主从节点之间通过高速串行总线进行互联。          ...
卷积神经网络的FPGA硬件加速——QMJ
微核实控的博客
02-22 2863
一、本次实验目的： 1）学习卷积神经网络； 2）FPGA中实现卷积运算。 二、脉动阵列结构设计    PE（处理单元）的内部结构设计如图1所示，PE完成乘加运算得到的结果保存于寄存器中，其中IN和W表示1个字节（8bit）的输入和权重数据传递到PE。 图1 PE结构   图2中，W、IN和OUT表示阵列的输入和输出，脉动阵列一行4个PE单元组成一组。   脉动阵列PE得到的结果传递到out1，out2，out3的步骤：   第一步：PE单元完成4次乘加运算；   第二步：将一组PE得到的结果，从左到
基于FPGA的车道线检测：SEU-Xilinx暑期学校项目设计
2019年SEU-Xilinx国际暑期学校的团队项目设计文档详细介绍了基于FPGA的车道线检测系统的设计。该项目由洪刚、王琳、刘海军和白鹏程四位组员合作完成，他们利用PYNQ-Z2平台的强大功能，实现了视频车道线检测功能。...
【MATLAB与FPGA加速自动驾驶】：硬件加速与视觉识别技术的结合
[【MATLAB与FPGA加速自动驾驶】：硬件加速与视觉识别技术的结合](https://img-blog.csdnimg.cn/img_convert/22da07b6543b5ed8eeff5b7f02b99052.png) # 1. MATLAB与FPGA在自动驾驶中的基础应用 MATLAB（Matrix ...
基于ARM FPGA的心电信号疲劳驾驶检测方法的研究.pdf
07-13
FPGA在信号处理、数据存储和实时系统中使用普遍，因其能够实现高性能的并行处理而被广泛应用于各种硬件加速器中。 3. 心电信号（ECG）：心电信号是心脏电活动的生物电位变化的记录，通常通过在人体皮肤上放置电极来...
【硬件加速的秘密】：如何用FPGA实现图像处理的革命性突破
![【硬件加速的秘密】：如何用FPGA实现图像处理的革命性突破]...# 1. FPGA与图像处理的革命性突破 ## 1.1 前言：FPGA在现代技术中的崛起 在数字世界的快速发展中，现场可编程门阵列（F
FPGA硬件加速平台
12-16
本平台是一款基于PCI Express总线架构的高性能FPGA算法加速卡，该FPGA加速板卡采用Xilinx的高性能7系列FPGA作为运算节点。
如何用FPGA实现算法的硬件加速
08-03
当设计者试图从算法中获得最佳性能但软件方法已无计可施时，可以尝试通过硬件/软件重新划分来进行加速。FPGA易于实现软件模块和硬件模块的相互交换，且不必改变处理器或进行板级变动。本文阐述如何用FPGA来实现算法的硬件加速。
手机ADAS：基于OpenCL的车道线检测应用评估
10-16
【OpenCL】：Open Computing Language (OpenCL) 是一种跨平台的并行编程框架，它允许程序员利用各种设备，如GPU、CPU、FPGA等，进行高性能计算。OpenCL 提供了一个标准化的接口，使开发者能够编写能够在不同硬件上...
【人工智能】FPGA实现人工智能算法硬件加速学习笔记
最新发布
qq_45355365的博客
08-03 2476
【人工智能】FPGA实现人工智能算法硬件加速学习笔记
FPGA开发：实现高性能硬件加速
CqpFsharp的博客
09-22 387
在将HDL代码烧录到FPGA之前，您可以进行仿真和综合来验证设计的功能和性能。仿真可用于检查逻辑电路的正确性，而综合则将HDL代码转换为FPGA上的物理布局和逻辑电路。请注意，以上代码示例仅用于说明目的，实际的FPGA开发可能涉及更复杂的设计和算法。在进行FPGA开发时，建议参考相关的文档、教程和示例代码，以获得更详细和全面的指导。FPGA（现场可编程门阵列）是一种灵活且可重构的硬件平台，广泛应用于高性能计算、信号处理和加速任务。将该文件烧录到FPGA开发板上，以在硬件上运行您的设计。
fpga的硬件加速原理
僧哥叨叨叨
10-21 3314
fpga做硬件加速
软件算法在FPGA上实现边缘计算的硬件加速：综合工具
JhzDev的博客
09-19 552
通过使用综合工具，可以将软件算法转换为FPGA上的硬件描述语言，并最终实现在FPGA上的硬件加速。当然，上述示例只是一个简单的演示，实际应用中可能涉及更复杂的算法和更多的优化技术。但总体而言，软件算法在FPGA上实现硬件加速的综合工具为边缘计算提供了强大的计算能力和灵活性，为各种应用场景带来了巨大的潜力和机会。然后，将生成的比特流文件加载到FPGA芯片中，配置FPGA以实现向量加法的硬件加速。下面是一个简单的示例，展示了如何使用综合工具将一个简单的软件算法转换为FPGA上的硬件加速模块。
评论
添加红包

请填写红包祝福语或标题

个

红包个数最小为10个

元

红包金额最低5元

当前余额3.43元 前往充值 >
需支付：10.00元
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付元
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明：

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

余额充值