GDAL 的java调用问题总结

文章主要是对http://blog.csdn.net/liminlu0314/article/details/8395622该博文的一些问题的总结

首先，想用java编写GDAL的人不一定想从源码上编译GDAL，只是想用GDAL比较成熟的类库，所以编译环节其实可以省略，直接下载编译好的java版GDAL，例如下载地址：http://download.csdn.net/detail/liminlu0314/5016578，这个库的上传者与文章作者同一个人，是GDAL方面的大牛

其次，gdal.jar是需要导入工程的，但是其他的由jni生成的的dll，其实并不是添加jar引用，因为dll并不是jar包，会报不是zip的错误，我们只需要把几个dll文件放在工程根目录即可。如下图，用IntelliJ IDEA创建的工程，导入dll，注意，除了java文件夹里文件名带jni的dll，其中别忘了还得把GDAL\bin里的gdal19.dll拷过来，不然会报找不到类库的错误。

再次，参看http://blog.csdn.net/liminlu0314/article/details/8828907这篇文章，其中读写代码都有一些错误，

以下指出一些错误，对于一些newFeature等单词没拆分的错误就不一一列举了，只是对能想起来的标红以下

package gdal_test;


import org.gdal.gdal.gdal;
import org.gdal.ogr.DataSource;
import org.gdal.ogr.Feature;
import org.gdal.ogr.FeatureDefn;
import org.gdal.ogr.FieldDefn;
import org.gdal.ogr.Geometry;
import org.gdal.ogr.Layer;
import org.gdal.ogr.ogr;


