Никоноров И.Д. ЛР№3 Линейная фильтрация изображений (горизонтальное разбиение). Ядро Гаусса 3x3. TBB

1706-4/nikonorov_id/TBB/TBB.sln

@@ -0,0 +1,31 @@
+
+Microsoft Visual Studio Solution File, Format Version 12.00
+# Visual Studio Version 16
+VisualStudioVersion = 16.0.29519.181
+MinimumVisualStudioVersion = 10.0.40219.1
+Project("{8BC9CEB8-8B4A-11D0-8D11-00A0C91BC942}") = "TBB", "TBB\TBB.vcxproj", "{01ABCEF3-C146-4904-B53C-D353C62CC69D}"
+EndProject
+Global
+	GlobalSection(SolutionConfigurationPlatforms) = preSolution
+		Debug|x64 = Debug|x64
+		Debug|x86 = Debug|x86
+		Release|x64 = Release|x64
+		Release|x86 = Release|x86
+	EndGlobalSection
+	GlobalSection(ProjectConfigurationPlatforms) = postSolution
+		{01ABCEF3-C146-4904-B53C-D353C62CC69D}.Debug|x64.ActiveCfg = Debug|x64
+		{01ABCEF3-C146-4904-B53C-D353C62CC69D}.Debug|x64.Build.0 = Debug|x64
+		{01ABCEF3-C146-4904-B53C-D353C62CC69D}.Debug|x86.ActiveCfg = Debug|Win32
+		{01ABCEF3-C146-4904-B53C-D353C62CC69D}.Debug|x86.Build.0 = Debug|Win32
+		{01ABCEF3-C146-4904-B53C-D353C62CC69D}.Release|x64.ActiveCfg = Release|x64
+		{01ABCEF3-C146-4904-B53C-D353C62CC69D}.Release|x64.Build.0 = Release|x64
+		{01ABCEF3-C146-4904-B53C-D353C62CC69D}.Release|x86.ActiveCfg = Release|Win32
+		{01ABCEF3-C146-4904-B53C-D353C62CC69D}.Release|x86.Build.0 = Release|Win32
+	EndGlobalSection
+	GlobalSection(SolutionProperties) = preSolution
+		HideSolutionNode = FALSE
+	EndGlobalSection
+	GlobalSection(ExtensibilityGlobals) = postSolution
+		SolutionGuid = {4CC8E9E7-4985-4064-B28E-C630931FE0C0}
+	EndGlobalSection
+EndGlobal

1706-4/nikonorov_id/TBB/TBB/Gauss_tbb.cpp

