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++; } } } }