1.直接在demos上修改，需要准备Sigsbee模型的二进制文件，可以在本博客的Research links中下载，若是Segy格式先去掉卷头和道头，用Matlab读入二进制并用imresize函数调整模型大小。本文使用的模型大小为300*900,dx=dz=10

2.在su/src/demos/Migration_Inversion/Pre_stack/Premigpspi中直接修改。

 1.将parfile1.par和parfile2.par中的模型参数改成Sigsbee模型，其余参数可不要管，如下图：

 2.修改MakeData,将unif2注掉，并把所有的vfile改成已近准备好的Sigsbee二进制文件。按照自己的意愿修改观测系统参数NSHOTS,FHSHOT,DSHOT,如下图：

对切除初至的命令参数进行修改：

同理对Migration.sh中的相应模型参数、初至速度、观测系统进行修改即可。

./MakeData正演得到19炮记录：

执行./Migration.sh完成PSPI叠前偏移，结果如下图：