一种遥感影像瓦片生成及实时部署的方法技术

本发明专利技术涉及一种遥感影像瓦片生成及实时部署的方法，包括以下步骤：解析命令行，获取瓦片配置参数；创建日志文件，检查配置参数；获取所运行计算机硬件信息，创建线程池；确定瓦片切片标准，获取影像参数，计算影像的最大与最小层级；根据所运行计算机的硬件信息，对遥感影像动态分块，遍历每一块影像，将分块影像导入内存，计算分块影像的最小与最大层级以及仿射变换参数；遍历每一层级，进行瓦片索引编码；线程池并行切片；每一层级瓦片处理完之后，当前块切片结束；每一块影像处理完之后，切片结束。本发明专利技术切片速度快，效率高，不仅能通过网络发送数据保存到服务器，也可保存到本地。

A method of tile generation and real-time deployment of remote sensing image

【技术实现步骤摘要】
一种遥感影像瓦片生成及实时部署的方法
本专利技术属于影像切片
，特别涉及一种遥感影像瓦片生成及实时部署的方法。
技术介绍
目前，切片是网络地理信息系统中重要的功能。但是随着影像数据量的增长，影像尺寸的增大，网络地理信息系统的实时性面临挑战，影像的切片速度严重影响网络地理信息系统的使用。同时，在进行影像切片时，边界的瓦片只存在部分区域，然而目前针对边界瓦片的处理方法比较少。目前如何实现影像数据的高效切片，如何实现超大尺寸影像切片，如何实现影像边界部分区域瓦片处理，成为亟待解决的问题。目前常用的切片方法主要有两种，一种是串行切片算法，此方法速度慢，无法实时处理GB级以上的数据；一种是集群切片算法，此方法成本高，部署复杂。针对超大尺寸影像，中国专利文献公开了一种超大尺寸影像切片方法，此方法虽然提高了切片效率，但是实际运用中程序复杂且容易出错。
技术实现思路
为了解决现有技术硬件成本高、切片效率低、部署复杂以及边界问题，本专利技术提供一种高效且易于实现和部署的遥感影像瓦片生成及实时部署的方法。本专利技术提供一种高效且易于实现和部署的遥感影像的遥感影像瓦片生成及实时部署的方法，包括以下步骤：步骤1，命令行解析，判断配置参数解析方法，获取瓦片配置参数。解析方法包括静态模式和动态模式，静态模式读取指定路径的XML文件，动态模式解析调用程序时的命令行参数；步骤2，创建日志文件，检查各种配置参数，如果某一配置参数为空，则程序退出。包括影像路径、坐标系类型、瓦片规范、瓦片格式等；步骤3，获取所运行计算机硬件信息，所运行计算机硬件信息包括处理器的核数、线程个数、当前可用内存大小读写能力等。创建线程池，线程的个数为所运行计算机的核数，线程池包括：任务队列、空闲线程队列、忙碌线程队列、线程管理器；步骤4，确定瓦片标准，即瓦片规范与坐标系类型相互组合：wmts/webmercator、wmts/wgs84、tms/webmercator、tms/wgs84。获取遥感影像参数，包括：影像大小、波段数、位深以及仿射变化参数。根据全球经度范围和影像的经度范围计算最小层级，当全球经度小于影像经度范围时为最小层级，根据全球经度范围和单个瓦片经度范围计算最大层级，当全球经度范围小于单个瓦片的经度范围时为最大层级；步骤5，根据所运行计算机硬件信息，对遥感影像进行动态分块，影像的动态分块是根据所运行计算机当前内存大小的一半进行分块。遍历每一块影像，将分块导入内存，计算分块影像的最小与最大层级以及仿射变换参数；步骤6，起始层级为原始影像的最小层级，遍历每一层级，判断当前层级是否大于最大层级，如果是，则切片任务结束，如果否，则计算当前层级瓦片索引编码，并将切图任务添加到线程池进行并行切片；其中，每一层级的影像瓦片，按照从左到右，从上到下的顺序，依据级别、行号、列号、位置进行索引编码；步骤7，每一层级瓦片处理完之后，当前块切片结束；步骤8，每一块影像处理完之后，切片结束；进一步的，步骤6中基于CPU的并行影像切片，其中单个线程的流程包括：步骤601：空间经纬度坐标转换为物理像素坐标，创建空瓦片；步骤602：判断当前层级是否小于等于分块影像的最小层级，小于等于分块影像的最小层级，执行步骤603，否则执行步骤604；步骤603：瓦片层级缺失处理，即根据瓦片编号和影像大小、影像仿射变化参数计算瓦片的起点和终点经纬度，然后根据经纬度计算像素的起点和终点坐标。根据像素的起点坐标，判断瓦片的位置，计算影像的起点像素坐标和偏移量，瓦片的起点坐标和偏移量，对待写入区域降采样处理，将数据写入瓦片对应区域，共存在15种可能存在的情况，然后执行步骤605；步骤604：根据瓦片编号和影像大小、影像仿射变化参数计算瓦片的起点和终点经纬度，然后根据经纬度计算像素的起点和终点坐标。根据像素的起点坐标，判断瓦片的位置，计算影像的起点像素坐标和偏移量，瓦片的起点坐标和偏移量，对待写入区域降采样处理，将数据写入瓦片对应区域，共存在9种情况；步骤605：判断当前瓦片每个点的像素值，如果为0则是无值瓦片，如果是无值瓦片，执行步骤601，否则执行步骤606；步骤606：判断当前瓦片是否处于影像四周边界，如果是，则是边界瓦片，如果是边界瓦片则执行步骤607，否则执行步骤608；步骤607：通过网络从服务器获取相同编码的瓦片进行融合，遍历待融合瓦片，如果某一点的像素为0，则从服务器瓦片获取对应位置的像素进行填充；步骤608：添加瓦片格式，根据配置参数URL地址，通过网络发送到目的地；步骤609：根据配置参数，判断当前瓦片是否需要保存到本地，如果需要则执行步骤610，否则执行步骤611；步骤610：根据配置参数，将当前瓦片保存到执行文件夹；步骤611：当前线程结束；本专利技术提供一种遥感影像瓦片生成及实时部署的方法有益效果包括：a)遥感影像瓦片生成及实时部署的方法，具有切片速度快、硬件成本低、部署容易以及支持数量大、尺寸大的影像数据等特点；b)对影像数据动态分块，是按照所运行计算机可用内存大小的一半进行分块，尽可能多的导入内存，又能够保证计算应用程序的运行；c)对影像瓦片进行无值检测，可以解决具有黑边的遥感影像；d)判断瓦片是否是包含部分区域的瓦片，如果是，则对边界包含部分区域的瓦片进行融合处理。e)软件不需要环境配置，程序轻量，部署快速容易；附图说明图1为本专利技术的总体流程示意图。图2为瓦片的索引编码图例。图3为线程池原理示意图。图4为降采样原理示意图。图5为正常瓦片分块示意图。图6为层级缺失分块示意图。具体实施方式为了更好的说明本专利技术，下面将结合附图和实施例对本专利技术做进一步详细描述。如图1所示一种遥感影像瓦片生成及实时部署的方法，采用金字塔模型和影像处理算法对影像数据进行切片，包括以下步骤：步骤1，命令行解析，判断配置参数解析方法，获取瓦片配置参数。解析方法包括静态模式和动态模式，静态模式读取指定路径的XML文件，动态模式解析调用程序时的命令行参数；步骤2，创建日志文件，检查各种配置参数，如果某一配置参数为空程序退出。包括影像路径、坐标系类型、瓦片规范、瓦片格式等；步骤3，获取所运行计算机硬件信息，所运行计算机硬件信息包括处理器的核数、线程个数、当前可用内存大小读写能力等。创建线程池，线程的个数为所运行计算机的核数，线程池包括：任务队列、空闲线程队列、忙碌线程队列，线程管理器；步骤4，确定瓦片标准，是瓦片规范与坐标系类型相互组合：wmts/webmercator、wmts/wgs84、tms/webmercator、tms/wgs84。获取遥感影像参数，包括：影像大小、波段数、位深以及仿射变化参数。根据全球经度范围和影像的经度范围计算最小层级，当全球经度小于影像经度范围时为最大层级，根据全球经度范围和单个瓦片经度范围计算最大层级，当全球经度范围小于单个瓦片的经度范围时为最大层级；步骤5，根据所运行计算机硬件信息，对遥感影像进行动态分块，影像的动态分块是根据所运行计算机当前内存大小的一半进行分块。遍历每一块影像，将分块导入内存，计算分块影像的最小与最大层级以及仿射变换参数；步骤6，起始层级为原始影像的最小层级，遍历每一层级，判断当前层级是否大于最大层级，如果是，则切片任务结束，如果否，则计算当前层本文档来自技高网...

