Java用GDAL读写shapefile的方法示例

GDAL介绍

GDAL(Geospatial Data Abstraction Library)是一个在X/MIT许可协议下的开源栅格空间数据转换库。它利用抽象数据模型来表达所支持的各种文件格式。它还有一系列命令行工具来进行数据转换和处理。

GDAL官方网址:http://www.gdal.org/,它能支持当前流行的各种地图数据格式,包括栅格和矢量地图,具体参考官方网站。该库使用C/C++开发,在Java中使用需要自己编译,具体编译过程这里就不说了,下面来看看本文的主要内容吧。

Java使用GDAL读写shapefile的方法示例

读取shp文件,并把它转化为json

import org.gdal.ogr.*;
import org.gdal.ogr.Driver;
import org.gdal.gdal.*;

public class GdalShpTest {
 public static void main(String[] args) {
  // 注册所有的驱动
  ogr.RegisterAll();
  // 为了支持中文路径,请添加下面这句代码
  gdal.SetConfigOption("GDAL_FILENAME_IS_UTF8","YES");
  // 为了使属性表字段支持中文,请添加下面这句
  gdal.SetConfigOption("SHAPE_ENCODING","");

  String strVectorFile = "D:\\test\\NODE.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("打开驱动成功!" );
  dv.CopyDataSource(ds, "D:\\test\\node.json");
  System.out.println("转换成功!" );
 }
}

写shp文件

import org.gdal.ogr.*;
import org.gdal.gdal.*;
class writeShp2 {
 public static void main(String[] args) {
  writeShp2 readshpObj = new writeShp2();
  readshpObj.WriteVectorFile();
 }

 static void WriteVectorFile() {
  String strVectorFile = "D:\\test\\test.shp";

  ogr.RegisterAll();
  gdal.SetConfigOption("GDAL_FILENAME_IS_UTF8", "NO");
  gdal.SetConfigOption("SHAPE_ENCODING", "CP936");

  String strDriverName = "ESRI Shapefile";
  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);

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

  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);

  // 创建五角形要素
  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);

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

得到test.dbf, test.shp, test.shx。

test.dbf如下:

打开形状如下

总结

以上就是这篇文章的全部内容了,希望本文的内容对大家的学习或者工作能带来一定的帮助,如果有疑问大家可以留言交流。

声明:本文内容来源于网络,版权归原作者所有,内容由互联网用户自发贡献自行上传,本网站不拥有所有权,未作人工编辑处理,也不承担相关法律责任。如果您发现有涉嫌版权的内容,欢迎发送邮件至:notice#cainiaojc.com(发邮件时,请将#更换为@)进行举报,并提供相关证据,一经查实,本站将立刻删除涉嫌侵权内容。