public class MyTest {


public static void main(String[] args) {
//ReadVectorFile();
WriteVectorFile();


}

static void ReadVectorFile()
{  
  String strVectorFile = "D:\\s57polygon.shp";
  
        // 注册所有的驱动  
        ogr.RegisterAll();  
  
        // 为了支持中文路径，请添加下面这句代码  
        gdal.SetConfigOption("GDAL_FILENAME_IS_UTF8","YES");  
        // 为了使属性表字段支持中文，请添加下面这句  
        gdal.SetConfigOption("SHAPE_ENCODING","");  
  
        //打开数据  
        DataSource ds = ogr.Open(strVectorFile,0);  
       
        if (ds == null)  
        {  
                  System.out.println("打开文件【"+ strVectorFile + "】失败！" );  
                  return;  
        }  
         
        System.out.println(ds.GetDriver().getName());
        System.out.println("打开文件【"+ strVectorFile + "】成功！" );  
  
        // 获取该数据源中的图层个数，一般shp数据图层只有一个，如果是mdb、dxf等图层就会有多个  
        int iLayerCount = ds.GetLayerCount();  
  
        // 获取第一个图层  
        Layer oLayer = ds.GetLayerByIndex(0);  
        if (oLayer == null)  
        {  
                  System.out.println("获取第0个图层失败！\n");  
                  return;  
        }  
  
        // 对图层进行初始化，如果对图层进行了过滤操作，执行这句后，之前的过滤全部清空  
        oLayer.ResetReading();  
  
        // 通过属性表的SQL语句对图层中的要素进行筛选，这部分详细参考SQL查询章节内容  
        //oLayer.SetAttributeFilter("\"NAME99\"LIKE \"北京市市辖区\"");  
  
        // 通过指定的几何对象对图层中的要素进行筛选  
        //oLayer.SetSpatialFilter();  
  
        // 通过指定的四至范围对图层中的要素进行筛选  
        //oLayer.SetSpatialFilterRect();  
  
        // 获取图层中的属性表表头并输出  
        System.out.println("属性表结构信息：");  
        FeatureDefn oDefn =oLayer.GetLayerDefn();  
        int iFieldCount =oDefn.GetFieldCount();  
        for (int iAttr = 0; iAttr <iFieldCount; iAttr++)  
        {  
                  FieldDefn oField =oDefn.GetFieldDefn(iAttr);  
  
                  System.out.println(oField.GetNameRef()+ ": " +  
                                    oField.GetFieldTypeName(oField.GetFieldType())+ "(" +  
                           oField.GetWidth()+"."+ oField.GetPrecision() + ")");  
        }  
  
        // 输出图层中的要素个数  
        System.out.println("要素个数 = " + oLayer.GetFeatureCount(0));  
  
        Feature oFeature = null;  
        // 下面开始遍历图层中的要素  
        while ((oFeature =oLayer.GetNextFeature()) != null)  
        {  
                  System.out.println("当前处理第" + oFeature.GetFID() + "个:\n属性值：");  
                  // 获取要素中的属性表内容  
                  for (int iField = 0; iField< iFieldCount; iField++)  
                  {  
                           FieldDefn oFieldDefn= oDefn.GetFieldDefn(iField);  
                           int type =oFieldDefn.GetFieldType();  
  
                           switch (type)  
                           {  
                                    case ogr.OFTString: //这个case及以下的case都连写了
                                              System.out.println(oFeature.GetFieldAsString(iField)+ "\t");  
                                              break;  
                                    case ogr.OFTReal:  
                                              System.out.println(oFeature.GetFieldAsDouble(iField)+ "\t");  
                                              break;  
                                    case ogr.OFTInteger:  
                                              System.out.println(oFeature.GetFieldAsInteger(iField)+ "\t");  
                                              break;  
                                    default:  
                                              System.out.println(oFeature.GetFieldAsString(iField)+ "\t");  
                                              break;  
                           }  
                  }  
  
                  // 获取要素中的几何体  
                  Geometry oGeometry =oFeature.GetGeometryRef();  
  
                  // 为了演示，只输出一个要素信息  
                  break;  
        }  
  
        System.out.println("数据集关闭！");  
}  


static void WriteVectorFile()  
{  
        String strVectorFile ="E:\\TestPolygon.shp";  
  
        // 注册所有的驱动  
        ogr.RegisterAll();  
         
        // 为了支持中文路径，请添加下面这句代码  
        gdal.SetConfigOption("GDAL_FILENAME_IS_UTF8","NO");  
        // 为了使属性表字段支持中文，请添加下面这句  
        gdal.SetConfigOption("SHAPE_ENCODING","");  
  
        //创建数据，这里以创建ESRI的shp文件为例  
        String strDriverName = "ESRI Shapefile"; //这里有非常致命的错误，原文这里两个单词也连写了，后来总是报错“驱动不可用”，我读取了一个shp文件并用getDriver()和getName()获取了shp文件的名称，其实是有空格的。
        org.gdal.ogr.Driver oDriver =ogr.GetDriverByName(strDriverName);  
        if (oDriver == null)  
        {  
                  System.out.println(strVectorFile+ " 驱动不可用！\n");  
                  return;  
        }  
  
        // 创建数据源  
        DataSource oDS = oDriver.CreateDataSource(strVectorFile,null);  
        if (oDS == null)  
        {  
                  System.out.println("创建矢量文件【"+ strVectorFile +"】失败！\n" );  
                  return;  
        }  
  
        // 创建图层，创建一个多边形图层，这里没有指定空间参考，如果需要的话，需要在这里进行指定  
        Layer oLayer =oDS.CreateLayer("TestPolygon", null, ogr.wkbPolygon, null);  
        if (oLayer == null)  
        {  
                  System.out.println("图层创建失败！\n");  
                  return;  
        }  
  
        // 下面创建属性表  
        // 先创建一个叫FieldID的整型属性  
        FieldDefn oFieldID = new FieldDefn("FieldID", ogr.OFTInteger);  
        oLayer.CreateField(oFieldID, 1);  
  
        // 再创建一个叫FeatureName的字符型属性，字符长度为50  
        FieldDefn oFieldName = new FieldDefn("FieldName", ogr.OFTString);  //原文这种地方好像都连写了
        oFieldName.SetWidth(100);  
        oLayer.CreateField(oFieldName, 1);  
  
        FeatureDefn oDefn =oLayer.GetLayerDefn();  
  
        // 创建三角形要素  
        Feature oFeatureTriangle = new Feature(oDefn);  
        oFeatureTriangle.SetField(0, 0);  
        oFeatureTriangle.SetField(1, "三角形");  
        Geometry geomTriangle =Geometry.CreateFromWkt("POLYGON ((0 0,20 0,10 15,0 0))");  
        oFeatureTriangle.SetGeometry(geomTriangle);  
  
        oLayer.CreateFeature(oFeatureTriangle);  
  
        // 创建矩形要素  
        Feature oFeatureRectangle = new Feature(oDefn);  
        oFeatureRectangle.SetField(0, 1);  
        oFeatureRectangle.SetField(1, "矩形");  
        Geometry geomRectangle =Geometry.CreateFromWkt("POLYGON ((30 0,60 0,60 30,30 30,30 0))");  
        oFeatureRectangle.SetGeometry(geomRectangle);  
  
        oLayer.CreateFeature(oFeatureRectangle);  

         //这里是原文没有的，不知道什么原因，创建的shp存在于磁盘上，但是用arcgis打开却没有任何记录，我当时看代码就觉得肯定是只画到内存里了，其实并没有关闭数据流，没有保存到磁盘，所以随便点了点，发现有个同步到磁盘的方法，就试了试，结果就行了
        oLayer.SyncToDisk();
        oDS.SyncToDisk();
  
        // 创建五角形要素  //记得原文创建五角星这里会出现“corrup data”错误（调试的时候能发现错误出在这里），具体为什么没研究，只是直接把这注释掉了。
        /*
        Feature oFeaturePentagon = new Feature(oDefn);  
        oFeaturePentagon.SetField(0, 2);  
        oFeaturePentagon.SetField(1, "五角形");  
        Geometry geomPentagon =Geometry.CreateFromWkt("POLYGON ((70 0,85 0,90 15,80 30,65 15,700))");  
        oFeaturePentagon.SetGeometry(geomPentagon);  
  
        oLayer.CreateFeature(oFeaturePentagon);  
        */
  
        System.out.println("\n数据集创建完成！\n");  
}  


}

用arcgis打开发现有字符编码问题，我个郁闷啊，查了好多资料不知道为啥，调试了一下，发现就是UTF-8编码，后来我狠了狠心，换了个编码

gdal.SetConfigOption("SHAPE_ENCODING","GBK");  发现用GBK编码就没事。。。。天理何在啊