Test files preload and som stress function changes

MES_Wind/bin/Debug/MES_test/u_test.prj

+PROJCS["unnamed",GEOGCS["unnamed ellipse",DATUM["unknown",SPHEROID["unnamed",6378137,0],EXTENSION["PROJ4_GRIDS","@null"]],PRIMEM["Greenwich",0],UNIT["degree",0.0174532925199433]],PROJECTION["Mercator_2SP"],PARAMETER["standard_parallel_1",0],PARAMETER["central_meridian",0],PARAMETER["false_easting",0],PARAMETER["false_northing",0],UNIT["Meter",1],EXTENSION["PROJ4","+proj=merc +a=6378137 +b=6378137 +lat_ts=0.0 +lon_0=0.0 +x_0=0.0 +y_0=0 +k=1.0 +units=m +nadgrids=@null +wktext  +no_defs"]]
\ No newline at end of file

MES_Wind/bin/Debug/MES_test/v_test.asc

+ncols	20
+nrows	20
+xllcorner	3749431.702
+yllcorner	5557489.554
+cellsize	526.316
+NODATA_value	-9999.000
+0.000 2.101 4.179 6.211 8.174 10.047 11.808 13.439 14.921 16.238 17.375 18.319 19.061 19.592 19.906 19.999 19.872 19.524 18.960 18.186
+0.000 1.991 3.960 5.884 7.744 9.518 11.187 12.732 14.136 15.383 16.460 17.355 18.058 18.561 18.858 18.947 18.826 18.496 17.962 17.229
+0.000 1.880 3.740 5.558 7.314 8.989 10.565 12.024 13.350 14.528 15.546 16.391 17.055 17.530 17.810 17.894 17.780 17.469 16.964 16.272
+0.000 1.770 3.520 5.231 6.884 8.461 9.944 11.317 12.565 13.674 14.631 15.427 16.051 16.498 16.763 16.842 16.734 16.441 15.966 15.314
+0.000 1.659 3.300 4.904 6.453 7.932 9.322 10.610 11.780 12.819 13.717 14.463 15.048 15.467 15.715 15.789 15.688 15.413 14.968 14.357
+0.000 1.548 3.080 4.577 6.023 7.403 8.701 9.902 10.994 11.965 12.802 13.498 14.045 14.436 14.667 14.736 14.642 14.386 13.970 13.400
+0.000 1.438 2.860 4.250 5.593 6.874 8.079 9.195 10.209 11.110 11.888 12.534 13.042 13.405 13.620 13.684 13.596 13.358 12.972 12.443
+0.000 1.327 2.640 3.923 5.163 6.345 7.458 8.488 9.424 10.255 10.973 11.570 12.039 12.374 12.572 12.631 12.550 12.331 11.975 11.486
+0.000 1.217 2.420 3.596 4.733 5.817 6.836 7.781 8.638 9.401 10.059 10.606 11.035 11.343 11.524 11.579 11.505 11.303 10.977 10.529
+0.000 1.106 2.200 3.269 4.302 5.288 6.215 7.073 7.853 8.546 9.145 9.642 10.032 10.311 10.477 10.526 10.459 10.276 9.979 9.572
+0.000 0.995 1.980 2.942 3.872 4.759 5.593 6.366 7.068 7.692 8.230 8.678 9.029 9.280 9.429 9.473 9.413 9.248 8.981 8.614
+0.000 0.885 1.760 2.615 3.442 4.230 4.972 5.659 6.282 6.837 7.316 7.713 8.026 8.249 8.381 8.421 8.367 8.221 7.983 7.657
+0.000 0.774 1.540 2.288 3.012 3.702 4.350 4.951 5.497 5.982 6.401 6.749 7.022 7.218 7.334 7.368 7.321 7.193 6.985 6.700
+0.000 0.664 1.320 1.961 2.581 3.173 3.729 4.244 4.712 5.128 5.487 5.785 6.019 6.187 6.286 6.316 6.275 6.165 5.987 5.743
+0.000 0.553 1.100 1.635 2.151 2.644 3.107 3.537 3.927 4.273 4.572 4.821 5.016 5.156 5.238 5.263 5.229 5.138 4.989 4.786
+0.000 0.442 0.880 1.308 1.721 2.115 2.486 2.829 3.141 3.418 3.658 3.857 4.013 4.125 4.191 4.210 4.183 4.110 3.992 3.829
+0.000 0.332 0.660 0.981 1.291 1.586 1.864 2.122 2.356 2.564 2.743 2.893 3.010 3.093 3.143 3.158 3.138 3.083 2.994 2.871
+0.000 0.221 0.440 0.654 0.860 1.058 1.243 1.415 1.571 1.709 1.829 1.928 2.006 2.062 2.095 2.105 2.092 2.055 1.996 1.914
+0.000 0.111 0.220 0.327 0.430 0.529 0.621 0.707 0.785 0.855 0.914 0.964 1.003 1.031 1.048 1.053 1.046 1.028 0.998 0.957
+0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000

