首页 > 其他分享 >ENVI5.6.2中北京三号卫星A星数据处理方法

ENVI5.6.2中北京三号卫星A星数据处理方法

时间:2022-11-01 16:36:12浏览次数:75  
标签:选择 校正 全色 NNDiffuse 三号 ENVI5.6 控制点 图像 数据处理

北京三号卫星A星于2021年6月11日,搭载了0.5米全色、2米多光谱(蓝、绿、红、近红外波段)相机,成像宽幅优于23km。

ENVI5.6.2中完全支持北京三号卫星A星图像数据处理,包括数据打开、辐射定标、大气校正、正射校正、图像融合等处理。本文介绍利用北京三号卫星A星生产DOM产品,包括正射校正和图像融合处理。其他版本类似的处理过程。

ENVI5.6.2软件的试用请浏览:envi.geoscene.cn/envi_license

1、数据打开

请使用易智瑞公司遥感事业部开发的“中国国产卫星支持工具”打开北京三号卫星A星数据。

建议使用ENVI App Store进行ENVI扩展工具的安装与管理。

ENVI App Store下载地址:envi.geoscene.cn/appstore

  • 在App Store找到中国国产卫星支持工具,点击安装进行安装。
  • 重启 ENVI即可。

或手动安装扩展工具: envi.geoscene.cn/china

  • 启动ENVI。选择主菜单 File > Open As > China Satellites >BJ3A。选择XML文件打开数据。
  • 在Layer Manager中,右键选择菜单View Metadata,选择各个选项可以看到自动读取的包括RPC、定标参数、拍摄参数、成像时间等元数据信息。

2、正射校正

正射校正可以选择无控制点/手动选择控制点正射校正和参考影像自动选择控制点正射校正两种工具。

2.1无控制点/手动选择控制点正射校正

在Toolbox中选择/Geometric Correction/Orthorectification/ RPC Orthorectification Workflow工具,该工具是流程化操作方式操作比较简单,不选择控制点将以无控制点正射校正输出。

在RPC Orthorectification Workflow面板中,输出像元大小(Output Pixel Size)默认自动计算整个图像的平均值,可以手动修改为全色0.5米、多光谱2米。

注:NNDiffuse Pan Sharpening图像融合算法要求全色图像分辨率是多光谱图像的整数倍。

如下图为无控制点的正射校正结果叠加在天地图影像底图上,天地图影像放大最高级切片,局部地方有几米的几何偏差。在精度要求很高的应用场景还是需要利用控制点进行正射校正。

图:无控制点正射校正精度对比

2.2参考影像自动选择控制点正射校正

参考影像应选择与待正射校正接近分辨率的影像,如0.5~2米。

在Toolbox中选择/Geometric Correction/Orthorectification/RPC Orthorectification Using Reference Image,打开基于基准影像的RPC正射校正工具;在面板上进行数据选择和参数设置:

  • Input Raster:选择前面打开的北京三号卫星A星全色或者多光谱图像;
  • Input Reference Raster:具有准确地理位置信息的参考影像;
  • Input DEM Raster:覆盖待校正影像区域的DEM数据;
  • DEM Is Height Above Ellipsoid:DEM数据是否是椭球高;
  • Requested Number of GCPs:需要采集的控制点最大数量;
  • Search Windows Size搜索窗口大小,窗口越大找点精度越高,时间越长;
  • Output Coordinate System:可选项,如果不设置,则输出的正射校正结果是UTM WGS84坐标系的;
  • Output Pixel Size:输出像元大小,北京三号卫星A星全色为5米,多光谱为2米。

注:输出像元大小为可选项,如果不设置,则输出的像元大小为根据RPC自动计算的平均像元大小,建议手动设置输出数据的分辨率。

  • Image Resampling:重采样方法,默认双线性内插Bilinear;
  • Grid Spacing:像素栅格间距。保持默认10;这个值越大,校正速度越快,精度越低。
  • Output Raster:正射校正结果输出路径及文件名;
  • Output GCPs:控制点文件输出路径及文件名;

