C#+GDAL读写文件

读取shp文件：

  private void btnBrower_Click(object sender, EventArgs e)
         {
             OpenFileDialog dlg = new OpenFileDialog();
             dlg.Title = "打开ShapeFile数据";
             dlg.Filter = "ShapeFile数据(*.shp)|*.shp";
             if (dlg.ShowDialog() == DialogResult.OK)
             {
                 Ogr.RegisterAll();

                 string strVectorFile = dlg.FileName;
                 textBox1.Text = strVectorFile;
                 //打开数据
                 DataSource ds = Ogr.Open(strVectorFile, );
                 if (ds == null)
                 {
                     listBox1.Items.Add(string.Format("打开文件【{0}】失败！", strVectorFile));
                     return;
                 }
                 listBox1.Items.Add(string.Format("打开文件【{0}】成功！", strVectorFile));

                 // 获取该数据源中的图层个数，一般shp数据图层只有一个，如果是mdb、dxf等图层就会有多个
                 int iLayerCount = ds.GetLayerCount();

                 // 获取第一个图层
                 Layer oLayer = ds.GetLayerByIndex();
                 if (oLayer == null)
                 {
                     listBox1.Items.Add(string.Format("获取第{0}个图层失败！\n", ));
                     return;
                 }

                 // 对图层进行初始化，如果对图层进行了过滤操作，执行这句后，之前的过滤全部清空
                 oLayer.ResetReading();

                 // 通过属性表的SQL语句对图层中的要素进行筛选，这部分详细参考SQL查询章节内容
                 //oLayer.SetAttributeFilter("\"NAME99\"LIKE \"北京市市辖区\"");

                 // 通过指定的几何对象对图层中的要素进行筛选
                 //oLayer.SetSpatialFilter();  

                 // 通过指定的四至范围对图层中的要素进行筛选
                 //oLayer.SetSpatialFilterRect();  

                 // 获取图层中的属性表表头并输出
                 listBox1.Items.Add("属性表结构信息：");
                 FeatureDefn oDefn = oLayer.GetLayerDefn();
                 int iFieldCount = oDefn.GetFieldCount();
                 for (int iAttr = ; iAttr < iFieldCount; iAttr++)
                 {
                     FieldDefn oField = oDefn.GetFieldDefn(iAttr);

                     listBox1.Items.Add(string.Format("{0}:{1} ({2}.{3})", oField.GetNameRef(),
                              oField.GetFieldTypeName(oField.GetFieldType()),
                              oField.GetWidth(), oField.GetPrecision()));
                 }
                 // 输出图层中的要素个数
                 listBox1.Items.Add(string.Format("要素个数 = {0}", oLayer.GetFeatureCount()));
                 Feature oFeature = null;
                 // 下面开始遍历图层中的要素
                 while ((oFeature = oLayer.GetNextFeature()) != null)
                 {
                     Geometry geo = oFeature.GetGeometryRef();
                     wkbGeometryType wkb = geo.GetGeometryType();
                     listBox1.Items.Add(string.Format("当前处理第要素值：{0}", wkb.ToString()));
                     string strGml=geo.ExportToGML();
                     listBox1.Items.Add(strGml);
                     listBox1.Items.Add(string.Format("当前处理第{0}个: \n属性值：", oFeature.GetFID()));

                     // 获取要素中的属性表内容
                     for (int iField = ; iField < iFieldCount; iField++)
                     {
                         FieldDefn oFieldDefn = oDefn.GetFieldDefn(iField);
                         FieldType type = oFieldDefn.GetFieldType();
                         switch (type)
                         {
                             case FieldType.OFTString:
                                 listBox1.Items.Add(string.Format("{0}\t", oFeature.GetFieldAsString(iField)));
                                 break;
                             case FieldType.OFTReal:
                                 listBox1.Items.Add(string.Format("{0}\t", oFeature.GetFieldAsDouble(iField)));
                                 break;
                             case FieldType.OFTInteger:
                                 listBox1.Items.Add(string.Format("{0}\t", oFeature.GetFieldAsInteger(iField)));
                                 break;
                             default:
                                 listBox1.Items.Add(string.Format("{0}\t", oFeature.GetFieldAsString(iField)));
                                 break;
                         }
                     }
                     // 获取要素中的几何体
                     Geometry oGeometry = oFeature.GetGeometryRef();
                     // 为了演示，只输出一个要素信息
                     break;
                 }
                 listBox1.Items.Add("数据集关闭！");
             }
         }

