GDAL-java环境配置

gdal在java环境下使用，需要将gdal编译为java可用的版本，具体编译步骤，以后再补，今天先使用别人编译好的java版本gdal，下载地址：http://www.gisinternals.com/，根据自己的需要下载。

配置步骤：

1、下载已编译好的gdal，如我下载的 x64 位的 release-1900-x64-gdal-2-3-2-mapserver-7-2-1， 解压之后，会得到如下图所示文件：

image.png

2、bin目录下文件列表：

image.png

3、配置gdal dll路径

有2种方式测试成功：
第1种：将bin目录下dll扩展名文件（不包含bin子目录下文件夹内的dll），拷贝到 java 工程所用 jdk 目录下 bin 文件内，如jdk路径是：C:\Program Files\Java\jdk1.7.0_79，则需要将gdal bin目录下的dll（不包含bin子目录下文件夹内的dll）拷贝到C:\Program Files\Java\jdk1.7.0_79\bin目录；

第2种：将bin目录下dll扩展名文件（不包含bin子目录下文件夹内的dll）拷贝到java工程根目录下。

image.png

4、新建java工程

需要引入解压gdal后得到的bin目录下gdal\java目录下的 gdal.jar

image.png image.png

WriteShp.java代码：

package Vector;

import org.gdal.gdal.gdal;
import org.gdal.ogr.*;
import org.gdal.osr.SpatialReference;

public class WriteShp {
    public static void main(String[] args) {
        //指定文件的名字和路径
        String strVectorFile ="E:\\test\\test.shp";

        // 注册所有的驱动
        ogr.RegisterAll();

        //这个可以看到支持哪些格式的驱动，不用的时候可以注释掉
        int count =ogr.GetDriverCount();
        for(int i=0;i<count;i++){
            System.out.println(ogr.GetDriver(i).GetName());
        }

        // 为了支持中文路径，请添加下面这句代码
        gdal.SetConfigOption("GDAL_FILENAME_IS_UTF8","YES");

        // 为了使属性表字段支持中文，请添加下面这句
        gdal.SetConfigOption("SHAPE_ENCODING","CP936");

        //创建数据，这里以创建ESRI的shp文件为例
        String strDriverName = "ESRI Shapefile";

        //创建一个文件，根据strDriverName扩展名自动判断是创建shp文件或其他文件
        Driver oDriver =ogr.GetDriverByName(strDriverName);
        if (oDriver == null) {
            System.out.println(strDriverName+ " 驱动不可用！\n");
            return;
        }

        // 创建数据源
        DataSource oDS = oDriver.CreateDataSource(strVectorFile,null);
        if (oDS == null) {
            System.out.println("创建矢量文件【"+ strVectorFile +"】失败！\n" );
            return;
        }

        // 创建图层，创建一个多边形图层，这里没有指定空间参考，如果需要的话，需要在这里进行指定
        //如果是mif或者tab，其实是可以多元素混合存放的，shp每个图层只能存放点线面中的一种。
        //所以tab创建图层的时候，选择什么都不影响后面的操作和存储结果

        SpatialReference spatialReference = new SpatialReference();
        spatialReference.SetWellKnownGeogCS("WGS84");

        Layer oLayer =oDS.CreateLayer("TestPolygon", spatialReference, ogr.wkbPolygon, null);
        if (oLayer == null) {
            System.out.println("图层创建失败！\n");
            return;
        }

        // 下面创建属性表
        // 先创建一个叫FieldID的整型属性
        FieldDefn oFieldID = new FieldDefn("FieldID", ogr.OFTInteger);
        oLayer.CreateField(oFieldID, 1);

        // 再创建一个叫名称的字符型属性，字符长度为50
        FieldDefn oFieldName = new FieldDefn("名称", ogr.OFTString);
        oFieldName.SetWidth(100);
        oLayer.CreateField(oFieldName, 1);

        FeatureDefn oDefn =oLayer.GetLayerDefn();

        // 创建三角形要素
        //使用wkb方式创建几何对象
        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);

        // 创建五角形要素
        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,70 0))");
        oFeaturePentagon.SetGeometry(geomPentagon);
        oLayer.CreateFeature(oFeaturePentagon);

        //写入文件
        oLayer.SyncToDisk();
        oDS.SyncToDisk();

        System.out.println("\n数据集创建完成！\n");
    }
}

ReadTiff.java代码：

package Raster;

import org.gdal.gdal.Band;
import org.gdal.gdal.Dataset;
import org.gdal.gdal.Driver;
import org.gdal.gdal.gdal;
import org.gdal.gdalconst.gdalconstConstants;

public class ReadTiff {
    public static void main(String[] args) {
        gdal.AllRegister();

        String rasterFilePath = "E:\\Test\\TestOpenLayers\\data\\GF1_PMS2_E121.7_N39.9_20170609_L1A0002429674_Fusion.tiff";

        Dataset dataset = gdal.Open(rasterFilePath,
                gdalconstConstants.GA_ReadOnly);
        if (dataset == null) {
            System.out.println("GDAL read error: " + gdal.GetLastErrorMsg());
        }

        Driver driver = dataset.GetDriver();
        System.out.println("driver short name: " + driver.getShortName());
        System.out.println("driver long name: " + driver.getLongName());
        System.out.println("metadata list: " + driver.GetMetadata_List());

        int xsize = dataset.getRasterXSize();
        int ysize = dataset.getRasterYSize();
        int count = dataset.getRasterCount();
        String proj = dataset.GetProjection();
        Band band = dataset.GetRasterBand(1);

        // 左上角点坐标 lon lat: transform[0]、transform[3]
        // 像素分辨率 x、y方向 : transform[1]、transform[5]
        // 旋转角度: transform[2]、transform[4])
        double[] transform = dataset.GetGeoTransform();
        for(int i=0; i<transform.length; i++){
            System.out.println("transform: " + transform[i]);
        }
    }
}

image.png