基于geotools实现geotif切瓦片

最新推荐文章于 2024-08-19 10:12:02 发布
最新推荐文章于 2024-08-19 10:12:02 发布
阅读量4.6k 收藏 35
点赞数 7
分类专栏： 地理信息 文章标签： geotools geotif 瓦片 webp
版权声明：本文为博主原创文章，遵循 CC 4.0 BY-SA 版权协议，转载请附上原文出处链接和本声明。
本文链接：https://blog.csdn.net/oYinHeZhiGuang/article/details/121595708
版权

1. 前言

本篇文章主要提供了geotools实现瓦片的切割。程序默认的瓦片大小是256✖256。另外，提供了保存png格式和webp格式的瓦片文件。

2. pom依赖

<?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">
    <modelVersion>4.0.0</modelVersion>

    <groupId>cn.surpass.geotools</groupId>
    <artifactId>geotools-tile</artifactId>
    <version>1.0-SNAPSHOT</version>

    <name>geotools-tile</name>
    <url>http://www.example.com</url>

    <properties>
        <project.build.sourceEncoding>UTF-8</project.build.sourceEncoding>
        <maven.compiler.source>1.8</maven.compiler.source>
        <maven.compiler.target>1.8</maven.compiler.target>
        <geotools.version>19.2</geotools.version>
    </properties>

    <dependencies>

        <!--日志依赖包-->
        <dependency>
            <groupId>org.slf4j</groupId>
            <artifactId>slf4j-log4j12</artifactId>
            <version>1.6.1</version>
        </dependency>
        <dependency>
            <groupId>org.slf4j</groupId>
            <artifactId>slf4j-log4j12</artifactId>
            <version>1.6.1</version>
        </dependency>

        <!--webp依赖包-->
        <dependency>
            <groupId>org.sejda.imageio</groupId>
            <artifactId>webp-imageio</artifactId>
            <version>0.1.6</version>
        </dependency>

        <!--geotools依赖包-->
        <dependency>
            <groupId>org.geotools</groupId>
            <artifactId>gt-geotiff</artifactId>
            <version>${geotools.version}</version>
        </dependency>
        <dependency>
            <groupId>org.geotools</groupId>
            <artifactId>gt-image</artifactId>
            <version>${geotools.version}</version>
        </dependency>
        <dependency>
            <groupId>org.geotools</groupId>
            <artifactId>gt-imagemosaic</artifactId>
            <version>${geotools.version}</version>
        </dependency>
        <dependency>
            <groupId>org.geotools</groupId>
            <artifactId>gt-render</artifactId>
            <version>${geotools.version}</version>
        </dependency>
        <dependency>
            <groupId>org.geotools</groupId>
            <artifactId>gt-imageio-ext-gdal</artifactId>
            <version>${geotools.version}</version>
        </dependency>

        <!--单元测试包-->
        <dependency>
            <groupId>junit</groupId>
            <artifactId>junit</artifactId>
            <version>4.11</version>
            <scope>test</scope>
        </dependency>
    </dependencies>

    <repositories>
        <repository>
            <id>osgeo</id>
            <name>OSGeo Release Repository</name>
            <url>https://repo.osgeo.org/repository/release/</url>
            <snapshots>
                <enabled>false</enabled>
            </snapshots>
            <releases>
                <enabled>false</enabled>
            </releases>
        </repository>
        <repository>
            <id>osgeo-snapshot</id>
            <name>OSGeo Snapshot Repository</name>
            <url>https://repo.osgeo.org/repository/snapshot/</url>
            <snapshots>
                <enabled>false</enabled>
            </snapshots>
            <releases>
                <enabled>false</enabled>
            </releases>
        </repository>
    </repositories>

    <build>
        <pluginManagement>
            <plugins>
                <plugin>
                    <artifactId>maven-clean-plugin</artifactId>
                    <version>3.1.0</version>
                </plugin>
                <plugin>
                    <artifactId>maven-resources-plugin</artifactId>
                    <version>3.0.2</version>
                </plugin>
                <plugin>
                    <artifactId>maven-compiler-plugin</artifactId>
                    <version>3.8.0</version>
                </plugin>
                <plugin>
                    <artifactId>maven-surefire-plugin</artifactId>
                    <version>2.22.1</version>
                </plugin>
                <plugin>
                    <artifactId>maven-jar-plugin</artifactId>
                    <version>3.0.2</version>
                </plugin>
                <plugin>
                    <artifactId>maven-install-plugin</artifactId>
                    <version>2.5.2</version>
                </plugin>
                <plugin>
                    <artifactId>maven-deploy-plugin</artifactId>
                    <version>2.8.2</version>
                </plugin>
                <plugin>
                    <artifactId>maven-site-plugin</artifactId>
                    <version>3.7.1</version>
                </plugin>
                <plugin>
                    <artifactId>maven-project-info-reports-plugin</artifactId>
                    <version>3.0.0</version>
                </plugin>
            </plugins>
        </pluginManagement>
    </build>
