使用cpxfiddle画图，例如phase，magnitude等
2024-05-30 15:07:56
1.cpxfiddle使用说明：
输入复数数据转为标准输出，输入可以为强度，相位，实部或者虚部。输出可以多视，裁剪，镜像，比例处理等。
函数说明：
-w 线段长度（宽度，方位向像元数，x方向）-e [1.0] 指数： 输出=比例*输出^指数-s [1.0] 比例： 输出 = 比例 * 输出 ^ 指数-l [1] 第一（方位向azimuth Y）line-L [all] 最后（方位向azimuth Y）line-p [1] 第一（距离向range X）pixel-P [all] 最后（距离向range X）pixel-M [1/1] 多视指数在X/Y方向（range/azuimth）-S [1/1] 裁剪指数在X/Y方向（range/azuimth）-q [mag] 输出：NORMAL | MAG | PHASE | MIXED | REAL | IMAGnormal = (real,imag)magnitude = sqrt(real^ 2 + imag^2)phase = atan2(imag,real)mixd = phase overlay, only with -o sunrasterreal = line [2*j],imag = line [2*j +1]-f [cr4] Input format identifier:CC1 | CUC1 | CI2 | CI4 | CR4 | CR8for complex 2x1B char, unsigned char, short integer, integer, float, double(major row order pixel interleaved complex data.)-o [ascii] Output format identifier (to stdout):UCHAR | SUNRASTER | FLOAT | SHORT | ASCIIuchar, short, sunraster, and float optionswrite binary to stdout!-c [gray] Colormap for "-o sunraster" option.FILENAME | GRAY | JET | HOT | COOL | BERT | HSVIf cmap is a filename, then a ascii file with256 lines, containing r g b [0:255] is assumed.If cmap is a identifier, this colormap is used.-m code Flag to mirror file content in X or Y direction:X | Y | XY | YXIf mirroring Y then the first output line isline number P (default last in file). This meansthat if -Y is specified, the last line notnecessarely is line -l (firstline default 1)The same is true for mirroring in X direction.-r [rmin/rmax] uses given minimum and maximum range on data as scalling parameters. Basically, it skips data sampling step for statistics (min,max,mean) computationand defines (min,max,mean) values from given parameters.It only effects magnitude, normal and mixed mode outputs.For example, use it to scale coherence maps to 0-1 range.Tip: no need to use --ignorenan when you are using -r option.--ignorenan ignores nan values during calculation of min, max and mean for scaling of sunraster outputs. It only effects magnitude, normal and mixed mode output types.Additionally, when binary output is done for magnitude images 0 is set instead of nan values. --fullstat compute min, max and mean evaluating every pixel for scaling of sunraster outputs. -B b|s Swap bytes.if -Bs then swap to host order.(call ntohs or ntohl, thus swapping from big endianto host type, which is small endian for X86 PC's.)if -Bb then swap to network order.(call htons or htonl, thus swapping to big endianfrom host type, which is small endian for X86 PC's.)On a big endian platform both these functions aredefined as NULL macros, and do nothing.-b Add a scalebar in lower right corner for -o sunrasterOnly for -q phase and mixed.-H bytes Skip Header bytes, e.g., for GENESIS SUNraster floats.-V Verbose.-h[elp] This help.
2.函数自带案例：
Check first few values in complex short file (width 100, height 200):cpxfiddle -w100 -fci2 -qnormal -oascii -L5 -P3 -V -- cpxfileTo generate a grd file for use with GMT, without creating a largetmpfile, use something like:cpxfiddle -w100 -fci2 -qphase -ofloat -- cpxfile | \xyz2grd -Gfile.grd -R1/100/1/200 -ZfTo crop a binary complex float file of width 781 and height 501:(redirect, or your screen will be flooded with binary output!)cpxfiddle -w781 -fcr4 -qnormal -ofloat -l101 -L200 \-p51 -P750 cpxfile > cpxfile.croppedTo convert the phase of the same file to SUNraster format:cpxfiddle -w781 -fcr4 -qphase -osunraster \-ccmap.rainbow -- cpxfile > phase.rasdisplay phase.rasWhere the colortable file is generated by (using GMT):makecpt -Crainbow -T0/256/1 | \awk '{if (NR>3&&NR<260){print $2, $3, $4}}' > cmap.rainbowTo create a sort of 8 bit pixmap, visualize with ImageMagick'sdisplay as gray format:cpxfiddle -qphase -fcr4 -ouchar -w1 -- cpxfile > phase.ucdisplay -size 781x501 gray:phase.ucThe file phase.uc could also be converted to another format, e.g.:rawtopgm 781 501 phase.uc > phase.pgmxv phase.pgmBut easier would have been to use cpxfiddle to generate theSUNraster file. This can be done by:cpxfiddle -w100 -qmag -osunraster -e0.3 -s1.1 file.cr4 > file.rasWhere file.cr4 is complex real4, -e is used to equalize the histogram, and -s to threshold the data (larger scale meansmore white). To overlay the phase as layer over the magnitude:cpxfiddle -w100 -qmixed -osunraster -e0.3 -- file.cr4 > file.rasxv file.rasTo read the header of a (non complex!) sunraster file (8 int32 values):cpxfiddle -w1 -L4 -qnormal -fci4 -- file.rasThis program can be tricked to crop float files (non complex)This program can be used for conversions of files, trick with -w1To crop a bigendian file on a linux PC, swap and crop like:cpxfiddle -w100 -fci2 -qnormal -oshort -Bs -- cpxfile > newfileTo read the content of a (non complex!) float vector with 32 byte headercpxfiddle -w1 -qnormal -fr4 -H32 -oascii -- file.rasTo scale the content of a float file: coherence map to [ 0 1 ] range.cpxfiddle -w 1000 -q normal -o sunraster -b -c gray -f r4 cohmap.r4 -r 0,1 > cohmap.ras BUGS:-M with -m: First mirrored then multilooked.-B: Not tested properly, but simply calls functions.all: If cpxfile is a symbolic link, filesize cannot bedetermined, and therefor the number of lines.Solution: Use the full path or use a hard link.(lot of options, binary to screen, but these are features.)SEE ALSO:cpx2ps, display, rawtopgm, makecpt (GMT)Please send your comments to: TUDelft, doris_users@tudelft.nl(BTW: view this help by piping stderr: cpxfiddle -h |& more in cshor cpxfiddle -h 2>&1 | less in bash.)
3. 实例：
1. 画相关性图
cpxfiddle -w 5000 -q normal -o sunraster -b -c gray -M 2/2 -r 0.0/1.0 -f r4 -l1 -p1 -P5000 ccoherence.raw > ccoherence.ras
2. 画高度图
cpxfiddle -w 5000 -q normal -o sunraster -b -c jet -M 2/2 -f r4 -l1 -p1 -P5000 dem_radar.raw > comprefdem_height.ras
3. 画相位图
cpxfiddle -w 5000 -q normal -o sunraster -b -c jet -M 2/2 -f r4 -l1 -p1 -P5000 refdem.raw > comprefdem_phase.ras
或者
cpxfiddle -w 5000 -q phase -o sunraster -b -c jet -M 2/2 -f cr4 -l1 -p1 -P5000 cint.0.5 filtered > interferogram_filt_pha.ras
4. 画强度图
cpxfiddle -w 5000 -e 0.3 -s 1.0 -q mag -o sunraster -b -c gray -M 2/2 -f cr4 -l1 -p1 -P5000 cint.0.5filtered > interferogram_filt_mag.ras
5. 画强度和相位混合图
cpxfiddle -w 5000 -e 0.3 -s 1.2 -q mixed -o sunraster -b -c hsv -M 2/2 -f cr4 -l1 -p1 -P5000 cint.0.5filtered > interferogram_filt_mix.ras
使用cpxfiddle画图，例如phase，magnitude等相关推荐
- 用Python做信号处理
用Python做信号处理 声明:本文中设计的知识和代码大部分来自:芥末的无奈的博客_CSDN博客-音频处理,c++,keras领域博主 以及 凌逆战 - 博客园 (cnblogs.com) 两位大神所 ...
- 【开源电机驱动】H 桥驱动-软件篇
原文地址:http://www.modularcircuits.com/blog/articles/h-bridge-secrets/h-bridge-control/ 本文为作者翻译校正稿件,含个人 ...
- MorphSVGPlugin from GreenSock 的源码注释分析
不管是ActionScript(flash)还是JavaScript,GreenSock的GSAP几乎是公认的最强动画库. 常用的就是其中的TweenLite.TweenMax.TimelineLit ...
- 深度学习中的语音信号处理基础
文章目录 音频处理流程 常用谱:幅度谱.梅尔谱 时域 --> 频域 分帧 窗长 帧移 语音信号特征获取流程 梅尔谱 使用 librosa 提取梅尔谱 使用 tacotron 获取梅尔谱(推荐) ...
- 相位变化对于真实信号的影响是什么？
一个信号经过傅里叶分解后变成Phase+Magnitude两部分 形象的说, Phase代表信号的形状,Magnitude代表对应形状的大小比例 一维信号的确不太直观,我们来看二维的信号(图像)吧,上 ...
- python rgb转hsv_Python colors.hsv_to_rgb方法代码示例
本文整理汇总了Python中matplotlib.colors.hsv_to_rgb方法的典型用法代码示例.如果您正苦于以下问题:Python colors.hsv_to_rgb方法的具体用法?Pyt ...
- matlab 画图像频谱图
matlab 画图像频谱图 需要注意的地方 图像格式 代码 处理单个图像 批量处理 需要注意的地方 图像格式 JPG: 有损压缩,24bit真彩色,不支持动画,不支持透明色. PNG 无损压缩,PNG ...
- python matplotlib画图的几个实例--latex，坐标系等
文章目录 实例1 学会使用tex/latex 实例2 学会画坐标轴 2.1过程 2.2 典型例子 2.2.1 一条带箭头的竖线 2.2.2 坐标系 2.2.3 坐标系上画三角函数 实例1 学会使用te ...
- MATLAB中plot()画图的颜色线型和希腊字母参数设置
转载自:https://www.cnblogs.com/sddai/p/5399462.html y 黄色 · 点线 m 粉 ...
最新文章
- OSS.Core基于Dapper封装(表达式解析+Emit)仓储层的构思及实现
- Android各层推荐开发书籍及参考资料
- mysql --The MEMORY Storage Engine--官方文档
- 怎么样得到平台相关的换行符？
- 小熊的人生回忆（三）
- canvas综合应用绘制哆啦A梦
- 上海消保委评饿了么“多等5分钟”功能：逻辑上有问题
- C# 泛型2---排序
- 基于matlab的am调制与仿真,基于MATLAB的AM调制及解调系统仿真分解
- 明日复明日,明日何其多,我生待明日,万事成蹉跎(分享)
- 1320. 拯救奶牛
- 初步学习网站分析工具
- 【题目集02丨PTA】PTA基础编程题目集（自己做的）
- 好用的电视盒子软件推荐：无广告看电视我选这两款
- 数据库介绍与数据库安装
- QT获取微秒级时间戳
- 计算机检测维修与数据恢复招标,计算机检测维修与数据恢复实训室企业招标文件.doc...
- XiaoHu日志 7/29~8/9
- 捻深红尚透,谁信道、花亦通灵。多态
- 浅谈：基于物联网的固定资产管理技术