GDAL几何修复，java几何拓扑修复、修复几何图形自相交工具类，支持shp、gdb几何数据格式，shp拓扑修复，gdb拓扑修复

package util; import org.gdal.gdal.gdal; import org.gdal.ogr.*; import java.io.IOException; import java.math.BigDecimal; /** * gdal几何修复工具类 * */ public class GdalMakeValidUtil { static { gdal.AllRegister(); gdal.SetConfigOption("GDAL_FILENAME_IS_UTF8","YES"); gdal.SetConfigOption("SHAPE_ENCODING","CP936"); } public static void main(String args[]) throws IOException { //打开数据 String shpPath = "D:\\新建文件夹\\拓扑错误示例.shp"; DataSource dataSource = openSource(shpPath, true); if (dataSource == null){ return; } //修复 makeValid(dataSource.GetLayer(0)); //数据存盘 dataSource.FlushCache(); dataSource.delete(); } /** * 几何修复 * @param layer 图层 * 对原图层进行修复，并增加“修复情况”字段 */ public static void makeValid(Layer layer){ if (layer == null){ return; } FeatureDefn featureDefn = layer.GetLayerDefn(); //图层增加修复情况字段 String fieldName = "修复情况"; int fieldNum = 0; while (featureDefn.GetFieldIndex(fieldName) != -1 || fieldNum > 9){ fieldName = "修复情况" + fieldNum; fieldNum += 1; } FieldDefn fieldDefn = new FieldDefn(); fieldDefn.SetName(fieldName); fieldDefn.SetType(ogr.OFTString); fieldDefn.SetWidth(200); layer.CreateField(fieldDefn); //修复 layer.ResetReading(); long featureCount = layer.GetFeatureCount(); for (long i=0;i<featureCount;i++){ Feature feature = layer.GetNextFeature(); Geometry geometry = feature.GetGeometryRef(); if (!geometry.IsValid()){ double oldArea = MathCalculateUtils.round(geometry.GetArea(),4); Geometry geometryNew = geometry.MakeValid(); if (!geometryNew.IsValid()){ feature.SetField(fieldName,"几何修复失败"); }else { double newArea = MathCalculateUtils.round(geometryNew.GetArea(),4); BigDecimal areaSubtract = new BigDecimal(String.valueOf(oldArea)).subtract(new BigDecimal(String.valueOf(newArea))); if (areaSubtract.compareTo(new BigDecimal("0.0")) != 0){ feature.SetField(fieldName,"几何修复完成，修复前后面积差" + areaSubtract.toString() + "平方米，注意检查几何图形是否正确"); }else { feature.SetField(fieldName,"几何修复完成"); } feature.SetGeometry(geometryNew); } }else { feature.SetField(fieldName,"几何正常，无需修复"); } layer.SetFeature(feature); } } /** * 打开数据源 * @param path 数据路径 * @param update 是否更新源数据 * @return */ public static DataSource openSource(String path,boolean update){ if (StringUtil.isEmpty(path)){ return null; } if (update){ return ogr.Open(path,1); }else { return ogr.Open(path,0); } } }