</project>
  • 这里需要有几个地方说明一下，引用以下的以来实现webp的输出：
      <dependency>
    <groupId>org.sejda.imageio</groupId>
    <artifactId>webp-imageio</artifactId>
    <version>0.1.6</version>
  </dependency>
  • maven中心库没有geotools的依赖，所以要加以下的仓库配置。
    <repositories>
  <repository>
    <id>osgeo</id>
    <name>OSGeo Release Repository</name>
    <url>https://repo.osgeo.org/repository/release/</url>
    <snapshots>
      <enabled>false</enabled>
    </snapshots>
    <releases>
       <enabled>false</enabled>
    </releases>
  </repository>
  <repository>
    <id>osgeo-snapshot</id>
    <name>OSGeo Snapshot Repository</name>
    <url>https://repo.osgeo.org/repository/snapshot/</url>
    <snapshots>
      <enabled>false</enabled>
    </snapshots>
    <releases>
      <enabled>false</enabled>
    </releases>
  </repository>
</repositories>

3.代码

package cn.surpass.geotools;

import org.geotools.coverage.grid.GridCoverage2D;
import org.geotools.coverage.grid.io.AbstractGridCoverage2DReader;
import org.geotools.coverage.grid.io.AbstractGridFormat;
import org.geotools.coverage.grid.io.GridFormatFinder;
import org.geotools.coverage.processing.CoverageProcessor;
import org.geotools.factory.Hints;
import org.geotools.gce.geotiff.GeoTiffFormat;
import org.geotools.geometry.Envelope2D;
import org.geotools.geometry.jts.ReferencedEnvelope;
import org.opengis.geometry.Envelope;
import org.opengis.parameter.ParameterValueGroup;
import org.opengis.referencing.crs.CoordinateReferenceSystem;
import org.slf4j.Logger;
import org.slf4j.LoggerFactory;

import javax.imageio.ImageIO;
import javax.media.jai.PlanarImage;
import java.awt.image.BufferedImage;
import java.awt.image.RenderedImage;
import java.io.File;
import java.io.FileOutputStream;
import java.io.IOException;

/**
 * @author surpassliang
 */
public class ImageTiles {

    private final static Logger LOGGER = LoggerFactory.getLogger(ImageTiles.class);

    private static final Integer TILE_SIZE = 256;

    public static void handlerImage2Tiles(String geoTifPath, String saveDir, ImageFormat imageFormat) {
        LOGGER.info(String.format("传入参数==>geoTifPath:%s;saveDir:%s;imageFormat:%s", geoTifPath, saveDir, imageFormat.name()));
        //获取GridCoverage2D 对象
        GridCoverage2D gridCoverage = GetGridCoverage2D(geoTifPath);

        //创建目标文件夹
        File saveFolder = new File(saveDir);
        if (saveFolder.exists()) {
            saveFolder.delete();
        }
        saveFolder.mkdirs();
        //影像切片
        handlerTile(gridCoverage, saveFolder, imageFormat);
        LOGGER.info("切片完成");
    }