读取个人地理数据库mdb：

  private void btnBrower_Click(object sender, EventArgs e)
         {
             OpenFileDialog dlg = new OpenFileDialog();
             dlg.Title = "打开Personal Geodatabase";
             dlg.Filter = "Personal Geodatabase数据(*.mdb)|*.mdb";
             if (dlg.ShowDialog() == DialogResult.OK)
             {
                 Ogr.RegisterAll();

                 string strVectorFile = dlg.FileName;
                 textBox1.Text = strVectorFile;
                 string strInfo = string.Empty;
                 //打开数据
                 DataSource ds = Ogr.Open(strVectorFile, );
                 if (ds == null)
                 {
                     strInfo += string.Format("打开文件【{0}】失败！\n", strVectorFile);
                     return;
                 }
                 strInfo += string.Format("打开文件【{0}】成功！\n", strVectorFile);

                 // 获取该数据源中的图层个数，一般shp数据图层只有一个，如果是mdb、dxf等图层就会有多个
                 int iLayerCount = ds.GetLayerCount();
                 strInfo += string.Format("个人地理数据库共包含{0}个图层！\n", iLayerCount);
                 // 获取第一个图层
                 Layer oLayer = ds.GetLayerByIndex();
                 if (oLayer == null)
                 {
                     strInfo += string.Format("获取第{0}个图层失败！\n", );
                     return;
                 }
                 strInfo += string.Format("第0个图层名称：{0}\n", oLayer.GetName());
                 // 对图层进行初始化，如果对图层进行了过滤操作，执行这句后，之前的过滤全部清空
                 oLayer.ResetReading();

                 // 通过属性表的SQL语句对图层中的要素进行筛选，这部分详细参考SQL查询章节内容
                 //oLayer.SetAttributeFilter("\"NAME99\"LIKE \"北京市市辖区\"");

                 // 通过指定的几何对象对图层中的要素进行筛选
                 //oLayer.SetSpatialFilter();  

                 // 通过指定的四至范围对图层中的要素进行筛选
                 //oLayer.SetSpatialFilterRect();  

                 // 获取图层中的属性表表头并输出
                 strInfo += "属性表结构信息：\n";
                 FeatureDefn oDefn = oLayer.GetLayerDefn();
                 int iFieldCount = oDefn.GetFieldCount();
                 for (int iAttr = ; iAttr < iFieldCount; iAttr++)
                 {
                     FieldDefn oField = oDefn.GetFieldDefn(iAttr);

                     strInfo += string.Format("{0}:{1} ({2}.{3})\n", oField.GetNameRef(),
                              oField.GetFieldTypeName(oField.GetFieldType()),
                              oField.GetWidth(), oField.GetPrecision());
                 }
                 // 输出图层中的要素个数
                 strInfo += string.Format("要素个数 = {0}\n", oLayer.GetFeatureCount());
                 Feature oFeature = null;
                 // 下面开始遍历图层中的要素
                 while ((oFeature = oLayer.GetNextFeature()) != null)
                 {
                     Geometry geo = oFeature.GetGeometryRef();
                     if (geo==null)
                     {
                         break;
                     }
                     wkbGeometryType wkb = geo.GetGeometryType();
                     strInfo += string.Format("当前处理第{0}要素值：{1}\n", oFeature.GetFID(), wkb.ToString());
                     string strGml = geo.ExportToGML();
                     strInfo += strGml;
                     strInfo += string.Format("\n当前处理第{0}个: \n属性值：", oFeature.GetFID());

                     // 获取要素中的属性表内容
                     for (int iField = ; iField < iFieldCount; iField++)
                     {
                         FieldDefn oFieldDefn = oDefn.GetFieldDefn(iField);
                         FieldType type = oFieldDefn.GetFieldType();
                         switch (type)
                         {
                             case FieldType.OFTString:
                                 strInfo += string.Format("{0}\t", oFeature.GetFieldAsString(iField));
                                 break;
                             case FieldType.OFTReal:
                                 strInfo += string.Format("{0}\t", oFeature.GetFieldAsDouble(iField));
                                 break;
                             case FieldType.OFTInteger:
                                 strInfo += string.Format("{0}\t", oFeature.GetFieldAsInteger(iField));
                                 break;
                             default:
                                 strInfo += string.Format("{0}\t", oFeature.GetFieldAsString(iField));
                                 break;
                         }
                     }
                     // 获取要素中的几何体
                     Geometry oGeometry = oFeature.GetGeometryRef();
                     // 为了演示，只输出一个要素信息
                     break;
                 }
                 strInfo += "\n数据集关闭！";
                 richTextBox1.Text = strInfo;
             }
         }