public static List <Minutia> ExtractMinutiae(double[,] imgBytes)
{
var lsEnhanced = SymmetryHelper.EstimateLS(imgBytes, sigma1, sigma2);
var psEnhanced = SymmetryHelper.EstimatePS(imgBytes, sigma1, sigma2);
//ImageHelper.SaveComplexArrayAsHSV(lsEnhanced,"C:\\temp\\lsenh.png");
//ImageHelper.SaveArray(NormalizeArray(psEnhanced.Select2D(x=>x.Magnitude)), "C:\\temp\\psenh.png");
var psi = KernelHelper.Zip2D(psEnhanced,
lsEnhanced.Select2D(x => x.Magnitude), (x, y) =>
x * (1.0d - y));
return(SearchMinutiae(psi, lsEnhanced, psEnhanced));
}
/// <summary>
/// Loads a tileset.
/// </summary>
/// <param name="sender"></param>
/// <param name="e"></param>
private void buttonLoadTileset_Click(object sender, EventArgs e)
{
if (openFileDialogTileset.ShowDialog() != DialogResult.OK)
{
return;
}
try
{
var document = new XmlDocument();
document.Load(openFileDialogTileset.FileName);
var root = document.DocumentElement;
_Tiles.Clear();
_Neighbors.Clear();
var tilemapDirectory = new FileInfo(openFileDialogTileset.FileName).DirectoryName;
// Load tiles
uint id = 0;
var filenameToId = new Dictionary <string, uint>();
foreach (XmlNode tile in root.SelectNodes("tiles/tile"))
{
var filename = tile.Attributes["name"].Value;
var symmetry = SymmetryHelper.Parse(tile.Attributes["symmetry"].Value);
double weight = 1.0;
if (tile.Attributes["weight"] != null)
{
weight = double.Parse(tile.Attributes["weight"].Value);
}
var nbOfSyms = SymmetryHelper.NumberOfPossibleOrientations(symmetry);
if (File.Exists(Path.Combine(tilemapDirectory, filename + " 0.png")))
{
var data = new Array2D <int> [nbOfSyms];
for (uint i = 0; i < data.Length; i++)
{
data[i] = FetchImageData(Path.Combine(tilemapDirectory, filename + " " + i + ".png"));
}
_Tiles.Add(new Tile <int>(data, symmetry, weight));
}
else
{
var imageData = FetchImageData(Path.Combine(tilemapDirectory, filename + ".png"));
_Tiles.Add(new Tile <int>(imageData, symmetry, weight));
}
filenameToId.Add(filename, id++);
}
// Load neighbors
foreach (XmlNode neighbor in root.SelectNodes("neighbors/neighbor"))
{
var tile1 = neighbor.Attributes["left"].Value;
var tile2 = neighbor.Attributes["right"].Value;
var sep1 = tile1.LastIndexOf(' ');
uint orien1 = 0;
if (sep1 > 0 && uint.TryParse(tile1.Substring(sep1 + 1), out orien1))
{
tile1 = tile1.Substring(0, sep1);
}
var sep2 = tile2.LastIndexOf(' ');
uint orien2 = 0;
if (sep2 > 0 && uint.TryParse(tile2.Substring(sep2 + 1), out orien2))
{
tile2 = tile2.Substring(0, sep2);
}
_Neighbors.Add(new TilingNeighbor(
new TileWithOrientation(filenameToId[tile1], orien1),
new TileWithOrientation(filenameToId[tile2], orien2)));
}
Logger.LogSuccess($"Loaded {_Tiles.Count} tiles and {_Neighbors.Count} neighbors!");
}
catch (Exception ex)
{
Logger.LogFailure(ex.Message);
MessageBox.Show(ex.ToString(), "Unexpected error", MessageBoxButtons.OK, MessageBoxIcon.Error);
}
}
public static double[,] EnhanceImage(double[,] imgBytes)
{
Stopwatch sw = new Stopwatch();
sw.Start();
var g1 = ImageSizeHelper.Reduce2(imgBytes, 1.7d);
var g2 = ImageSizeHelper.Reduce2(g1, 1.21d);
var g3 = ImageSizeHelper.Reduce2(g2, K);
var g4 = ImageSizeHelper.Reduce2(g3, K);
var p3 = ImageSizeHelper.Expand2(g4, K, new Size(g3.GetLength(0), g3.GetLength(1)));
var p2 = ImageSizeHelper.Expand2(g3, K, new Size(g2.GetLength(0), g2.GetLength(1)));
var p1 = ImageSizeHelper.Expand2(g2, 1.21d, new Size(g1.GetLength(0), g1.GetLength(1)));
var l3 = ContrastEnhancement(KernelHelper.Subtract(g3, p3));
var l2 = ContrastEnhancement(KernelHelper.Subtract(g2, p2));
var l1 = ContrastEnhancement(KernelHelper.Subtract(g1, p1));
//SaveArray(l3, "C:\\temp\\l3.png");
//SaveArray(l1, "C:\\temp\\l1.png");
//SaveArray(l2, "C:\\temp\\l2.png");
var ls1 = SymmetryHelper.EstimateLS(l1, sigma1, sigma2);
var ls2 = SymmetryHelper.EstimateLS(l2, sigma1, sigma2);
var ls3 = SymmetryHelper.EstimateLS(l3, sigma1, sigma2);
//SaveComplexArrayAsHSV(ls1, "C:\\temp\\ls1.png");
//SaveComplexArrayAsHSV(ls2, "C:\\temp\\ls2.png");
//SaveComplexArrayAsHSV(ls3, "C:\\temp\\ls3.png");
var ls2Scaled =
KernelHelper.MakeComplexFromDouble(
ImageSizeHelper.Expand2(ls2.Select2D(x => x.Real), K, new Size(l1.GetLength(0), l1.GetLength(1))),
ImageSizeHelper.Expand2(ls2.Select2D(x => x.Imaginary), K, new Size(l1.GetLength(0), l1.GetLength(1))));
var multiplier = KernelHelper.Subtract(ls1.Select2D(x => x.Phase), ls2Scaled.Select2D(x => x.Phase));
for (int x = 0; x < ls1.GetLength(0); x++)
{
for (int y = 0; y < ls1.GetLength(1); y++)
{
ls1[x, y] *= Math.Abs(Math.Cos(multiplier[x, y]));
}
}
DirectionFiltering(l1, ls1, tau1, tau2);
DirectionFiltering(l2, ls2, tau1, tau2);
DirectionFiltering(l3, ls3, tau1, tau2);
var ll2 = ImageSizeHelper.Expand2(l3, K, new Size(l2.GetLength(0), l2.GetLength(1)));
l2 = KernelHelper.Add(ll2, l2);
var ll1 = ImageSizeHelper.Expand2(l2, 1.21d, new Size(l1.GetLength(0), l1.GetLength(1)));
l1 = KernelHelper.Add(ll1, l1);
var ll0 = ImageSizeHelper.Expand2(l1, 1.7d, new Size(imgBytes.GetLength(0), imgBytes.GetLength(1)));
ll0 = ContrastEnhancement(ll0);
sw.Stop();
var enhanced = RearrangeArray(ll0, 0, 255);
return(enhanced);
}
