C#用GDAL/OGR库创建与写Shape文件,测试通过[原创]

2009年7月19日　周日11:24

GDAL/OGR是非常著名的开源GIS库，GDAL是对栅格数据进行操作，而OGR是对矢量数据进行操作，它们相当于一个通用数据访问库。甚至ESRI的产品中都用了此库。

现在，GDAL/OGR也提供了对.NET的编译支持，具体的编译与使用过程大家GOOGLE一下，网上资料非常详尽。其官网上有一些C，C++，Python的示例代码。 但没有C#的示例代码，但程序之间大同小异，下面我试的一个C#创建Shape文件，并往其中添加Point的示例代码，希望对各位有用。

 

 1            //注册Ogr库
 2            string pszDriverName = "ESRI Shapefile";
 3            OSGeo.OGR.Ogr.RegisterAll();
 4
 5            //调用对Shape文件读写的Driver接口
 6            OSGeo.OGR.Driver poDriver = OSGeo.OGR.Ogr.GetDriverByName(pszDriverName);
 7            if (poDriver == null)
 8                MessageBox.Show("Driver Error");
 9
10           //用此Driver创建Shape文件
11            OSGeo.OGR.DataSource poDS;
12            poDS = poDriver.CreateDataSource("point_out.shp", null);
13            if (poDS == null)
14                MessageBox.Show("DataSource Creation Error");
15
16            //创建层Layer
17            OSGeo.OGR.Layer poLayer;
18            poLayer = poDS.CreateLayer("point_out", null, OSGeo.OGR.wkbGeometryType.wkbPoint, null);
19            if (poLayer == null)
20                MessageBox.Show("Layer Creation Failed");
21            
22            //创建属性列两列
23            OSGeo.OGR.FieldDefn oField = new OSGeo.OGR.FieldDefn("名称", OSGeo.OGR.FieldType.OFTString);
24            oField.SetWidth(16);
25            OSGeo.OGR.FieldDefn oField2 = new OSGeo.OGR.FieldDefn("高度", OSGeo.OGR.FieldType.OFTInteger);            
26            poLayer.CreateField(oField, 1);
27            poLayer.CreateField(oField2, 0);
28
29            //创建一个Feature,一个Point
30            OSGeo.OGR.Feature poFeature = new Feature(poLayer.GetLayerDefn());
31            OSGeo.OGR.Geometry pt = new Geometry(OSGeo.OGR.wkbGeometryType.wkbPoint);
32            //往创建的Shape中写入二十个点,
33            for (int i = 100; i < 120; i++)
34            poDS.Dispose();

这里还没有涉及到定义投影坐标系的问题，后面将会继续完善起来。

此文链接：http://www.cnblogs.com/webgis8/archive/2009/07/19/1526425.html

From：giswei之修道院