GeoTools初体验之读取SHP文件

大名鼎鼎的GeoTools早有耳闻，但是一直没有亲自Lu过。只知道著名的GeoServer是在GeoTools上开发的。嗯，爱屋及乌，对GeoTools也有莫名的好感。前段时间不那么忙了，想看看这个库能不能够实现一个功能：我有几百个栅格文件，需要根据其数值为每个栅格创建SLD文件。文件数目比较多所以考虑使用GIS的库或工具自己来实现批量处理。能够马上想到两个方案，第一使用GDAL读取raster文件，但是不确定GDAL能不能够生成SLD文件。这个可能不是很重要，因为大概率的情况下能够找到另外的库来做。这个方案可以用C++或Python来实现，问题应该不大。第二个方案就是尝尝鲜用没有Lu过的例如GeoTools。

 

在好奇心的驱使下，开始GeoTools的吐血之旅。这个库根据官网文档和前浪的说法，应该是很强大很好用的。虽说也看到有同学的博文说是文档中的代码没有一个能够跑通的，开始我不信，后来。。。。。。觉得自己怎么说也是GIS老程序员了，应该能够轻松搞定。结果花了几炷香的时间连raster文件都没有读出来，心里顿时哇凉哇凉地。

 

今天晚上再次收拾情绪上路，想着先参考官方文档从最简单的做起 - 那就先读取SHP文件吧。花了一个多小时，终于能够让代码运行成功。下面是JAVA代码和mvn文件，供同学们参考。有问题可以到QQ群“开源GIS技术交流群”讨论交流。

import java.io.File;
import java.util.HashMap;
import java.util.Map;
import java.net.URL;

import org.geotools.data.DataStore;
import org.geotools.data.FeatureSource;
import org.geotools.data.shapefile.ShapefileDataStoreFactory;
import org.geotools.data.shapefile.ShapefileDataStore;
import org.geotools.feature.FeatureCollection;
import org.geotools.feature.FeatureIterator;
import org.opengis.feature.simple.SimpleFeature;
import org.opengis.feature.simple.SimpleFeatureType;


public class Main {

    public static void main(String[] args) throws Exception {
        if (args.length >= 1) {
            String shpFilePath = args[0];
            File file = new File(shpFilePath);
            URL url = file.toURI().toURL();

            ShapefileDataStoreFactory factory = new ShapefileDataStoreFactory();
            Map<String, Object> params = new HashMap<>();
            params.put("url", url);
            DataStore dataStore = factory.createDataStore( params );
            ShapefileDataStore shpDataStore = (ShapefileDataStore)dataStore;

            FeatureSource<SimpleFeatureType, SimpleFeature> source = shpDataStore.getFeatureSource();
            
            FeatureCollection<SimpleFeatureType, SimpleFeature> collection = source.getFeatures();

            try (FeatureIterator<SimpleFeature> features = collection.features()) {
                while (features.hasNext()) {
                    SimpleFeature feature = features.next();
                    System.out.print(feature.getID());
                    System.out.print(": ");
                    System.out.println(feature.getDefaultGeometryProperty().getValue());
                }
            }

            System.out.println("=========== SUCCESS ============");
        }
    }
}

<?xml version="1.0" encoding="UTF-8"?>
<project xmlns="http://maven.apache.org/POM/4.0.0"
         xmlns:xsi="http://www.w3.org/2001/XMLSchema-instance"
         xsi:schemaLocation="http://maven.apache.org/POM/4.0.0 http://maven.apache.org/xsd/maven-4.0.0.xsd">
    <parent>
        <artifactId>sldgen</artifactId>
        <groupId>net.brucejia</groupId>
        <version>1.0-SNAPSHOT</version>
    </parent>
    <modelVersion>4.0.0</modelVersion>
    <properties>
        <maven.compiler.source>1.8</maven.compiler.source>
        <maven.compiler.target>1.8</maven.compiler.target>
    </properties>

    <artifactId>shp2txt</artifactId>
    <dependencies>
        <!-- https://mvnrepository.com/artifact/org.geotools/gt-main -->
        <dependency>
            <groupId>org.geotools</groupId>
            <artifactId>gt-main</artifactId>
            <version>${geotools.version}</version>
        </dependency>
        <dependency>
            <groupId>org.geotools</groupId>
            <artifactId>gt-shapefile</artifactId>
            <version>${geotools.version}</version>
        </dependency>
    </dependencies>
</project>

下一次准备Lu一下读取raster文件，敬请期待。