for (int i = 0; i < src->rows; i++) {
			for (int j = 0; j < src->cols; j++) {
				if (src->at<Vec3b>(i, j) == src->at<Vec3b>(y, x)) {
					nom += (i - rb)*(j - cb);
					den += (j - cb)*(j - cb) - (i - rb)*(i - rb);
					if (src->at<Vec3b>(i, j) != src->at<Vec3b>(i - 1, j - 1))	{ perimeter++; continue; }
					if (src->at<Vec3b>(i, j) != src->at<Vec3b>(i - 1, j))		{ perimeter++; continue; }
					if (src->at<Vec3b>(i, j) != src->at<Vec3b>(i - 1, j + 1))	{ perimeter++; continue; }
					if (src->at<Vec3b>(i, j) != src->at<Vec3b>(i, j - 1))		{ perimeter++; continue; }
					if (src->at<Vec3b>(i, j) != src->at<Vec3b>(i, j + 1))		{ perimeter++; continue; }
					if (src->at<Vec3b>(i, j) != src->at<Vec3b>(i + 1, j - 1))	{ perimeter++; continue; }
					if (src->at<Vec3b>(i, j) != src->at<Vec3b>(i + 1, j))		{ perimeter++; continue; }
					if (src->at<Vec3b>(i, j) != src->at<Vec3b>(i + 1, j + 1))	{ perimeter++; continue; }

				}
			}
		}