GDAL数据批量导出gdb、shapefile等,Java使用GDAL实现通用的数据导出工具类,gdb导出,shp导出,GDAL导出,几何数据导出,shapefile导出,java导出

7 篇文章 1 订阅
5 篇文章 2 订阅

目录

一、实现代码

1.字段注解GdalExportField

2.创建数据源

3.wkt转GDAL几何

4.设置要素核心方法

5.实体对象LandDTO

6.调用示例

二、数据导出效果

1.导出gdb效果

 2.导出shapefile效果

三、使用方法


最近的项目遇到一个需要导出点、线、面、表导出gdb的业务需求,经过思考比较选择了使用GDAL,并且为了便于后续使用,基于注解反射编写了通用的数据导出工具类,只需要将自定义的注解添加在实体字段上便可实现对字段的导出数设置,使用非常便捷。对GIS中常用的数据格式gdb和shapefile进行了验证,二者导出的使用方式相同,并且完美实现字段名、字段类型、字段长度、别名、是否允许空值等属性的控制,下面来看导出效果和工具类代码。

工具类源代码下载:https://download.csdn.net/download/tylkhx/85196838

一、实现代码

1.字段注解GdalExportField

通过注解来指定输出图层的字段信息


@Target(ElementType.FIELD)
@Retention(RetentionPolicy.RUNTIME)
public @interface GdalExportField {
    String name() default "";  
    int type() default ogr.OFTString;
    int length() default 100;
    String alternativeName() default "";
    boolean nullAble() default true; 
    String defaultValue() default "";
    int subType() default 0;
    boolean isGeom() default false;
}

2.创建数据源

创建输出数据的数据源

  public static DataSource createDataSource(String savePath, String driverName){
    if (StringUtil.isEmpty(savePath) || StringUtil.isEmpty(driverName)){
      return null;
    }
    Driver driverOut = ogr.GetDriverByName(driverName);
    return driverOut.CreateDataSource(savePath);
  }

3.wkt转GDAL几何

将java内的几何由wkt格式转换成GDAL几何对象

 public static Geometry createGeometry(String wkt){
    if (StringUtil.isEmpty(wkt)){
      return null;
    }
    return Geometry.CreateFromWkt(wkt);
  }

4.设置要素核心方法

将要素添加到输出图层的核心方法

  public static void setFeature(Layer layer, List<Object> dataList, Class dataClass){
    if (layer == null || dataClass == null){
      return;
    }

    Field[] fields = dataClass.getDeclaredFields();
    if (fields == null || fields.length == 0){
      return;
    }
    for (int i=0;i<fields.length;i++){
      Field fieldClass = fields[i];
      fieldClass.setAccessible(true);
      if (fieldClass.isAnnotationPresent(GdalExportField.class)){
        GdalExportField gdalExportField = fieldClass.getAnnotation(GdalExportField.class);
        if (gdalExportField.isGeom()){
          continue;
        }
        String fieldClassName = fieldClass.getName();
        if (StringUtil.notEmpty(gdalExportField.name())){
          fieldClassName = gdalExportField.name();
        }
        FieldDefn fieldDefn = new FieldDefn();
        fieldDefn.SetName(fieldClassName);
        fieldDefn.SetType(gdalExportField.type());
        fieldDefn.SetWidth(gdalExportField.length());
        String alternativeName = gdalExportField.alternativeName();
        if (StringUtil.isEmpty(alternativeName)){
          alternativeName = fieldClassName;
        }
        fieldDefn.SetAlternativeName(alternativeName);
        if (!gdalExportField.nullAble()){
          fieldDefn.SetNullable(0);
        }
        if (StringUtil.notEmpty(gdalExportField.defaultValue())){
          fieldDefn.SetDefault(gdalExportField.defaultValue());
        }
        fieldDefn.SetSubType(gdalExportField.subType());
        layer.CreateField(fieldDefn);
      }
    }

    if (dataList == null || dataList.isEmpty()){
      return;
    }
    FeatureDefn featureDefn = layer.GetLayerDefn();
    for (Object o:dataList){
      Feature feature = new Feature(featureDefn);
      for (int k=0;k<fields.length;k++){
        Field field = fields[k];
        if (field.isAnnotationPresent(GdalExportField.class)){
          GdalExportField gdalField = field.getAnnotation(GdalExportField.class);
          String fieldName = field.getName();
          if (StringUtil.notEmpty(gdalField.name())){
            fieldName = gdalField.name();
          }
          int filedType = gdalField.type();
          Object fieldValue = null;
          try {
            fieldValue = field.get(o);
          }catch (Exception e){
            e.printStackTrace();
          }
          if (fieldValue == null){
            continue;
          }
          if (gdalField.isGeom()){
            if (layer.GetGeomType() == ogr.wkbNone){ 
              continue;
            }
            Geometry geometry = null;
            if (TYPE_GADL_Geometry.equals(fieldValue.getClass().getName())){
              geometry = (Geometry) fieldValue;
            }else {
              geometry = createGeometry(fieldValue.toString());
            }
            if (geometry == null){
              continue;
            }
            feature.SetGeometry(geometry);
          }else if (filedType == ogr.OFTDate || filedType == ogr.OFTDateTime){
            Date date = (Date) fieldValue;
            Calendar calendar = Calendar.getInstance();
            calendar.setTime(date);
            int year = calendar.get(Calendar.YEAR);
            int month = calendar.get(Calendar.MONTH)+1;
            int day = calendar.get(Calendar.DATE);
            int hour = calendar.get(Calendar.HOUR_OF_DAY);
            int minute = calendar.get(Calendar.MINUTE);
            float second = calendar.get(Calendar.SECOND);
            feature.SetField(fieldName,year,month,day,hour,minute,second,0);
          }else if (filedType == ogr.OFTInteger){
            feature.SetField(fieldName,String.valueOf(fieldValue));
          }else if (filedType == ogr.OFTReal){
            feature.SetField(fieldName,Double.valueOf(String.valueOf(fieldValue)));
          }else if (filedType == ogr.OFTString){
            feature.SetField(fieldName,String.valueOf(fieldValue));
          }else if (filedType == ogr.OFTBinary){
            feature.SetField(fieldName,String.valueOf(fieldValue));
          }
        }
      }
      layer.CreateFeature(feature);
    }
  }

