Java-GDAL读取ArcGIS FileGDB

之前的博客中我们已经介绍了关于ArcGIS FileGDB数据是什么,在这里我们就暂时不做重复了

1. 安装GDAL,在官网上下载安装文件(记得与电脑版本、JDK版本保持一致)

https://gdal.org/index.html

链接:https://pan.baidu.com/s/1WFbx9ap3ixukysanWRujJg
提取码:zaqz
–来自百度网盘超级会员V4的分享

2. Java 使用GDAL 读取FileGDB

2.1 Maven配置

<dependencys>
    <dependency>
        <groupId>org.gdal</groupId>
        <artifactId>gdal</artifactId>
        <version>1.11.2</version>
    </dependency>
</dependencys>

具体可参考:https://blog.csdn.net/weixin_42504649/article/details/114217332

2.2 代码片段

2.2.1 读取GDB


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

import java.sql.Date;
import java.sql.Time;
import java.sql.Timestamp;
import java.time.LocalDate;
import java.time.LocalDateTime;
import java.time.LocalTime;
/**
 * @Classname GDALHelper
 * @Description GDAL读取GDB
 * @Date 2021/4/19 16:27
 * @Created by xiaocai
 */
public class GDALHelper {

    static {
        gdal.AllRegister();
    }

    public static void main(String[] args) {
        Driver driver = ogr.GetDriverByName("OpenFileGDB");
        DataSource dataSource = driver.Open("C:\\Users\\xiaocai\\Desktop\\XYBASE.gdb", 0);
        for (int i = 0; i < dataSource.GetLayerCount(); i++) {
            //获取图层
            Layer layer = dataSource.GetLayer(i);
            System.out.println(i + "\t" + layer.GetName());
            do {//获取图层下的要素
                Feature feature = layer.GetNextFeature();
                if (null == feature) {
                    break;
                }
                //获取样式,这里是空,待确定
                System.out.println(feature.GetStyleString());
                //获取geometry
                System.out.println(feature.GetGeometryRef().ExportToWkt());//也可以ExportToWkb() gml等等
                //获取属性
                for (int p = 0; p < feature.GetFieldCount(); p++) {
                    System.out.println(feature.GetFieldDefnRef(p).GetName()
                            + "\t" +
                            getProperty(feature, p));
                }
                System.out.println("---------");
            } while (true);
        }
    }

    private static Object getProperty(Feature feature, int index) {
        int type = feature.GetFieldType(index);
        PropertyGetter propertyGetter;
        if (type < 0 || type >= propertyGetters.length) {
            propertyGetter = stringPropertyGetter;
        } else {
            propertyGetter = propertyGetters[type];
        }
        try {
            return propertyGetter.get(feature, index);
        } catch (Exception e) {
            e.printStackTrace();
            return null;
        }
    }

    /**
     * 属性获取器
     */
    @FunctionalInterface
    private interface PropertyGetter {
        Object get(Feature feature, int index);
    }

    private static final PropertyGetter stringPropertyGetter = (feature, index) -> feature.GetFieldAsString(index);

    /**
     * feature.GetFieldType(index)得到一个属性类型的int值,该值对应具体类型
     */
    private static final PropertyGetter[] propertyGetters = new PropertyGetter[]{
            (feature, index) -> feature.GetFieldAsInteger(index),//0	Integer
            (feature, index) -> feature.GetFieldAsIntegerList(index),//1	IntegerList
            (feature, index) -> feature.GetFieldAsDouble(index),//2	Real
            (feature, index) -> feature.GetFieldAsDoubleList(index),//3	RealList
            stringPropertyGetter,//4	String
            (feature, index) -> feature.GetFieldAsStringList(index),//5	StringList
            stringPropertyGetter,//6	(unknown)
            stringPropertyGetter,//7	(unknown)
            (feature, index) -> feature.GetFieldAsBinary(index),//8	Binary
            (feature, index) -> {
                int[] pnYear = new int[1];
                int[] pnMonth = new int[1];
                int[] pnDay = new int[1];
                int[] pnHour = new int[1];
                int[] pnMinute = new int[1];
                float[] pfSecond = new float[1];
                int[] pnTZFlag = new int[1];
                feature.GetFieldAsDateTime(index, pnYear, pnMonth, pnDay, pnHour, pnMinute, pfSecond, pnTZFlag);
                Date date = Date.valueOf(LocalDate.of(pnYear[0], pnMonth[0], pnDay[0]));
                return date;
            },//9	Date
            (feature, index) -> {
                int[] pnYear = new int[1];
                int[] pnMonth = new int[1];
                int[] pnDay = new int[1];
                int[] pnHour = new int[1];
                int[] pnMinute = new int[1];
                float[] pfSecond = new float[1];
                int[] pnTZFlag = new int[1];
                feature.GetFieldAsDateTime(index, pnYear, pnMonth, pnDay, pnHour, pnMinute, pfSecond, pnTZFlag);
                float fSecond = pfSecond[0];
                int s = (int) fSecond;
                int ns = (int) (1000000000 * fSecond - s);
                Time time = Time.valueOf(LocalTime.of(pnHour[0], pnMinute[0], s, ns));
                return time;
            },// 10	Time
            (feature, index) -> {
                int[] pnYear = new int[1];
                int[] pnMonth = new int[1];
                int[] pnDay = new int[1];
                int[] pnHour = new int[1];
                int[] pnMinute = new int[1];
                float[] pfSecond = new float[1];
                int[] pnTZFlag = new int[1];
                feature.GetFieldAsDateTime(index, pnYear, pnMonth, pnDay, pnHour, pnMinute, pfSecond, pnTZFlag);
                float fSecond = pfSecond[0];
                int s = (int) fSecond;
                int ns = (int) (1000000000 * fSecond - s);
                LocalDateTime localDateTime = LocalDateTime.of(
                        LocalDate.of(pnYear[0], pnMonth[0], pnDay[0]),
                        LocalTime.of(pnHour[0], pnMinute[0], s, ns)
                );
                Timestamp timestamp = Timestamp.valueOf(localDateTime);
                return timestamp;
            },//11	DateTime
            (feature, index) -> feature.GetFieldAsInteger64(index),//12	Integer64
            (feature, index) -> feature.GetFieldAsIntegerList(index),//13 Integer64List
            //>=14	(unknown)
    };
}

