haoNote

原创 用python生成动图

【代码】用python生成透明背景动图。

转载 矢量算符运算法则

来源：

原创 2023中国科学院文献情报中心期刊分区表（GEOCHEMISTRY & GEOPHYSICS 地球化学与地球物理）

原创 MATLAB计算一元二次方程ax^2+bx+c=d并保留5位小数

【代码】MATLAB计算一元二次方程ax^2+bx+c=d并保留5位小数。

原创 zip分卷压缩与解压

【代码】zip分卷压缩与解压。

原创 Madagascar v3.2 绘图: shell脚本示例

【代码】Madagascar v3.2 绘图shell脚本示例。

原创 MATLAB绘制波数谱

【代码】MATLAB绘制波数谱。

原创 Linux解压rar文件

【代码】Linux解压rar文件。

原创 使用SU画波数谱

suaddhead < pp.dat n1=234 n2=663 |sushw key=dt a=1000| suspeck1k2 | suximage cmap=rgb1 legend=1 perc=99

原创 使用SU做波数域滤波

# 滤出低波数suaddhead < pp_lap_pow.dat n1=234 n2=663 | sukfilter k=0.,0.3,0.7,1 amps=1,1,0,0 verbose=1 |sushw key=dt a=1000 | sustrip> pp_lap_pow_lk.dat# 滤出中波数suaddhead < pp_lap_pow.dat n1=234 n2=663 | sukfilter k=0.,0.3,0.7,1 amps=0,1,1,0 verbose=

原创 用SU保存矢量图

用SU保存矢量图img=$1psimage < ${img}.dat n1=240 \label1="Depth(km)" label2="Trace" labelfont=Times-Roman labelsize=24 \f1=0 d1=0.055 f2=0 d2=1 perc=99 \legend=0 lstyle=vertright lbeg=-0.1 lend=0.1 ldnum=0.05 lgrid=dot \height=6.0 width=10.0 \verbose=

原创 SU应用场景3：地震记录抽稀

# 从in.su每隔20道抽出一道存到out.susuwind < in.su key=tracl j=20 >out.su

原创 SU应用场景2：地震记录重采样

SU应用场景2：地震记录插值问题描述命令问题描述现有地震记录640道(无道头)，每道256个采样点，采样间隔5ms.现在需要将采样点插值成1280个采样点命令suaddhead < input.dat ns=256 | sushw key=dt a=5000 >output1.susuresamp <output1.su nt=1280 dt=1000 >output2.susustrip < output2.su >output2.dat...

原创 SU应用场景1：从多个二进制炮记录各抽一道拼在一起最后画频谱

suaddhead < file1.dat ns=100 | suwind dt=0.001 key=tracl min=1 max=1 >f1.susuaddhead < file2.dat ns=100 | suwind dt=0.001 key=tracl min=2 max=2 >f2.susuaddhead < file3.dat ns=100 | suwind dt=0.001 key=tracl min=3 max=3 >f3.sucat f1.su

原创 使用Madagascar实现三维地震数据体的动态展示

#!/bin/shplot(){ segyFile=3D.segy segyDir=. imgDir=. title=3D ny= number Of Inline nx= number Of Crossline crl0= the first Crossline to show inl0= the first Inline to show gain=98 inlshow= the default inline to show

原创 使用Madagascar绘制剖面图

目录1.查看要绘制的segy文件的总道数，以确定绘制范围2.快速预览从第a道起的b道剖面（注意:a+b1.查看要绘制的segy文件的总道数，以确定绘制范围sfsegyread < 文件名.segy endian=1 | sfinn1 : 每一道的采样点数n2 : 总道数2.快速预览从第a道起的b道剖面（注意:a+b<n2）sfsegyread < 文件名.segy endian=1 | sfput d1=0.01 unit1=km label1=Depth title= 标题

原创 基于Ubuntu16.04的Anconda3安装及虚拟环境配置（python3.7 + cuda10.1 + cuDNN7.6 + tensorflow2.1.0）

一、安装Anconda在清华大学开源镜像站或者Anconda官网选择适合你机器的版本下载。我这里下载的安装脚本文件名为Anaconda3-2020.11-Linux-x86_64.sh运行安装脚本sh Anaconda3-2020.11-Linux-x86_64.sh输入配置环境变量sudo gedit /etc/profile打开文本五后在最后一行加上：export PATH=~/anaconda3/bin:$PATH一、查看显卡信息lspci | grep -i vga我手上

原创 使用python3读取segy文件道头信息

# code for read headers of each traceimport matplotlib.pyplot as mpimport numpy as np import sysimport structimport binasciifileName = "BAOFU_amp.segy"nTrace = 4nSample = 3001fSegy = open(fileName,"rb")fSegyHeader = {'tracr':[5,8] ,'ep':[17,2

原创 使用python3读IBM格式segy文件并绘制变密度剖面图

# code for read IBM segy format seismic data files import matplotlib.pyplot as mpimport numpy as np import sysimport structimport binasciidef ibm2ieee(ibm): if ibm == 0: return 0.0 sign = ibm >> 31 & 0x01 exponent = ibm

