GDAL(Geospatial Data Abstraction Library)是一个开源栅格空间数据转换库。它利用抽象数据模型来表达所支持的各种文件格式。它还有一系列命令行工具来进行数据转换和处理。
参考网址：http://www.gisinternals.com/sdk/
1，GDAL-java写入shp文件
WriteShp.java代码：
package Vector; import org.gdal.gdal.gdal; import org.gdal.ogr.*; import org.gdal.osr.SpatialReference; //将空间数据写入shp文件 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_DATA路径（gdal根目录下的bin\gdal-data） gdal.SetConfigOption("GDAL_DATA","F:\\GDAL学习文件夹\\release-1900-x64-gdal-2-3-2-mapserver-7-2-1\\bin\\gdal-data"); // 为了支持中文路径，请添加下面这句代码 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.ImportFromEPSG(4326); 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"); } }
2，GDAL-java读取shp文件信息
ReadShp.java代码:
package Vector; import org.gdal.gdal.gdal; import org.gdal.ogr.*; import org.gdal.osr.SpatialReference; import java.util.HashMap; import java.util.Map; //读取某个磁盘路径下的shp文件信息（范围、属性等） public class ReadShp { public static void main(String[] args) throws Exception { //指定文件的名字和路径 String strVectorFile ="E:\\test\\test.shp"; // 注册所有的驱动 ogr.RegisterAll(); //配置GDAL_DATA路径（gdal根目录下的bin\gdal-data） gdal.SetConfigOption("GDAL_DATA","F:\\GDAL学习文件夹\\release-1900-x64-gdal-2-3-2-mapserver-7-2-1\\bin\\gdal-data"); // 为了支持中文路径，请添加下面这句代码 gdal.SetConfigOption("GDAL_FILENAME_IS_UTF8","YES"); // 为了使属性表字段支持中文，请添加下面这句 gdal.SetConfigOption("SHAPE_ENCODING","CP936"); //读取数据，这里以ESRI的shp文件为例 String strDriverName = "ESRI Shapefile"; //创建一个文件，根据strDriverName扩展名自动判断驱动类型 org.gdal.ogr.Driver oDriver =ogr.GetDriverByName(strDriverName); if (oDriver == null) { System.out.println(strDriverName+ " 驱动不可用！\n"); return; } DataSource dataSource = oDriver.Open(strVectorFile); Layer layer = dataSource.GetLayer("test"); String layerName = layer.GetName(); System.out.println("图层名称："+layerName); SpatialReference spatialReference = layer.GetSpatialRef(); //System.out.println(spatialReference); System.out.println("空间参考坐标系："+spatialReference.GetAttrValue("AUTHORITY",0)+spatialReference.GetAttrValue("AUTHORITY",1)); double[] layerExtent = layer.GetExtent(); System.out.println("图层范围：minx:"+layerExtent[0]+",maxx:"+layerExtent[1]+",miny:"+layerExtent[2]+",maxy:"+layerExtent[3]); FeatureDefn featureDefn = layer.GetLayerDefn(); int fieldCount = featureDefn.GetFieldCount(); Map<String,Object> fieldMap = new HashMap(); for(int i=0; i<fieldCount; i++){ FieldDefn fieldDefn = featureDefn.GetFieldDefn(i); //得到属性字段类型 int fieldType = fieldDefn.GetFieldType(); String fieldTypeName = fieldDefn.GetFieldTypeName(fieldType); //得到属性字段名称 String fieldName = fieldDefn.GetName(); fieldMap.put(fieldTypeName,fieldName); } long featureCount = layer.GetFeatureCount(); System.out.println("图层要素个数："+featureCount); for(int i=0; i<featureCount; i++){ Feature feature = layer.GetFeature(i); Object[] arr = fieldMap.values().toArray(); for(int k=0; k<arr.length; k++){ String fvalue = feature.GetFieldAsString(arr[k].toString()); System.out.print(" 属性名称:"+arr[k].toString()+",属性值:"+fvalue); } System.out.println(); } } }
3，GDAL-java读取tiff文件信息
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]); } } }
4，shp转换geojson
package gdal; import org.gdal.*; import org.gdal.ogr.DataSource; import org.gdal.ogr.Driver; import org.gdal.ogr.ogr; import org.gdal.gdal.*; public class Test { public static void main(String[] args) { // 注册所有的驱动 ogr.RegisterAll(); // 为了支持中文路径，请添加下面这句代码 gdal.SetConfigOption("GDAL_FILENAME_IS_UTF8","YES"); // 为了使属性表字段支持中文，请添加下面这句 gdal.SetConfigOption("SHAPE_ENCODING",""); //shp文件所在的位置 String strVectorFile = "D:\\sichuan\\sichuanPointALL.shp"; //打开数据 DataSource ds = ogr.Open(strVectorFile,0); if (ds == null) { System.out.println("打开文件失败！" ); return; } System.out.println("打开文件成功！" ); Driver dv = ogr.GetDriverByName("GeoJSON"); if (dv == null) { System.out.println("打开驱动失败！" ); return; } System.out.println("打开驱动成功！" ); //输出geojson的位置及文件名 dv.CopyDataSource(ds, "D:\\data\\sichuan.geojson"); System.out.println("转换成功！" ); } }
【参考】
GDAL项目官方网站：http://www.gdal.org
GDAL在Java中的使用：http://trac.osgeo.org/gdal/wiki/GdalOgrInJava
