C# 中SharpMap的简单使用实例详解

using BruTile.Predefined;using GeoAPI.CoordinateSystems.Transformations;using ProjNet.CoordinateSystems;using ProjNet.CoordinateSystems.Transformations;using SharpMap;using SharpMap.Data.Providers;using SharpMap.Layers;using SharpMap.Rendering;using SharpMap.Rendering.Thematics;using SharpMap.Styles;using System;using System.Collections.Generic;using System.Data;using System.Data.OleDb;using System.Drawing;using System.Drawing.Drawing2D;using System.Drawing.Text;using System.Linq;using System.Text;using Point = GeoAPI.Geometries.Coordinate;namespace DemoSharpMap{ public class SharpMapHelper { private const string XlsConnectionString = "Provider={2};Data Source={0}//{1};Extended Properties=/"Excel 8.0;HDR=Yes;IMEX=1/""; public static Map InitializeMap(MapType tt,float angle) {  Map map = null;  switch (tt)  {  case MapType.RunLine:   map = InitializeMapOsmWithXls(angle);   break;  case MapType.MapInfo:   map = InitializeMapinfo(angle);   break;  case MapType.ShapeFile:   map = InitializeMapOrig(angle);   break;  case MapType.Static:   map = InitializeMapOsmWithXls2(angle);   break;  default:   map = InitializeMapOsmWithXls(angle);   break;  }  return map; } /// <summary> /// MapInfo格式的地图文件 /// </summary> /// <param name="angle"></param> /// <returns></returns> private static Map InitializeMapinfo(float angle) {  //Initialize a new map of size 'imagesize'  Map map = new Map();  //Set up the countries layer  VectorLayer layCountries = new VectorLayer("Countries");  //Set the datasource to a shapefile in the App_data folder  try  {  layCountries.DataSource = new Ogr("GeoData/MapInfo/countriesMapInfo.tab");  }  catch (TypeInitializationException ex)  {  if (ex.Message == "The type initializer for 'OSGeo.OGR.Ogr' threw an exception.")  {   throw new Exception(   String.Format(    "The application threw a PINVOKE exception. You probably need to copy the unmanaged dll's to your bin directory. They are a part of fwtools {0}. You can download it from: http://home.gdal.org/fwtools/",    GdalRasterLayer.FWToolsVersion));  }  throw;  }  //Set fill-style to green  layCountries.Style.Fill = new SolidBrush(Color.Green);  //Set the polygons to have a black outline  layCountries.Style.Outline = Pens.Black;  layCountries.Style.EnableOutline = true;  layCountries.SRID = 4326;  //Set up a river layer  VectorLayer layRivers = new VectorLayer("Rivers");  //Set the datasource to a shapefile in the App_data folder  layRivers.DataSource = new Ogr("GeoData/MapInfo/riversMapInfo.tab");  //Define a blue 1px wide pen  layRivers.Style.Line = new Pen(Color.Blue, 1);  layRivers.SRID = 4326;  //Set up a river layer  VectorLayer layCities = new VectorLayer("Cities");  //Set the datasource to a shapefile in the App_data folder  layCities.DataSource = new Ogr("GeoData/MapInfo/citiesMapInfo.tab");  layCities.Style.SymbolScale = 0.8f;  layCities.MaxVisible = 40;  layCities.SRID = 4326;  //Set up a country label layer  LabelLayer layLabel = new LabelLayer("Country labels");  layLabel.DataSource = layCountries.DataSource;  layLabel.Enabled = true;  layLabel.LabelColumn = "Name";  layLabel.Style = new LabelStyle();  layLabel.Style.ForeColor = Color.White;  layLabel.Style.Font = new Font(FontFamily.GenericSerif, 12);  layLabel.Style.BackColor = new SolidBrush(Color.FromArgb(128, 255, 0, 0));  layLabel.MaxVisible = 90;  layLabel.MinVisible = 30;  layLabel.Style.HorizontalAlignment = LabelStyle.HorizontalAlignmentEnum.Center;  layLabel.SRID = 4326;  layLabel.MultipartGeometryBehaviour = LabelLayer.MultipartGeometryBehaviourEnum.Largest;  //Set up a city label layer  LabelLayer layCityLabel = new LabelLayer("City labels");  layCityLabel.DataSource = layCities.DataSource;  layCityLabel.Enabled = true;  layCityLabel.LabelColumn = "Name";  layCityLabel.Style = new LabelStyle();  layCityLabel.Style.ForeColor = Color.Black;  layCityLabel.Style.Font = new Font(FontFamily.GenericSerif, 11);  layCityLabel.MaxVisible = layLabel.MinVisible;  layCityLabel.Style.HorizontalAlignment = LabelStyle.HorizontalAlignmentEnum.Left;  layCityLabel.Style.VerticalAlignment = LabelStyle.VerticalAlignmentEnum.Bottom;  layCityLabel.Style.Offset = new PointF(3, 3);  layCityLabel.Style.Halo = new Pen(Color.Yellow, 2);  layCityLabel.TextRenderingHint = TextRenderingHint.AntiAlias;  layCityLabel.SmoothingMode = SmoothingMode.AntiAlias;  layCityLabel.SRID = 4326;  layCityLabel.LabelFilter = LabelCollisionDetection.ThoroughCollisionDetection;  layCityLabel.Style.CollisionDetection = true;  //Add the layers to the map object.  //The order we add them in are the order they are drawn, so we add the rivers last to put them on top  map.Layers.Add(layCountries);  map.Layers.Add(layRivers);  map.Layers.Add(layCities);  map.Layers.Add(layLabel);  map.Layers.Add(layCityLabel);  //增加Layers  var xlsPath = string.Format(XlsConnectionString, System.IO.Directory.GetCurrentDirectory(), "GeoData//Henan.xls", Properties.Settings.Default.OleDbProvider);  var ds = GetDataFromExcel(xlsPath, "Cities");  //var ct = GetCoordinateTransformation();  //TransCoordinate(ds, ct);  string columeName = "Rotation";  //Add Rotation Column  AddColumeToDataSet(ds, columeName, -angle);  var xlsLayer = GetLayerFromDataSet2(ds, Color.GreenYellow);//Set up provider  map.Layers.Add(xlsLayer); //Add layer to map  map.Center = xlsLayer.Envelope.Centre;// new Point(0, 0);  map.MapScale = 350;  //map.Center = new Point(0, 0);  //_ogrSampleDataset = "MapInfo";  //Matrix mat = new Matrix();  //mat.RotateAt(angle, map.WorldToImage(map.Center));  //map.MapTransform = mat;  //map.ZoomToBox(xlsLayer.Envelope);  return map; } /// <summary> /// ShapeFile /// </summary> /// <param name="angle"></param> /// <returns></returns> private static Map InitializeMapOrig(float angle) {  //Initialize a new map of size 'imagesize'  Map map = new Map();  //Set up the countries layer  VectorLayer layCountries = new VectorLayer("Countries");  //Set the datasource to a shapefile in the App_data folder  layCountries.DataSource = new ShapeFile("GeoData/World/countries.shp", true);  //Set fill-style to green  layCountries.Style.Fill = new SolidBrush(Color.FromArgb(64, Color.Green));  //Set the polygons to have a black outline  layCountries.Style.Outline = Pens.Black;  layCountries.Style.EnableOutline = true;  layCountries.SRID = 4326;  //Set up a river layer  VectorLayer layRivers = new VectorLayer("Rivers");  //Set the datasource to a shapefile in the App_data folder  layRivers.DataSource = new ShapeFile("GeoData/World/rivers.shp", true);  //Define a blue 1px wide pen  layRivers.Style.Line = new Pen(Color.Blue, 1);  layRivers.SRID = 4326;  //Set up a cities layer  VectorLayer layCities = new VectorLayer("Cities");  //Set the datasource to a shapefile in the App_data folder  layCities.DataSource = new ShapeFile("GeoData/World/cities.shp", true);  layCities.Style.SymbolScale = 0.8f;  layCities.MaxVisible = 40;  layCities.SRID = 4326;  //Set up a country label layer  LabelLayer layLabel = new LabelLayer("Country labels");  layLabel.DataSource = layCountries.DataSource;  layLabel.Enabled = true;  layLabel.LabelColumn = "Name";  layLabel.Style = new LabelStyle();  layLabel.Style.ForeColor = Color.White;  layLabel.Style.Font = new Font(FontFamily.GenericSerif, 12);  layLabel.Style.BackColor = new SolidBrush(Color.FromArgb(128, 255, 0, 0));  layLabel.MaxVisible = 90;  layLabel.MinVisible = 30;  layLabel.Style.HorizontalAlignment = LabelStyle.HorizontalAlignmentEnum.Center;  layLabel.SRID = 4326;  layLabel.MultipartGeometryBehaviour = LabelLayer.MultipartGeometryBehaviourEnum.Largest;  layLabel.LabelFilter = LabelCollisionDetection.ThoroughCollisionDetection;  layLabel.Style.CollisionDetection = true;  layLabel.LabelPositionDelegate = fdr => fdr.Geometry.InteriorPoint.Coordinate;  layLabel.PriorityColumn = "POPDENS";  //Set up a city label layer  LabelLayer layCityLabel = new LabelLayer("City labels");  layCityLabel.DataSource = layCities.DataSource;  layCityLabel.Enabled = true;  layCityLabel.LabelColumn = "Name";  layCityLabel.Style = new LabelStyle();  layCityLabel.Style.ForeColor = Color.Black;  layCityLabel.Style.Font = new Font(FontFamily.GenericSerif, 11);  layCityLabel.MaxVisible = layLabel.MinVisible;  layCityLabel.Style.HorizontalAlignment = LabelStyle.HorizontalAlignmentEnum.Left;  layCityLabel.Style.VerticalAlignment = LabelStyle.VerticalAlignmentEnum.Bottom;  layCityLabel.Style.Offset = new PointF(3, 3);  layCityLabel.Style.Halo = new Pen(Color.Yellow, 2);  layCityLabel.TextRenderingHint = TextRenderingHint.AntiAlias;  layCityLabel.SmoothingMode = SmoothingMode.AntiAlias;  layCityLabel.SRID = 4326;  layCityLabel.LabelFilter = LabelCollisionDetection.ThoroughCollisionDetection;  layCityLabel.Style.CollisionDetection = true;  layCityLabel.PriorityColumn = "POPULATION";  layCityLabel.Theme = new GradientTheme(layCityLabel.PriorityColumn, 250000, 5000000,  new LabelStyle  {   MaxVisible = 10,   CollisionBuffer = new Size(0, 0),   CollisionDetection = true,   Enabled = true,   ForeColor = Color.LightSlateGray,   Halo = new Pen(Color.Silver, 1),   HorizontalAlignment = LabelStyle.HorizontalAlignmentEnum.Center,   VerticalAlignment = LabelStyle.VerticalAlignmentEnum.Middle,   Font = new Font(GenericFontFamilies.SansSerif.ToString(), 8f, FontStyle.Regular)  },  new LabelStyle  {   MaxVisible = layLabel.MinVisible,   CollisionBuffer = new Size(3, 3),   CollisionDetection = true,   Enabled = true,   ForeColor = Color.LightSlateGray,   Halo = new Pen(Color.Silver, 5),   HorizontalAlignment = LabelStyle.HorizontalAlignmentEnum.Center,   VerticalAlignment = LabelStyle.VerticalAlignmentEnum.Middle,   Font = new Font(GenericFontFamilies.SansSerif.ToString(), 16f, FontStyle.Bold)  });  bool ignoreLength = false;  var layRiverLabel = new LabelLayer("River labels")  {  DataSource = layRivers.DataSource,  Enabled = true,  LabelColumn = "Name",  TextRenderingHint = TextRenderingHint.AntiAlias,  SmoothingMode = SmoothingMode.AntiAlias,  SRID = 4326,  LabelFilter = LabelCollisionDetection.ThoroughCollisionDetection,  MultipartGeometryBehaviour = LabelLayer.MultipartGeometryBehaviourEnum.CommonCenter,  Style =      new LabelStyle      {      ForeColor = Color.DarkBlue,      Font = new Font(FontFamily.GenericSansSerif, 11),      HorizontalAlignment = LabelStyle.HorizontalAlignmentEnum.Center,      VerticalAlignment = LabelStyle.VerticalAlignmentEnum.Middle,      //CollisionDetection = true,      Halo = new Pen(Color.Azure, 2),      IgnoreLength = ignoreLength,      Offset = new PointF(0, -10)      },  };  //Add the layers to the map object.  //The order we add them in are the order they are drawn, so we add the rivers last to put them on top  //map.BackgroundLayer.Add(AsyncLayerProxyLayer.Create(layCountries));  map.Layers.Add(layCountries);  map.Layers.Add(layRivers);  map.Layers.Add(layCities);  map.Layers.Add(layLabel);  map.Layers.Add(layCityLabel);  map.Layers.Add(layRiverLabel);  //增加Layers  var xlsPath = string.Format(XlsConnectionString, System.IO.Directory.GetCurrentDirectory(), "GeoData//Henan.xls", Properties.Settings.Default.OleDbProvider);  var ds = GetDataFromExcel(xlsPath, "Cities");  //var ct = GetCoordinateTransformation();  //TransCoordinate(ds, ct);  string columeName = "Rotation";  //Add Rotation Column  AddColumeToDataSet(ds, columeName, -angle);  var xlsLayer = GetLayerFromDataSet2(ds, Color.GreenYellow);//Set up provider  map.Layers.Add(xlsLayer); //Add layer to map     //limit the zoom to 360 degrees width     //map.MaximumZoom = 360;     //map.BackColor = Color.LightBlue;  //map.Zoom = 360;  map.Center = xlsLayer.Envelope.Centre;// new Point(0, 0);  map.MapScale = 350;  //Matrix mat = new Matrix();  //mat.RotateAt(angle, map.WorldToImage(map.Center));  //map.MapTransform = mat;  //map.ZoomToBox(xlsLayer.Envelope);  return map; } /// <summary> /// 在线显示，圆点显示轨迹 /// </summary> /// <param name="angle"></param> /// <returns></returns> private static Map InitializeMapOsmWithXls(float angle) {  var map = new Map();  var tileLayer = new TileAsyncLayer(  KnownTileSources.Create(KnownTileSource.OpenStreetMap), "TileLayer - OSM with XLS");  tileLayer.SRID = 4326;  map.BackgroundLayer.Add(tileLayer);  //Get data from excel  var xlsPath = string.Format(XlsConnectionString, System.IO.Directory.GetCurrentDirectory(), "GeoData//Cities.xls", Properties.Settings.Default.OleDbProvider);  var ds = GetDataFromExcel(xlsPath, "Cities");  var ds1 = GetDataFromExcel(xlsPath, "Cities2");  var ct = GetCoordinateTransformation();  TransCoordinate(ds, ct);  TransCoordinate(ds1, ct);  string columeName = "Rotation";  //Add Rotation Column  AddColumeToDataSet(ds, columeName, -angle);  AddColumeToDataSet(ds1, columeName, -angle);  var xlsLayer = GetLayerFromDataSet(ds, Color.GreenYellow);//Set up provider  map.Layers.Add(xlsLayer); //Add layer to map  var xlsLayer1 = GetLayerFromDataSet(ds1, Color.Red);  map.Layers.Add(xlsLayer1);  var xlsLabelLayer = GetLabelLayerByVectorLayer(xlsLayer, "XLSLabel");  xlsLabelLayer.Theme = new SharpMap.Rendering.Thematics.FontSizeTheme(xlsLabelLayer, map) { FontSizeScale = 1000f };  map.Layers.Add(xlsLabelLayer);  map.ZoomToBox(xlsLayer.Envelope.ExpandedBy(xlsLayer1.Envelope));  return map; } /// <summary> /// 在线显示，图标显示轨迹 /// </summary> /// <param name="angle"></param> /// <returns></returns> private static Map InitializeMapOsmWithXls2(float angle) {  var map = new Map();  var tileLayer = new TileAsyncLayer(  KnownTileSources.Create(KnownTileSource.OpenStreetMap), "TileLayer - OSM with XLS");  tileLayer.SRID = 4326;  map.BackgroundLayer.Add(tileLayer);  //Get data from excel  var xlsPath = string.Format(XlsConnectionString, System.IO.Directory.GetCurrentDirectory(), "GeoData//Henan.xls", Properties.Settings.Default.OleDbProvider);  var ds = GetDataFromExcel(xlsPath, "Cities");  var ct = GetCoordinateTransformation();  TransCoordinate(ds, ct);  string columeName = "Rotation";  //Add Rotation Column  AddColumeToDataSet(ds, columeName, -angle);  var xlsLayer = GetLayerFromDataSet2(ds, Color.GreenYellow);//Set up provider  map.Layers.Add(xlsLayer); //Add layer to map  var xlsLabelLayer = GetLabelLayerByVectorLayer(xlsLayer, "XLSLabel");  xlsLabelLayer.Theme = new FontSizeTheme(xlsLabelLayer, map) { FontSizeScale = 1000f };  map.Layers.Add(xlsLabelLayer);  map.ZoomToBox(xlsLayer.Envelope);  return map; } /// <summary> /// 从Excel中读取数据 /// </summary> private static DataSet GetDataFromExcel(string xlsPath, string sheetName) {  DataSet ds = new DataSet("XLS");  string sql = string.Format("SELECT * FROM [{0}$];", sheetName);  using (var cn = new OleDbConnection(xlsPath))  {  cn.Open();  using (var da = new OleDbDataAdapter(new OleDbCommand(sql, cn)))  {   da.Fill(ds);  }  }  return ds; } /// <summary> /// 获取坐标转换对象 /// </summary> /// <returns></returns> private static ICoordinateTransformation GetCoordinateTransformation() {  //The SRS for this datasource is EPSG:4326, therefore we need to transfrom it to OSM projection  var ctf = new CoordinateTransformationFactory();  var cf = new CoordinateSystemFactory();  var epsg4326 = cf.CreateFromWkt("GEOGCS[/"WGS 84/",DATUM[/"WGS_1984/",SPHEROID[/"WGS 84/",6378137,298.257223563,AUTHORITY[/"EPSG/",/"7030/"]],AUTHORITY[/"EPSG/",/"6326/"]],PRIMEM[/"Greenwich/",0,AUTHORITY[/"EPSG/",/"8901/"]],UNIT[/"degree/",0.01745329251994328,AUTHORITY[/"EPSG/",/"9122/"]],AUTHORITY[/"EPSG/",/"4326/"]]");  var epsg3857 = cf.CreateFromWkt("PROJCS[/"Popular Visualisation CRS / Mercator/", GEOGCS[/"Popular Visualisation CRS/", DATUM[/"Popular Visualisation Datum/", SPHEROID[/"Popular Visualisation Sphere/", 6378137, 0, AUTHORITY[/"EPSG/",/"7059/"]], TOWGS84[0, 0, 0, 0, 0, 0, 0], AUTHORITY[/"EPSG/",/"6055/"]],PRIMEM[/"Greenwich/", 0, AUTHORITY[/"EPSG/", /"8901/"]], UNIT[/"degree/", 0.0174532925199433, AUTHORITY[/"EPSG/", /"9102/"]], AXIS[/"E/", EAST], AXIS[/"N/", NORTH], AUTHORITY[/"EPSG/",/"4055/"]], PROJECTION[/"Mercator/"], PARAMETER[/"False_Easting/", 0], PARAMETER[/"False_Northing/", 0], PARAMETER[/"Central_Meridian/", 0], PARAMETER[/"Latitude_of_origin/", 0], UNIT[/"metre/", 1, AUTHORITY[/"EPSG/", /"9001/"]], AXIS[/"East/", EAST], AXIS[/"North/", NORTH], AUTHORITY[/"EPSG/",/"3857/"]]");  var ct = ctf.CreateFromCoordinateSystems(epsg4326, epsg3857);  return ct; } /// <summary> /// 转换地球经纬度到坐标 /// </summary> /// <param name="ds"></param> /// <param name="ct"></param> private static void TransCoordinate(DataSet ds, ICoordinateTransformation ct) {  foreach (System.Data.DataRow row in ds.Tables[0].Rows)  {  if (row["X"] == DBNull.Value || row["Y"] == DBNull.Value) continue;  var coords = new[] { Convert.ToDouble(row["X"]), Convert.ToDouble(row["Y"]) };  coords = ct.MathTransform.Transform(coords);  row["X"] = coords[0];  row["Y"] = coords[1];  } } /// <summary> /// 增加列 /// </summary> /// <param name="ds"></param> /// <param name="columeName"></param> /// <param name="columeValue"></param> private static void AddColumeToDataSet(DataSet ds, string columeName, float columeValue) {  ds.Tables[0].Columns.Add(columeName, typeof(float));  foreach (System.Data.DataRow row in ds.Tables[0].Rows)  {  row["Rotation"] = -columeValue;  } } /// <summary> /// 轨迹用点表示 /// </summary> /// <param name="ds"></param> /// <param name="c"></param> /// <returns></returns> private static VectorLayer GetLayerFromDataSet(DataSet ds, Color c) {  var xlsProvider = new DataTablePoint(ds.Tables[0], "OID", "X", "Y");  var xlsLayer = new VectorLayer("XLS", xlsProvider)  { Style = new VectorStyle() { PointColor = new SolidBrush(c) } };  return xlsLayer; } /// <summary> /// 获取带图标的图层 /// </summary> /// <param name="ds"></param> /// <param name="c"></param> /// <returns></returns> private static VectorLayer GetLayerFromDataSet2(DataSet ds, Color c) {  var xlsProvider = new DataTablePoint(ds.Tables[0], "OID", "X", "Y");  var xlsLayer = new VectorLayer("XLS", xlsProvider)  { Style = { Symbol=Properties.Resources.redflag} };  return xlsLayer; } private static LabelLayer GetLabelLayerByVectorLayer(VectorLayer xlsLayer, string layerName) {  var xlsLabelLayer = new LabelLayer(layerName)  {  DataSource = xlsLayer.DataSource,  LabelColumn = "NAME",  //PriorityColumn = "Population",  Style =   {   CollisionBuffer = new System.Drawing.SizeF(2f, 2f),   CollisionDetection = true   },  LabelFilter = LabelCollisionDetection.ThoroughCollisionDetection  };  return xlsLabelLayer; } } public enum MapType { ShapeFile = 0, MapInfo = 1, RunLine = 2,//运行轨迹 Static = 3 //定点数据 }}