public Histogram CreateDipAzimuthHistogram(int numberOfClasses)
{
Histogram histogram = new Histogram(numberOfClasses, -Math.PI, +Math.PI);
IRegularGrid2D dZdX = CreatePartialDerivativeIGrid();
IRegularGrid2D dZdY = CreatePartialDerivativeJGrid();
for (int i = 0; i < SizeI; ++i)
{
for (int j = 0; j < SizeJ; ++j)
{
double?x = dZdX[i, j];
double?y = dZdY[i, j];
if (x.HasValue && y.HasValue)
{
if (x != 0.0 && y != 0.0)
{
double dipAzimuth = Math.Atan2(-y.Value, -x.Value); // Negate x and y to give down-dip rather than up-dip
histogram.Accumulate(dipAzimuth);
}
}
}
}
return(histogram);
}
/// <summary>
/// Generate histograms for the image. One histogram is generated for each color channel.
/// You will need to call the Refresh function to do the painting afterward.
/// </summary>
/// <param name="image">The image to retrieve histogram from</param>
/// <param name="numberOfBins">The number of bins for each histogram</param>
public void GenerateHistograms(IImage image, int numberOfBins)
{
IImage[] channels = image.Split();
Type imageType = Toolbox.GetBaseType(image.GetType(), "Image`2");
IColor typeOfColor = Activator.CreateInstance(imageType.GetGenericArguments()[0]) as IColor;
String[] channelNames = Reflection.ReflectColorType.GetNamesOfChannels(typeOfColor);
System.Drawing.Color[] colors = Reflection.ReflectColorType.GetDisplayColorOfChannels(typeOfColor);
float minVal, maxVal;
#region Get the maximum and minimum color intensity values
System.Type typeOfDepth = imageType.GetGenericArguments()[1];
if (typeOfDepth == typeof(Byte))
{
minVal = 0.0f;
maxVal = 256.0f;
}
else
{
#region obtain the maximum and minimum color value
double[] minValues, maxValues;
System.Drawing.Point[] minLocations, maxLocations;
image.MinMax(out minValues, out maxValues, out minLocations, out maxLocations);
double min = minValues[0], max = maxValues[0];
for (int i = 1; i < minValues.Length; i++)
{
if (minValues[i] < min) min = minValues[i];
if (maxValues[i] > max) max = maxValues[i];
}
#endregion
minVal = (float)min;
maxVal = (float)max;
}
#endregion
for (int i = 0; i < channels.Length; i++)
using (Histogram hist = new Histogram(numberOfBins, new RangeF(minVal, maxVal)))
{
hist.Accumulate(new IImage[1] { channels[i] });
AddHistogram(channelNames[i], colors[i], hist);
}
}