MES_Wind/bin/Debug/MES_test/v_test.asc.aux.xml

+<PAMDataset>
+  <PAMRasterBand band="1">
+    <Metadata>
+      <MDI key="STATISTICS_MAXIMUM">19.999000549316</MDI>
+      <MDI key="STATISTICS_MEAN">6.9478149894997</MDI>
+      <MDI key="STATISTICS_MINIMUM">0</MDI>
+      <MDI key="STATISTICS_STDDEV">5,63996171951294</MDI>
+    </Metadata>
+  </PAMRasterBand>
+</PAMDataset>

MES_Wind/bin/Debug/MES_test/v_test.prj

+PROJCS["unnamed",GEOGCS["unnamed ellipse",DATUM["unknown",SPHEROID["unnamed",6378137,0],EXTENSION["PROJ4_GRIDS","@null"]],PRIMEM["Greenwich",0],UNIT["degree",0.0174532925199433]],PROJECTION["Mercator_2SP"],PARAMETER["standard_parallel_1",0],PARAMETER["central_meridian",0],PARAMETER["false_easting",0],PARAMETER["false_northing",0],UNIT["Meter",1],EXTENSION["PROJ4","+proj=merc +a=6378137 +b=6378137 +lat_ts=0.0 +lon_0=0.0 +x_0=0.0 +y_0=0 +k=1.0 +units=m +nadgrids=@null +wktext  +no_defs"]]
\ No newline at end of file

MES_Wind/frmMain.Designer.cs

@@ -214,7 +214,7 @@
             this.btnLoadWindY.Name = "btnLoadWindY";
             this.btnLoadWindY.Size = new System.Drawing.Size(96, 46);
             this.btnLoadWindY.TabIndex = 3;
-            this.btnLoadWindY.Text = "Load Ywind";
+            this.btnLoadWindY.Text = "Load test files";
             this.btnLoadWindY.UseVisualStyleBackColor = true;
             this.btnLoadWindY.Click += new System.EventHandler(this.btnLoadWindY_Click);
             // 

MES_Wind/frmMain.cs

@@ -56,6 +56,8 @@ namespace MES_Wind
             double vwind = 0;
             double umod = 0;
             double sinwind = 0;
+            double anglewind = 0;
+            double angleline = Math.Atan2((endY - startY),(endX - startX));
             double climwind = 0;
             double distance = Math.Sqrt((endX - startX) * (endX - startX) + (endY - startY) * (endY - startY));
             double distpropD = distance / dThreshold;
@@ -80,14 +82,12 @@ namespace MES_Wind
                         curY = curY + constYdiff;
                     }
                     Coordinate coords = new Coordinate(curX,curY);