【技术保护点】
1.一种遥感影像瓦片生成及实时部署的方法，其特征在于，包括如下步骤：步骤1，命令行解析，判断配置参数解析方法，获取瓦片配置参数；步骤2，创建日志文件，检查配置参数，如果某一配置参数为空，则退出；步骤3，获取所运行计算机硬件信息，并创建线程池；步骤4，确定瓦片标准，获取遥感影像参数，根据全球经度范围和影像的经度范围计算最小层级，当全球经度小于影像经度范围时为最小层级，根据全球经度范围和单个瓦片经度范围计算最大层级，当全球经度范围小于单个瓦片的经度范围时为最大层级；步骤5，根据所运行计算机硬件信息，对遥感影像进行动态分块，遍历每一分块影像，将分块影像导入内存，计算分块影像的最小与最大层级以及仿射变换参数；步骤6，起始层级为原始影像的最小层级，遍历每一层级，判断当前层级是否大于最大层级，如果是，则切片任务结束，如果否，则计算当前层级瓦片索引编码，并将切图任务添加到线程池进行并行切片；其中每一层级的影像瓦片，按照从左到右，从上到下的顺序，依据级别、行号、列号、位置进行索引编码；步骤7，每一层级瓦片处理完之后，当前块切片结束；步骤8，每一块影像处理完之后，切片结束。

