public override Tuple <IImageRaster <IRaster2DInteger, int>, IImageRaster <IRaster2DInteger, float>, IList <int[]>, IList <float[]> > Initialize(IRaster2DInteger raster) { IImageRaster2D <int> region_image = new ImageRaster2D <int>(raster); IImageRaster2D <float> distance_image = new ImageRaster2D <float>(raster, Single.MaxValue); Tuple <IList <int[]>, IList <float[]> > Centroids = InitializeCentroids(raster); IList <int[]> cluster_spatial_centroids = Centroids.Item1; IList <float[]> cluster_feature_centroids = Centroids.Item2; int[] neigbourhood_dimensions = d_cluster_dimensions_double; int neigbourhood_element_count = ((neigbourhood_dimensions[0] * 2) + 1) * ((neigbourhood_dimensions[1] * 2) + 1); int[] element_index_indexes = new int[neigbourhood_element_count]; for (int centroid_index = 0; centroid_index < cluster_spatial_centroids.Count; centroid_index++) { raster.GetNeigbourhoodElementIndexesRBA(cluster_spatial_centroids[centroid_index], neigbourhood_dimensions, element_index_indexes); for (int element_index_index = 0; element_index_index < neigbourhood_element_count; element_index_index++) { int element_index = element_index_indexes[element_index_index]; if (element_index != -1) { int[] element_coordinates = raster.GetElementCoordinates(element_index); float distance = FunctionDistanceEuclidean.ComputeStatic(element_coordinates, cluster_spatial_centroids[centroid_index]); if (distance < distance_image.GetElementValue(element_index)) { distance_image.SetElementValue(element_index, distance); region_image.SetElementValue(element_index, centroid_index); } } } } return(new Tuple <IImageRaster <IRaster2DInteger, int>, IImageRaster <IRaster2DInteger, float>, IList <int[]>, IList <float[]> >(region_image, distance_image, cluster_spatial_centroids, cluster_feature_centroids)); }
private void Start() { Complex[] block; while (is_running) { if (blocker.OutputQueue.TryDequeue(out block)) { // do fft MathNet.Numerics.IntegralTransforms.Fourier.Forward(block); //TODO Shift image? float sum = 0; for (int y_index = 0; y_index < image.Raster.Size1; y_index++) { //image.SetElementValue(x_index, y_index, (float)random.NextDouble()); float value = 0; for (int f_index = 0; f_index < f_rate; f_index++) { if ((y_index * f_rate) + f_index < (block.Length)) { value += (float)(block[(y_index * f_rate) + f_index].Magnitude * ((float)((y_index * f_rate) + f_index))); } } value /= f_rate; sum += value; image.SetElementValue(x_index, y_index, value); } sum /= image.Raster.Size1 * 0.5f; for (int y_index = 0; y_index < image.Raster.Size1; y_index++) { image.SetElementValue(x_index, y_index, image.GetElementValue(x_index, y_index) / sum); } x_index = (x_index + 1) % image.Raster.Size0; } //Render image Render(); } }
private ImageRaster2D <Color> IntializeOutputImage(IImageRaster2D <Color> input_image, int output_width, int output_height) { ImageRaster2D <Color> output_image = new ImageRaster2D <Color>(output_width, output_height); Random random = new Random(); /* randomize the last few columns from the input image */ for (int element_index = 0; element_index < output_image.Raster.ElementCount; element_index++) { int random_x = random.Next((input_image.Raster.Size0 - 1) - (scale * 2)) + scale; int random_y = random.Next((input_image.Raster.Size1 - 1) - (scale * 2)) + scale; output_image.SetElementValue(element_index, input_image.GetElementValue(random_x, random_y)); } return(output_image); }
public IImageRaster <IRaster2DInteger, bool> Render(Tuple <List <int>, IRaster3DInteger> render_source) { List <int> elements_true = render_source.Item1; IRaster3DInteger source_raster = render_source.Item2; float[,] coordinates_image = new float[elements_true.Count, 3]; float[,] elements_projection = new float[elements_true.Count, 3]; int [] coordinates = new int [3]; for (int element_index_index = 0; element_index_index < elements_true.Count; element_index_index++) { source_raster.GetElementCoordinatesRBA(elements_true[element_index_index], coordinates); coordinates_image[element_index_index, 0] = coordinates[0]; coordinates_image[element_index_index, 1] = coordinates[1]; coordinates_image[element_index_index, 2] = coordinates[2]; } for (int element_index = 0; element_index < elements_true.Count; element_index++) { elements_projection[element_index, 0] = ((coordinates_image[element_index, 0] - image_focus[0]) * projection_vector_x[0]) + ((coordinates_image[element_index, 1] - image_focus[1]) * projection_vector_x[1]) + ((coordinates_image[element_index, 2] - image_focus[2]) * projection_vector_x[2]); elements_projection[element_index, 1] = ((coordinates_image[element_index, 0] - image_focus[0]) * projection_vector_y[0]) + ((coordinates_image[element_index, 1] - image_focus[1]) * projection_vector_y[1]) + ((coordinates_image[element_index, 2] - image_focus[2]) * projection_vector_y[2]); elements_projection[element_index, 2] = ((coordinates_image[element_index, 0] - image_focus[0]) * projection_vector_z[0]) + ((coordinates_image[element_index, 1] - image_focus[1]) * projection_vector_z[1]) + ((coordinates_image[element_index, 2] - image_focus[2]) * projection_vector_z[2]); } Raster2DInteger bitmap_raster = new Raster2DInteger(bitmap_size_x, bitmap_size_y); ImageRaster2D <bool> destination_image = new ImageRaster2D <bool>(bitmap_size_x, bitmap_size_y); for (int element_index = 0; element_index < elements_true.Count; element_index++) { int x_target = (int)elements_projection[element_index, 0] + (bitmap_raster.Size0 / 2); int y_target = (int)elements_projection[element_index, 1] + (bitmap_raster.Size1 / 2); if (bitmap_raster.ContainsCoordinates(x_target, y_target)) { destination_image.SetElementValue(x_target, y_target, true); } } return(destination_image); }
public IImageRaster2D <float> Render(IImageSpace <float, float> source_image) { IImageRaster2D <float> destination_image = new ImageRaster2D <float>(); float[] coordinates = new float [render_origen.Length]; for (int index_y = 0; index_y < this.resolution_y; index_y++) { for (int index_x = 0; index_x < this.resolution_x; index_x++) { ToolsCollection.CopyRBA(this.render_origen, coordinates); for (int index_dimension = 0; index_dimension < this.dimension_count; index_dimension++) { coordinates[index_dimension] += (render_stride_x[index_dimension] * index_x) + (render_stride_y[index_dimension] * index_y); } destination_image.SetElementValue(index_x, index_y, source_image.GetLocationValue(coordinates)); } } return(destination_image); }
public static void DistanceTransform2DSlowRBA(ImageRaster2D <bool> mask_image, float[] voxel_size, ImageRaster2D <float> distance_image) { IRaster2DInteger raster = mask_image.Raster; int size_0 = mask_image.Raster.Size0; int size_1 = mask_image.Raster.Size1; // Apply the transform in the x-direction //for (int plane_index = 0; plane_index < size_z * size_y; plane_index++) //{ Parallel.For(0, size_0, index_0 => { for (int index_1 = 0; index_1 < size_1; index_1++) { float best_distance = float.MaxValue; if (mask_image.GetElementValue(index_0, index_1)) { best_distance = 0; } else { for (int index_0_1 = 0; index_0_1 < size_0; index_0_1++) { for (int index_1_1 = 0; index_1_1 < size_1; index_1_1++) { if (mask_image.GetElementValue(index_0_1, index_1_1)) { float distance = ToolsMath.Sqrt(ToolsMath.Sqr((index_0_1 - index_0) * voxel_size[0]) + ToolsMath.Sqr((index_1_1 - index_1) * voxel_size[0])); if (distance < best_distance) { best_distance = distance; } } } } } distance_image.SetElementValue(index_0, index_1, best_distance); } }); }