@@ -0,0 +1,198 @@
+#include <iostream>
+#include <stdio.h>
+#include <ctime>
+#include "tbb/parallel_for.h"
+#include "tbb/blocked_range2d.h"
+#include "tbb/task_scheduler_init.h"
+#include <omp.h>
+#include <opencv2/opencv.hpp>
+#include "opencv2/core/core.hpp"
+#include "opencv2/highgui/highgui.hpp"
+#include "opencv2/imgproc/imgproc.hpp"
+#include <chrono>
+
+using namespace cv;
+using namespace tbb;
+using namespace std;
+
+#define M_PI 3.14159265358979323846
+
+int Clamp(int value, int min, int max)
+{
+	if (value < min)
+		return min;
+
+	if (value > max)
+		return max;
+
+	return value;
+}
+
+void deletekernel(double* kernel, int k) {
+
+	delete[] kernel;
+}
+
+
+//-------------------------------------------------------------//
+//-------------------------------------------------------------//
+//-------------------------------------------------------------//
+//-------------------------------------------------------------//
+//-------------------------------------------------------------//
+
+double* createKernel(int rows, int cols) {
+	int tmp = rows * cols;
+	double* mas;
+	mas = new double[tmp];
+	return mas;
+}
+void InitKern(double* kernel, int k, double sigma)
+{
+	int radius = (k / 2);
+	double sum = 0.0;
+	double tmp = 1 / (2 * M_PI * sigma * sigma);
+	int step = radius;
+
+	for (int i = -radius; i <= radius; i++)
+	{
+		for (int j = -radius; j <= radius; j++)
+		{
+			kernel[j + step] = tmp * exp(-(i * i + j * j) / (2 * sigma * sigma));
+			sum += kernel[j + step];
+		}
+		step += k;
+	}
+
+	for (int i = 0; i < k * k; i++)
+	{
+		kernel[i] /= sum;
+	}
+}
+
+void mat_cols_mult(double* kernel, const Mat& input, Mat output, int i, int radius)
+{
+	for (size_t j = 1; j < input.cols; j++)
+	{
+		double tmp = 0;
+		for (int l = -radius; l <= radius; l++)
+			for (int w = -radius; w <= radius; w++)
+			{
+				int idX = Clamp(i + l, 0, input.rows - 1);
+				int idY = Clamp(j + w, 0, input.cols - 1);
+				tmp += double(input.at<uchar>(idX, idY)) * kernel[l + w + 2 * (l + 2)];
+			}
+		output.at<uchar>(i, j) = int(tmp);
+	}
+}
+
+
+void parallel_matrix_multiply(double* kernel, const Mat& input, Mat& output, int k)
+{
+	int radius = k / 2;
+	parallel_for(blocked_range<size_t>(1, input.rows), [=](const blocked_range<size_t>& r)
+		{
+			for (size_t i = r.begin(); i != r.end(); ++i)
+			{
+				mat_cols_mult(kernel, input, output, i, radius);
+			}
+		});
+}
+
+void Sequential_Gauss(const Mat& input, Mat& output, double* kernel, int kern_size)
+{
+	int radius = int(kern_size / 2);
+	for (int i = 1; i < input.rows; i++) {
+		for (int j = 1; j < input.cols; j++) {
+			double tmp = 0.0;
+			for (int l = -radius; l <= radius; l++) {
+				for (int w = -radius; w <= radius; w++) {
+					int idX = Clamp(i + l, 0, input.rows - 1);
+					int idY = Clamp(j + w, 0, input.cols - 1);
+					tmp += double(input.at<uchar>(idX, idY)) * kernel[l + w + 2 * (l + 2)];
+				}
+			}
+			output.at<uchar>(i, j) = int(tmp);
+		}
+	}
+}
+
+bool check(const Mat& image1, const Mat& image2)
+{
+	Mat res;
+	bitwise_xor(image1, image2, res);
+	if (countNonZero(res) > 0)
+	{
+		cout << res;
+		return false;
+	}
+	else
+		return true;
+}
+
+int main(int argc, char** argv) {
+
+	string orig_image_path;
+	string result_image_path;
+	int thread_nums;
+
+	if (argc < 4)
+	{
+		orig_image_path += "pic.jpg";
+		result_image_path += "result.jpg";
+		thread_nums = 6;
+	}
+	else
+	{
+		orig_image_path += argv[1];
+		result_image_path += argv[2];
+	}
+
+	Mat original_image;
+	original_image = imread(orig_image_path, IMREAD_COLOR);
+
+	if (original_image.empty())
+	{
+		cout << "No image data" << endl;
+		return -1;
+	}
+
+	Mat gray_image;
+	Mat seq_result;
+	Mat TBB_result;
+	cvtColor(original_image, gray_image, COLOR_BGR2GRAY);
+	cvtColor(original_image, seq_result, COLOR_BGR2GRAY);
+	cvtColor(original_image, TBB_result, COLOR_BGR2GRAY);
+
+	int k = 3;
+
+	double* kernel;
+	kernel = createKernel(k, k);
+	InitKern(kernel, k, 7);
+
+	auto begin = std::chrono::steady_clock::now();
+	task_scheduler_init init(thread_nums);
+	parallel_matrix_multiply(kernel, gray_image, TBB_result, k);
+	auto end = std::chrono::steady_clock::now();
+	auto elapsed_ms = std::chrono::duration_cast<std::chrono::milliseconds>(end - begin);
+	cout << "The TBB time is: " << elapsed_ms.count() << " ms\n";
+
+	auto s_begin = std::chrono::steady_clock::now();
+	Sequential_Gauss(gray_image, seq_result, kernel, k);
+	auto s_end = std::chrono::steady_clock::now();
+	auto s_elapsed_ms = std::chrono::duration_cast<std::chrono::milliseconds>(s_end - s_begin);
+	cout << "The sequential time is: " << s_elapsed_ms.count() << " ms" << endl;
+
+	if (check(seq_result, TBB_result))
+	{
+		imwrite(result_image_path, TBB_result);
+		namedWindow("my_result", WINDOW_AUTOSIZE);
+		namedWindow("Gray image", WINDOW_AUTOSIZE);
+		imshow("my_result", TBB_result);
+		imshow("Gray image", gray_image);
+		waitKey(0);
+	}
+
+	deletekernel(kernel, k);
+
+	return 0;
+}

1706-4/nikonorov_id/TBB/TBB/TBB.vcxproj

