版权声明:本文为博主原创文章,转载请注明原文出处!

作者:阿振

写作时间:2020年4月14日 周二


使用SNAP JAVA API处理Sentinel-1数据

SNAP软件使用Java语言开发,提供了Python接口snappy,官方教程中也多以Python接口进行示范。但是我在使用Python接口过程中,发现并不是很好用,你必须要同时懂Java语言才能很好地使用Python接口,在IDEA中使用Python接口的代码基本上没有提示,报错了也是Java的错误提示。而且,Java本来是运行在虚拟机上的语言,效率不高,再用Python包一层,更加降低了运行效率。

所以,对于我来说,SNAP的Python接口就是鸡肋,没有使用的必要。直接使用Java接口,方便程序直接进行调试,有问题可以直接去看源代码,解决了官方资料不足的问题。

SNAP GPF的使用范式

SNAP推荐使用GPF(Graph Processing Framework)进行数据处理。GPF的使用也很简单,首先基于Operation创建Product,然后写入Product。在写入过程中会自动执行Operation完成你想要的数据处理流程,处理的算法是封装在Operation中的。

使用GPF的好处是你可以进行多个Operation的流程处理,前一个处理结果直接进入后一个处理,不需要写入磁盘,可以减少磁盘IO带来的时间消耗。

对Sentinel-1 GRD数据的处理案例

下面以对Sentinel-1 GRD数据进行一系列预处理流程演示SNAP Java接口的使用。

预处理的流程包括:首先对GRD各个极化波进行段辐射校正输出后向散射系数(Sigma nought),然后进行数据裁剪,最后进行地形校正的。

下面是使用Better Java (Kotlin)进行entinel-1 GRD数据处理的源码。其中,GPF.createProduct方法需要传入Operation的名称和参数,这些信息都可以通过查看Java源代码的方式找到。

代码运行过程中需要下载DEM数据,可能会比较耗时。如果不对数据裁剪,由于原始数据太大,可能会造成内存溢出。我的笔记本根本跑不动。

