DorisInSAR处理基本步骤和基本参数设定
doris使用手册2
三、 input.s_initial 批处理PROCESS s_readfiles // reading of slc filesPROCESS s_crop // convert data to internal format PROCESS s_porbits // get precise orbits 此步骤和第一步是一样的处理过程。
是对副图像处理。
四、 input.coregistration 批处理 PROCESS coarseorb通过轨道参数粗略计算主图像和副图像之间的平移量,精度在30个象元左右。
这是一个快速的得到初始偏移量的方法,在精匹配之前必须要做粗匹配。
这里的偏移量是这样定义的:(,)(,)(,)s m P l p P l p offset l p =+这个算法步骤可以分为3步:1、利用主图像的中心像元(l,p ),计算出对应的在参考椭球体上的点P (x,y,z )。
2、根据多普勒方程计算出副图像的卫星位置,得到其P 点,然后再计算得到其副图像的中心像元(l,p )。
3、根据主副图像之间的中心像元(l,p )差值来得出偏移量。
这里要用到3个方程:多普勒方程E1,距离方程E2,椭圆体方程E3。
1222223222ˆ:0:()0:10ss light range dx x x E xdx E dx dx v t x y z E a a b=-⋅=⋅-=++-=PROCESS coarsecorr 粗校正,利用相关技术计算主图像和副图像之间的象元级的方位向和距离向平移量。
幅值图像在这里需要用到,还要计算空间域或者谱域的相关系数。
主副图像之间具有最高相关系数的偏移量点就被作为该点的偏移量,也就是最相似的点。
CC_METHOD magfft /magspace有2种方法可以选择,magfft是最大估计法,magspace是相干系数法。
这里的2个方法都是在幅值图像中计算相关系数。
magfft这种方法是基于干涉条纹频谱的最大信噪比进行图像配准的算法,也称最大谱估计法,由Gabriel等(1988)在处理SIR-B数据时提出,算法也比较复杂,但是匹配精度很高。
doris insert语句
doris insert语句Doris是一个开源的分布式列式存储系统,支持实时的OLAP查询和分析。
它采用了分布式架构,可以横向扩展以处理大规模数据集。
在使用Doris时,我们可以通过INSERT语句将数据插入到Doris中,下面是一些关于Doris INSERT语句的示例。
1. 插入单行数据INSERT INTO table_name (column1, column2, column3) VALUES (value1, value2, value3);这是最基本的INSERT语句,用于向指定表中插入一行数据。
column1, column2, column3是表中的列名,value1, value2, value3是要插入的数据值。
2. 插入多行数据INSERT INTO table_name (column1, column2, column3)VALUES (value1, value2, value3), (value4, value5, value6), (value7, value8, value9);如果需要同时插入多行数据,可以在VALUES子句中使用逗号分隔的多个value列表。
3. 插入查询结果INSERT INTO table_name (column1, column2, column3)SELECT column1, column2, column3 FROM another_table WHERE condition;除了直接插入值之外,还可以通过SELECT语句查询出结果后再插入到目标表中。
这对于批量数据导入非常有用。
4. 插入部分列数据INSERT INTO table_name (column1, column2)VALUES (value1, value2);如果只需要插入表中的部分列数据,可以在INSERT语句中指定要插入的列名和对应的值。
5. 插入默认值INSERT INTO table_name DEFAULT VALUES;如果在插入数据时,某些列使用默认值,可以使用DEFAULT VALUES 关键字。
starrocks使用手册
StarRocks,也被称为Apache Doris,是一款高性能、高并发、低延迟的MPP(大规模并行处理)数据库。
以下是StarRocks的使用手册:1. 安装与配置:安装StarRocks之前,需要先准备一个合适的操作系统环境,如Linux。
然后,可以从StarRocks的官方网站或GitHub仓库下载安装包,按照说明进行安装。
安装完成后,需要进行配置。
主要配置文件为`conf/server.conf`,需要设置数据库的基本信息,如端口号、内存大小等。
2. 创建数据库:在StarRocks中,数据库是存储数据的基本单位。
可以使用以下SQL 语句创建一个新的数据库:```sqlCREATE DATABASE database_name;```3. 创建表:在数据库中,表是存储数据的基本单位。
可以使用以下SQL语句创建一个新的表:```sqlCREATE TABLE table_name (column1 data_type,column2 data_type,...);```其中,`column1`, `column2`等是列名,`data_type`是对应的数据类型。
4. 插入数据:可以使用以下SQL语句向表中插入数据:```sqlINSERT INTO table_name (column1, column2, ...) VALUES (value1, value2, ...);```5. 查询数据:可以使用以下SQL语句查询表中的数据:```sqlSELECT * FROM table_name;```这将返回表中的所有数据。
你也可以指定查询的列,使用条件语句进行筛选等。
6. 其他操作:StarRocks还支持更新、删除、排序等操作,具体使用方式可以参考StarRocks的官方文档或使用手册。
doris 数据写入流程
doris 数据写入流程下载温馨提示:该文档是我店铺精心编制而成,希望大家下载以后,能够帮助大家解决实际的问题。
文档下载后可定制随意修改,请根据实际需要进行相应的调整和使用,谢谢!并且,本店铺为大家提供各种各样类型的实用资料,如教育随笔、日记赏析、句子摘抄、古诗大全、经典美文、话题作文、工作总结、词语解析、文案摘录、其他资料等等,如想了解不同资料格式和写法,敬请关注!Download tips: This document is carefully compiled by theeditor. I hope that after you download them,they can help yousolve practical problems. The document can be customized andmodified after downloading,please adjust and use it according toactual needs, thank you!In addition, our shop provides you with various types ofpractical materials,such as educational essays, diaryappreciation,sentence excerpts,ancient poems,classic articles,topic composition,work summary,word parsing,copy excerpts,other materials and so on,want to know different data formats andwriting methods,please pay attention!Doris 数据写入流程一、准备工作阶段。
在进行 Doris 数据写入之前,需要进行充分的准备。
Doris处理过程及参数说明
1,命令介绍.生成模板:run –g编辑模板:编辑第一步模板:run –e1编辑第二步模板:run –e2编辑第三步模板:run –e3编辑第四步模板:run –e4编辑第五步模板:run –e5编辑第六步模板:run –e6编辑第七步模板:run –e7编辑第八步模板:run –e9(每行前加 C 表示屏蔽该行。
打开模板文件后,按“i“表示可以对进行插入操作,如果删去时,直接用delete或后退键。
按ESC键时则退出当前的编辑状态。
若要退出文件时,则先按ESC,然后:q!,或:wq,其中:q!表示退出不保存修改,:wq为保存修改。
)运行处理:运行第一步:run –s1运行第二步:run –s2运行第三步:run –s3运行第四步:run –s4运行第五步:run –s5运行第六步:run –s6运行第七步:run –s7运行第八步:run –s92,参数解说及设置。
2.1,第一步:读取主影像信息。
(1)需要处理步骤。
必需:m_readfiles(读取主影像头文件信息).必需:m_crop(读取主影像数据体信息).不需:m_ovs(主影像过采样步骤):作干涉处理时,一般不运行。
可选:m_prorbits(读取主影像精轨文件):若要使用提供的精密轨道数据时,就执行该行,否则直接C掉。
(2)小步骤参数说明及设置。
############################################################INSAR通用参数设置,每步都有的。
###########################################################名称:LOGFILE用处:用以记载DORIS处理过程中所有信息。
文件名:doris.log.名称:M_RESFILE用处:用以记载主影像文件所有信息。
文件名:master.res名称:S_RESFILE用处:用以记载辅影像文件所有信息。
doris 接入流程 标准
Doris 接入流程标准一、需求分析在进行Doris 接入之前,首先要明确业务需求和数据量大小,评估Doris 的性能和扩展性是否满足要求。
同时,还需要考虑现有系统与Doris 的兼容性和对接方式。
二、接入准备1. 安装和配置Doris:根据官方文档,安装和配置Doris,确保其正常运行。
2. 确定数据源和目标:确定需要迁移的数据源和Doris 的目标数据库,确保数据可以正确地导入到Doris 中。
三、数据迁移1. 数据导出:从数据源中导出需要迁移的数据,通常可以使用工具或编写脚本实现。
2. 数据转换:对导出的数据进行必要的转换和清洗,以满足Doris 的数据格式和要求。
3. 数据导入:将转换后的数据导入到Doris 的目标数据库中,通常可以使用Doris 自带的导入工具或编写脚本实现。
四、系统配置1. 调整Doris 参数:根据实际情况,调整Doris 的参数,如内存、CPU、磁盘等,以优化性能和扩展性。
2. 配置数据源连接:配置Doris 与数据源的连接,确保数据可以正确地导入到Doris 中。
3. 配置安全策略:根据实际需求,配置Doris 的安全策略,如访问控制、数据加密等。
五、测试与验证1. 功能测试:测试Doris 系统的各项功能是否正常,如查询、插入、更新、删除等操作。
2. 性能测试:测试Doris 系统的性能是否满足要求,如查询速度、并发处理能力等。
3. 数据验证:验证迁移后的数据是否正确,如数据完整性、一致性等。
六、监控与优化1. 监控Doris 系统运行状态:监控Doris 系统的各项指标,如CPU、内存、磁盘使用情况等,确保系统正常运行。
2. 优化性能:根据监控结果,优化Doris 系统的性能,如调整参数、优化查询语句等。
3. 故障恢复:在出现故障时,及时进行故障恢复,确保系统正常运行。
七、文档整理1. 编写操作手册:编写Doris 系统的操作手册,包括安装、配置、使用等方面的说明。
doris读写流程
doris读写流程
在Doris中,数据的读写流程大致如下:
在写流程中,ColumnWriter 将各自数据、索引生成的 Page 顺序写入到文件中。
具体流程可能包括:数据经过压缩、编码、序列化等过程,最后以分页为单位,将每一页写入到底层的存储文件中。
在这个过程中,还会进行一些列的存储格式的优化,以提高数据查询效率。
至于读流程,一个名为OlapScanner的组件会对一个tablet数据读取操作进行整体的封装。
主要内容如下:
1. 根据查询指定的version版本查找出需要读取的RowsetReader。
2. 设置查询信息,包括tablet、读取类型reader_type=READER_QUERY、是否进行聚合、version。
3. 设置查询条件信息,包括filter过滤字段、is_nulls字段。
4. 置返回列信息。
5. 设置查询的key_ranges范围。
6. 初始化Reader对象。
在初始化Reader对象后,会进行一系列的处理,包括对读取参数的处理以及对不同模型读取方式的调整。
然后,会初始化查询条件对象。
至此,整个读流程就完成了。
以上是Doris中数据的读写流程的大致介绍,更具体和详细的内容可能需要查阅相关文档或者与专业的技术人员交流才能获取。
INSAR技术原理及方案
INSAR沉降监测技术方案简介一.技术方案1.1 主要功能INSAR系统平台集成了功能强大的合成孔径雷达(Synthetic Aperture Radar)影像数据处理系统,整合了多种先进高效的处理算法,并支持用户嵌入开发自主算法模型。
基于并行多线程处理技术,能够快速处理分析海量雷达影像数据。
除了具有传统InSAR算法模型之外,INSAR系统平台还具有以下技术特点:-大气分析技术;-季节性形变分析技术;-非线性形变分析技术;-城市区域PS点层析功能。
- 全图形用户界面;●软件界面全图形化;●处理进度监控管理全图形化;●中间处理结果、报告图形化显示;●最终结果图形化显示输出;- 采用并行多线程处理技术,高效处理分析InSAR数据集;- 高度的兼容性:●跨平台处理系统:支持Windows、Linux和Mac OS;●支持所有在轨SAR传感器数据模式,包括最新发射的ESASentinel, PAL-SAR2等;●兼容其他软件系统文件格式:GAMMA,DORIS,ROI_PAC等;- 工程级别处理精度:●识别区分建筑物和地面上的PS点,区分建筑物形变和地面形变;●精确反演影像采集时间的大气相位图,剔除大气相位残值;●精确反演由温度变化导致的季节性热形变;精确监测非线性形变;对于城市区域或强反射硬目标形变速率监测精度可达3-5毫米。
1.2工作流程雷达影像数据的地面沉降监测服务基于人造多普勒雷达的相位干涉测量技术(InSAR),通过精确测量监测数据集中雷达影像的电磁波相位变化情况,去除大气和噪声干扰,获取地表的真实形变信息。
雷达干涉测量技术自问世以来就被广泛应用于测量和形变监测领域并形成了系统成熟的工程化解决方案,根据监测数据集情况、监测技术手段和测量结果的不同可以将地面沉降监测项目分成3个生产阶段:DInSAR、SBAS-InSAR、PSInSAR。
1.2.1 DInSAR地面沉降监测工作流程在工作启动初期,由于监测数据集影像数量较少,可以采用DInSAR(差分干涉)技术对监测数据进行定性处理分析,快速探测监测区域内形变中心分布情况及其影响范围,粗略估算形变速率及预测发展趋势。
doris with as 用法
doris with as 用法Doris With As 是一种在Python编程语言中常用的语法结构,用于将一个变量赋值给另一个变量,实现两个变量指向同一个对象的效果。
本文将详细介绍Doris With As 的使用方式以及相关注意事项。
一、Doris With As 的基本用法Doris With As 语法的基本格式如下:```with 表达式 as 变量:# 执行代码块,使用变量进行相关操作```在这个语法结构中,表达式指的是一个需要执行的代码块或函数,变量则是用于接收表达式返回结果的变量。
使用Doris With As 的主要目的是确保在代码执行完毕后,相关资源可以被正确释放,例如文件的关闭、数据库连接的断开等。
Doris With As 可以自动处理这些资源释放的过程,使得代码更加简洁,同时避免了手动释放资源可能导致的遗漏或错误。
二、Doris With As 的示例为了更好地理解Doris With As 的使用方式,以下是一个具体的示例:```with open('data.txt', 'r') as file:content = file.read()print(content)```在这个示例中,`open('data.txt', 'r')` 是一个表达式,用于打开名为'data.txt' 的文件,并赋值给变量 `file`。
`as` 关键字后的 `file` 是用于接收打开文件对象的变量。
使用 `with` 结构后,可以在代码块内部执行相关操作,而不需要手动关闭文件。
三、Doris With As 的注意事项在使用Doris With As 时,需要注意以下几个问题:1. Doris With As 只能用于支持上下文管理器协议的对象。
上下文管理器是指实现了 `__enter__` 和 `__exit__` 方法的对象。
用InSAR 技术进行形变监测的研究
第31卷第3期2006年5月测绘科学Science of Surveying and M app ingVol131No13M ay收稿日期:2006Ο02Ο14用I nSAR技术进行形变监测的研究成英燕①,贾有良②,党亚民①,章传银①(①中国测绘科学研究院,北京 100039;②天津市测绘院,天津 300041)【摘 要】卫星合成孔径雷达干涉技术((I nS AR)通过对地面同一地区两次或多次平行观测成像得到的复图象对进行处理可获得数字地面模型,用以监测地表面变化。
本文简单叙述了(I nS AR)基本原理、处理流程及方法,并通过天津S AR数据的处理,详细介绍了S AR数据处理整个过程,之后给出了基于ERSΟ2天津地区原始数据的实验结果形成的干涉图及数字高程图。
【关键词】合成孔径雷达;差分干涉;数字高程【中图分类号】P228 【文献标识码】A 【文章编号】1009Ο2307(2006)03Ο0056Ο031 引言I nS AR(I nterfer ometry Synthetic Aperture Radar)技术主要用于探测地表形变、地面沉降,有越来越多的国家和地区利用I nS AR来探测诸如由于矿山开采、地震、火山运动、地下水开采等引起的地面沉降现象。
干涉雷达优点较多:具有全天候工作能力,发射的微波对地物有一定的穿透能力,能提供光学遥感所不能提供的信息,且是主动式工作方式。
因ERSΟ1已处于瘫痪,目前用于I nS AR技术研究的数据来源主要有:ERSΟ2、SI R2C/XS AR、RADRS ATΟ1、ERS21、T OPS AR和SE AS AT等。
对于欧洲雷达卫星ERSΟ2和加拿大雷达卫星RADRS ATΟ1,采用干涉技术来产生DE M,监测地面位移变化,精度可以达到mm量级。
因此,该技术手段特别适于解决大面积的滑坡、崩塌、泥石流以及地裂缝、地面沉降等地质灾害的监测预报,是一项快速、经济的空间探测高新技术。
DorisInSAR处理基本步骤和基本参数设定
1,Doris处理可以分步处理,也可以整体处理。
Doris处理要先生成模板文件,doris每步处理都对应一个模板文件,这些模板文件包含了doris每步处理所需要的参数,因此用户可以直接打开模板文件,用以编辑更改然后保存,也可以在linux终端中采用Vi命令,或nedit(相当于window下的写字板)来进行编辑保存。
2,查看版本信息,终端键入:doris 即可查看。
3,生成模板。
命令:run –g说明:会在当前工作目录下生成模板文件,它包含三个文件夹目录:Inputfiles,Outdata ,Outinfo。
Outdata用于存放每步生成的结果数据文件。
Outinfo用于存放每步处理的log文件,记录中间处理过程。
Inputfiles中存放的就是模板文件。
它包含的文件:对这些模板进行编辑时,文件中每行前面加c和空格,表示程序处理时不读取该行参数,否则都会执行到的,编辑参数选项时直接更改后面的的数值即可。
4,doris处理步骤总体介绍。
Doris的INSAR处理共分八步:每步又可以分为许多小步,可以通过编辑来选择执行或不执行某些小步骤(小步骤在模板文件中的标识符为process)第一步run –s1:读入主影像的头文件及数据体信息以及轨道信息(可以使用粗轨),其对应的模板文件 input.m_initial.,这步包含:读取主影像头文件信息,读取主影像数据体信息,主影像过采样,读取主影像精轨文件。
(如下:加c 表示屏蔽该行,以下都同)第二步run –s2:读入辅影像的头文件及数据体信息以及轨道信息(可以使用粗轨),其对应的模板文件:input.s_initial。
这步包含:读取辅影像头文件信息,读取辅影像数据体信息,辅影像过采样,读取辅影像精轨文件。
第三步run –s3:影像配准。
对应的模板文件:input.coregistration.包含步骤有:一级轨道级配准,二级像素级粗配准,三级子像素级精配准,主影像方位向滤波,辅影像方位向滤波,以及影像距离向滤波。
doris使用手册1
一、前言Doris软件全称为Delft Object-oriented Radar Interferometric Software,是荷兰Delft大学的Kampes等人用面向对象的C++语言编写的。
该软件主要开发目的是为研究地表的三维地形及地表变形,为完全免费软件,其源代码是开放的,利于研究者使用。
该软件的主要运行环境为Unix系统,但其在Linux及Cygwin 环境下运行效率也很高。
Doris数据处理大致可以分为四个模块:原始数据读取模块、配准和参考椭球相位计算模块、复相位图像和相干图计算模块、DEM和形变干涉图生成模块。
每个模块中有独立的小模块,在进行数据处理时,每个模块进行独立的计算,处理完后有输出纪录处理的结果及数据,下一模块在执行时则先读取所需的参数及数据,这样每个模块都可以单独进行测试,同时也可以选择不同的算法和软件来实现其中的任意流程,以了解其运算的情形来选择最合适的参数值。
由于Doris 处理的对象是图像矩阵,所以,为了便于再开发和再研究,doris各模块运行产生的中间数据皆为裸矩阵(无头无尾文件)。
Doris的文件命令是由多个批处理文件组成,文件里先规定了要进程(PROCESS)哪些步骤,然后在每个进程(PROCESS)的步骤后,添加执行的关键字(CARD)。
各主要PROCESS如下表所示:上面已经概要地介绍了每个进程的作用,下面的章节将对具体的进程和关键字做介绍并附加具体实验结果。
为了方便起见,本文将一整套Doris工作流程分为七个步骤,每个步骤由一个批处理文件进行。
这七个批处理文件是:1.input.m_initial2.input.s_initial3.input.coregistration4.input.resample5.input.products6.input.filter_unwrap7.input.s2h_geocode以上七个文件有着相同的规格,包括题头、一般性选项、I/O文件和进程四部分。
DORIS软件应用于D-InSAR技术的研究
可用 于地 壳变形 、 震变形 、 面沉降 、 地 地 山体 滑坡 、
火 山变形 和冰川 移 动等较 大 面积 的地 表 变形 测量 和大 型工 程 ( 如水 库 、 坝 、 梁 和 管 线 ) 大 桥 的形 变
监测 , 尤其 在 突发 性 的 形 变监 测 中可 获 得 较 好 的
V0. No 3 J21 . S t 2 08 ep . 0
D R S软 件 应 用 于 D—IS R技 术 的 研 究 O I nA
边 磊 王 猛2 ,
,. 海大学 土木工程 学院 , 苏 南 京 1河 江 209 ; 1o 8
\. 2 山东省淄博市公路管理局 青莱高速公路路政大队, 山东 淄博
干涉 图 或 精 确 的 D M。这 种 方 法 能 感 受 到 m E m
常规遥 感 测量 技术 不但 需 要 地 面控 制 点 , 而
且 受 时间 、 间及 气 象 条 件 的 限制 。新 近 发 展 的 空 星载合 成孔径 雷 达干 涉 测量 (nA 技术 主要 应 IS R) 用 于 常规光学 遥 感 图像难 以获 得 的地 区进 行 地形
借 助 于 D R S软 件 进行 的 D—IS R 处理 可 以获得 m 级 的形 变量 。 O I nA m 关键词 : 面 变形 ; O I 地 D R S软件 ; D—IS n AR技 术 ; a 古城 Bm
中图分类 号 : P 5 . T 7 11
文 献标 识码 : A
文 章编 号 :6 1— 3 2 20 ) 3— 0 5— 4 17 5 2 ( 0 8 0 0 0 0
软件 是使 用 面向对 象 的 C++ 言 编 写 的 , 若 干 语 以 公 共 开 源 软 件 ( Fv , E O B, N P U, 如 fw G T R S A H r G T等 ) 辅 助 , 要 运 行 环 境 为 U i 作 系 M 为 主 nx操
基于时序InSAR技术的蔚县地面形变监测
基于时序InSAR技术的蔚县地面形变监测仝云霄;高井祥;陈宇【摘要】合成孔径雷达干涉测量(InSAR)技术在城市地表形变监测中具有高精度、高分辨率、低成本、空间连续监测等优势.以河北省张家口市蔚县为研究区域,获取了20景TerraSAR-X影像,分别采用PS-InSAR和SBAS-In-SAR两种时序技术得到了蔚县城区2015年6月9日到2016年1月4日期间地面平均沉降速率,结果表明,两种时序InSAR技术监测沉降结果具有很高一致性,而且相关性较高,达到0.9以上,城区中心发生了比较严重的地表沉降,形成了一个明显的沉降漏斗,沉降速率达到了30 mm/a,地下水的严重开采以及建筑物的加速构建是引起该区域地面沉降的主要原因,该研究可以为城市地面沉降治理、整体规划等提供参考意义.【期刊名称】《现代测绘》【年(卷),期】2019(042)002【总页数】6页(P1-6)【关键词】PS-InSAR;SBAS-InSAR;城市地面沉降;TerraSAR-X【作者】仝云霄;高井祥;陈宇【作者单位】中国矿业大学,江苏省资源环境信息工程重点实验室,江苏徐州221116;中国矿业大学,江苏省资源环境信息工程重点实验室,江苏徐州 221116;中国矿业大学,江苏省资源环境信息工程重点实验室,江苏徐州 221116【正文语种】中文【中图分类】TP722.60 引言地面沉降是指在一定的地表面积内不断加大对地下资源开发而产生的地面形变现象,是一种缓慢性地质灾害,可导致不可补偿的永久性环境资源损失。
地面沉降的产生与人类的各种经济活动密切相关,近年来,对地下水、石油、天然气、矿产资源的过度开采,以及隧道、地铁等地下工程的建设导致我国很多城市和地区出现了比较严重的沉降现象。
据统计,目前我国发生地面沉降的城市和地区大约有100多个,其中华北平原、长江三角洲以及汾渭盆地沉降较为严重[1],因此有必要加强对城市地面沉降的监测。
传统的形变监测方法主要有大地水准测量、GNSS测量和全站仪测量等,这些方法监测精度较高,在形变监测领域发挥了重要的作用,但是也存在着成本高、工作周期长、工作量大、无法获取空间连续信息等缺陷。
用doris进行INSAR的操作过程
以下摘自:/blog/static/839782482009599 583826/Xiaoyan@hp-47c03a5dcc20/cygdrive/e/insar/doristest_Data/Your_work$ cd e:Xiaoyan@hp-47c03a5dcc20/cygdrive/e$ cd insar/Xiaoyan@hp-47c03a5dcc20/cygdrive/e/insar$ cd testXiaoyan@hp-47c03a5dcc20/cygdrive/e/insar/test$ run -genvironment variable PAGER not set, please set this variable to specify your prefered viewer. For example in csh shell type: setenv PAGER moresetenv PAGER less (adviced)Missing filename ("less --help" for help)using PAGER moreenvironment variable EDITOR not set, please set this variable to specify your prefered editor. For example in csh shell type:setenv EDITOR visetenv EDITOR neditsetenv EDITOR xemacsNEdit: Can't open displayusing EDITOR virun: ***INFO*** generation of inputfiles finished.run: ***INFO*** inputdir: Inputfilesrun: ***INFO*** outputdir (stout): Outinforun: ***INFO*** outputdir (matrices): Outdatarun: ***INFO*** Using executable in path: dorisrun: ***INFO*** or (if -d option present): doris.debug run: ***TIP*** run -e1; run -s1Xiaoyan@hp-47c03a5dcc20/cygdrive/e/insar/test$ run -s1environment variable PAGER not set, please set this variable to specify your prefered viewer. For example in csh shell type: setenv PAGER moresetenv PAGER less (adviced)Missing filename ("less --help" for help)using PAGER moreenvironment variable EDITOR not set, please set this variableto specify your prefered editor. For example in csh shell type: setenv EDITOR visetenv EDITOR neditsetenv EDITOR xemacsNEdit: Can't open displayusing EDITOR vidoris Inputfiles/input.m_initial >>Outinfo/stdout.input.m_initialPROGRESS: Interpretation inputoptionsfile finished.total cpu: 0 min 0.015 secPROGRESS: Finished initializationPROGRESS: Start M_READFILES.PROGRESS: readvolume finished.WARNING : c4missingrawlines: 2PROGRESS: readleader finished.PROGRESS: readdat (header info) finished.PROGRESS: Finished M_READFILES.total cpu: 0 min 0.046 secPROGRESS: Start M_PORBITS.(getorb: reading arclist E:/insar/zhengjiang/orbit/ers1/arclist) (getorb: orbit file E:/insar/zhengjiang/orbit/ers1/ODR.467 ...leapsecond: 960101)PROGRESS: getorb: program finished ok.PROGRESS: Finished M_PORBITS.total cpu: 0 min 0.34 secPROGRESS: Orbit: interpolation coefficients computed. PROGRESS: Start M_CROP.PROGRESS: Start cropping slc data.WARNING : writeslc: length of record 1 = "16012"; expected "19612" for ERS SLC (CEOS, full scene).PROGRESS: WRITESLC: 0%PROGRESS: WRITESLC: 10%PROGRESS: WRITESLC: 20%PROGRESS: WRITESLC: 30%PROGRESS: WRITESLC: 40%PROGRESS: WRITESLC: 50%PROGRESS: WRITESLC: 60%PROGRESS: WRITESLC: 70%PROGRESS: WRITESLC: 80%PROGRESS: WRITESLC: 90%PROGRESS: WRITESLC: 100%PROGRESS: Finished M_CROP.total cpu: 0 min 12.497 secPROGRESS: calling preview for cropped masterPROGRESS: Start PREVIEW generation.PROGRESS: SUNraster file created of: Outdata/master.raw (see also file: ./master_mag.ras.sh)total cpu: 0 min 12.588 sec...An expert is a man who has made all the mistakes, which can be made,in a very narrow field.--- WARNING SUMMARY ---There were 2 messages:1: c4missingrawlines: 22: writeslc: length of record 1 = "16012"; expected "19612" for ERS SLC (CEOS, full scene).Redirected output in file:Outinfo/stdout.input.m_initialview with: run -v1NOTE: before run -s1, you can do run -ql to get a quicklook processing.which will run step 1 and 2 if not done already.Xiaoyan@hp-47c03a5dcc20/cygdrive/e/insar/test$ @(#)Doris software, $Revision: 3.11 $, $Author: kampes $ cpxfiddle: INFO: Number of output pixels: 500Xiaoyan@hp-47c03a5dcc20/cygdrive/e/insar/test$ run -s2environment variable PAGER not set, please set this variable to specify your prefered viewer. For example in csh shell type: setenv PAGER moresetenv PAGER less (adviced)Missing filename ("less --help" for help)using PAGER moreenvironment variable EDITOR not set, please set this variable to specify your prefered editor. For example in csh shell type: setenv EDITOR visetenv EDITOR neditsetenv EDITOR xemacsNEdit: Can't open displayusing EDITOR vidoris Inputfiles/input.s_initial >>Outinfo/stdout.input.s_initialPROGRESS: Interpretation inputoptionsfile finished.total cpu: 0 min 0.015 secPROGRESS: Finished initializationPROGRESS: Orbit: interpolation coefficients computed. PROGRESS: Start S_READFILES.WARNING : Found identifier, but not so desired word!, exiting routine...PROGRESS: readvolume finished.WARNING : c4SWSTchange: 1PROGRESS: readleader finished.PROGRESS: readdat (header info) finished.PROGRESS: Finished S_READFILES.total cpu: 0 min 0.015 secPROGRESS: Start S_PORBITS.(getorb: reading arclist E:/insar/zhengjiang/orbit/ERS2/arclist) (getorb: orbit file E:/insar/zhengjiang/orbit/ERS2/ODR.074) PROGRESS: getorb: program finished ok.PROGRESS: Finished S_PORBITS.total cpu: 0 min 0.31 secPROGRESS: Orbit: interpolation coefficients computed. PROGRESS: Start S_CROP.PROGRESS: Start cropping slc data.ERROR : [[2205]]: Unable to open input file:/cdrom/scene1/dat_01.001 in source file: at line: 2915Caught error of EXCEPTION class!It is: file errortotal cpu: 0 min 0.607 sec--- WARNING SUMMARY ---There were 2 messages:1: Found identifier, but not so desired word!, exiting routine... 2: c4SWSTchange: 1Redirected output in file:Outinfo/stdout.input.s_initialview with: run -v2Xiaoyan@hp-47c03a5dcc20/cygdrive/e/insar/test$ run -s2environment variable PAGER not set, please set this variable to specify your prefered viewer. For example in csh shell type: setenv PAGER moresetenv PAGER less (adviced)Missing filename ("less --help" for help)using PAGER moreenvironment variable EDITOR not set, please set this variable to specify your prefered editor. For example in csh shell type: setenv EDITOR visetenv EDITOR neditsetenv EDITOR xemacsNEdit: Can't open displayusing EDITOR virun: ***WARNING*** (redirected) outputfile:Outinfo/stdout.input.s_initial existed, I renamed it to Outinfo/stdout.input.s_initial.3852doris Inputfiles/input.s_initial >>Outinfo/stdout.input.s_initialPROGRESS: Interpretation inputoptionsfile finished.total cpu: 0 min 0.015 secPROGRESS: Finished initializationPROGRESS: Orbit: interpolation coefficients computed. PROGRESS: Start S_READFILES.WARNING : Found identifier, but not so desired word!, exiting routine...PROGRESS: readvolume finished.WARNING : c4SWSTchange: 1PROGRESS: readleader finished.PROGRESS: readdat (header info) finished.PROGRESS: Finished S_READFILES.total cpu: 0 min 0.078 secPROGRESS: Start S_PORBITS.(getorb: reading arclist E:/insar/zhengjiang/orbit/ERS2/arclist) (getorb: orbit file E:/insar/zhengjiang/orbit/ERS2/ODR.074) PROGRESS: getorb: program finished ok.PROGRESS: Finished S_PORBITS.total cpu: 0 min 0.308 secPROGRESS: Orbit: interpolation coefficients computed. PROGRESS: Start S_CROP.PROGRESS: Start cropping slc data.WARNING : writeslc: length of record 1 = "16012"; expected "19612" for ERS SLC (CEOS, full scene).PROGRESS: WRITESLC: 0%PROGRESS: WRITESLC: 10%PROGRESS: WRITESLC: 20%PROGRESS: WRITESLC: 30%PROGRESS: WRITESLC: 40%PROGRESS: WRITESLC: 50%PROGRESS: WRITESLC: 60%PROGRESS: WRITESLC: 70%PROGRESS: WRITESLC: 80%PROGRESS: WRITESLC: 90%PROGRESS: WRITESLC: 100%PROGRESS: Finished S_CROP.total cpu: 0 min 15.417 secPROGRESS: calling preview for cropped slavePROGRESS: Start PREVIEW generation.PROGRESS: SUNraster file created of: Outdata/slave.raw (see also file: ./slave_mag.ras.sh)total cpu: 0 min 15.524 sec...tip: The "run" script can be used to generate template input files.--- WARNING SUMMARY ---There were 3 messages:1: Found identifier, but not so desired word!, exiting routine... 2: c4SWSTchange: 13: writeslc: length of record 1 = "16012"; expected "19612" for ERS SLC (CEOS, full scene).Redirected output in file:Outinfo/stdout.input.s_initialview with: run -v2Xiaoyan@hp-47c03a5dcc20/cygdrive/e/insar/test$ @(#)Doris software, $Revision: 3.11 $, $Author: kampes $ cpxfiddle: INFO: Number of output pixels: 500Xiaoyan@hp-47c03a5dcc20/cygdrive/e/insar/test$ run -s3environment variable PAGER not set, please set this variable to specify your prefered viewer. For example in csh shell type: setenv PAGER moresetenv PAGER less (adviced)Missing filename ("less --help" for help)using PAGER moreenvironment variable EDITOR not set, please set this variable to specify your prefered editor. For example in csh shell type: setenv EDITOR visetenv EDITOR neditsetenv EDITOR xemacsNEdit: Can't open displayusing EDITOR vidoris Inputfiles/input.coregistration >>Outinfo/stdout.input.coregistrationPROGRESS: Interpretation inputoptionsfile finished.total cpu: 0 min 0.031 secPROGRESS: Finished initializationPROGRESS: Orbit: interpolation coefficients computed. PROGRESS: Orbit: interpolation coefficients computed. PROGRESS: Start COARSE_ORB.PROGRESS: Coarse precise orbits coregistration finished. PROGRESS: Finished COARSE_ORB.total cpu: 0 min 0.312 secPROGRESS: Start COARSE_CORR.PROGRESS: COARSE_CORR: 0%PROGRESS: COARSE_CORR: 10%PROGRESS: COARSE_CORR: 20%PROGRESS: COARSE_CORR: 30%PROGRESS: COARSE_CORR: 40%PROGRESS: COARSE_CORR: 50%PROGRESS: COARSE_CORR: 60%PROGRESS: COARSE_CORR: 70%PROGRESS: COARSE_CORR: 80%PROGRESS: COARSE_CORR: 90%PROGRESS: COARSE_CORR: 100%PROGRESS: Coarse coregistration based on correlation finished. PROGRESS: Finished COARSE_CORR.total cpu: 0 min 12.796 secPROGRESS: Start FILTAZI (master).PROGRESS: FILTAZI: 0%PROGRESS: Finished FILTAZI for master.total cpu: 0 min 20.061 secPROGRESS: Start FILTAZI (slave).PROGRESS: FILTAZI: 0%PROGRESS: FILTAZI: 10%PROGRESS: FILTAZI: 20%PROGRESS: FILTAZI: 30%PROGRESS: FILTAZI: 40%PROGRESS: FILTAZI: 50%PROGRESS: FILTAZI: 60%PROGRESS: Finished FILTAZI slave.total cpu: 0 min 28.484 secPROGRESS: Start FINE.PROGRESS: FINE: 0%PROGRESS: FINE: 10%PROGRESS: FINE: 20%PROGRESS: FINE: 30%PROGRESS: FINE: 40%PROGRESS: FINE: 50%PROGRESS: FINE: 60%PROGRESS: FINE: 70%PROGRESS: FINE: 80%PROGRESS: FINE: 90%PROGRESS: FINE: 100%PROGRESS: Fine coregistration finished. PROGRESS: Finished FINE.total cpu: 2 min 41.359 sectotal cpu: 2 min 41.453 sec...Circular definition: see Definition, circular.--- WARNING SUMMARY ---There were no messages.Redirected output in file:Outinfo/stdout.input.coregistrationview with: run -v3Xiaoyan@hp-47c03a5dcc20/cygdrive/e/insar/test$ gmtset: Command not found.psxy: Command not found.psxy: Command not found.psxy: Command not found.pstext: Command not found.pstext: Command not found.pstext: Command not found.psxy: Command not found.pstext: Command not found.pstext: Command not found.Xiaoyan@hp-47c03a5dcc20/cygdrive/e/insar/test$ run -s4environment variable PAGER not set, please set this variable to specify your prefered viewer. For example in csh shell type: setenv PAGER moresetenv PAGER less (adviced)Missing filename ("less --help" for help)using PAGER moreenvironment variable EDITOR not set, please set this variable to specify your prefered editor. For example in csh shell type:setenv EDITOR visetenv EDITOR neditsetenv EDITOR xemacsNEdit: Can't open displayusing EDITOR vidoris Inputfiles/input.resample >>Outinfo/stdout.input.resamplePROGRESS: Interpretation inputoptionsfile finished. total cpu: 0 min 0.031 secPROGRESS: Finished initializationPROGRESS: Orbit: interpolation coefficients computed. PROGRESS: Orbit: interpolation coefficients computed. PROGRESS: Start COREGPM.PROGRESS: Start iteration 0PROGRESS: Start iteration 1PROGRESS: Start iteration 2PROGRESS: Start iteration 3PROGRESS: Start iteration 4PROGRESS: Start iteration 5PROGRESS: Start iteration 6PROGRESS: Start iteration 7PROGRESS: Start iteration 8PROGRESS: Start iteration 9PROGRESS: Start iteration 10PROGRESS: Start iteration 11PROGRESS: Start iteration 12PROGRESS: Start iteration 13PROGRESS: Start iteration 14PROGRESS: Start iteration 15PROGRESS: Start iteration 16PROGRESS: Start iteration 17PROGRESS: Start iteration 18PROGRESS: Start iteration 19PROGRESS: Start iteration 20PROGRESS: finished computation of coregistration parameters. PROGRESS: Finished COREGPM.total cpu: 0 min 2.782 secPROGRESS: Start RESAMPLE.PROGRESS: Resample: normalized lookup table created (kernel and axis).PROGRESS: RESAMPLE: 0%gmtset: Command not found.psxy: Command not found.pstext: Command not found. psxy: Command not found. pstext: Command not found. psxy: Command not found.psxy: Command not found.psxy: Command not found.psxy: Command not found.psxy: Command not found. pstext: Command not found. pshistogram: Command not found. pshistogram: Command not found. gv: Unable to open the display. gv: Unable to open the display. gv: Unable to open the display. PROGRESS: RESAMPLE: 10% PROGRESS: RESAMPLE: 20% PROGRESS: RESAMPLE: 30% PROGRESS: RESAMPLE: 40% PROGRESS: RESAMPLE: 50% PROGRESS: RESAMPLE: 60% PROGRESS: RESAMPLE: 70% PROGRESS: RESAMPLE: 80%PROGRESS: RESAMPLE: 90%PROGRESS: Finished RESAMPLE.total cpu: 1 min 18.797 secPROGRESS: calling preview for resampled slavePROGRESS: Start PREVIEW generation.PROGRESS: SUNraster file created of: Outdata/slave.resampled (see also file: ./slave_rs_mag.ras.sh)@(#)Doris software, $Revision: 3.11 $, $Author: kampes $ cpxfiddle: INFO: Number of output pixels: 454total cpu: 1 min 18.95 sec...A Mathematician is a machine for turning coffee into theorems.--- WARNING SUMMARY ---There were no messages.Redirected output in file: Outinfo/stdout.input.resample view with: run -v4Xiaoyan@hp-47c03a5dcc20/cygdrive/e/insar/test$run -s5environment variable PAGER not set, please set this variableto specify your prefered viewer. For example in csh shell type: setenv PAGER moresetenv PAGER less (adviced)Missing filename ("less --help" for help)using PAGER moreenvironment variable EDITOR not set, please set this variable to specify your prefered editor. For example in csh shell type: setenv EDITOR visetenv EDITOR neditsetenv EDITOR xemacsNEdit: Can't open displayusing EDITOR vidoris Inputfiles/input.products >>Outinfo/stdout.input.productsPROGRESS: Interpretation inputoptionsfile finished.total cpu: 0 min 0.047 secPROGRESS: Finished initializationPROGRESS: Orbit: interpolation coefficients computed. PROGRESS: Orbit: interpolation coefficients computed. PROGRESS: Start INTERFERO.PROGRESS: INTERFERO: 0%PROGRESS: finished compinterfero.PROGRESS: Finished INTERFERO.total cpu: 0 min 1.688 secPROGRESS: calling preview for complex interferogram PROGRESS: Start PREVIEW generation.PROGRESS: SUNraster file created of: Outdata/master_slave.cint (see also file: ./interferogram_mag.ras.sh)PROGRESS: Start PREVIEW generation.@(#)Doris software, $Revision: 3.11 $, $Author: kampes $ cpxfiddle: INFO: Number of output pixels: 454PROGRESS: SUNraster file created of: Outdata/master_slave.cint (see also file: ./interferogram_pha.ras.sh)PROGRESS: Start PREVIEW generation.@(#)Doris software, $Revision: 3.11 $, $Author: kampes $ cpxfiddle: INFO: Number of output pixels: 454PROGRESS: SUNraster file created of: Outdata/master_slave.cint (see also file: ./interferogram_mix.ras.sh)PROGRESS: Start COMPREFPHA.WARNING : Deviation quite large. Decrease degree or number of points?@(#)Doris software, $Revision: 3.11 $, $Author: kampes $ cpxfiddle: INFO: Number of output pixels: 454PROGRESS: Finished COMPREFPHA.total cpu: 0 min 2.194 secPROGRESS: Start SUBTRREFPHA.PROGRESS: SUBTRREFPHA: buffer: 1; line: 127 : 1100PROGRESS: Finished SUBTRREFPHA.total cpu: 0 min 3.193 secPROGRESS: calling preview for refpha corrected image PROGRESS: Start PREVIEW generation.PROGRESS: SUNraster file created of: Outdata/master_slave.srp (see also file: ./interferogram_srp_mag.ras.sh) PROGRESS: Start PREVIEW generation.@(#)Doris software, $Revision: 3.11 $, $Author: kampes $ cpxfiddle: INFO: Number of output pixels: 227PROGRESS: SUNraster file created of: Outdata/master_slave.srp (see also file: ./interferogram_srp_pha.ras.sh) PROGRESS: Start PREVIEW generation.@(#)Doris software, $Revision: 3.11 $, $Author: kampes $ cpxfiddle: INFO: Number of output pixels: 227PROGRESS: SUNraster file created of: Outdata/master_slave.srp (see also file: ./interferogram_srp_mix.ras.sh) PROGRESS: Start COHERENCE.PROGRESS: COHERENCE: 0%@(#)Doris software, $Revision: 3.11 $, $Author: kampes $ cpxfiddle: INFO: Number of output pixels: 227PROGRESS: finished compcoherence.PROGRESS: Finished COHERENCE.total cpu: 0 min 12.093 secPROGRESS: calling preview for coherencePROGRESS: Start PREVIEW generation.PROGRESS: Trying to trick cpxfiddle for REAL4 format (should work) PROGRESS: SUNraster file created of: Outdata/master_slave.coh (see also file: ./coherence.ras.sh)@(#)Doris software, $Revision: 3.11 $, $Author: kampes $ cpxfiddle: WARNING: using real input, width in non-complex pixels total cpu: 0 min 12.214 sec...This time it will surely run.--- WARNING SUMMARY ---There was 1 message:1: Deviation quite large. Decrease degree or number of points?Redirected output in file: Outinfo/stdout.input.products view with: run -v5Xiaoyan@hp-47c03a5dcc20/cygdrive/e/insar/test$ run -s6environment variable PAGER not set, please set this variable to specify your prefered viewer. For example in csh shell type: setenv PAGER moresetenv PAGER less (adviced)Missing filename ("less --help" for help)using PAGER moreenvironment variable EDITOR not set, please set this variable to specify your prefered editor. For example in csh shell type: setenv EDITOR visetenv EDITOR neditsetenv EDITOR xemacsNEdit: Can't open displayusing EDITOR vidoris Inputfiles/input.filter_unwrap >>Outinfo/stdout.input.filter_unwrapPROGRESS: Interpretation inputoptionsfile finished.total cpu: 0 min 0.03 secPROGRESS: Finished initializationPROGRESS: Orbit: interpolation coefficients computed. PROGRESS: Orbit: interpolation coefficients computed. PROGRESS: Start UNWRAP.PROGRESS: Configuration file for snaphu: snaphu.conf created. PROGRESS: Finished snaphu_unwrap.PROGRESS: Finished UNWRAP.total cpu: 0 min 1.357 secPROGRESS: calling preview for unwrapped interferogram PROGRESS: Start PREVIEW generation.PROGRESS: Trying to trick cpxfiddle for HGT format (should work) PROGRESS: SUNraster file created of: Outdata/master_slave.uint (see also file: ./unwrapped_interferogram.ras.sh)@(#)Doris software, $Revision: 3.11 $, $Author: kampes $total cpu: 0 min 1.462 sec...What we know is not much. What we do not know is immense.--- WARNING SUMMARY ---There were no messages.Redirected output in file:Outinfo/stdout.input.filter_unwrapview with: run -v6Xiaoyan@hp-47c03a5dcc20/cygdrive/e/insar/test$ run -s7environment variable PAGER not set, please set this variable to specify your prefered viewer. For example in csh shell type: setenv PAGER moresetenv PAGER less (adviced)Missing filename ("less --help" for help)using PAGER moreenvironment variable EDITOR not set, please set this variable to specify your prefered editor. For example in csh shell type: setenv EDITOR visetenv EDITOR neditsetenv EDITOR xemacsNEdit: Can't open displayusing EDITOR vidoris Inputfiles/input.s2h_geocode >>Outinfo/stdout.input.s2h_geocodePROGRESS: Interpretation inputoptionsfile finished.total cpu: 0 min 0.016 secPROGRESS: Finished initializationPROGRESS: Orbit: interpolation coefficients computed. PROGRESS: Orbit: interpolation coefficients computed. PROGRESS: Start SLANT2H.PROGRESS: S2H: schwabisch: STEP1: compute reference phase for Nheights.PROGRESS: S2H: schwabisch: STEP2: estimate coefficients 1d polynomial.PROGRESS: S2H: schwabisch: STEP3: estimate coefficients for 2d polynomial.PROGRESS: S2H: schwabisch: STEP4: compute height for all pixels. PROGRESS: SLANT2H: 0%PROGRESS: finished slant2hschwabisch.PROGRESS: Finished SLANT2H.total cpu: 0 min 1.172 secPROGRESS: calling preview for slant to heightPROGRESS: Start PREVIEW generation.PROGRESS: Trying to trick cpxfiddle for REAL4 format (should work) PROGRESS: SUNraster file created of: Outdata/heights.raw (see also file: ./height_in_radarcoordinates.ras.sh) PROGRESS: Start GEOCODE.PROGRESS: GEOCODE: 0%@(#)Doris software, $Revision: 3.11 $, $Author: kampes $ cpxfiddle: WARNING: using real input, width in non-complex pixels PROGRESS: finished geocode.PROGRESS: Finished GEOCODE.total cpu: 0 min 12.934 sectotal cpu: 0 min 12.996 sec...On a Korean kitchen knife:"WARNING KEEP OUT OF CHILDREN."--- WARNING SUMMARY ---There were no messages.Redirected output in file:Outinfo/stdout.input.s2h_geocodeview with: run -v7Xiaoyan@hp-47c03a5dcc20/cygdrive/e/insar/test$。
doris使用手册3
六、 input.products 批处理PROCESS interfero计算得到干涉图,输出文件是master_slave.cint ,输出格式是complex_real4。
这里干涉图的定义为:**I M S R =⋅⋅*表示取共扼,.表示逐点相乘,I 表示干涉图,M 表示主图像,S 表示副图像,R 表示相位(振幅为1)。
其干涉图的相位主值表示为2arctan (.)imag real I I φ=这里2arctan 是代表取第四象限的值。
φ表示相位图,I 表示干涉图。
M S R φφφφ=--INT_METHOD OVERSAMPLE在计算干涉条纹数据之前,可以进行过采样,采样因子至少为2。
过采样实际上是做内插,采样因子为2即在原数据之间内插2个点。
这种采样最好在频域内进行,先对原始数据做FFT 变化,然后做过采样,再由逆变换后的数据计算干涉条纹。
INT_OUT_CINT filename需要输入前面输出的.cint 格式的文件。
INT_MULTILOOK 5 1为了对干涉条纹进行显示和评估,在过采样之后进行多视处理。
它实际上是一种压缩显示方法,显示的结果与已有图件形态差异较小一些。
干涉条纹的多视主力中有采取方位向为20,距离向为4的方案,其比例为5:1。
PROCESS comprefphaSAR 图像中任一点可视为雷达波束与地球表面的交点。
它由地面椭球模型、SAR 多普勒方程和卫星与目标点之间的距离方程三个关系式来确定。
这里选用WGS84参考椭球模型作为高度为零的参考基准面,去除由此椭球面引入的平地相位。
S1图3-5 去平地效应示意图所以,干涉场景区域的参考基准面上的每一点都就满足椭球方程:122222=++b P a P P z y x式中,a 、b 分别为椭球长半轴和短半轴。
由SAR 成像原理可知,目标点P 满足距离方程:R S P =-||式中,R=R 0+col*spc ,R 0是图像第一列的卫星的斜距,spc 表示距离向分辨率。
doris使用手册1
一、前言Doris软件全称为Delft Object-oriented Radar Interferometric Software,是荷兰Delft大学的Kampes等人用面向对象的C++语言编写的。
该软件主要开发目的是为研究地表的三维地形及地表变形,为完全免费软件,其源代码是开放的,利于研究者使用。
该软件的主要运行环境为Unix系统,但其在Linux及Cygwin 环境下运行效率也很高。
Doris数据处理大致可以分为四个模块:原始数据读取模块、配准和参考椭球相位计算模块、复相位图像和相干图计算模块、DEM和形变干涉图生成模块。
每个模块中有独立的小模块,在进行数据处理时,每个模块进行独立的计算,处理完后有输出纪录处理的结果及数据,下一模块在执行时则先读取所需的参数及数据,这样每个模块都可以单独进行测试,同时也可以选择不同的算法和软件来实现其中的任意流程,以了解其运算的情形来选择最合适的参数值。
由于Doris 处理的对象是图像矩阵,所以,为了便于再开发和再研究,doris各模块运行产生的中间数据皆为裸矩阵(无头无尾文件)。
Doris的文件命令是由多个批处理文件组成,文件里先规定了要进程(PROCESS)哪些步骤,然后在每个进程(PROCESS)的步骤后,添加执行的关键字(CARD)。
各主要PROCESS如下表所示:上面已经概要地介绍了每个进程的作用,下面的章节将对具体的进程和关键字做介绍并附加具体实验结果。
为了方便起见,本文将一整套Doris工作流程分为七个步骤,每个步骤由一个批处理文件进行。
这七个批处理文件是:1.input.m_initial2.input.s_initial3.input.coregistration4.input.resample5.input.products6.input.filter_unwrap7.input.s2h_geocode以上七个文件有着相同的规格,包括题头、一般性选项、I/O文件和进程四部分。
doris 压缩参数
doris 压缩参数杜瑞斯(Doris)压缩算法是一种广泛应用于数据压缩领域的算法,其优势在于能够在保持高压缩比的同时,实现较快的压缩速度。
为了充分发挥杜瑞斯(Doris)压缩算法的性能,了解和合理设置压缩参数至关重要。
一、杜瑞斯(Doris)压缩算法简介杜瑞斯(Doris)压缩算法是一种基于块的压缩算法,其主要思想是将数据分成若干块,并对每个块进行独立压缩。
与其他压缩算法相比,杜瑞斯(Doris)算法在压缩过程中充分考虑了数据的相关性,从而实现了更高的压缩比和更快的压缩速度。
二、杜瑞斯(Doris)压缩参数概述杜瑞斯(Doris)压缩算法包含多个压缩参数,这些参数可以根据具体需求进行调整,以达到最佳压缩效果。
主要压缩参数如下:1.块大小:设定数据分块的大小,通常根据数据特点和压缩设备性能进行选择。
2.预测长度:预测下一个块的数据长度,以调整压缩策略。
3.阈值:设定一个阈值,当数据块长度超过该阈值时,采用更积极的压缩策略。
4.模式选择:根据数据特点,选择合适的压缩模式,如哈夫曼编码、算术编码等。
5.迭代次数:设定压缩过程中的最大迭代次数,以平衡压缩比和压缩速度。
三、杜瑞斯(Doris)压缩参数详解1.块大小:合理设置块大小可以提高压缩效果。
一般来说,块越大,压缩比越高,但压缩速度可能降低;块越小,压缩速度越高,但压缩比可能降低。
可以根据实际需求和设备性能进行调整。
2.预测长度:预测长度的设置会影响压缩算法的预判能力。
预测长度越长,压缩效果越好,但计算复杂度也会增加。
通常根据数据特点进行选择。
3.阈值:合理设置阈值可以提高压缩效率。
当数据块长度超过阈值时,说明数据具有一定的冗余度,可以采用更积极的压缩策略。
4.模式选择:根据数据特点选择合适的压缩模式。
哈夫曼编码适用于频度较高的数据,算术编码适用于频度较低的数据。
5.迭代次数:迭代次数会影响压缩效果和速度。
一般来说,迭代次数越多,压缩效果越好,但压缩速度会降低。
doris监控指标及阈值
doris监控指标及阈值摘要:一、引言二、doris 监控指标概述1.CPU 使用率2.内存使用率3.磁盘I/O4.网络I/O5.连接数6.线程池使用情况三、doris 监控阈值设置1.正常阈值范围2.警告阈值3.严重阈值四、doris 监控指标处理建议1.CPU 使用率过高2.内存使用率过高3.磁盘I/O过高4.网络I/O过高5.连接数过多6.线程池使用情况不佳五、总结正文:一、引言Doris 是一款高性能的分布式关系型数据库,为了确保其稳定运行,我们需要密切关注一些关键性能指标。
本文将详细介绍Doris 的监控指标及阈值设置,并提供相应的处理建议。
二、doris 监控指标概述Doris 提供了丰富的监控指标,主要包括以下几个方面:1.CPU 使用率:衡量数据库服务器处理能力的指标,一般情况下,我们希望其使用率在50%-70% 之间。
2.内存使用率:反映数据库服务器内存资源的使用情况,理想情况下应保持在一个较低水平,如60% 左右。
3.磁盘I/O:衡量数据库服务器磁盘读写性能,过高可能会导致数据读写延迟。
正常情况下,其值应保持在较低水平。
4.网络I/O:反映数据库服务器网络带宽的使用情况,过高可能会导致数据传输速度下降。
理想情况下,其值应保持在一个较低水平。
5.连接数:表示数据库服务器当前正在处理的连接数量,过高可能会导致服务器压力过大。
一般情况下,我们希望其值控制在1000-2000 之间。
6.线程池使用情况:线程池是Doris 用于处理并发任务的机制,如果线程池使用率过高,可能会导致任务处理速度下降。
理想情况下,线程池使用率应保持在70%-80% 之间。
三、doris 监控阈值设置针对上述监控指标,我们需要设置合适的阈值,以便在出现异常情况时能够及时发现并处理。
以下是一些建议的阈值设置:1.正常阈值范围:一般为指标的平均值±1-2 倍标准差。
2.警告阈值:正常阈值的1.5-2 倍。
doris数据写入流程
doris数据写入流程我们需要创建一个Doris的表,定义表的结构和数据类型。
在创建表的过程中,我们可以指定表的分区方式、排序方式等属性,以便更好地组织和存储数据。
接下来,我们需要准备要写入的数据。
可以将数据存储在本地文件或者其他数据库中,然后通过导入工具将数据导入到Doris中。
导入工具提供了多种导入方式,如从本地文件导入、从其他数据库导入等。
我们可以根据实际情况选择合适的导入方式。
在数据导入之前,我们需要确保数据的格式和类型与表的定义相匹配。
如果数据不符合表的定义,可能会导致写入失败或者数据被截断。
因此,在导入数据之前,我们需要进行数据的预处理,将数据转换为符合Doris表结构的格式。
接下来,我们可以使用Doris提供的写入API将数据写入到表中。
写入API提供了多种方式,如单条写入、批量写入等。
我们可以根据实际情况选择合适的写入方式。
在写入数据时,我们需要注意设置合适的写入参数,如写入模式、写入超时时间等。
在写入数据的过程中,Doris会对数据进行压缩和编码,以减少数据的存储空间和提高查询性能。
同时,Doris还会进行数据的分片和分布式存储,以实现数据的高可用和高并发访问。
数据写入完成后,我们可以通过Doris提供的查询API对数据进行查询和分析。
查询API提供了多种查询方式,如单表查询、多表关联查询等。
我们可以根据实际需求进行查询,并通过结果进行统计和分析。
除了使用API进行数据写入,我们还可以使用Doris提供的可视化工具进行数据写入。
可视化工具提供了图形界面,可以方便地进行数据导入、数据预览和数据写入等操作。
通过可视化工具,我们可以更加直观地了解数据写入的过程和结果。
在数据写入过程中,我们还需要注意数据的一致性和完整性。
Doris 提供了事务支持,可以保证数据的一致性和完整性。
在写入数据时,我们可以使用事务机制来保证多个写操作的原子性和隔离性。
总结起来,Doris数据写入流程包括创建表、准备数据、导入数据、写入数据和查询数据等步骤。
- 1、下载文档前请自行甄别文档内容的完整性,平台不提供额外的编辑、内容补充、找答案等附加服务。
- 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
- 3、如文档侵犯您的权益,请联系客服反馈,我们会尽快为您处理(人工客服工作时间:9:00-18:30)。
1,Doris处理可以分步处理,也可以整体处理。
Doris处理要先生成模板文件,doris每步处理都对应一个模板文件,这些模板文件包含了doris每步处理所需要的参数,因此用户可以直接打开模板文件,用以编辑更改然后保存,也可以在linux终端中采用Vi命令,或nedit(相当于window下的写字板)来进行编辑保存。
2,查看版本信息,终端键入:doris 即可查看。
3,生成模板。
命令:run –g
说明:会在当前工作目录下生成模板文件,它包含三个文件夹目录:Inputfiles,Outdata ,Outinfo。
Outdata用于存放每步生成的结果数据文件。
Outinfo用于存放每步处理的log文件,记录中间处理过程。
Inputfiles中存放的就是模板文件。
它包含的文件:
对这些模板进行编辑时,文件中每行前面加c和空格,表示程序处理时不读取该行参数,否则都会执行到的,编辑参数选项时直接更改后面的的数值即可。
4,doris处理步骤总体介绍。
Doris的INSAR处理共分八步:每步又可以分为许多小步,可以通过编辑来选择执行或不执行某些小步骤(小步骤在模板文件中的标识符为process)
第一步run –s1:读入主影像的头文件及数据体信息以及轨道信息(可以使用粗轨),其对应的模板文件 input.m_initial.,这步包含:读取主影像头文件信息,读取主影像数据体信息,主影像过采样,读取主影像精轨文件。
(如下:加c 表示屏蔽该行,以下都同)
第二步run –s2:读入辅影像的头文件及数据体信息以及轨道信息(可以使用粗轨),其对应的模板文件:input.s_initial。
这步包含:读取辅影像头文件信息,读取辅影像数据体信息,辅影像过采样,读取辅影像精轨文件。
第三步run –s3:影像配准。
对应的模板文件:input.coregistration.包含步骤有:一级轨道级配准,二级像素级粗配准,三级子像素级精配准,主影像方位向滤波,辅影像方位向滤波,以及影像距离向滤波。
第四步run –s4:影像重采样。
对应的模板文件:input.resample。
包含步骤有:估计配准多项式系数,影像重采样,以及距离向滤波。
第五步run –s5:生成干涉图。
对应的模板文件:input.products.。
包含步骤有:干涉图生成,除平地效应,以及生成相干系数图。
其中去平地效应有两种方式:
一种是基于轨道去平地效应(如果在第一二步提供精轨,那么就是基于精轨去平地效应,否则是基于粗轨去平地效应)。
其中comprefpha步骤为用轨道来计算平地相位。
Subtrrefpha 步骤为减去平地相位。
第二种是基于粗DEM去平地效应。
其中comprefdem 为利用提供的粗DEM来反算求得平地相位,subtrrefdem为减去由粗DEM反算得到的平地相位。
第六步run –s6:干涉图滤波及解缠。
对应的模板文件:input.filter_unwrap.包含的步骤有:干涉图滤波,相位解缠,以及差分干涉处理(针对Dinsar)
注意的是doris没有提供强有力的解缠程序,它只提供一种简单的枝切法解缠,这对于复杂地形下的干涉图解缠会存在很大误差,因此它一般调用比较有名的snaphu解缠算法来完成解缠过程。
Snaphu解缠也可以单独拿出来进行处理,
第七步run –s7:DEM生成。
对应的模板文件:input.s2h_geocode. 。
它包含的步骤有:相位高程转换,地理编码
第八步run –s9(整体的步骤包括第三至第七):doris分步处理前七步即完成,这里的第八步实际上是整体处理,它直接包含第三至第七步的全过程,对参数设置进行了简化,它对应的模板设置文件:input.quicklook.。
该模板文件是将前三到七步模板文件中认为对insar处理比较关键的参数进行了一个集中设置。
因此若要进行整体处理时,必须也先要运行第一步和第二步,然后再运行第八步,即可。
1,命令介绍.
生成模板:run –g
编辑模板:
编辑第一步模板:run –e1
编辑第二步模板:run –e2
编辑第三步模板:run –e3
编辑第四步模板:run –e4
编辑第五步模板:run –e5
编辑第六步模板:run –e6
编辑第七步模板:run –e7
编辑第八步模板:run –e9
(每行前加 C 表示屏蔽该行。
打开模板文件后,按“i“表示可以对进行插入操作,如果删去时,直接用delete或后退键。
按ESC键时则退出当前的编辑状态。
若要退出文件时,则先按ESC,然后:q!,或:wq,其中:q!表示退出不保存修改,:wq为保存修改。
)运行处理:
运行第一步:run –s1
运行第二步:run –s2
运行第三步:run –s3
运行第四步:run –s4
运行第五步:run –s5
运行第六步:run –s6
运行第七步:run –s7
运行第八步:run –s9
2,参数解说及设置。
2.1,第一步:读取主影像信息。
(1)需要处理步骤。
必需:m_readfiles(读取主影像头文件信息).
必需:m_crop(读取主影像数据体信息).
不需:m_ovs(主影像过采样步骤):作干涉处理时,一般不运行。
可选:m_prorbits(读取主影像精轨文件):若要使用提供的精密轨道数据时,就执行该行,否则直接C掉。
(2)小步骤参数说明及设置。
############################################################
INSAR通用参数设置,每步都有的。
###########################################################
名称:LOGFILE
用处:用以记载DORIS处理过程中所有信息。
文件名:doris.log.
名称:M_RESFILE
用处:用以记载主影像文件所有信息。