public void LoadDBFFile(string FileName, string IDColumn, string NameColumn, string LengthColumn, string LRColumn)
{
NewMessage("Reading: " + FileName);
using (DBFReader dbr = new DBFReader(FileName))
{
while (!dbr.EndOfData)
{
var data = dbr.ReadNext();
int catchid = int.Parse(data[IDColumn].ToString());
string type = data[NameColumn].ToString();
double length = double.Parse(data[LengthColumn].ToString());
double lengthReductionFactor = double.Parse(data[LRColumn].ToString());
StreamClassification streamclass;
if (StreamClasses.TryGetValue(type, out streamclass))
{
Tuple<double, double> previous;
if (!ReductionFactors.TryGetValue(catchid, out previous))
{
previous = new Tuple<double, double>(0, 0);
ReductionFactors.Add(catchid, previous);
}
DataRow datarow = Data.Rows.Find(new object[] { catchid });
if (datarow == null)
{
datarow = Data.NewRow();
datarow[0] = catchid;
Data.Rows.Add(datarow);
}
double summer = 0;
double winter = 0;
if (length > 0)
{
summer = 0.01 * MultiplicationFactor * Math.Pow(streamclass.StreamDepthSummer / (length * lengthReductionFactor / (streamclass.StreamVelocitySummer * 365.0 * 86400.0)), Exponent);
winter = 0.01 * MultiplicationFactor * Math.Pow(streamclass.StreamDepthWinter / (length * lengthReductionFactor / (streamclass.StreamVelocityWinter * 365.0 * 86400.0)), Exponent);
}
int i = StreamClasses.IndexOfKey(streamclass.StreamType);
datarow[i + 1] = length;
datarow[i + 4] = winter;
datarow[i + 8] = summer;
var newfact = new Tuple<double, double>(1.0 - (1.0 - previous.Item1) * (1 - summer), 1.0 - (1.0 - previous.Item2) * (1 - winter));
datarow[7] = newfact.Item2;
datarow[11] = newfact.Item1;
ReductionFactors[catchid] = newfact;
}
}
}
}
public ShapeReader(string FileName)
: base()
{
_fileName = FileName;
if (!File.Exists(FileName))
throw new FileNotFoundException("Could not find the file: " + Path.GetFullPath(FileName));
_data = new DBFReader(FileName);
// Open shapefile
_shapePointer = ShapeLib.SHPOpen(FileName, "rb");
int nbentities = 0;
double[] arr1 = new double[4];
double[] arr2 = new double[4];
ShapeLib.SHPGetInfo(_shapePointer, ref nbentities, ref type, arr1, arr2);
NoOfEntries = nbentities;
}
public ShapeReader(string FileName)
: base()
{
_fileName = FileName;
if (!File.Exists(FileName))
{
throw new FileNotFoundException("Could not find the file: " + Path.GetFullPath(FileName));
}
_data = new DBFReader(FileName);
// Open shapefile
_shapePointer = ShapeLib.SHPOpen(FileName, "rb");
int nbentities = 0;
double[] arr1 = new double[4];
double[] arr2 = new double[4];
ShapeLib.SHPGetInfo(_shapePointer, ref nbentities, ref type, arr1, arr2);
NoOfEntries = nbentities;
}
public override void Initialize(DateTime Start, DateTime End, IEnumerable<Catchment> Catchments)
{
base.Initialize(Start, End, Catchments);
//Source, Catchment
Dictionary<string, int> Sources = new Dictionary<string,int>();
Dictionary<string, XYPoint> PointSources = new Dictionary<string,XYPoint>();
List<XYPoint> points = new List<XYPoint>();
//Read in the points
using (ShapeReader sr = new ShapeReader(ShapeFile.FileName))
{
foreach (var gd in sr.GeoData)
{
XYPoint xp = gd.Geometry as XYPoint;
if (!PointSources.ContainsKey(gd.Data[ShapeFile.ColumnNames[0]].ToString())) //Some sources are listed multiple times
if(gd.Data[ShapeFile.ColumnNames[1]].ToString()=="land")
PointSources.Add(gd.Data[ShapeFile.ColumnNames[0]].ToString(), gd.Geometry as XYPoint);
}
}
NewMessage("Distributing sources in catchments");
//Distribute points on catchments
Parallel.ForEach(PointSources, (p) =>
{
foreach (var c in Catchments)
{
if (c.Geometry.Contains(p.Value.X, p.Value.Y))
{
//if (c.CoastalZones.Any(co=>co.Contains(p.Value.X, p.Value.Y)))
lock (Lock)
Sources.Add(p.Key, c.ID);
break;
}
}
});
NewMessage(PointSources.Count +" point sources distributed on " + Sources.Values.Distinct().Count().ToString() + " catchments");
Dictionary<string, int> entries = new Dictionary<string, int>();
foreach (var source in Sources.Keys)
entries.Add(source, 0);
NewMessage("Reading outlet data");
//Read source data and distrubute on catchments
using (DBFReader dbf = new DBFReader(DBFFile.FileName))
{
int k = 0;
for (int i = 0; i < dbf.NoOfEntries; i++)
{
string id = dbf.ReadString(i, DBFFile.ColumnNames[0]).Trim();
int year = dbf.ReadInt(i, DBFFile.ColumnNames[1]);
double value = dbf.ReadDouble(i, DBFFile.ColumnNames[2]);
int catchid;
if (Sources.TryGetValue(id, out catchid))
{
entries[id] += 1;
Dictionary<int, double> timevalues;
if (!YearlyData.TryGetValue(catchid, out timevalues))
{
timevalues = new Dictionary<int, double>();
YearlyData.Add(catchid, timevalues);
}
if (timevalues.ContainsKey(year))
timevalues[year] += value;
else
{
if (year >= Start.Year) //This also removes some -99 years
timevalues.Add(year, value);
}
}
else
{
//if we have outlet data but no source placed
k++;
}
}
NewMessage(k + " time series entries have no points");
}
foreach (var c in YearlyData.Values)
foreach (var val in c.ToList())
c[val.Key] = val.Value/ ((DateTime.DaysInMonth(val.Key,2) + 337.0) * 86400.0);
NewMessage(entries.Values.Count(v => v == 0) + " points have no time series data");
NewMessage("Initialized");
}
/// <summary>
/// Reads the particles and distributes them on the dictionaries
/// </summary>
/// <param name="shapefilename"></param>
public IEnumerable<Particle> ReadParticleFile(string shapefilename, Func<DBFReader, int, SinkType,bool> ExcludeThis, bool DrainToBoundaryOption)
{
List<int> RedoxedParticles = new List<int>();
Dictionary<int, Particle> NonRedoxedParticles = new Dictionary<int, Particle>();
NewMessage("Reading particles from: " + shapefilename);
int k = 0;
using (DBFReader sr = new DBFReader(shapefilename))
{
try
{
for (int i = 0; i < sr.NoOfEntries; i++)
{
int id = sr.ReadInt(i, "ID");
SinkType sinktype = (SinkType)Enum.Parse(typeof(SinkType), sr.ReadString(i, "SinkType"));
if (sinktype == SinkType.Active_cell) //Do not create a new particle
RedoxedParticles.Add(id);
else if (ExcludeThis != null && ExcludeThis(sr, i, sinktype))
{
k++;
}
else
{
Particle p = new Particle();
p.Registration = sr.ReadInt(i, "Registrati");
p.XStart = sr.ReadDouble(i, "X-Birth");
p.YStart = sr.ReadDouble(i, "Y-Birth");
p.ZStart = sr.ReadDouble(i, "Z-Birth");
p.X = sr.ReadDouble(i, "X-Reg");
p.Y = sr.ReadDouble(i, "Y-Reg");
p.Z = sr.ReadDouble(i, "Z-Reg");
p.TravelTime = sr.ReadDouble(i, "TravelTime");
p.SinkType = sinktype;
NonRedoxedParticles.Add(id, p);
if(DrainToBoundaryOption && p.SinkType == SinkType.Drain_to_Boundary)
{
p.X = p.XStart;
p.Y = p.YStart;
p.Z = p.ZStart;
p.Registration = 0;
}
}
}
}
catch (Exception e)
{
NewMessage("Error reading data");
}
}
int nosink = 0;
//Set the registration on all particles that have be redoxed
foreach (var pid in RedoxedParticles)
if (NonRedoxedParticles.ContainsKey(pid))
NonRedoxedParticles[pid].Registration = 1;
else
nosink++;
NewMessage("Excluded: " + k + " particles.");
return NonRedoxedParticles.Values;
}
