// B1. Làm trơn Gaussian
					GaussianBlur(sourceImage, sourceImage, Size(5, 5), 0, 0);

					// B2. Tính Edge Strength và Edge Direction bằng cách lọc Sobel theo x và y
					// Edge Strength = Gx + Gy; Edge Direction = artan(Gy/Gx);
					EdgeDetector ed;
					Mat EdgeStrength;

					float H_Sob[9] = { 1, 0, -1, 2, 0, -2, 1, 0, -1 };
					float V_Sob[9] = { -1, -2, -1, 0, 0, 0, 1, 2, 1 };
					vector<float> H_kernel;
					vector<float> V_kernel;
					float t;

					H_kernel.assign(H_Sob, H_Sob + 9);
					V_kernel.assign(V_Sob, V_Sob + 9);
					t = 1.0f / 4;
					Mat horizon;	//Wx
					Mat vertical;	//Wy

					// tích chập chiều ngang
					_Convo.SetKernel(H_kernel, 3, 3);

					_Convo.DoConvolution(sourceImage, horizon);

					//Tích chập chiều dọc
					_Convo.SetKernel(V_kernel, 3, 3);
					_Convo.DoConvolution(sourceImage, vertical);

					// Edge Strength
					EdgeStrength = abs(horizon) + abs(vertical);


					//Edge Direction
					Mat flagPixelDel = Mat::zeros(EdgeStrength.size(), CV_8U);

					// Đi lần lượt từng pixel của ảnh nguồn, bỏ pixel ở rìa,  và áp dụng kĩ thuật nonmaximum suppression:
					for (int i = 0; i < EdgeStrength.rows; i ++)
					{

						for (int j = 0; j < EdgeStrength.cols; j ++)
						{
							// Nếu pixel ở rìa ảnh thì xét không phải biên cạnh.
							if (i == 0 || j == 0 || i == EdgeStrength.rows - 1 || j == EdgeStrength.cols - 1)
								flagPixelDel.at<uchar>(i, j) = 255;
							else
							{
								float res = (atan2(vertical.at<short>(i, j), horizon.at<short>(i, j)) * 180) / PI;
								// Xét pixel kề theo hướng cạnh. Nếu không phải cực đại cục bộ thì loại, không phải là biên cạnh.
								if (res <= 22.5 || res > 157.5)
								{
									if (EdgeStrength.at<short>(i + 1, j) > EdgeStrength.at<short>(i, j) || EdgeStrength.at<short>(i - 1, j) > EdgeStrength.at<short>(i, j))
										flagPixelDel.at<uchar>(i, j) = 255;

								}
								else
									if (res > 22.5 && res <= 67.5)
									{
										if (EdgeStrength.at<short>(i + 1, j + 1) > EdgeStrength.at<short>(i, j) || EdgeStrength.at<short>(i - 1, j - 1) > EdgeStrength.at<short>(i, j))

											flagPixelDel.at<uchar>(i, j) = 255;
									}
									else
										if (res > 67.5 && res <= 112.5)
										{
											if (EdgeStrength.at<short>(i, j + 1) > EdgeStrength.at<short>(i, j) || EdgeStrength.at<short>(i, j - 1) > EdgeStrength.at<short>(i, j))

												flagPixelDel.at<uchar>(i, j) = 255;

										}
										else
										{
											if (EdgeStrength.at<short>(i + 1, j - 1) > EdgeStrength.at<short>(i, j) || EdgeStrength.at<short>(i - 1, j + 1) > EdgeStrength.at<short>(i, j))
												flagPixelDel.at<uchar>(i, j) = 255;
										}
							}
						}
					}


					horizon.release();
					vertical.release();

					//Áp dụng kĩ thuật hysteresis để loại pixel thừa.

					double min, max;
					minMaxLoc(EdgeStrength, &min, &max);
					// Ngưỡng cho Canny.
					double T0 = min + 80;
					double T1 = min + 100;


					int neighborI[8] = {-1, -1, -1, 0, 1, 1, 1, 0};
					int neighborJ[8] = { -1, 0, 1, 1, 1, 0, -1, -1};
					destinationImage.release();
					destinationImage = Mat::zeros(EdgeStrength.size(), CV_8U) ;
					// Nếu giá trị độ lớn đạo hàm lớn hơn T1 thì kết nạp là hướng cạnh.
					for (int i = 0; i < EdgeStrength.rows; i ++)
						for (int j = 0; j < EdgeStrength.cols; j ++)
						{
							if (flagPixelDel.at<uchar>(i, j) == 0)
								if (EdgeStrength.at<short>(i, j) >= T1)
									destinationImage.at<uchar>(i, j) = 255;
						}

						cout << "jehe" << endl;

						Mat flagTemp = Mat::zeros(EdgeStrength.size(), CV_8U);
						// Nếu giá trị độ lớn đạo hàm những pixel thoả mãn điều kiện cực đại cục bộ, nhưng không lớn hơn T1, mà lớn hơn T2 và pixel kề là biên cạnh thì kết nạp.
						for (int i = 0; i < EdgeStrength.rows; i ++)
							for (int j = 0; j < EdgeStrength.cols; j ++)
							{
								if (flagPixelDel.at<uchar>(i, j) == 0)
									if (destinationImage.at<uchar>(i, j) == 0)
										if (EdgeStrength.at<short>(i, j) >= T0)
										{
											// Nếu pixel ở rìa ảnh thì không xét.
											if (i == 0 || j == 0 || i == EdgeStrength.rows - 1 || j == EdgeStrength.cols - 1)
												continue;
											else
											{
												for (int ij = 0; ij < 8; ij ++)
													if (destinationImage.at<uchar>(i + neighborI[ij], j + neighborJ[ij]) == 255)
													{
														flagTemp.at<uchar>(i, j) = 255;
														break;
													}

											}
										}
							}
							destinationImage = destinationImage + flagTemp;

							flagPixelDel.release();
							flagTemp.release();
				}