-                    RcIndex rowColumnU = Uwind_raster.DataSet.Bounds.ProjToCell(coords);
-                    uwind = Uwind_raster.DataSet.Value[rowColumnU.Row, rowColumnU.Column];
-                    RcIndex rowColumnV = Vwind_raster.DataSet.Bounds.ProjToCell(coords);
-                    vwind = Vwind_raster.DataSet.Value[rowColumnV.Row, rowColumnV.Column];
-                    RcIndex rowColumnC = Vwind_raster.DataSet.Bounds.ProjToCell(coords);
-                    climwind = Vwind_raster.DataSet.Value[rowColumnC.Row, rowColumnC.Column];
+                    uwind = interpol(coords, Uwind_raster);
+                    vwind = interpol(coords, Vwind_raster);
+                    climwind = interpol(coords, clim_layer);
                     umod = Math.Sqrt(uwind*uwind + vwind*vwind);
-                    sinwind = Math.Sin(Math.Atan2(uwind,vwind));
+                    anglewind = Math.Atan2(vwind,uwind) - angleline;
+                    sinwind = Math.Sin(anglewind);
                     double C_height = 1.0;
                     if (umod < 20)
                     { //wind is too low 
@@ -117,10 +117,26 @@ namespace MES_Wind
                    
                     lineCheckPoint.Add(chkpnt);
                 }
+            }
+            else
+            {
+                 
             }
             return lineCheckPoint;
         }