【技术特征摘要】
1.一种遥感影像瓦片生成及实时部署的方法，其特征在于，包括如下步骤：步骤1，命令行解析，判断配置参数解析方法，获取瓦片配置参数；步骤2，创建日志文件，检查配置参数，如果某一配置参数为空，则退出；步骤3，获取所运行计算机硬件信息，并创建线程池；步骤4，确定瓦片标准，获取遥感影像参数，根据全球经度范围和影像的经度范围计算最小层级，当全球经度小于影像经度范围时为最小层级，根据全球经度范围和单个瓦片经度范围计算最大层级，当全球经度范围小于单个瓦片的经度范围时为最大层级；步骤5，根据所运行计算机硬件信息，对遥感影像进行动态分块，遍历每一分块影像，将分块影像导入内存，计算分块影像的最小与最大层级以及仿射变换参数；步骤6，起始层级为原始影像的最小层级，遍历每一层级，判断当前层级是否大于最大层级，如果是，则切片任务结束，如果否，则计算当前层级瓦片索引编码，并将切图任务添加到线程池进行并行切片；其中每一层级的影像瓦片，按照从左到右，从上到下的顺序，依据级别、行号、列号、位置进行索引编码；步骤7，每一层级瓦片处理完之后，当前块切片结束；步骤8，每一块影像处理完之后，切片结束。2.如权利要求1所述的一种遥感影像瓦片生成及实时部署的方法，其特征在于：步骤6中线程池中单个线程进行切片的流程如下，步骤601，空间经纬度坐标转换为物理像素坐标，创建空瓦片；步骤602，判断当前层级是否小于等于分块影像的最小层级，小于等于分块影像的最小层级，执行步骤603，否则执行步骤604；步骤603，瓦片层级缺失处理，即根据瓦片编号和影像大小、影像仿射变化参数计算瓦片的起点和终点经纬度，然后根据经纬度计算像素的起点和终点坐标，根据像素的起点坐标，判断瓦片的位置，计算影像的起点像素坐标和偏移量，瓦片的起点坐标和偏移量，对待写入区域降采样处理，将数据写入瓦片对应区域，然后执行步骤605；步骤604，根据瓦片编号和影像大小、影像仿射变化参数计算瓦片的...

【专利技术属性】
技术研发人员：李治江，刘亚旗，乐晓峰，杨振雷，张过，
申请(专利权)人：武汉大学，
类型：发明
国别省市：湖北,42