2.2.2 在PostGIS中创建要素类

    System.out.println(String.format("%s   %s   ", layer.GetName(), layer.GetGeomType()));
    if (layer.GetGeomType() == ogr.wkbPoint||layer.GetGeomType()==ogr.wkbPointM||layerGetGeomType()==ogr.wkbPointZM) {
        pCreateLayer = pgDataSource.CreateLayer(layer.GetName(), layer.GetSpatialRef(), ogr.wkbPoint, null);
    } else if (layer.GetGeomType() ==ogr.wkbLineString||layer.GetGeomType() ==ogrwkbMultiLineString||layer.GetGeomType() ==ogr.wkbMultiLineStringM||layer.GetGeomTyp() ==ogr.wkbMultiLineStringZM) {
        pCreateLayer = pgDataSource.CreateLayer(layer.GetName(), layer.GetSpatialRef(), ogr.wkbMultiLineString, null);
    } else {
        pCreateLayer = pgDataSource.CreateLayer(layer.GetName(), layer.GetSpatialRef(), ogr.wkbMultiPolygon, null);
    }
    FeatureDefn pFeatureDefn = layer.GetLayerDefn();
    int iFieldCount = pFeatureDefn.GetFieldCount();
    for (int iAttr = 0; iAttr < iFieldCount; iAttr++) {
        FieldDefn oField = pFeatureDefn.GetFieldDefn(iAttr);

        if (!oField.GetName().toUpperCase().equals("SHAPE_LENGTH") && !oField.GetName().toUpperCase().equals("SHAPE_AREA")) {
            pCreateLayer.CreateField(oField);
        }
    }

    pgDataSource.FlushCache();

pgDataSource.FlushCache();比较关键,需要同步缓存到数据库,否则在变量释放前是找不到数据的

2.2.3. 插入要素类到PostGIS中

public void inputGDBFeaClass2PostGIS(Layer layer, Layer pPGLayer, Connection connection) {

        PreparedStatement preparedStatement = null;
        try{
            int index = 0;
            Map<String, Object> mapValues = new LinkedHashMap<>();
            String strInsertSql = getPgClassInsertSql(layer.GetName(), layer, mapValues,pPGLayer);
            preparedStatement = connection.prepareStatement(strInsertSql);

            Feature oFeature = null;
//                FeatureDefn pCreateFeatureDefn =pPGLayer.GetLayerDefn();
            // 下面开始遍历图层中的要素
            while ((oFeature = layer.GetNextFeature()) != null) {
                //获取属性
                for (int p = 0; p < oFeature.GetFieldCount(); p++) {
                    if(mapValues.containsKey(oFeature.GetFieldDefnRef(p).GetName().toLowerCase()))
                    {
                        mapValues.put(oFeature.GetFieldDefnRef(p).GetName().toLowerCase(),getProperty(oFeature, p));
                    }
                }

                // 带高程信息转为二维平面数据
                Geometry pGeometry=oFeature.GetGeometryRef();
                pGeometry.FlattenTo2D();
                mapValues.put("wkb_geometry",pGeometry.ExportToWkt());


                Integer MapIndex = 1;
                for (String key : mapValues.keySet()) {
                    if(key.toString().equals("wkb_geometry"))
                    {
                        preparedStatement.setString(MapIndex++, String.format("SRID=4528;%s",mapValues.get(key).toString()));
                    }
                    else
                    {
                        preparedStatement.setObject(MapIndex++, mapValues.get(key));
                    }
                }

                preparedStatement.addBatch();
                index++;
                if (index % 10000 == 0) {
                    System.out.println(index);
                    preparedStatement.executeBatch();
                    preparedStatement.clearBatch();
                }
            }
            preparedStatement.executeBatch();
        }
        catch (Exception ex)
        {
            ex.printStackTrace();
        }

    }

