def METHOD_NAME(image: np.ndarray, x: float, y: float, w: float, h: float) -> Stats:
	"""
	Given (x, y, w, h) denoting a resolution-independent region, return the color stats of the region.
	Return the average color of the region and the std dev of the region.
	"""
	start_x, start_y, end_x, end_y = _get_region(image, x, y, w, h)
	crop = image[start_y:end_y, start_x:end_x].astype(float)
	num_pixels = crop.shape[0] * crop.shape[1]
	crop_sum = np.sum(crop, (0, 1))
	crop_avg = crop_sum / num_pixels
	avg_sum = np.sum(crop_avg)
	crop_color_ratio = [1/3., 1/3., 1/3.] if avg_sum <= 1e-6 else crop_avg / avg_sum
	crop_color_ratio = np.array([crop_color_ratio[2], crop_color_ratio[1], crop_color_ratio[0]])
	crop_sqr_sum = np.sum(np.square(crop), (0, 1))
	crop_stddev = np.sqrt((crop_sqr_sum - (np.square(crop_sum) / num_pixels)) / (num_pixels-1))
	crop_stddev_sum = np.sum(crop_stddev)
	return Stats(
		avg = np.flip(crop_avg), # from BGR to RGB
		avg_sum = avg_sum,
		stddev = np.flip(crop_stddev), # from BGR to RGB
		stddev_sum = crop_stddev_sum,
		color_ratio = crop_color_ratio,
		crop_size=(end_x - start_x, end_y - start_y),
	)