Ejemplo n.º 1
0
    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;
          }
        }
      }
    }
Ejemplo n.º 2
0
    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;
     
    }
Ejemplo n.º 3
0
        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;
        }
Ejemplo n.º 4
0
    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");
    }
Ejemplo n.º 5
0
    /// <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;

    }