for (int i = wlo; i < whi; i++){
				for (int j = hlo; j < hhi; j++){
					Gaussian g = new Gaussian();
					double[][] filter = g.gaussianFilter(radius, sigma);

					// iterate over the blur matrix
					for (int k = 0; k < filter.length; k++){
						for (int l = 0; l < filter.length; l++){

							int x = i + k;
							int y = j + l;

							// check to see if we go out of range of the image
							if (x >= width){
								x = width - 1;
							}
							if (y >= height){
								y = height - 1;
							}
							// sum up all the contributions
							pixels[width*i + j].R += filter[k][l] * original[width*x + y].R;
							pixels[width*i + j].G += filter[k][l] * original[width*x + y].G;
							pixels[width*i + j].B += filter[k][l] * original[width*x + y].B;
						}
					}

				}
			}