-
+        public IFeatureSet brokenpoints(List<CheckPoint> chklist)
+        {
+            IFeatureSet points = new FeatureSet(FeatureType.Point);
+            Coordinate chkcords = new Coordinate();
+            return points;
+        }
+        public double interpol(Coordinate coords, IMapRasterLayer raster)
+        {
+         double rval = 0;
+            RcIndex rc= raster.DataSet.Bounds.ProjToCell(coords);
+            rval = raster.DataSet.Value[rc.Row, rc.Column];
+         return  rval;
+        }
 
         private void bntLoadWindX_Click(object sender, EventArgs e)
         {
@@ -131,7 +147,15 @@ namespace MES_Wind
 
         private void btnLoadWindY_Click(object sender, EventArgs e)
         {
-            map1.AddLayer();
+            string curDir = Environment.CurrentDirectory;
+            string path_to_tests = "\\MES_test\\";
+            string file_path = curDir + path_to_tests;
+
+            map1.AddLayer(file_path+"u_test.asc");
+            map1.AddLayer(file_path + "v_test.asc");
+            map1.AddLayer(file_path + "clim5_test.asc");
+            map1.AddLayer(file_path + "clim10_test.asc");
+            map1.AddLayer(file_path + "powerlines.shp");
             map1.ZoomToMaxExtent();
         }
 
@@ -200,8 +224,8 @@ namespace MES_Wind
                 //copy line layer FeatureSet
                 IFeatureSet pwlineSet = pwlLayer.DataSet;
                 // new FeatureSet for resulting broken powerlines
-                IFeatureSet brklineSet = new FeatureSet(FeatureType.Line);
-                DataTable dt = pwlineSet.DataTable;
+                //IFeatureSet brklineSet = new FeatureSet(FeatureType.Line);
+                //DataTable dt = pwlineSet.DataTable;
                 List<CheckPoint> fullCheckList = new List<CheckPoint>();
                 foreach (IFeature feature in pwlineSet.Features)
                 {
@@ -210,14 +234,14 @@ namespace MES_Wind
                     DataRow featureData = feature.DataRow;
                     int id = int.Parse(featureData["PW_ID"].ToString());
                     int year = int.Parse(featureData["Year"].ToString());
-                    double height = double.Parse(featureData["height"].ToString());
+                    double height = double.Parse(featureData["height_m"].ToString());
                     int power = int.Parse(featureData["Power"].ToString());
                     LineString linestr = feature.BasicGeometry as LineString;
                     if (linestr != null)
                     { // case if powerline consists of one line
                         // get coordinates list
                         IList<Coordinate> points = linestr.Coordinates;
-                        IFeature brklineFeature = brklineSet.AddFeature(linestr);
+                        
 
                         // cycle throw all points in line
                         for (int i=1; i< points.Count; i++)
@@ -239,17 +263,38 @@ namespace MES_Wind
                         }
 
                         fullCheckList.AddRange(lineCheckList);
-                        brklineFeature.CopyAttributes(feature);
+                        
                     }
                     else
                     {//case if powerline is multiline
                         MultiLineString multiline = feature.BasicGeometry as MultiLineString;
-                        if ( multiline != null){ 
+                        if ( multiline != null){
+                            
                             foreach (IGeometry line in multiline.Geometries)
                             {
                                 IList<Coordinate> points = line.Coordinates;
+                                for (int i = 1; i < points.Count; i++)
+                                {
+                                    List<CheckPoint> segmentCheckList = new List<CheckPoint>();
+                                    double x1 = points[i - 1].X;
+                                    double y1 = points[i - 1].Y;
+                                    double x2 = points[i].X;
+                                    double y2 = points[i].Y;
+                                    if (power > 5 && power < 330)
+                                    {
+                                        segmentCheckList = CalcBrkPoint(x1, y1, x2, y2, distThreshold, u_rasterLayer, v_rasterLayer, clim10_rasterLayer, height);
+                                    }
+                                    else
+                                    {
+                                        segmentCheckList = CalcBrkPoint(x1, y1, x2, y2, distThreshold, u_rasterLayer, v_rasterLayer, clim5_rasterLayer, height);
+                                    }
+                                    lineCheckList.AddRange(segmentCheckList);
+                                }
+
+                                fullCheckList.AddRange(lineCheckList);
+
                             }
-                            IFeature brklineFeature = brklineSet.AddFeature(multiline);
+                            
                             MessageBox.Show("Works");
                         }
 

MES_Wind/obj/Debug/MES_Wind.csproj.FileListAbsolute.txt

@@ -188,3 +188,89 @@ C:\Users\Geophyslab-laptop\Documents\Mes_wind\MES_Wind\obj\Debug\MES_Wind.frmGra
 C:\Users\Geophyslab-laptop\Documents\Mes_wind\MES_Wind\obj\Debug\MES_Wind.frmMain.resources
 C:\Users\Geophyslab-laptop\Documents\Mes_wind\MES_Wind\obj\Debug\MES_Wind.Properties.Resources.resources
 C:\Users\Geophyslab-laptop\Documents\Mes_wind\MES_Wind\obj\Debug\MES_Wind.csproj.GenerateResource.Cache
+C:\MES_Wind\MES_Wind\obj\Debug\MES_Wind.exe
+C:\MES_Wind\MES_Wind\obj\Debug\MES_Wind.pdb
+C:\MES_Wind\MES_Wind\bin\Debug\MES_Wind.exe
+C:\MES_Wind\MES_Wind\bin\Debug\MES_Wind.pdb
+C:\MES_Wind\MES_Wind\bin\Debug\DotSpatial.Analysis.pdb
+C:\MES_Wind\MES_Wind\bin\Debug\DotSpatial.Analysis.xml
+C:\MES_Wind\MES_Wind\bin\Debug\DotSpatial.Compatibility.pdb
+C:\MES_Wind\MES_Wind\bin\Debug\DotSpatial.Compatibility.xml
+C:\MES_Wind\MES_Wind\bin\Debug\DotSpatial.Controls.pdb
+C:\MES_Wind\MES_Wind\bin\Debug\DotSpatial.Controls.xml
+C:\MES_Wind\MES_Wind\bin\Debug\DotSpatial.Data.pdb
+C:\MES_Wind\MES_Wind\bin\Debug\DotSpatial.Data.xml
+C:\MES_Wind\MES_Wind\bin\Debug\DotSpatial.Data.Forms.pdb
+C:\MES_Wind\MES_Wind\bin\Debug\DotSpatial.Data.Forms.xml
+C:\MES_Wind\MES_Wind\bin\Debug\DotSpatial.Extensions.pdb
+C:\MES_Wind\MES_Wind\bin\Debug\DotSpatial.Extensions.xml
+C:\MES_Wind\MES_Wind\bin\Debug\DotSpatial.Modeling.Forms.pdb
+C:\MES_Wind\MES_Wind\bin\Debug\DotSpatial.Modeling.Forms.xml
+C:\MES_Wind\MES_Wind\bin\Debug\DotSpatial.Mono.pdb
+C:\MES_Wind\MES_Wind\bin\Debug\DotSpatial.Mono.xml
+C:\MES_Wind\MES_Wind\bin\Debug\DotSpatial.Positioning.pdb
+C:\MES_Wind\MES_Wind\bin\Debug\DotSpatial.Positioning.xml
+C:\MES_Wind\MES_Wind\bin\Debug\DotSpatial.Positioning.Design.pdb
+C:\MES_Wind\MES_Wind\bin\Debug\DotSpatial.Positioning.Design.xml
+C:\MES_Wind\MES_Wind\bin\Debug\DotSpatial.Positioning.Forms.pdb
+C:\MES_Wind\MES_Wind\bin\Debug\DotSpatial.Positioning.Forms.xml
+C:\MES_Wind\MES_Wind\bin\Debug\DotSpatial.Projections.pdb
+C:\MES_Wind\MES_Wind\bin\Debug\DotSpatial.Projections.xml
+C:\MES_Wind\MES_Wind\bin\Debug\DotSpatial.Projections.Forms.pdb
+C:\MES_Wind\MES_Wind\bin\Debug\DotSpatial.Projections.Forms.xml
+C:\MES_Wind\MES_Wind\bin\Debug\DotSpatial.Serialization.pdb
+C:\MES_Wind\MES_Wind\bin\Debug\DotSpatial.Serialization.xml
+C:\MES_Wind\MES_Wind\bin\Debug\DotSpatial.Symbology.pdb
+C:\MES_Wind\MES_Wind\bin\Debug\DotSpatial.Symbology.xml
+C:\MES_Wind\MES_Wind\bin\Debug\DotSpatial.Symbology.Forms.pdb
+C:\MES_Wind\MES_Wind\bin\Debug\DotSpatial.Symbology.Forms.xml
+C:\MES_Wind\MES_Wind\bin\Debug\DotSpatial.Topology.pdb
+C:\MES_Wind\MES_Wind\bin\Debug\DotSpatial.Topology.xml
+C:\MES_Wind\MES_Wind\bin\Debug\ZedGraph.xml
+C:\MES_Wind\MES_Wind\bin\Debug\de-DE\DotSpatial.Controls.resources.dll
+C:\MES_Wind\MES_Wind\bin\Debug\it-IT\DotSpatial.Controls.resources.dll
+C:\MES_Wind\MES_Wind\bin\Debug\pt-BR\DotSpatial.Controls.resources.dll
+C:\MES_Wind\MES_Wind\bin\Debug\zh-CN\DotSpatial.Controls.resources.dll
+C:\MES_Wind\MES_Wind\bin\Debug\zh-CN\DotSpatial.Data.Forms.resources.dll
+C:\MES_Wind\MES_Wind\bin\Debug\zh-CN\DotSpatial.Modeling.Forms.resources.dll
+C:\MES_Wind\MES_Wind\bin\Debug\zh-MO\DotSpatial.Modeling.Forms.resources.dll
+C:\MES_Wind\MES_Wind\bin\Debug\cs\DotSpatial.Projections.Forms.resources.dll
+C:\MES_Wind\MES_Wind\bin\Debug\de\DotSpatial.Projections.Forms.resources.dll
+C:\MES_Wind\MES_Wind\bin\Debug\de-DE\DotSpatial.Projections.Forms.resources.dll
+C:\MES_Wind\MES_Wind\bin\Debug\el\DotSpatial.Projections.Forms.resources.dll
+C:\MES_Wind\MES_Wind\bin\Debug\it-IT\DotSpatial.Projections.Forms.resources.dll
+C:\MES_Wind\MES_Wind\bin\Debug\ja-JP\DotSpatial.Projections.Forms.resources.dll
+C:\MES_Wind\MES_Wind\bin\Debug\lt\DotSpatial.Projections.Forms.resources.dll
+C:\MES_Wind\MES_Wind\bin\Debug\nl\DotSpatial.Projections.Forms.resources.dll
+C:\MES_Wind\MES_Wind\bin\Debug\pt\DotSpatial.Projections.Forms.resources.dll
+C:\MES_Wind\MES_Wind\bin\Debug\pt-BR\DotSpatial.Projections.Forms.resources.dll
+C:\MES_Wind\MES_Wind\bin\Debug\ro\DotSpatial.Projections.Forms.resources.dll
+C:\MES_Wind\MES_Wind\bin\Debug\th\DotSpatial.Projections.Forms.resources.dll
+C:\MES_Wind\MES_Wind\bin\Debug\vi-VN\DotSpatial.Projections.Forms.resources.dll
+C:\MES_Wind\MES_Wind\bin\Debug\zh-CN\DotSpatial.Projections.Forms.resources.dll
+C:\MES_Wind\MES_Wind\bin\Debug\de-DE\DotSpatial.Symbology.resources.dll
+C:\MES_Wind\MES_Wind\bin\Debug\cs\DotSpatial.Symbology.Forms.resources.dll
+C:\MES_Wind\MES_Wind\bin\Debug\de-DE\DotSpatial.Symbology.Forms.resources.dll
+C:\MES_Wind\MES_Wind\bin\Debug\en-US\DotSpatial.Symbology.Forms.resources.dll
+C:\MES_Wind\MES_Wind\bin\Debug\it\DotSpatial.Symbology.Forms.resources.dll
+C:\MES_Wind\MES_Wind\bin\Debug\it-IT\DotSpatial.Symbology.Forms.resources.dll
+C:\MES_Wind\MES_Wind\bin\Debug\pt-BR\DotSpatial.Symbology.Forms.resources.dll
+C:\MES_Wind\MES_Wind\bin\Debug\zh-CN\DotSpatial.Symbology.Forms.resources.dll
+C:\MES_Wind\MES_Wind\bin\Debug\de\ZedGraph.resources.dll
+C:\MES_Wind\MES_Wind\bin\Debug\es\ZedGraph.resources.dll
+C:\MES_Wind\MES_Wind\bin\Debug\fr\ZedGraph.resources.dll
+C:\MES_Wind\MES_Wind\bin\Debug\hu\ZedGraph.resources.dll
+C:\MES_Wind\MES_Wind\bin\Debug\it\ZedGraph.resources.dll
+C:\MES_Wind\MES_Wind\bin\Debug\ja\ZedGraph.resources.dll
+C:\MES_Wind\MES_Wind\bin\Debug\pt\ZedGraph.resources.dll
+C:\MES_Wind\MES_Wind\bin\Debug\ru\ZedGraph.resources.dll
+C:\MES_Wind\MES_Wind\bin\Debug\sk\ZedGraph.resources.dll
+C:\MES_Wind\MES_Wind\bin\Debug\sv\ZedGraph.resources.dll
+C:\MES_Wind\MES_Wind\bin\Debug\tr\ZedGraph.resources.dll
+C:\MES_Wind\MES_Wind\bin\Debug\zh-CN\ZedGraph.resources.dll
+C:\MES_Wind\MES_Wind\bin\Debug\zh-tw\ZedGraph.resources.dll
+C:\MES_Wind\MES_Wind\obj\Debug\MES_Wind.csprojResolveAssemblyReference.cache
+C:\MES_Wind\MES_Wind\obj\Debug\MES_Wind.frmGraph.resources
+C:\MES_Wind\MES_Wind\obj\Debug\MES_Wind.frmMain.resources
+C:\MES_Wind\MES_Wind\obj\Debug\MES_Wind.Properties.Resources.resources
+C:\MES_Wind\MES_Wind\obj\Debug\MES_Wind.csproj.GenerateResource.Cache