public void GetPolygonsTest()
{
var m = new HydroNumerics.MikeSheTools.Core.Model(@"E:\dhi\data\dkm\dk2\result\DK2_v3_gvf_PT_100p_24hr.she");
var precip = new HydroNumerics.MikeSheTools.DFS.DFS2(m.Input.MIKESHE_FLOWMODEL.Climate.PrecipitationRate.FULLY_DISTRIBUTED.DFS_2D_DATA_FILE.FILE_NAME);
var actual = XYPolygon.GetPolygons(precip);
System.Data.DataTable dt = new System.Data.DataTable();
dt.Columns.Add("ID", typeof(int));
int k = 0;
using (HydroNumerics.Geometry.Shapes.ShapeWriter sw = new Geometry.Shapes.ShapeWriter(@"d:\temp\out.shp"))
{
for (int i = 0; i < precip.NumberOfColumns; i++)
{
for (int j = 0; j < precip.NumberOfRows; j++)
{
var drow = dt.NewRow();
drow[0] = k;
sw.Write(new GeoRefData()
{
Geometry = actual[i, j], Data = drow
});
k++;
}
}
}
}