    private static void handlerTile(GridCoverage2D gridCoverage, File saveFolder, ImageFormat imageFormat) {
        //获取gridCoverage边界
        Envelope2D coverageEnvelope = gridCoverage.getEnvelope2D();
        double minX = coverageEnvelope.getMinX();
        double maxX = coverageEnvelope.getMaxX();
        double minY = coverageEnvelope.getMinY();
        double maxY = coverageEnvelope.getMaxY();
        LOGGER.info(String.format("图像的四至范围:(%f,%f)~(%f,%f)",minX,minY,maxX,maxY));

        //获取像素大小
        RenderedImage sourceImage = gridCoverage.getRenderableImage(0, 1).createDefaultRendering();
        int rasterXSize = sourceImage.getWidth();
        int rasterYSize = sourceImage.getHeight();
        LOGGER.info("图像像素大小:X = {};Y = {}",rasterXSize,rasterYSize);

        //获取geoTif坐标系统
        CoordinateReferenceSystem targetCRS = gridCoverage.getCoordinateReferenceSystem();
        LOGGER.info("图像的坐标系统为:{}",targetCRS.toString());

        //计算瓦片数量
        int htc = (int) Math.ceil(rasterXSize * 1.0 / TILE_SIZE);
        int vtc = (int) Math.ceil(rasterYSize * 1.0 / TILE_SIZE);

        //地理瓦片长宽
        double geoTileWidth = (maxX - minX) / (double) htc;
        double geoTileHeight = (maxY - minY) / (double) vtc;

        for (int hIndex = 0; hIndex < htc; hIndex++) {
            for (int vIndex = 0; vIndex < vtc; vIndex++) {
                LOGGER.info("Processing tile at indices i: {} and j: {}", vIndex, hIndex);

                //创建瓦片envelope
                double startX = (hIndex * geoTileWidth) + minX;
                double startY = (vIndex * geoTileHeight) + minY;
                Envelope envelope = new ReferencedEnvelope(startX, startX + geoTileWidth, startY, startY + geoTileHeight, targetCRS);

                //根据范围切割瓦片
                GridCoverage2D finalCoverage = cropCoverage(gridCoverage, envelope);

                //保存文件
                File tileFile = new File(saveFolder,String.format("%d_%d.%s",vIndex,hIndex,imageFormat.getSuffixName()));

                //保存文件
                BufferedImage tempImage = PlanarImage.wrapRenderedImage(finalCoverage.getRenderedImage()).getAsBufferedImage();
                saveTile(tempImage, tileFile, imageFormat);
            }
        }
    }

    private static void saveTile(BufferedImage image, File tileFile, ImageFormat imageFormat){
        try (FileOutputStream fos = new FileOutputStream(tileFile)) {
            ImageIO.write(image, imageFormat.getFormatName(), fos);
            fos.flush();
        } catch (IOException e) {
            throw new RuntimeException(e.getMessage(),e);
        }
    }

    private static GridCoverage2D cropCoverage(GridCoverage2D gridCoverage, Envelope envelope) {
        CoverageProcessor processor = CoverageProcessor.getInstance();
        ParameterValueGroup coverageCrop = processor.getOperation("CoverageCrop").getParameters();
        coverageCrop.parameter("Source").setValue(gridCoverage);
        coverageCrop.parameter("Envelope").setValue(envelope);
        return (GridCoverage2D) processor.doOperation(coverageCrop);
    }

    private static GridCoverage2D GetGridCoverage2D(String geoTifPath) {
        File geoTifFile = new File(geoTifPath);
        if (!geoTifFile.exists()) {
            throw new RuntimeException("geoTif文件不存在");
        }
        AbstractGridFormat format = GridFormatFinder.findFormat(geoTifFile);

        if (!(format instanceof GeoTiffFormat)) {
            throw new RuntimeException(String.format("传入文件\"%s\"不是geoTif格式文件", geoTifPath));
        }
        Hints hints = new Hints(Hints.FORCE_LONGITUDE_FIRST_AXIS_ORDER, Boolean.TRUE);

        AbstractGridCoverage2DReader geoTifReader = format.getReader(geoTifFile, hints);
        try {
            return geoTifReader.read(null);
        } catch (Exception e) {
            throw new RuntimeException(e.getMessage(), e);
        }
    }

    public enum ImageFormat {
        FORMAT_WEBP("webp", "webp"),
        FORMAT_PNG("png", "png");
        private final String formatName;
        private final String suffixName;
        ImageFormat(String formatName, String suffixName) {
            this.formatName = formatName;
            this.suffixName = suffixName;
        }
        public String getFormatName() {
            return formatName;
        }
        public String getSuffixName() {
            return suffixName;
        }
    }
}

这里需要说明的一下，private static GridCoverage2D cropCoverage(GridCoverage2D gridCoverage, Envelope envelope)方法也可以有如下的实现：

private static GridCoverage2D cropCoverage(GridCoverage2D gridCoverage, Envelope envelope) {
  Operations ops = new Operations(null);
  return (GridCoverage2D) ops.crop(gridCoverage, envelope);
}

4. 方法调用