原创 使用python3读IEEE格式segy文件并绘制变密度剖面图

# code for read IEEE segy format seismic data files import matplotlib.pyplot as mpimport numpy as np import sysimport structimport binasciifileName = "2267x1401.segy"nTrace = 2207nSample = 1401fSegy = open(fileName,"rb")data = np.zeros((nSample,

原创 判断segy文件的编码格式是IBM还是IEEE

# code for judge a segy file IEEE or IBM format = 'nothing'fileName = "vp_marmousi-ii.segy"fSegy = open(fileName,"rb")fSegy.seek(3225,0) tempValue = fSegy.read(1)hexValue = tempValue.hex()decValue = int(hexValue,16)if decValue == 1: format = 'I

原创 如何基于ubuntu16.04完整安装Madagascar3.1.1

基于ubbuntu16.04安装Madagascar3.1.11.准备库2.下载安装包3.解压安装包3.安装4.测试1.准备库在你的终端执行以下命令：sudo apt-get install libxaw7-devsudo apt-get install libnetpbm10-devsudo apt-get install libtiff5-devsudo apt-get install libgd-devsudo apt-get install libplplot-devsudo apt

转载 Linux操作系统：如何实现vim的跨文件复制粘贴

用vim打开一个文件，例如：a.cpp在普通模式下，输入：":sp"（不含引号）横向切分一个窗口，或者":vsp"纵向切分一个窗口，敲入命令后，你将看到两个窗口打开的是同一个文件在普通模式下，输入：":e b.cpp"，在其中一个窗口里打开另一个文件切换到含有源文件（a.cpp）的窗口，在普通模式下，把光标移到你需要复制内容的起始行，然后输入你想复制的行的数量（从光标所在行往下计算），在行数后面接着输入yy，这样就将内容复制到临时寄存器里了（在 普通模式下ctrl+w，再按一下w，可以在两个窗口之间.

原创 Linux操作系统：如何解决cannot open shared object file: No such file or directory一类问题

最近在linux操作系统运行程序时出现了libiomp5.so: cannot open shared object file: No such file or directory这种报错，比较简单的解决步骤是：定位所需库文件的位置loacte libiomp5.so2.拷贝该库文件如果能在本地找到该库文件，则将其拷贝到运行程序的路径如果本地没有该库文件，则从其他计算机拷贝...

原创 用SeismicUnix（SU）读segy文件画地震剖面

1、读.segy生成.susegyread tape=文件名.segy verbose=1 ns=采样点数 conv=1>文件名.su2、用.su画变密度剖面图（红白蓝）suximage < 文件名.su legend=1 d2=tr.d2 d1=tr.d1 title=“图名” cmap=rgb1 &3.自定义时窗时：suximage < 文件名.su legend=1 d2=tr.d2 d1=tr.d1 title=“图名” cmap=rgb1 x1beg=纵坐标起点

原创 c++ 输出保留两位小数怎样实现

题目描述：班上有学生若干名，给出每名学生的年龄（整数），求班上所有学生的平均年龄，保留到小数点后两位。输入描述:第一行有一个整数n（1<= n <= 100），表示学生的人数。其后n行每行有1个整数，取值为15到25。输出描述:可能有多组测试数据，对于每组数据，输出一行，该行包含一个浮点数，为要求的平均年龄，保留到小数点后两位。示例1输入21817输出17....

原创 STL—map的使用

参考资料题目描述：给出一个01字符串（长度不超过100），求其每一个子串出现的次数。输入描述:输入包含多行，每行一个字符串。输出描述:对每个字符串，输出它所有出现次数在1次以上的子串和这个子串出现的次数，输出按字典序排序。示例1：输入10101输出0 201 21 310 2101 2#include <iostream>#include <...

原创 吃糖果

题目描述:名名的妈妈从外地出差回来，带了一盒好吃又精美的巧克力给名名（盒内共有 N 块巧克力，20 > N >0）。 妈妈告诉名名每天可以吃一块或者两块巧克力。 假设名名每天都吃巧克力，问名名共有多少种不同的吃完巧克力的方案。 例如： 如果N=1，则名名第1天就吃掉它，共有1种方案； 如果N=2，则名名可以第1天吃1块，第2天吃1块，也可以第1天吃2块，共有2种方案； 如果N=3，则...

原创 c++ 自学笔记：多态

c++ 自学笔记：多态1.什么是多态2.动态多态的原理1.什么是多态C++面向对象三大特征：封装、继承、多态多态分两类：静态多态：包括 函数重载、运算符重载，复用函数名动态多态：包括 派生类、虚函数实现运行时多态静态多态与动态多态区别：静态多态的函数地址早绑定——编译阶段确定函数地址动态多态的函数地址晚绑定——运行阶段确定函数地址动态多态满足条件：有继承关系子类重写...

原创 子数组异或查询

