Aaru/DiscImageChef.Core/Entropy.cs

// /***************************************************************************
// The Disc Image Chef
// ----------------------------------------------------------------------------
//
// Filename       : Entropy.cs
// Author(s)      : Natalia Portillo <claunia@claunia.com>
//
// Component      : Core algorithms.
//
// --[ Description ] ----------------------------------------------------------
//
//     Calculates the entropy of an image
//
// --[ License ] --------------------------------------------------------------
//
//     This program is free software: you can redistribute it and/or modify
//     it under the terms of the GNU General Public License as
//     published by the Free Software Foundation, either version 3 of the
//     License, or (at your option) any later version.
//
//     This program is distributed in the hope that it will be useful,
//     but WITHOUT ANY WARRANTY; without even the implied warranty of
//     MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the
//     GNU General Public License for more details.
//
//     You should have received a copy of the GNU General Public License
//     along with this program.  If not, see <http://www.gnu.org/licenses/>.
//
// ----------------------------------------------------------------------------
// Copyright © 2011-2019 Natalia Portillo
// ****************************************************************************/

using System;
using System.Collections.Generic;
using System.Linq;
using DiscImageChef.Checksums;
using DiscImageChef.CommonTypes.Interfaces;
using DiscImageChef.CommonTypes.Structs;
using DiscImageChef.Console;

namespace DiscImageChef.Core
{
    public class Entropy
    {
        bool        debug;
        IMediaImage inputFormat;
        bool        verbose;

        public Entropy(bool debug, bool verbose, IMediaImage inputFormat)
        {
            this.debug       = debug;
            this.verbose     = verbose;
            this.inputFormat = inputFormat;
        }

        public event InitProgressHandler   InitProgressEvent;
        public event UpdateProgressHandler UpdateProgressEvent;
        public event EndProgressHandler    EndProgressEvent;
        public event InitProgressHandler   InitProgress2Event;
        public event UpdateProgressHandler UpdateProgress2Event;
        public event EndProgressHandler    EndProgress2Event;

        public EntropyResults[] CalculateTracksEntropy(bool duplicatedSectors)
        {
            List<EntropyResults> entropyResultses = new List<EntropyResults>();

            try
            {
                List<Track> inputTracks = inputFormat.Tracks;

                InitProgressEvent?.Invoke();

                foreach(Track currentTrack in inputTracks)
                {
                    EntropyResults trackEntropy = new EntropyResults {Track = currentTrack.TrackSequence, Entropy = 0};
                    UpdateProgressEvent
                      ?.Invoke($"Entropying track {currentTrack.TrackSequence} of {inputTracks.Max(t => t.TrackSequence)}",
                               currentTrack.TrackSequence, inputTracks.Max(t => t.TrackSequence));

                    ulong[]      entTable              = new ulong[256];
                    ulong        trackSize             = 0;
                    List<string> uniqueSectorsPerTrack = new List<string>();

                    trackEntropy.Sectors = currentTrack.TrackEndSector - currentTrack.TrackStartSector + 1;
                    DicConsole.VerboseWriteLine("Track {0} has {1} sectors", currentTrack.TrackSequence,
                                                trackEntropy.Sectors);

                    InitProgress2Event?.Invoke();

                    for(ulong i = currentTrack.TrackStartSector; i <= currentTrack.TrackEndSector; i++)
                    {
                        UpdateProgress2Event
                          ?.Invoke($"Entropying sector {i             + 1} of track {currentTrack.TrackSequence}",
                                   (long)(currentTrack.TrackEndSector - (i + 1)),
                                   (long)trackEntropy.Sectors);
                        byte[] sector = inputFormat.ReadSector(i, currentTrack.TrackSequence);

                        if(duplicatedSectors)
                        {
                            string sectorHash = Sha1Context.Data(sector, out _);
                            if(!uniqueSectorsPerTrack.Contains(sectorHash)) uniqueSectorsPerTrack.Add(sectorHash);
                        }

                        foreach(byte b in sector) entTable[b]++;

                        trackSize += (ulong)sector.LongLength;
                    }

                    EndProgress2Event?.Invoke();

                    trackEntropy.Entropy += entTable.Select(l => (double)l / (double)trackSize)
                                                    .Select(frequency => -(frequency * Math.Log(frequency, 2))).Sum();

                    if(duplicatedSectors) trackEntropy.UniqueSectors = uniqueSectorsPerTrack.Count;

                    entropyResultses.Add(trackEntropy);
                }

                EndProgressEvent?.Invoke();
            }
            catch(Exception ex)
            {
                if(debug) DicConsole.DebugWriteLine("Could not get tracks because {0}", ex.Message);
                else DicConsole.ErrorWriteLine("Unable to get separate tracks, not calculating their entropy");
            }

            return entropyResultses.ToArray();
        }

        public EntropyResults CalculateMediaEntropy(bool duplicatedSectors)
        {
            EntropyResults entropy       = new EntropyResults {Entropy = 0};
            ulong[]        entTable      = new ulong[256];
            ulong          diskSize      = 0;
            List<string>   uniqueSectors = new List<string>();

            entropy.Sectors = inputFormat.Info.Sectors;
            DicConsole.WriteLine("Sectors {0}", entropy.Sectors);
            InitProgressEvent?.Invoke();
            for(ulong i = 0; i < entropy.Sectors; i++)
            {
                UpdateProgressEvent?.Invoke($"Entropying sector {i + 1}", (long)(i + 1), (long)entropy.Sectors);
                byte[] sector = inputFormat.ReadSector(i);

                if(duplicatedSectors)
                {
                    string sectorHash = Sha1Context.Data(sector, out _);
                    if(!uniqueSectors.Contains(sectorHash)) uniqueSectors.Add(sectorHash);
                }

                foreach(byte b in sector) entTable[b]++;

                diskSize += (ulong)sector.LongLength;
            }

            EndProgressEvent?.Invoke();

            entropy.Entropy += entTable.Select(l => (double)l / (double)diskSize)
                                       .Select(frequency => -(frequency * Math.Log(frequency, 2))).Sum();

            if(duplicatedSectors) entropy.UniqueSectors = uniqueSectors.Count;

            return entropy;
        }
    }

    public struct EntropyResults
    {
        public uint   Track;
        public double Entropy;
        public int?   UniqueSectors;
        public ulong  Sectors;
    }
}