public XLayer ReadShapefile(string shpfilename)
{
//打开文件和读取工具
FileStream fsr = new FileStream(shpfilename, FileMode.Open);
BinaryReader br = new BinaryReader(fsr);
//读取文件头
ShapefileHeader sfh = (ShapefileHeader)
XTools.FromBytes(br, typeof(ShapefileHeader));
//获得空间对象类型
SHAPETYPE ShapeType = (SHAPETYPE)Enum.Parse(typeof(SHAPETYPE), sfh.ShapeType.ToString());
//获得空间范围
XExtent extent = new XExtent(new XVertex(sfh.Xmin, sfh.Ymin),
new XVertex(sfh.Xmax, sfh.Ymax));
//读dbf
string dbfFileName = shpfilename.Replace(".shp", ".dbf");
List <XField> fields = new List <XField>();
List <XAttribute> attributes = new List <XAttribute>();
ReadDBFFile(dbfFileName, fields, attributes);
//构建图层
XLayer layer = new XLayer(shpfilename, ShapeType, extent, fields);
int index = 0;
while (br.PeekChar() != -1)
{
//读记录头
RecordHeader rh = (RecordHeader)
XTools.FromBytes(br, typeof(RecordHeader));
int RecordLength = FromBigToLittle(rh.RecordLength) * 2 - 4;
byte[] RecordContent = br.ReadBytes(RecordLength);
//开始读实际的空间数据
if (ShapeType == SHAPETYPE.point)
{
XPointSpatial point = ReadPoint(RecordContent);
layer.AddFeature(new XFeature(point, attributes[index]));
}
else if (ShapeType == SHAPETYPE.line)
{
List <XLineSpatial> lines = ReadLines(RecordContent);
foreach (XLineSpatial line in lines)
{
layer.AddFeature(new XFeature(line, new XAttribute(attributes[index])));
}
}
else if (ShapeType == SHAPETYPE.polygon)
{
List <XPolygonSpatial> polygons = ReadPolygons(RecordContent);
foreach (XPolygonSpatial polygon in polygons)
{
layer.AddFeature(new XFeature(polygon, new XAttribute(attributes[index])));
}
}
index++;
}
//关闭读取工具和文件
br.Close();
fsr.Close();
return(layer);
}
public XLayer ReadFile(string filename)
{
FileStream fsr = new FileStream(filename, FileMode.Open);
BinaryReader br = new BinaryReader(fsr);
//读文件头
MyFileHeader mfh = (MyFileHeader)XTools.FromBytes(br, typeof(MyFileHeader));
//读图层名称
string name = XTools.ReadString(br);
//读属性字段结构
List <XField> fields = ReadFields(mfh.FieldCount, br);
//定义图层
SHAPETYPE ShapeType =
(SHAPETYPE)Enum.Parse(typeof(SHAPETYPE),
mfh.ShapeType.ToString());
XExtent extent = new XExtent(
new XVertex(mfh.MinX, mfh.MinY),
new XVertex(mfh.MaxX, mfh.MaxY));
XLayer layer = new XLayer(name, ShapeType, extent, fields);
//读空间对象类型
for (int i = 0; i < mfh.FeatureCount; i++)
{
XSpatial spatial = ReadSpatial(ShapeType, br);
XAttribute attribute = ReadAttribute(br, fields);
layer.AddFeature(new XFeature(spatial, attribute));
}
//关闭文件并返回结果
br.Close();
fsr.Close();
return(layer);
}
public List <XLayer> ReadJSONFile(string filename)
{
List <XLayer> layers = new List <XLayer>();
StreamReader sr = new StreamReader(filename, Encoding.UTF8);
Newtonsoft.Json.Linq.JObject o = Newtonsoft.Json.Linq.JObject.Parse(sr.ReadToEnd()); //读入一个json文件
Newtonsoft.Json.Linq.JToken features = o["features"]; //读取所有空间对象数组
List <Newtonsoft.Json.Linq.JToken> flst = features.ToList();
List <XField> fields = new List <XField>(); /*ReadFields(mfh.FieldCount, br);*/
XLayer pointlayer = new XLayer(filename, SHAPETYPE.point, null, fields);
XLayer linelayer = new XLayer(filename, SHAPETYPE.line, null, fields);
XLayer polygonlayer = new XLayer(filename, SHAPETYPE.polygon, null, fields);
double maxx = Double.MinValue;
double minx = Double.MaxValue;
double maxy = Double.MinValue;
double miny = Double.MaxValue;
for (int index = 0; index < flst.Count; index++)
{
List <XVertex> _Vertexes = new List <XVertex>(); //记录feature的节点
string featuretype = (string)(flst[index]["geometry"]["type"]); //读取对象类型 包括Point LineString Polygon
Newtonsoft.Json.Linq.JToken properties = flst[index]["properties"]; //读取对象的所有字段与属性
//if (featuretype == "Polygon") flst[index]["geometry"]["coordinates"][0]
if (featuretype == "Point")//如果是点对象 则按点对象的数据结构进行处理
{
var location = flst[index]["geometry"]["coordinates"];
double lon = (double)location.First;
double lat = (double)location.Last;
Console.WriteLine(lon);
if (lon > maxx)
{
maxx = lon;
}
else if (lon < minx)
{
minx = lon;
}
if (lat > maxy)
{
maxy = lat;
}
else if (lat < miny)
{
miny = lat;
}
_Vertexes.Add(new XVertex(lon, lat));
ReadJSONField(properties, pointlayer);//右边返回一个Fields泛型赋给左边图层字段
}
else if (featuretype == "LineString")//如果是线对象 则按线对象的数据结构进行处理
{
foreach (var location in flst[index]["geometry"]["coordinates"])
{
double lon = (double)location.First;
double lat = (double)location.Last;
//Console.WriteLine(lon);
if (lon > maxx)
{
maxx = lon;
}
else if (lon < minx)
{
minx = lon;
}
if (lat > maxy)
{
maxy = lat;
}
else if (lat < miny)
{
miny = lat;
}
_Vertexes.Add(new XVertex(lon, lat));
ReadJSONField(properties, linelayer);//右边返回一个Fields泛型赋给左边图层字段
}
}
else if (featuretype == "Polygon")//如果是多边形对象 则按多边形对象的数据结构进行处理
{
foreach (var singlepoly in flst[index]["geometry"]["coordinates"])
{
foreach (var location in singlepoly)
{
double lon = (double)location.First;
double lat = (double)location.Last;
Console.WriteLine(lon);
if (lon > maxx)
{
maxx = lon;
}
else if (lon < minx)
{
minx = lon;
}
if (lat > maxy)
{
maxy = lat;
}
else if (lat < miny)
{
miny = lat;
}
_Vertexes.Add(new XVertex(lon, lat));
ReadJSONField(properties, polygonlayer);//右边返回一个Fields泛型赋给左边图层字段
}
}
}
if (featuretype == "Point")
{
XFeature onefeature = new XFeature(new XPointSpatial(_Vertexes[0]), null); //暂时默认无属性
pointlayer.AddFeature(onefeature); //给线图层添加对象
}
else if (featuretype == "LineString")
{
XFeature onefeature = new XFeature(new XLineSpatial(_Vertexes), null); //暂时默认无属性
linelayer.AddFeature(onefeature); //给线图层添加对象
}
else if (featuretype == "Polygon")
{
XFeature onefeature = new XFeature(new XPolygonSpatial(_Vertexes), null); //暂时默认无属性
polygonlayer.AddFeature(onefeature); //给线图层添加对象
}
}
//Console.WriteLine(minx.ToString(), miny.ToString(), maxx.ToString(), maxy.ToString());
XExtent layersextent = new XExtent(new XVertex(minx, miny), new XVertex(maxx, maxy));
layers.Add(polygonlayer);
layers.Add(linelayer);
layers.Add(pointlayer);
foreach (XLayer ll in layers)
{
ll.Extent = layersextent;//给三种layer赋上共同的extent
ll.DrawAttributeOrNot = false;
}
sr.Close();
return(layers);
}