ImageTiles.handlerImage2Tiles("C:\\Users\\lianghj\\Desktop\\TB.tif","C:\\Users\\Admin\\Desktop\\tile", ImageTiles.ImageFormat.FORMAT_WEBP);
ImageTiles.handlerImage2Tiles("C:\\Users\\lianghj\\Desktop\\TB.tif","C:\\Users\\Admin\\Desktop\\tile", ImageTiles.ImageFormat.FORMAT_PNG);

5. 效果展示

这里主要展示两种不同的格式数据量的大小，这里差距还是挺大的。所以，现在还是很流行的图片格式。
在这里插入图片描述

surpassLiang
关注
专栏目录
GeoTools实战指南: 切割矢量地图数据并导出为Google-XYZ标准的瓦片
孙霸天的专栏
05-11 1700
这个教程主要介绍了如何使用 GeoTools 库将矢量数据（shapefile）切割成不同缩放级别的切片。首先，它讨论了如何读取矢量数据、创建样式和过滤器，然后通过遍历切片来渲染它们。最后，保存切片为符合 Google Maps XYZ 标准的 PNG 格式。
GeoTools实战指南: 切割矢量地图数据并导出为TMS标准的瓦片
孙霸天的专栏
05-10 1046
在本教程中，我们将介绍如何使用Java和GeoTools库来切割矢量地图数据并导出为地图**TMS**标准的瓦片。这对于创建自定义地图服务以及生成本地地图缓存非常有用。你可以了解如何将矢量地图数据切割成适用于Web地图服务的瓦片。
推荐项目：tiles-to-tiff——将XYZ瓦片地图轻松转换为地理参考TIFF图像的神器
最新发布
gitblog_01063的博客
08-19 553
推荐项目：tiles-to-tiff——将XYZ瓦片地图轻松转换为地理参考TIFF图像的神器 tiles-to-tiffPython script for converting XYZ raster tiles for slippy maps to a georeferenced TIFF image 项目地址:https://gitcode.com/gh_mirrors/ti/tiles-to...
GeoTools剪切多边形tif
m0_58742205的博客
08-10 678
【代码】GeoTools剪切多边形tif。
geotools 影像切片 缩放
GIS摆渡人
02-04 3368
图像瓦片应用 ImageLab教程涵盖了加载和渲染;本教程将演示如何使用CoverageProcessor和friends在覆盖率上直接执行基本操作(如裁剪和缩放)，并使用Arguments工具简化命令行处理。我们将创建一个简单的实用程序应用程序来“平铺”一个覆盖范围(为了简单起见，只是简单地细分信封)并可选地缩放结果平铺 预先准备 其中许多可能已经在前面示例中的pom.xml中了，但至少要...
GeoTools加载瓦片并合成Tiff文件，附带经纬度坐标-Java
m0_37956733的博客
03-21 602
在日常工作中，需要根据某个经纬度范围获取当前的瓦片数据，再输出成Tiff文件数据，并追加经纬度信息。
GIS讲堂第二课-地图切片简介以及OL中的调用
weixin_30448603的博客
03-04 130
概述：在大家的支持与帮助下，“GIS讲堂”第二课已落下帷幕，在此对大家的支持与帮助表示衷心的感谢，同时呢，也给大家致歉，由于感冒的缘故，讲课的时候的各种毛病还请大家多多包涵，下面就本节课程的内容做一个公布于说明。课程内容：1、课程主题 地图切片简介以及OL中的调用2、课程主要内容 课程主要内容包括两部分：地图切片简介和OL中切片的调用的与展示。a、地图切片简介 包...
springboot项目-基于geotools的倾斜摄影大地坐标转自定义投影坐标源码.zip
08-31
springboot项目-基于geotools的倾斜摄影大地坐标转自定义投影坐标源码.zipspringboot项目-基于geotools的倾斜摄影大地坐标转自定义投影坐标源码.zipspringboot项目-基于geotools的倾斜摄影大地坐标转自定义投影坐标...
GeoTools实战指南: 首次尝试自定义静态矢量切片
孙霸天的专栏
05-03 672
本教程将介绍如何使用GeoTools库在Java中将一个矢量图层（shapefile 文件）切分成多个图像切片，并保存为 PNG 格式。
java基于geotools实现shapefile几何数据批量读取通用工具类，实现shp数据通用批量读取，shapefile读取
05-28
1.java实现shapefile几何数据批量读取，并封装为通用的工具类，实现shp数据的...3.1.本工具类基于geotools开发，开发使用geotools版本为19.2，几何JTS拓扑套件版本为jts-1.13 4.使用前请自行配置好geotools及jts的jar包
瓦片图切片工具
07-14
openlayers用它不需要另写路径算法
GeotiffParser（tiff格式数据解析器）
08-04
tiff格式数据解析为前端opengl可读取渲染格式的数据
wms源码java-gs-mvt:用于提供Mapnik矢量切片的GeoserverWMS扩展
06-05
wms源码 java 地理服务器 MVT 扩展 概述 Geoserver 的这个扩展增加了作为WMS或Slippy Map Tile请求的结果以协议缓冲区输出格式交付的可能性。 此版本 (0.3.X) 已经过开发和测试，但也可能适用于之前的版本。 从 Geoserver 2.14 及其升级到 Geotools 20.0 开始，由于使用的 JTS 库版本发生变化，此插件将无法工作。 对于 Geoserver 版本 2.14，需要从 jts 1.14 - jts 1.15+ 迁移。 此过程将在单独的分支中完成。 生成的矢量切片可以例如由 WebGL JS 客户端（如 或 ）呈现。 入门 要构建插件，需要 geoserver 源代码。 可在 或 上找到。 您可以使用 maven 构建gs-mvt-0.3.X.jar或直接从 . 为了开始使用gs-mvt-0.3.X.jar包和依赖库protobuf-java-2.6.1.jar必须复制到 geoserver 的 lib 目录geoserver/WEB-INF/lib 。 启动 GeoServer 后，格式application/x-proto
GEOTIFF java api
04-19
GEOTIFF java api
GeoTools实战指南: 理解和使用GeoTools来处理栅格数据创建TMS切片
孙霸天的专栏
05-14 1413
欢迎回到我的编程博客，今天我们要探讨一段Java代码，它使用GeoTools库读取栅格数据文件（例如GeoTIFF格式），并将其渲染成图片切片，这里的切片是TMS格式的。
java使用geotool完成geotiff文件的读写
hwk6155的博客
11-24 2255
geotools版本问题，尽量使用高版本（推荐28.0，不建议使用30.0 因为需要jdk版本大于11）低版本geotools生成的tiff文件会存在以下问题。3.读取tiff文件 并重新写入。2.生成tiff文件。
geotools实现wmts服务
weixin_43031156的博客
07-19 790
【代码】geotools实现wmts服务。
Java geotools实现热力图，生成tiff文件
ljy1998dsb的博客
09-16 2009
基于Java/geotools生成tiff格式的热力图，
geotools 图片切为瓦片
06-06
使用 GeoTools 将图片切为瓦片需要进行以下步骤： 1. 加载图片 使用 GeoTools 的 `ImageIO.read()` 方法加载图片。例如： ``` File file = new File("image.jpg"); BufferedImage image = ImageIO.read(file); ``` 2. 创建瓦片 使用 `TiledImage` 类创建瓦片。例如： ``` TiledImage tiledImage = new TiledImage(0, 0, tileWidth, tileHeight, 0, 0, image.getType()); ``` 其中 `tileWidth` 和 `tileHeight` 表示每个瓦片的宽度和高度。 3. 将图片切分为瓦片 使用 `ImageUtil` 类的 `crop()` 方法将图片切分为瓦片。例如： ``` int rows = image.getHeight() / tileHeight; int cols = image.getWidth() / tileWidth; for (int row = 0; row < rows; row++) { for (int col = 0; col < cols; col++) { Rectangle rect = new Rectangle(col * tileWidth, row * tileHeight, tileWidth, tileHeight); BufferedImage tileImage = ImageUtil.crop(image, rect); tiledImage.setTile(col, row, tileImage); } } ``` 4. 保存瓦片 使用 `ImageIO.write()` 方法保存瓦片。例如： ``` for (int row = 0; row < rows; row++) { for (int col = 0; col < cols; col++) { BufferedImage tileImage = tiledImage.getTile(col, row); File tileFile = new File("tile_" + row + "_" + col + ".jpg"); ImageIO.write(tileImage, "jpg", tileFile); } } ``` 这样就可以将图片切分为瓦片并保存到指定路径了。
评论 8
添加红包

请填写红包祝福语或标题

个

红包个数最小为10个

元

红包金额最低5元

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

打赏作者

surpassLiang

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

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

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

打赏作者

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

抵扣说明：

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

余额充值