5.实体对象LandDTO

在LandDTO的字段中要添加GdalExportField注解

public class LandDTO {
    @GdalExportField(length = 19)
    private String ZDDM;            // 宗地代码
    @GdalExportField(length = 50,name = "ZDMC",alternativeName = "宗地名称",nullAble = false)
    private String name;            // 宗地名称
    @GdalExportField(length = 50,type = 2)
    private Double ZDMJ;            // 宗地面积
    @GdalExportField(isGeom = true)
    private Geometry geom;

    public String getZDDM() {
        return ZDDM;
    }

    public void setZDDM(String ZDDM) {
        this.ZDDM = ZDDM;
    }

    public String getName() {
        return name;
    }

    public void setName(String name) {
        this.name = name;
    }

    public Double getZDMJ() {
        return ZDMJ;
    }

    public void setZDMJ(Double ZDMJ) {
        this.ZDMJ = ZDMJ;
    }

    public Geometry getGeom() {
        return geom;
    }

    public void setGeom(Geometry geom) {
        this.geom = geom;
    }
}

6.调用示例

  public static void main(String args[]) throws IOException {
    List<LandDTO> landDTOList = new ArrayList<>();
    LandDTO landDTO = new LandDTO();
    landDTO.setZDDM("542145210001");
    landDTO.setName("张三");
    landDTO.setZDMJ(135.24);
    landDTO.setGeom(createGeometry("Polygon((2 2,6 2,8 8,2 6,2 2))"));
    landDTOList.add(landDTO);

    LandDTO landDTO2 = new LandDTO();
    landDTO2.setZDDM("542145210002");
    landDTO2.setName("李四");
    landDTO2.setZDMJ(232.94);
    landDTO2.setGeom(createGeometry("Polygon((12 12,16 12,18 18,12 16,12 12))"));
    landDTOList.add(landDTO2);

    SpatialReference sr = new SpatialReference();
    sr.ImportFromEPSG(4326);
    String saveGdbPath = "C:/Users/dl/Desktop/新建文件夹/测试输出gdb.gdb";
    DataSource dataSourceGDB = GdalUtil.createDataSource(saveGdbPath, "FileGDB");
    Layer gdbLayer = dataSourceGDB.CreateLayer("ZDXX", sr, ogr.wkbPolygon, null);
    GdalUtil.setFeature(gdbLayer,(List) landDTOList,LandDTO.class);

    Layer gdbTable = dataSourceGDB.CreateLayer("ZDB", sr, ogr.wkbNone, null);
    GdalUtil.setFeature(gdbTable,(List) landDTOList,LandDTO.class);
    dataSourceGDB.SyncToDisk();
    dataSourceGDB.FlushCache();
    dataSourceGDB.delete();
    String saveShpPath = "C:/Users/dl/Desktop/新建文件夹/测试输出shp.shp";
    DataSource dataSourceShp = GdalUtil.createDataSource(saveShpPath, "ESRI Shapefile");
    Layer ShpLayer = dataSourceShp.CreateLayer("ZDXX", sr, ogr.wkbPolygon, null);
    GdalUtil.setFeature(ShpLayer,(List) landDTOList,LandDTO.class);
    dataSourceShp.SyncToDisk();
    dataSourceShp.FlushCache();
    dataSourceShp.delete();
  }

