コード例 #1
0
ファイル: PointLayer.cs プロジェクト: hanchao/DotSpatial
 private void Configure(IFeatureSet inFeatureSet)
 {
     if (inFeatureSet == null) throw new ArgumentNullException("inFeatureSet");
     if (inFeatureSet.FeatureType != FeatureType.Point &&
         inFeatureSet.FeatureType != FeatureType.MultiPoint) throw new PointFeatureTypeException();
     
     Symbology = new PointScheme();
     Symbology.SetParentItem(this);
 }
コード例 #2
0
        private PointScheme ColorBox(IEnumerable <Color> colors)
        {
            PointScheme ps = new PointScheme();

            ps.Categories.Clear();
            foreach (Color color in colors)
            {
                IColorable c = _template as IColorable;
                if (c != null)
                {
                    c.Color = color;
                }
                PointCategory pc = new PointCategory(_template);
                ps.Categories.Add(pc);
            }
            ps.Categories[0].FilterExpression = "[" + _classificationField + "] < ";
            return(ps);
        }
コード例 #3
0
ファイル: PointLayer.cs プロジェクト: yangkf1985/DotSpatial
        private void Configure(IFeatureSet inFeatureSet)
        {
            FeatureType ft = inFeatureSet.FeatureType;

            if (ft != FeatureType.Point && ft != FeatureType.MultiPoint && ft != FeatureType.Unspecified)
            {
                throw new PointFeatureTypeException();
            }
            if (inFeatureSet.NumRows() == 0)
            {
                MyExtent = new Extent(-180, -90, 180, 90);
            }
            if (inFeatureSet.NumRows() == 1)
            {
                MyExtent = inFeatureSet.Extent.Copy();
                MyExtent.ExpandBy(10, 10);
            }
            Symbology = new PointScheme();
        }
コード例 #4
0
 private PointScheme ColorBox(IEnumerable<Color> colors)
 {
     PointScheme ps = new PointScheme();
     ps.Categories.Clear();
     foreach (Color color in colors)
     {
         IColorable c = _template as IColorable;
         if (c != null)
         {
             c.Color = color;
         }
         PointCategory pc = new PointCategory(_template);
         ps.Categories.Add(pc);
     }
     ps.Categories[0].FilterExpression = "[" + _classificationField + "] < ";
     return ps;
 }
コード例 #5
0
ファイル: PointLayer.cs プロジェクト: DIVEROVIEDO/DotSpatial
 private void Configure(IFeatureSet inFeatureSet)
 {
     FeatureType ft = inFeatureSet.FeatureType;
     if (ft != FeatureType.Point && ft != FeatureType.MultiPoint && ft != FeatureType.Unspecified)
     {
         throw new PointFeatureTypeException();
     }
     if (inFeatureSet.NumRows() == 0)
     {
         MyExtent = new Extent(-180, -90, 180, 90);
     }
     if (inFeatureSet.NumRows() == 1)
     {
         MyExtent = inFeatureSet.Extent.Copy();
         MyExtent.ExpandBy(10, 10);
     }
     Symbology = new PointScheme();
 }
