geotools 转换 shpfile 坐标系

简要说明

在java利用geotools工具，将shpfile的坐标系转为指定坐标系。目前测试发现，国外坐标系转换有问题，转换处理的shpfile再去geoserver发布时，经纬度调转了，例如坐标系4547转3857；不过转国内坐标系4326/4490时，却正常。暂时还未定位到哪里问题，也欢迎留言指导一下。

maven依赖

<dependency>
            <groupId>org.geotools</groupId>
            <artifactId>gt-shapefile</artifactId>
            <version>27.0</version>
        </dependency>
        <dependency>
            <groupId>org.geotools.jdbc</groupId>
            <artifactId>gt-jdbc-postgis</artifactId>
            <version>27.0</version>
        </dependency>

        <dependency>
            <groupId>org.geotools</groupId>
            <artifactId>gt-referencing</artifactId>
            <version>27.0</version>
        </dependency>
        <dependency>
            <groupId>org.geotools</groupId>
            <artifactId>gt-render</artifactId>
            <version>27.0</version>
        </dependency>
        <dependency>
            <groupId>org.geotools</groupId>
            <artifactId>gt-coverage</artifactId>
            <version>27.0</version>
        </dependency>
        <dependency>
            <groupId>it.geosolutions.imageio-ext</groupId>
            <artifactId>imageio-ext-utilities</artifactId>
            <version>1.4.4</version>
        </dependency>
        <!-- https://mvnrepository.com/artifact/com.conversantmedia/disruptor -->
        <dependency>
            <groupId>com.conversantmedia</groupId>
            <artifactId>disruptor</artifactId>
            <version>1.2.15</version>
        </dependency>
        <!-- https://mvnrepository.com/artifact/org.geotools/gt-geojson -->
        <dependency>
            <groupId>org.geotools</groupId>
            <artifactId>gt-geojson</artifactId>
            <version>27.0</version>
        </dependency>
        <!-- https://mvnrepository.com/artifact/org.geotools/gt-epsg-hsql -->
        <dependency>
            <groupId>org.geotools</groupId>
            <artifactId>gt-epsg-hsql</artifactId>
            <version>27.0</version>
        </dependency>
        <dependency>
            <groupId>org.geotools</groupId>
            <artifactId>gt-main</artifactId>
            <version>27.0</version>
        </dependency>
        <dependency>
            <groupId>org.geotools</groupId>
            <artifactId>gt-opengis</artifactId>
            <version>27.0</version>
        </dependency>
        <dependency>
            <groupId>si.uom</groupId>
            <artifactId>si-quantity</artifactId>
            <version>2.0.1</version>
        </dependency>
        <dependency>
            <groupId>si.uom</groupId>
            <artifactId>si-units</artifactId>
            <version>2.0.1</version>
        </dependency>
        <dependency>
            <groupId>tech.uom.lib</groupId>
            <artifactId>uom-lib-common</artifactId>
            <version>2.0</version>
        </dependency>
        <dependency>
            <groupId>systems.uom</groupId>
            <artifactId>systems-common</artifactId>
            <version>2.0.2</version>
        </dependency>

样例代码

/**
*dataPath 源数据路径
*srid 指定的坐标系
* outputPath 输出的路径
**/
public void changeShpSrid(String dataPath, Integer srid, String outputPath) {
        ShapefileDataStore inputDataStore = null;
        ShapefileDataStore outputDataStore = null;
        try {
            File inputFile = new File(dataPath);
            //判断文件是否存在
            if(!inputFile.exists()) {
            	System.out.println("源文件不存在");
            	return;
            }
            Map<String, Serializable> inputParams = new HashMap<>();
            inputParams.put("url", inputFile.toURI().toURL());
            inputDataStore = (ShapefileDataStore)DataStoreFinder.getDataStore(inputParams);
            //一般shpfile的编码都是gbk
            inputDataStore.setCharset(Charset.forName("GBK"));

            // 获取输入Shapefile的FeatureSource和Schema
            String typeName = inputDataStore.getTypeNames()[0];
            SimpleFeatureType schema = inputDataStore.getSchema(typeName);

            // 创建输出Shapefile
            File outputFile = new File(outputPath);
            ShapefileDataStoreFactory dataStoreFactory = new ShapefileDataStoreFactory();
            Map<String, Serializable> outputParams = new HashMap<>();
            outputParams.put("url", outputFile.toURI().toURL());
            outputParams.put("create spatial index", Boolean.TRUE);
            outputDataStore = (ShapefileDataStore) dataStoreFactory.createNewDataStore(outputParams);
            outputDataStore.createSchema(schema);
            outputDataStore.setCharset(Charset.forName("GBK"));

            // 设置目标坐标系
            CoordinateReferenceSystem targetCRS = CRS.decode("EPSG:" + srid);
            outputDataStore.forceSchemaCRS(targetCRS);

            // 获取FeatureWriter来写入输出Shapefile
            String outputTypeName = outputDataStore.getTypeNames()[0];
            FeatureWriter<SimpleFeatureType, SimpleFeature> writer = outputDataStore.getFeatureWriterAppend(outputTypeName, Transaction.AUTO_COMMIT);

            CRSAuthorityFactory factory = CRS.getAuthorityFactory(true);
            SimpleFeatureIterator featureIterator = inputDataStore.getFeatureSource(typeName).getFeatures().features();

            Feature firstFeature = featureIterator.next();
            //拿到第一 feature的目的是为了拿到源文件的坐标系
            //第一条数据转坐标系
            Geometry firstGeometry = (Geometry) ((SimpleFeature) firstFeature).getDefaultGeometry();
            SimpleFeatureType sourceSimpleFeatureType = inputDataStore.getSchema();
            CoordinateReferenceSystem coordinateReferenceSystem = sourceSimpleFeatureType.getCoordinateReferenceSystem();
            changeNewSrid(factory, coordinateReferenceSystem, (SimpleFeature)firstFeature, firstGeometry, srid);


            // 读取输入Shapefile中的每个Feature，转换坐标系，并写入输出Shapefile
            while (featureIterator.hasNext()) {
                SimpleFeature feature = featureIterator.next();
                Geometry sourceGeometry = (Geometry) feature.getDefaultGeometry();
                SimpleFeature fNew = writer.next();
                fNew.setAttributes(feature.getAttributes());

                changeNewSrid(factory, coordinateReferenceSystem, fNew, sourceGeometry, srid);
            }

            featureIterator.close();
            // 关闭writer和datastore
            writer.close();
        } catch (IOException e) {
            System.out.println("文件流异常")
        } catch (FactoryException e) {
            System.out.println("获取geotools处理工厂异常");
        } finally {
            inputDataStore.dispose();
            outputDataStore.dispose();
        }
    }

private void changeNewSrid(CRSAuthorityFactory factory, CoordinateReferenceSystem fileSrid, SimpleFeature feature
            , Geometry sourceGeometry , Integer targetSrid) {
        CoordinateReferenceSystem targetCoordinate = null;
        try {
            //目标坐标系
            targetCoordinate = factory.createCoordinateReferenceSystem("EPSG:" + targetSrid);
            //建立转换
            MathTransform transform = CRS.findMathTransform(fileSrid, targetCoordinate, true);
            //转换
            Geometry res = JTS.transform(sourceGeometry, transform);
            feature.setDefaultGeometry(res);
        } catch (FactoryException e) {
            System.out.println("获取转换工厂异常");
        } catch (TransformException e) {
            System.out.println("shpfile转换坐标系失败");
        }
    }