二、数据导出效果

1.导出gdb效果

 2.导出shapefile效果

三、使用方法

1.在实体对象要导出的字段上添加GdalExportField注解,并在注解中指定字段名(如果不知道,则默认使用实体对象的字段名)、字段类型(默认字符串)、字段长度、别名、是否允许空值等信息。关注微行公众号:“GIS工具乐园”,持续更新开发技巧分享哦。

2.调用工具类的createDataSource()方法创建数据源(DataSource),创建数据源时指定要导出的数据格式:“FileGDB”为GDB,“ESRI Shapefile”为shapefile。

3.调用GDAL DataSource的自带方法CreateLayer()创建图层,在创建图层时指定图层名、图层类型等。

4.调用工具类的setFeature()方法将要导出的数据添加的图层中,该方法会创建要素,并反射获取GdalExportField注解的内容,设置到要素字段,并添加要素到图层中。

5.依次调用GDAL DataSource的自带方法SyncToDisk()、FlushCache()、 delete()方法将数据写出到磁盘

6.以上操作在工具类的main方法有示例

  • 5
    点赞
  • 7
    收藏
    觉得还不错? 一键收藏
  • 打赏
    打赏
  • 8
    评论
您可以使用Java GDAL库来读取Shapefile数据并获取最大经纬度和最小经纬度。下面是一个示例代码片段,展示了如何实现这一功能: ```java import org.gdal.gdal.Dataset; import org.gdal.ogr.*; public class ShapefileReader { public static void main(String[] args) { // 指定Shapefile文件路径 String shapefilePath = "path/to/your/shapefile.shp"; // 注册所有的GDAL驱动 gdal.AllRegister(); // 打开Shapefile数据集 Dataset dataset = ogr.Open(shapefilePath, 0); if (dataset == null) { System.out.println("无法打开Shapefile文件"); System.exit(1); } // 获取第一个图层 Layer layer = dataset.GetLayer(0); // 获取图层的空间参考 SpatialReference srs = layer.GetSpatialRef(); // 获取图层的要素数量 int featureCount = layer.GetFeatureCount(); // 定义最大经纬度和最小经纬度的初始值 double minX = Double.MAX_VALUE; double minY = Double.MAX_VALUE; double maxX = -Double.MAX_VALUE; double maxY = -Double.MAX_VALUE; // 遍历所有要素获取最大经纬度和最小经纬度 for (int i = 0; i < featureCount; i++) { Feature feature = layer.GetFeature(i); Geometry geometry = feature.GetGeometryRef(); // 获取要素的最小外包矩形 Envelope envelope = new Envelope(); geometry.GetEnvelope(envelope); // 更新最大经纬度和最小经纬度 minX = Math.min(minX, envelope.MinX); minY = Math.min(minY, envelope.MinY); maxX = Math.max(maxX, envelope.MaxX); maxY = Math.max(maxY, envelope.MaxY); // 释放要素和几何体 feature.delete(); geometry.delete(); } // 打印最大经纬度和最小经纬度 System.out.println("最小经度: " + minX); System.out.println("最小纬度: " + minY); System.out.println("最大经度: " + maxX); System.out.println("最大纬度: " + maxY); // 释放图层和数据集 layer.delete(); dataset.delete(); } } ``` 请确保您已经正确安装了Java GDAL库,并将您的Shapefile文件路径替换为实际路径。运行此代码片段后,您将获得最大经纬度和最小经纬度的输出结果。

“相关推荐”对你有帮助么?

  • 非常没帮助
  • 没帮助
  • 一般
  • 有帮助
  • 非常有帮助
提交
评论 8
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包

打赏作者

GIS工具开发

你的鼓励将是我创作的最大动力

¥1 ¥2 ¥4 ¥6 ¥10 ¥20
扫码支付:¥1
获取中
扫码支付

您的余额不足,请更换扫码支付或充值

打赏作者

实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值