SU（seismic Unix）专栏
关注
分享
扫一扫
文章平均质量分 63
Coder802
你的努力，终将成就不可替代的自己。
展开
-
数组在内存中存储方式
SU（Seismic Unix）中多维（≥2）数组的存储是列优先（columnmajor）的CUDA中多维数组（≥2）数组的 存储是行优先（row major）的用SU开的数组，将数据传给CUDA后还是列优先的。原创 2020-05-23 16:13:47 · 732 阅读 · 0 评论 -
Ubuntu16.04 下convert 命令 将eps转tif,jpg,png,pdf格式
在Ubuntu16.04 中用SU生成的图件一般是eps格式的，我们需要将其转换为tif,jpg,png or pdf等格式。用ImageMagick的convert命令可以实现这些格式转换操作：1、eps转tifconvert -density 600 hs_shot.eps -resize 100% hs_shot.tif2、eps转jpgconvert -density...原创 2019-10-27 16:14:24 · 4366 阅读 · 0 评论 -
Ubuntu16.04 下SU画图，批量和单个
SU画EPS图很方便，也很简单，当单个画图时用下面的命令： supsimage < shot_new/hs_shot$i.su \ d1num=0.5 d2num=100 \ wbox=30 hbox=6 \ threecolor=1 bps=24 perc=97.5 \ wrgb=1.0,0,0 grgb=1.0...原创 2019-10-27 16:00:08 · 573 阅读 · 0 评论 -
基于Madagascar的二维地震声波波动方程正演模拟
最近在将SU写的地震勘探的程序迁移到Madagascar上，初步尝试，写了一个二维声波方程正演程序，很简单，也很基本，只能说出波场快照，没有吸收边界条件，贴出来，供大家参考。代码和脚本如下：#include &lt;time.h&gt;#include "rsf.h"#define FSIZE sizeof(float)static float ricker (float t, floa...原创 2018-12-24 15:03:49 · 3153 阅读 · 3 评论 -
SU编程和Madagascar编程中的函数对比和区别--终端提示函数
之前一直用SU开发地震勘探的程序，在老板的要求下逐步转向Madagascar编程，做程序的迁移。在语言转换过程中，受限于SU的编程习惯，始终想用SU中的一下函数，但是在编译时会会有“warning”提示。所以我将SU和Madagascar中功能相同的函数进行对比分析，以备常用。1、终端信息提示函数在SU中，该函数为warn()，用法swarn(&amp;amp;quot; *****Program Complete! ...原创 2018-12-23 15:37:34 · 718 阅读 · 0 评论 -
CUDA、SU、MPI和Madagascar混合编程的Makefile文件配置
本人从事地震勘探的科研工作，在Ubuntu系统中用CUDA、SU、MPI和Madagascar做C语言的混合编程。多语言混合编程，关键之处在于Makefile文件的配置Makefile文件的配置。在此给出一个简单的示例，供大家参考。Makefile文件：# # Makefile for GPU program based on SU,Madagascar,MPI and CUDA# Hom...原创 2018-11-17 18:49:36 · 707 阅读 · 4 评论 -
SU操作说明--SU格式转SEGY格式
在使用SU命令将SU格式的地震数据转化成SEGY格式的地震数据时，有两个步骤需要执行。 （1）生成SEGY的 ASC II 或者二进制道头，为写SEGY数据的命令做准备：segyhdrs &lt; data1.su （2）...原创 2018-05-24 09:45:53 · 3964 阅读 · 4 评论 -
SU命令的功能及基本用法-- xmovie
xmovie：以动画形式查看一组数据，通常用于查看一组波场快照，其输入数据是多个波场快照数据组成的一个数据体。如下：首先，我们先将波场快照数据合成一个大的数据集，使用cat命令，部分脚本如下：for (( i = 5; i <=2500 ; i+=20 )); do echo $i cat wf_f_$i.out >> wf_f_all.bindone然后，用xmov...原创 2018-03-29 14:43:33 · 7417 阅读 · 0 评论 -
SU命令的功能及基本用法--psmerge
psmerge：将多个eps格式的图件堆叠在一起，形成组合图，请注意图件的顺序，因为在堆叠时是向前覆盖的，即最后一张放在最上面。基本用法：psmerge in=Fig1.eps in=Fig2.eps in=Fig3.eps > Fig_all.eps...原创 2018-03-04 11:32:27 · 1344 阅读 · 0 评论 -
基于SU的快速傅里叶变换（FFT）
在地震勘探的数据处理过程中，我们经常遇到FFT的问题，我从SU中抽去了一部分代码，自己写了一个FFT的程序，可以做单道、多道、平均谱、和Log谱等，现在将代码贴出来，供大家参考：/* Copyright (c) China University of Petroleum , 2018. *//* All rights reserved. ...原创 2018-03-17 15:30:16 · 1451 阅读 · 9 评论 -
SU数据新旧格式转换问题
最近在用SU做处理，用本机的SU程序对甲方给的数据进行处理时，总是报道头有问题。查了一下，发现甲方给的数据使用低版本的SU生成的，所以需要转换成新版本的数据格式。SU旧格式转新格式的命令为：suoldtonewsuoldtonew < old.su > new.su...原创 2018-03-02 19:22:37 · 1432 阅读 · 2 评论 -
SU（Seismic Unix）与CUDA的混合编程
当前，众多的学者在Linux下面用SU做地震勘探的编程，并与CUDA做融合，进行GPU编程开发。本人研究FWI（全波形反演），也是用CUDA+SU的方式编程C语言的程序。那么在编写MakeFile文件的时候，如何将CUDA的库文件与SU的库文件进行结合呢？下面我给大家一个例子，大家可以参考着进行自己的编写。SU+CUDA的MakeFile文件示例：# Makefile for GP原创 2018-01-10 10:36:16 · 632 阅读 · 0 评论 -
Ubuntu16.04上安装SU（Seismic Unix）的基本步骤
Seismic Unix（SU）是科罗拉多矿业学院的波场研究中心（CWP）开发的一款开源地震数据处理软件，当前很多学者在Linux下用SU做地震数据的处理与编程。我本人是研究FWI的，也在Linux下面用SU做编程来进行地震勘探的科研。分享一下在Linux下面安装SU的一点经验。操作系统：64位的Ubuntu16.04.3；机器：工作站（DELL）SU41R3：下载地址：链接：https://p...原创 2018-01-10 10:58:48 · 13762 阅读 · 20 评论