有一个正整数数组 arr，现给你一个对应的查询数组 queries，其中 queries[i] = [Li, Ri]。对于每个查询 i，请你计算从 Li 到 Ri 的 XOR 值（即 arr[Li] xor arr[Li+1] xor … xor arr[Ri]）作为本次查询的结果。并返回一个包含给定查询 queries 所有结果的数组。示例 1：输入：arr = [1,3,4,8], qu...

原创 解码字母到整数映射

给你一个字符串 s，它由数字（‘0’ - ‘9’）和 ‘#’ 组成。我们希望按下述规则将 s 映射为一些小写英文字符：字符（‘a’ - ‘i’）分别用（‘1’ - ‘9’）表示。字符（‘j’ - ‘z’）分别用（‘10#’ - ‘26#’）表示。返回映射之后形成的新字符串。题目数据保证映射始终唯一。示例 1：输入：10#11#12输出：jkab示例 2：输入：1326#输出：acz...

原创 内存分区模型

内存分区模型1.总览2.程序运行前3.程序运行后4.new操作符1.总览代码区：存放函数体的二进制代码，由操作系统进行管理的区域全局区：存放全局变量和静态变量以及常量的区域栈区 ：由编译器自动分配释放，存放函数的参数值、局部变量等堆区：由程序员分配和释放，若程序员不释放，程序结束时由操作系统回收不同区域存放的数据，生命周期不同2.程序运行前在程序编译后，生成了可执行程序（如Wi...

原创 C语言的typedef和#define

1.typedef用于给已有的类型名指定一个新的类型名typedef int INTEGER;typedef float REAL;INTEGER i,j; //等价于：int i,j;REAL x,y; //等价于：float x,y;/*========================================================*/typedef st...

原创 共用体、枚举类型

共用体定义与引用union 共用体名{ 共用体成员}变量表列;共用体名.共用体成员共用体变量所占内存长度等于最长的成员的长度同一内存段可以存放不同类型的成员，但在每一瞬间只能存放其中一种共用体变量起作用的成员是最后一次存放的成员共用体变量的地址和它各个成员的地址是同一个地址枚举类型声明方式enum weekday{sun,mon,tue,wed,thu,...

原创 单链表的建立及其结点的查找、删除、插入

目录一、静态链表二、动态链表1.动态列表相关函数2.动态链表的建立3.动态链表的遍历查找4.动态链表的删除5.动态链表的插入一、静态链表静态链表中，所有结点都是在程序中定义的，不是临时开辟的，也不能用完后释放。例子：#include <stdio.h>struct student{ char *name; float score; struct student *ne...

原创 C指针——详细笔记

目录概述指针变量的定义、引用指针变量作函数参数一维数组与指针一维数组名作函数参数多维数组与指针字符串与指针指向函数的指针结构体与指针指针数组、数组指针指向指针的指针指针运算概述在程序中一般是通过变量名对内存单元进行存取操作的。其实程序经编译后已经将变量名转换为变量的地址，对变量值的存取都是通过地址进行的 。指针是变量的地址，指针变量是只能存放地址的变量。指针变量的定义、引用1.定义指针变...

原创 位运算与位段

1. 位运算总结2. 位段在数据通信领域，控制信息往往只占一个字节中的一个或几个二进制位，常常在一个字节中存放几个信息。c语言允许在一个结构体中以位为单位来指定其成员所占内存长度，这种以位为单位的成员称为“位段”或“位域”。尤其要注意各位段的引用范围【注】1字节=8位1.位段举例struct packed_data{ // 尤其要注意各位段的...

原创 预处理命令

预处理命令不是c语言本身的组成部分，不能直接对它们进行编译。经过预处理后的程序不再包括预处理命令，由编译程序对预处理后的源程序进行通常的编译处理，得到可供执行的目标代码。c提供的预处理功能主要有3种：宏定义文件包含条件编译1.宏定义1.1不带参数的宏定义1.作用：用指定标识符来代替一个字符串。#define 标识符 字符串2.在预编译时将宏名替换成字符串的过程称为“宏展开”3...

原创 C语言文件操作——定位、出错检测

fseek函数1.作用：利用fseek函数通过对文件的位置指针的控制，可以实现对流式文件的随机读写。fseek一般用于二进制文件。2.调用形式fseek(文件指针,位移量,起始点);//位移量：以起始点为基点向前移动的字节数起始点被编码为：0，1，20：文件开始（SEEK_SET）1：当前位置(SEEK_CUR)2：文件末尾(SEEK_END)feof函数1.作用：判断文件是否...

原创 C语言文件操作——开、关、读、写

概述1.一个C文件是一个字节流或二进制流，它把数据看作是一连串的字符（字节），而不考虑记录的界限。在C语言中对文件的存取是以字符（字节）为单位的。（注：流式文件 —— 输出时不会自动增加回车换行符以作为记录结束的标志，输入时不以回车换行符作为记录的间隔）2.ANSI C标准使用缓冲文件系统处理文本文件和二进制文件。每个被使用的文件都在内存中开辟一个区，用来存放文件的有关信息，这些信息是保存在...