package cn.demoimport java.nio.file.Paths
import org.apache.commons.io.FilenameUtilsimport kotlin.collections.HashMapimport org.esa.snap.core.dataio.ProductIO
import org.esa.snap.core.gpf.GPF
import com.bc.ceres.core.PrintWriterConciseProgressMonitorimport org.locationtech.jts.io.WKTReaderfun main() {val srcPath =Paths.get("/Users/Demo/Desktop/S1A_IW_GRDH_1SDV_20200301T104455_20200301T104520_031481_03A00B_0A9F.zip")val srcProduct = ProductIO.readProduct(srcPath.toFile())val outDir = "/Users/Demo/Desktop"val baseName = FilenameUtils.getBaseName(srcPath.toString())GPF.getDefaultInstance().operatorSpiRegistry.loadOperatorSpis()for (polar in arrayOf("VV", "VH")) {// 首先进行辐射校正(CALIBRATION)// 对应是的org.esa.s1tbx.calibration.gpf.CalibrationOp类var parameters = HashMap<String, Any>()parameters["outputSigmaBand"] = trueparameters["selectedPolarisations"] = polarval caliProduct = GPF.createProduct("Calibration", parameters, srcProduct)// 然后进行裁剪,如果不裁剪,图像太大,容易OutOfMemory// 对应的是org.esa.snap.core.gpf.common.SubsetOp类val wktRect = "POLYGON((108.175 33.873,108.782 33.873,108.782 33.129,108.175 33.129,108.175 33.873))"parameters["geoRegion"] = WKTReader().read(wktRect)parameters["bandNames"] = "Sigma0_${polar}"val subsetProduct = GPF.createProduct("Subset", parameters, caliProduct)// 然后进行地形校正(TERRAIN CORRECTION)// 对应的是org.esa.s1tbx.sar.gpf.geometric.RangeDopplerGeocodingOpval corrPath = "${outDir}/${baseName}_Corrected_${polar}"parameters.clear()parameters["pixelSpacingInMeter"] = 10.0parameters["sourceBands"] = "Sigma0_${polar}"val corrProduct = GPF.createProduct("Terrain-Correction", parameters, subsetProduct)ProductIO.writeProduct(corrProduct, corrPath, "GeoTIFF",PrintWriterConciseProgressMonitor(System.out))// 最后进行对象销毁,释放内存空间caliProduct.dispose()subsetProduct.dispose()corrProduct.dispose()}srcProduct.dispose()
}

SNAP Java API处理Sentinel-1数据相关推荐

  1. 用java api 实现查询 Hive 数据

    版本:cdh5.4.7, hive1.1.0 pom文件 <dependency><groupId>org.apache.hive</groupId><art ...

  2. 2021年大数据Kafka(五):❤️Kafka的java API编写❤️

    全网最详细的大数据Kafka文章系列,强烈建议收藏加关注! 新文章都已经列出历史文章目录,帮助大家回顾前面的知识重点. 目录 系列历史文章 Kafka的java API编写 一.生产者代码 第一步: ...

  3. 2021年大数据ZooKeeper(五):ZooKeeper Java API操作

    目录 ZooKeeper Java API操作 引入maven坐标 节点的操作 ZooKeeper Java API操作 这里操作Zookeeper的JavaAPI使用的是一套zookeeper客户端 ...

  4. Java数据持久层框架 MyBatis之API学习八(Java API详解)

    对于MyBatis的学习而言,最好去MyBatis的官方文档:http://www.mybatis.org/mybatis-3/zh/index.html 对于语言的学习而言,马上上手去编程,多多练习 ...

  5. 大数据技术之_20_Elasticsearch学习_01_概述 + 快速入门 + Java API 操作 + 创建、删除索引 + 新建、搜索、更新删除文档 + 条件查询 + 映射操作

    大数据技术之_20_Elasticsearch学习_01 一 概述 1.1 什么是搜索? 1.2 如果用数据库做搜索会怎么样? 1.3 什么是全文检索和 Lucene? 1.4 什么是 Elastic ...

  6. 大数据技术之_20_Elasticsearch学习_01_概述 + 快速入门 + Java API 操作 + 创建、删除索引 + 新建、搜索、更新删除文档 + 条件查询 + 映射操作...

    一 概述1.1 什么是搜索?1.2 如果用数据库做搜索会怎么样?1.3 什么是全文检索和 Lucene?1.4 什么是 Elasticsearch?1.5 Elasticsearch 的适用场景1.6 ...

  7. Python中使用sentinelsat包自动下载Sentinel系列数据

    Python中使用sentinelsat包自动下载Sentinel系列数据 前言 普通下载 使用工具批量下载Sentinel 数据 使用SNAP批量下载 第三方软件或网站批量下载 Sentinel-1 ...

  8. 【ZooKeeper Notes 3】ZooKeeper Java API 使用样例

    查看PDF版本 转载请注明:@ni掌柜 nileader@gmail.com ZooKeeper是一个分布式的,开放源码的分布式应用程序协调服务框架,包含一组简单的原语集合.通过这些原语言的组合使用, ...

  9. hbase建索引java api_hbase java api样例(版本1.3.1,新API)

    hbase版本:1.3.1 目的:HBase新API的使用方法. 尝试并验证了如下几种java api的使用方法. 1.创建表 2.创建表(预分区) 3.单条插入 4.批量插入 5.批量插入(客户端缓 ...

  10. MyBatis——Java API

    Java API 既然你已经知道如何配置 MyBatis 和创建映射文件,你就已经准备好来提升技能了. MyBatis 的 Java API 就是你收获你所做的努力的地方.正如你即将看到的,和 JDB ...

最新文章

  1. Python Socket编程基础篇
  2. openLDAP的编译安装以及配置
  3. 如何查看oracle trace 文件
  4. c# 定位内存快速增长_c#如何避免内存分配瓶颈以提高多线程性能
  5. 简单的docker命令ubuntu系统
  6. IBM推出云深度学习图形芯片Tesla P100 GPU加速器
  7. 利用passssh,批量远程修改机器密码
  8. 回文数python_回文数 python
  9. OneData建设探索之路:SaaS收银运营数仓建设
  10. Android开发启动未注册的activity,Hook使用demo
  11. #把函数当作参数传给另一个函数
  12. Shiro【授权过滤器、与ehcache整合、验证码、记住我】
  13. 光学算法——经典枝切法(解包裹算法)
  14. 向量化编程:arrayfun及cellfun函数的使用
  15. 将xml类型的excel文件转换为标准的Excel格式文件
  16. 计算机系毕业论文ppt模板,计算机专业答辩PPT模板.ppt
  17. 2021-03-16PyCharm3.0默认快捷键(翻译的)PyCharm Default Keymap
  18. qt textbrowser的边界框怎样改变颜色_专访天使投资人续沛川:用深度思考打破人生边界,拥有张力一生...
  19. 大数据学习路线(包含全套视频教程)
  20. C/C++ 与 Python

热门文章

  1. c语言中用数学库函数求绝对值,c语言中求绝对值数学函数.doc
  2. WindowsServer2003 2008R2架设ASP网站实战
  3. java汽车租赁系统(java租赁系统java共享汽车java汽车租借系统)java汽车租赁管理系统java自行车租赁系统共
  4. 高等数学笔记:定积分换元谬误
  5. oppo小布机器人_oppo小布小布没反应
  6. java 套接字关联的通道_java.nio.channels.SocketChannel
  7. fabric 中 peer 和 couch 容器内部数据存放目录地址
  8. 48 款数据可视化分析工具大集合
  9. Spring中的 @Lazy注解简析
  10. Java 获取、删除Word文本框中的表格