@@ -0,0 +1,149 @@
+<?xml version="1.0" encoding="utf-8"?>
+<Project DefaultTargets="Build" xmlns="http://schemas.microsoft.com/developer/msbuild/2003">
+  <ItemGroup Label="ProjectConfigurations">
+    <ProjectConfiguration Include="Debug|Win32">
+      <Configuration>Debug</Configuration>
+      <Platform>Win32</Platform>
+    </ProjectConfiguration>
+    <ProjectConfiguration Include="Release|Win32">
+      <Configuration>Release</Configuration>
+      <Platform>Win32</Platform>
+    </ProjectConfiguration>
+    <ProjectConfiguration Include="Debug|x64">
+      <Configuration>Debug</Configuration>
+      <Platform>x64</Platform>
+    </ProjectConfiguration>
+    <ProjectConfiguration Include="Release|x64">
+      <Configuration>Release</Configuration>
+      <Platform>x64</Platform>
+    </ProjectConfiguration>
+  </ItemGroup>
+  <PropertyGroup Label="Globals">
+    <VCProjectVersion>16.0</VCProjectVersion>
+    <ProjectGuid>{01ABCEF3-C146-4904-B53C-D353C62CC69D}</ProjectGuid>
+    <RootNamespace>TBB</RootNamespace>
+    <WindowsTargetPlatformVersion>10.0</WindowsTargetPlatformVersion>
+  </PropertyGroup>
+  <Import Project="$(VCTargetsPath)\Microsoft.Cpp.Default.props" />
+  <PropertyGroup Condition="'$(Configuration)|$(Platform)'=='Debug|Win32'" Label="Configuration">
+    <ConfigurationType>Application</ConfigurationType>
+    <UseDebugLibraries>true</UseDebugLibraries>
+    <PlatformToolset>v142</PlatformToolset>
+    <CharacterSet>Unicode</CharacterSet>
+  </PropertyGroup>
+  <PropertyGroup Condition="'$(Configuration)|$(Platform)'=='Release|Win32'" Label="Configuration">
+    <ConfigurationType>Application</ConfigurationType>
+    <UseDebugLibraries>false</UseDebugLibraries>
+    <PlatformToolset>v142</PlatformToolset>
+    <WholeProgramOptimization>true</WholeProgramOptimization>
+    <CharacterSet>Unicode</CharacterSet>
+  </PropertyGroup>
+  <PropertyGroup Condition="'$(Configuration)|$(Platform)'=='Debug|x64'" Label="Configuration">
+    <ConfigurationType>Application</ConfigurationType>
+    <UseDebugLibraries>true</UseDebugLibraries>
+    <PlatformToolset>v142</PlatformToolset>
+    <CharacterSet>Unicode</CharacterSet>
+  </PropertyGroup>
+  <PropertyGroup Condition="'$(Configuration)|$(Platform)'=='Release|x64'" Label="Configuration">
+    <ConfigurationType>Application</ConfigurationType>
+    <UseDebugLibraries>false</UseDebugLibraries>
+    <PlatformToolset>v142</PlatformToolset>
+    <WholeProgramOptimization>true</WholeProgramOptimization>
+    <CharacterSet>Unicode</CharacterSet>
+  </PropertyGroup>
+  <Import Project="$(VCTargetsPath)\Microsoft.Cpp.props" />
+  <ImportGroup Label="ExtensionSettings">
+  </ImportGroup>
+  <ImportGroup Label="Shared">
+  </ImportGroup>
+  <ImportGroup Label="PropertySheets" Condition="'$(Configuration)|$(Platform)'=='Debug|Win32'">
+    <Import Project="$(UserRootDir)\Microsoft.Cpp.$(Platform).user.props" Condition="exists('$(UserRootDir)\Microsoft.Cpp.$(Platform).user.props')" Label="LocalAppDataPlatform" />
+  </ImportGroup>
+  <ImportGroup Label="PropertySheets" Condition="'$(Configuration)|$(Platform)'=='Release|Win32'">
+    <Import Project="$(UserRootDir)\Microsoft.Cpp.$(Platform).user.props" Condition="exists('$(UserRootDir)\Microsoft.Cpp.$(Platform).user.props')" Label="LocalAppDataPlatform" />
+  </ImportGroup>
+  <ImportGroup Label="PropertySheets" Condition="'$(Configuration)|$(Platform)'=='Debug|x64'">
+    <Import Project="$(UserRootDir)\Microsoft.Cpp.$(Platform).user.props" Condition="exists('$(UserRootDir)\Microsoft.Cpp.$(Platform).user.props')" Label="LocalAppDataPlatform" />
+  </ImportGroup>
+  <ImportGroup Label="PropertySheets" Condition="'$(Configuration)|$(Platform)'=='Release|x64'">
+    <Import Project="$(UserRootDir)\Microsoft.Cpp.$(Platform).user.props" Condition="exists('$(UserRootDir)\Microsoft.Cpp.$(Platform).user.props')" Label="LocalAppDataPlatform" />
+  </ImportGroup>
+  <PropertyGroup Label="UserMacros" />
+  <PropertyGroup Condition="'$(Configuration)|$(Platform)'=='Debug|Win32'">
+    <LinkIncremental>true</LinkIncremental>
+  </PropertyGroup>
+  <PropertyGroup Condition="'$(Configuration)|$(Platform)'=='Debug|x64'">
+    <LinkIncremental>true</LinkIncremental>
+  </PropertyGroup>
+  <PropertyGroup Condition="'$(Configuration)|$(Platform)'=='Release|Win32'">
+    <LinkIncremental>false</LinkIncremental>
+  </PropertyGroup>
+  <PropertyGroup Condition="'$(Configuration)|$(Platform)'=='Release|x64'">
+    <LinkIncremental>false</LinkIncremental>
+  </PropertyGroup>
+  <ItemDefinitionGroup Condition="'$(Configuration)|$(Platform)'=='Debug|Win32'">
+    <ClCompile>
+      <WarningLevel>Level3</WarningLevel>
+      <SDLCheck>true</SDLCheck>
+      <PreprocessorDefinitions>_DEBUG;_CONSOLE;%(PreprocessorDefinitions)</PreprocessorDefinitions>
+      <ConformanceMode>true</ConformanceMode>
+    </ClCompile>
+    <Link>
+      <SubSystem>Console</SubSystem>
+      <GenerateDebugInformation>true</GenerateDebugInformation>
+    </Link>
+  </ItemDefinitionGroup>
+  <ItemDefinitionGroup Condition="'$(Configuration)|$(Platform)'=='Debug|x64'">
+    <ClCompile>
+      <WarningLevel>Level3</WarningLevel>
+      <SDLCheck>true</SDLCheck>
+      <PreprocessorDefinitions>_DEBUG;_CONSOLE;%(PreprocessorDefinitions)</PreprocessorDefinitions>
+      <ConformanceMode>true</ConformanceMode>
+    </ClCompile>
+    <Link>
+      <SubSystem>Console</SubSystem>
+      <GenerateDebugInformation>true</GenerateDebugInformation>
+    </Link>
+  </ItemDefinitionGroup>
+  <ItemDefinitionGroup Condition="'$(Configuration)|$(Platform)'=='Release|Win32'">
+    <ClCompile>
+      <WarningLevel>Level3</WarningLevel>
+      <FunctionLevelLinking>true</FunctionLevelLinking>
+      <IntrinsicFunctions>true</IntrinsicFunctions>
+      <SDLCheck>true</SDLCheck>
+      <PreprocessorDefinitions>NDEBUG;_CONSOLE;%(PreprocessorDefinitions)</PreprocessorDefinitions>
+      <ConformanceMode>true</ConformanceMode>
+    </ClCompile>
+    <Link>
+      <SubSystem>Console</SubSystem>
+      <EnableCOMDATFolding>true</EnableCOMDATFolding>
+      <OptimizeReferences>true</OptimizeReferences>
+      <GenerateDebugInformation>true</GenerateDebugInformation>
+    </Link>
+  </ItemDefinitionGroup>
+  <ItemDefinitionGroup Condition="'$(Configuration)|$(Platform)'=='Release|x64'">
+    <ClCompile>
+      <WarningLevel>Level3</WarningLevel>
+      <FunctionLevelLinking>true</FunctionLevelLinking>
+      <IntrinsicFunctions>true</IntrinsicFunctions>
+      <SDLCheck>true</SDLCheck>
+      <PreprocessorDefinitions>NDEBUG;_CONSOLE;%(PreprocessorDefinitions)</PreprocessorDefinitions>
+      <ConformanceMode>true</ConformanceMode>
+      <AdditionalIncludeDirectories>C:\OpenCV\opencv\build\include;C:\tbb\include;%(AdditionalIncludeDirectories)</AdditionalIncludeDirectories>
+    </ClCompile>
+    <Link>
+      <SubSystem>Console</SubSystem>
+      <EnableCOMDATFolding>true</EnableCOMDATFolding>
+      <OptimizeReferences>true</OptimizeReferences>
+      <GenerateDebugInformation>true</GenerateDebugInformation>
+      <AdditionalLibraryDirectories>C:\OpenCV\opencv\build\x64\vc15\lib;C:\tbb\lib\intel64\vc14;%(AdditionalLibraryDirectories)</AdditionalLibraryDirectories>
+      <AdditionalDependencies>opencv_world420.lib;tbb.lib;%(AdditionalDependencies)</AdditionalDependencies>
+    </Link>
+  </ItemDefinitionGroup>
+  <ItemGroup>
+    <ClCompile Include="Gauss_tbb.cpp" />
+  </ItemGroup>
+  <Import Project="$(VCTargetsPath)\Microsoft.Cpp.targets" />
+  <ImportGroup Label="ExtensionTargets">
+  </ImportGroup>
+</Project>