private String getPgClassInsertSql(String strLayerName, Layer pSouceLayer, Map<String, Object> mapValues, Layer pPGTargetLayer) {
        List<String> listValues = new ArrayList<>();


        FeatureDefn pSourceFeatureDefn = pSouceLayer.GetLayerDefn();

        FeatureDefn pPGTargetFeatureDefn = pPGTargetLayer.GetLayerDefn();

        for (int i = 0; i < pSourceFeatureDefn.GetFieldCount(); i++) {
            FieldDefn pSourceFieldDefn = pSourceFeatureDefn.GetFieldDefn(i);
            if (pPGTargetFeatureDefn.GetFieldIndex(pSourceFieldDefn.GetName()) >= 0) {
                mapValues.put(pSourceFieldDefn.GetName().toLowerCase().toString(), null);
                listValues.add("?");
            }
        }
        mapValues.put("wkb_geometry", null);
        listValues.add("st_geomfromText(?)");

        String strInsertSql = String.format("insert into %s(%s) values(%s)", strLayerName, String.join(",", mapValues.keySet()), String.join(",", listValues));
        return strInsertSql;
    }
  • 3
    点赞
  • 6
    收藏
    觉得还不错? 一键收藏
  • 打赏
    打赏
  • 4
    评论
GDAL(Geospatial Data Abstraction Library)是一个开源的地理空间数据处理库,可以用于读取、写入和操作各种不同格式的地理空间数据文件。而GDB(Geodatabase)是ESRI公司的一种地理数据库格式,GDB文件可以存储各种地理空间数据,如矢量数据、栅格数据、拓扑数据等。 在C语言中使用GDAL读取GDB文件时,可以按照以下步骤进行: 1. 首先,需要在代码中包含GDAL的头文件。例如: ```c #include "gdal.h" ``` 2. 创建GDAL数据集对象,用于打开GDB文件并获取其中的数据。例如: ```c GDALDatasetH hDataset; // GDAL数据集对象指针 hDataset = GDALOpenEx("path_to_gdb_file", GDAL_OF_VECTOR, NULL, NULL, NULL); ``` 3. 获取数据集中的图层数量和每个图层的名称。例如: ```c int layerCount = GDALDatasetGetLayerCount(hDataset); // 获取图层数量 for (int i = 0; i < layerCount; i++) { OGRLayerH hLayer = GDALDatasetGetLayer(hDataset, i); // 获取第i个图层 const char* layerName = OGR_FD_GetName(GDAL_L_GetLayerDefn(hLayer)); // 获取图层名称 // 打印图层名称 printf("Layer %d: %s\n", i, layerName); } ``` 4. 读取图层中的要素数据。例如: ```c OGRLayerH hLayer = GDALDatasetGetLayer(hDataset, layerIndex); // 获取第layerIndex个图层 OGRFeatureH hFeature; OGR_L_ResetReading(hLayer); // 重置读取位置 while((hFeature = OGR_L_GetNextFeature(hLayer)) != NULL) { // 获取要素的属性值 OGRFeatureDefnH hFeatureDefn = OGR_L_GetLayerDefn(hLayer); int fieldCount = OGR_FD_GetFieldCount(hFeatureDefn); // 获取属性字段数量 for (int j = 0; j < fieldCount; j++) { OGRFieldDefnH hFieldDefn = OGR_FD_GetFieldDefn(hFeatureDefn, j); const char* fieldName = OGR_Fld_GetNameRef(hFieldDefn); // 获取字段名称 int fieldValue = OGR_F_GetFieldAsInteger(hFeature, j); // 获取字段值(整数类型) // 打印字段名称和值 printf("Field %s: %d\n", fieldName, fieldValue); } OGR_F_Destroy(hFeature); // 释放要素对象 } ``` 5. 最后,记得关闭GDAL数据集并释放资源。例如: ```c GDALClose(hDataset); // 关闭数据集 ``` 以上就是利用GDAL库在C语言中读取GDB文件的基本步骤。可根据实际需要进行进一步的数据处理和操作。
评论 4
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

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

打赏作者

自己的九又四分之三站台

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

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

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

打赏作者

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

抵扣说明:

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

余额充值