public IEnumerable <XYPolygon> DivideIntoGrid(int Factor)
{
ASCIIGrid ag = new ASCIIGrid();
ag.Data = new MathNet.Numerics.LinearAlgebra.Double.DenseMatrix(Factor, Factor);
ag.XOrigin = this.Points.Min(p => p.X);
ag.YOrigin = this.Points.Min(p => p.Y);
ag.GridSize = (Points.Max(p => p.X) - Points.Min(p => p.X)) / Factor;
var area = XYPolygon.GetPolygons(ag);
for (int i = 0; i < area.GetLength(0); i++)
{
for (int j = 0; j < area.GetLength(1); j++)
{
yield return(area[i, j]);
}
}
}
[STAThread]//Due to OpenFileDialog
static void Main(string[] args)
{
//Creates an open FileDialog
OpenFileDialog ofd = new OpenFileDialog();
ofd.Filter = "Known file types (*.asc)|*.asc"; //Only open .asc-files
ofd.Multiselect = false;
//Now show the dialog and continue if the user presses ok
if (ofd.ShowDialog() == DialogResult.OK)
{
//Prepare a dictiondary to hold the DFS2-Files
Dictionary<int, DFS2> _files = new Dictionary<int, DFS2>();
//Get the directory of the chosen file
string dir = Path.GetDirectoryName(ofd.FileName);
//Get all the file names in the directory sorted alphabetically
var AllAscFiles = Directory.GetFiles(dir,"*"+ Path.GetExtension(ofd.FileName), SearchOption.TopDirectoryOnly).OrderBy(var => var.ToString());
//Loop all the files
foreach (var file in AllAscFiles)
{
//Create an asciireader
ASCIIGrid asc = new ASCIIGrid();
asc.Load(file);
DFS2 dfs;
//Get the filename
string fileName =Path.GetFileNameWithoutExtension(file);
//The filekey is the last three digits in the filename
int FileKey = int.Parse(fileName.Substring(fileName.Length-3,3));
//Find the file in the dictionnary of files
if (!_files.TryGetValue(FileKey, out dfs))
{
//The file was not there
//Create a new DFS2
dfs = new DFS2( Path.Combine(Path.GetDirectoryName(file), fileName.Substring(11, fileName.Length - 11) + ".dfs2"), 1, asc);
//Set grid and geo info
dfs.NumberOfColumns = asc.NumberOfColumns;
dfs.NumberOfRows = asc.NumberOfRows;
dfs.XOrigin = asc.XOrigin;
dfs.YOrigin = asc.YOrigin;
dfs.GridSize = asc.GridSize;
//set time info
dfs.TimeOfFirstTimestep = DateTime.Parse(fileName.Substring(0, 10));
dfs.TimeStep = TimeSpan.FromDays(1);
//Set item info
dfs.FirstItem.EumItem = DHI.Generic.MikeZero.eumItem.eumIPrecipitationRate;
dfs.FirstItem.Name = "Radar precipitation";
dfs.FirstItem.EumUnit = DHI.Generic.MikeZero.eumUnit.eumUmillimeterPerDay;
//Add to dictionary
_files.Add(FileKey, dfs);
}
//Set the data of the next timestep
dfs.SetData(dfs.NumberOfTimeSteps, 1,asc.Data);
}
//All the files have been read
//Dispose the dfs-files
foreach (var dfs in _files.Values)
dfs.Dispose();
}
}
public IEnumerable<XYPolygon> DivideIntoGrid(int Factor)
{
ASCIIGrid ag = new ASCIIGrid();
ag.Data = new MathNet.Numerics.LinearAlgebra.Double.DenseMatrix(Factor, Factor);
ag.XOrigin = this.Points.Min(p => p.X);
ag.YOrigin = this.Points.Min(p => p.Y);
ag.GridSize = (Points.Max(p => p.X) -Points.Min(p => p.X)) / Factor;
var area = XYPolygon.GetPolygons(ag);
for (int i = 0; i < area.GetLength(0); i++)
for (int j = 0; j < area.GetLength(1); j++)
yield return area[i, j];
}