コード例 #6
0
        /// <summary>
        /// Creates a new raster with the specified cell size.  If the cell size
        /// is zero, this will default to the shorter of the width or height
        /// divided by 256.  If the cell size produces a raster that is greater
        /// than 8, 000 pixels in either dimension, it will be re-sized to
        /// create an 8, 000 length or width raster.
        /// </summary>
        /// <param name="fs">The featureset to convert to a raster.</param>
        /// <param name="extent">Force the raster to this specified extent.</param>
        /// <param name="cellSize">The double extent of the cell.</param>
        /// <param name="fieldName">The integer field index of the file.</param>
        /// <param name="outputFileName">The fileName of the raster to create.</param>
        /// <param name="driverCode">The optional GDAL driver code to use if using GDAL
        /// for a format that is not discernable from the file extension.  An empty string
        ///  is usually perfectly acceptable here.</param>
        /// <param name="options">For GDAL rasters, they can be created with optional parameters
        ///  passed in as a string array.  In most cases an empty string is perfectly acceptable.</param>
        /// <param name="progressHandler">An interface for handling the progress messages.</param>
        /// <returns>Generates a raster from the vectors.</returns>
        public static IRaster ToRaster(IFeatureSet fs, Extent extent, double cellSize, string fieldName,
                                       string outputFileName,
                                       string driverCode, string[] options, IProgressHandler progressHandler)
        {
            Extent env = extent;
            if (cellSize == 0)
            {
                if (env.Width < env.Height)
                {
                    cellSize = env.Width / 256;
                }
                else
                {
                    cellSize = env.Height / 256;
                }
            }
            int w = (int)Math.Ceiling(env.Width / cellSize);
            if (w > 8000)
            {
                w = 8000;
                cellSize = env.Width / 8000;
            }
            int h = (int)Math.Ceiling(env.Height / cellSize);
            if (h > 8000)
            {
                h = 8000;
            }
            Bitmap bmp = new Bitmap(w, h);
            Graphics g = Graphics.FromImage(bmp);
            g.Clear(Color.Transparent);
            g.SmoothingMode = SmoothingMode.None;
            g.TextRenderingHint = TextRenderingHint.SingleBitPerPixel;
            g.InterpolationMode = InterpolationMode.NearestNeighbor;
            Hashtable colorTable;
            MapArgs args = new MapArgs(new Rectangle(0, 0, w, h), env, g);

            switch (fs.FeatureType)
            {
                case FeatureType.Polygon:
                    {
                        MapPolygonLayer mpl = new MapPolygonLayer(fs);
                        PolygonScheme ps = new PolygonScheme();
                        colorTable = ps.GenerateUniqueColors(fs, fieldName);
                        mpl.Symbology = ps;
                        mpl.DrawRegions(args, new List<Extent> { env });
                    }
                    break;
                case FeatureType.Line:
                    {
                        MapLineLayer mpl = new MapLineLayer(fs);
                        LineScheme ps = new LineScheme();
                        colorTable = ps.GenerateUniqueColors(fs, fieldName);
                        mpl.Symbology = ps;
                        mpl.DrawRegions(args, new List<Extent> { env });
                    }
                    break;
                default:
                    {
                        MapPointLayer mpl = new MapPointLayer(fs);
                        PointScheme ps = new PointScheme();
                        colorTable = ps.GenerateUniqueColors(fs, fieldName);
                        mpl.Symbology = ps;
                        mpl.DrawRegions(args, new List<Extent> { env });
                    }
                    break;
            }
            Type tp = fieldName == "FID" ? typeof(int) : fs.DataTable.Columns[fieldName].DataType;
            // We will try to convert to double if it is a string
            if (tp == typeof(string))
            {
                tp = typeof(double);
            }
            InRamImageData image = new InRamImageData(bmp, env);
            ProgressMeter pm = new ProgressMeter(progressHandler, "Converting To Raster Cells", h);

            IRaster output;
            output = Raster.Create(outputFileName, driverCode, w, h, 1, tp, options);
            output.Bounds = new RasterBounds(h, w, env);

            double noDataValue = output.NoDataValue;

            if (fieldName != "FID")
            {
                // We can't use this method to calculate Max on a non-existent FID field.
                double dtMax = Convert.ToDouble(fs.DataTable.Compute("Max(" + fieldName + ")", ""));
                double dtMin = Convert.ToDouble(fs.DataTable.Compute("Min(" + fieldName + ")", ""));

                if (dtMin <= noDataValue && dtMax >= noDataValue)
                {
                    if (dtMax != GetFieldValue(tp, "MaxValue"))
                    {
                        output.NoDataValue = noDataValue;
                    }
                    else if (dtMin != GetFieldValue(tp, "MinValue"))
                    {
                        output.NoDataValue = noDataValue;
                    }
                }
            }

            List<RcIndex> locations = new List<RcIndex>();
            List<string> failureList = new List<string>();
            for (int row = 0; row < h; row++)
            {
                for (int col = 0; col < w; col++)
                {
                    Color c = image.GetColor(row, col);
                    if (c.A == 0)
                    {
                        output.Value[row, col] = output.NoDataValue;
                    }
                    else
                    {
                        if (colorTable.ContainsKey(c) == false)
                        {
                            if (c.A < 125)
                            {
                                output.Value[row, col] = output.NoDataValue;
                                continue;
                            }

                            // Use a color matching distance to pick the closest member
                            object val = GetCellValue(w, h, row, col, image, c, colorTable, locations);

                            output.Value[row, col] = GetDouble(val, failureList);
                        }
                        else
                        {
                            output.Value[row, col] = GetDouble(colorTable[c], failureList);
                        }
                    }
                }
                pm.CurrentValue = row;
            }
            const int maxIterations = 5;
            int iteration = 0;
            while (locations.Count > 0)
            {
                List<RcIndex> newLocations = new List<RcIndex>();
                foreach (RcIndex location in locations)
                {
                    object val = GetCellValue(w, h, location.Row, location.Column, image,
                                              image.GetColor(location.Row, location.Column), colorTable, newLocations);
                    output.Value[location.Row, location.Column] = GetDouble(val, failureList);
                }
                locations = newLocations;
                iteration++;
                if (iteration > maxIterations)
                {
                    break;
                }
            }

            pm.Reset();
            return output;
        }