参数设置好之后,点击OK开始执行。

3、图像融合

推荐使用NNDiffuse 图像融合方法,该方法能很好的保留纹理和光谱信息。由于NNDiffuse方法对储存顺序为BIP或者BIL的多光谱图像处理效率更高,处理效率大约为BSQ储存顺序2~3倍。因此我们可以利用/Raster Management/Convert Interleave工具对储存顺序为BSQ的多光谱图像转为BIP或者BIL。

在Toolbox中选择/Image Sharpening/NNDiffuse Pan Sharpening,在NNDiffuse Pan Sharpening参数面板中分别选择全色和多光谱正射校正结果,其他参数选择默认。

  • Pixel Size Ratio多光谱与全色图像的分辨率比值,不设置情况自动计算。NNDiffuse算法要求是整数倍。
  • Spatial Smoothness设置空间重采样插值核参数,不填写默认为分辨率比值62。
  • Intensity Smoothness平滑强度值,值越小得到的结果图像越清晰(类似高通滤波效果),可用于图像可视化。值越大得到的图像噪声越少,可用于图像分类或者定量反演。不填写默认为自动根据图像计算,可填写1~20之间的值。

图:NNDiffuse Pan Sharpening参数面板

如下图为真彩色和标准假彩色合成展示融合结果。

图:真彩色合成结果

图:标准假彩色合成结果

如下图为融合前后的效果对比。

图:耕地

图:清晰可见的电力塔和电力线

图:丘陵林地

图:村庄

标签:选择,校正,全色,NNDiffuse,三号,ENVI5.6,控制点,图像,数据处理
From: https://www.cnblogs.com/enviidl/p/16848150.html

相关文章

  • js数据处理4
    <template><ModuleBox2title="数据趋势图"bgType="2":cityShow="true"moduleType="sjqs":styleObj="styleObj"><divclass="echarts-box">......
  • 大数据处理架构Hadoop
    Hadoop可以支持多种编程语言HDFS分布式文件系统MapReduce分布式并行编程框架这两大核心组件用于解决怎么实现海量数据的分布式存储和怎么实现海量数据分布式处理Hiv......
  • python3使用libpcap库进行抓包及数据处理
    python版本:python3.9libpcap版本:1.11.0b7pythonlibpcap库是底层绑定c语言libpcap库的开发包,旨在提供python应用可访问的unixclibpcap库API(以及为win32系统提供的Npca......
  • 前端项目实战68-数据处理之一个数组和一个对象
    exportfunctionmyContact(target:any,source:any){for(const[key,value]ofObject.entries(source)){const[name,index]=key.split('-')//要......
  • js数据处理4
    <template><ModuleBox2title="数据趋势图"bgType="2":cityShow="true"moduleType="sjqs":styleObj="styleObj"><divclass="echarts......
  • python3使用mutagen进行音频元数据处理
    python版本:python3.9  mutagen版本:1.46.0mutagen是一个处理音频元数据的python模块,支持多种音频格式,是一个纯粹的python库,仅依赖python标准库,可在Python3.7及以上版......
  • 复杂场景数据处理的 OLTP 与 OLAP 融合实践
    本文首发于NebulaGraph公众号DagController介绍DagController是NebulaGraph企业版的系统,经过反复测试无误后进行了发布,它主要解决的是OLTP和OLAP的融合问题......
  • python 图形的数据处理 (折线图为例)
    1.通过json模块对数据进行处理ab173.com是懒人工具-json在线解析,可以通过他对json数据进行格式化的分析。"""演示可视化需求1:折线图开发"""importjsonfrompyec......
  • openmetadata 元数据处理标准
    openmetadata对于元数据的处理也是基于了schema主要是workflow这个定义,包含了不同的定义,基于自己定义的workflow进行执行数据的获取、处理以及写入参考workflow定义......
  • 基础09:数据处理之增删改
    一、插入数据1.1方式1:VALUES的方式添加使用这种语法一次只能向表中插入一条数据。情况1:为表的所有字段按默认顺序插入数据--语法格式如下:INSERTINTO表名VALU......