# FPOCS2013

**Repository Path**: nicklinyi/FPOCS2013

## Basic Information

- **Project Name**: FPOCS2013
- **Description**: 2013-Computers & Geosciences-On analysis-based two-step interpolation methods for randomly sampled seismic data
- **Primary Language**: Unknown
- **License**: Not specified
- **Default Branch**: master
- **Homepage**: None
- **GVP Project**: No

## Statistics

- **Stars**: 1
- **Forks**: 1
- **Created**: 2018-12-11
- **Last Updated**: 2022-03-19

## Categories & Tags

**Categories**: Uncategorized

**Tags**: None

## README

Name of the program: IST-POCS-FISTA-FPOCS comparison package for seismic data interpolation

Name  of the manuscript: On the analysis-based two-step interpolation methods for the randomly sampled seismic data

Author(s) details:  

	Pengliang Yang(1), Jinghuai Gao(2), Wenchao Chen(3)

	Address: 	Institute of Wave and Information, 
			School of Electronic and Information Engineering, 
			Xi'an Jiaotong University, 
			Xi'an,  P.R. China

	Email(s):	(1) ypl_2100@yahoo.com.cn
			(2) jhgao@mail.xjtu.edu.cn
			(3) wencchen@mail.xjtu.edu.cn


#### Instructions:	

Reproducible research work in the paper, just run

0. 	spectrum_illustration.m

1.	generate_testdata.m

2. 	demo_synthetic1

	demo_synthetic2

	demo_field

These demos will produce a txt file to show the  result details.

Additinal notes: 
(1)	This paper uses FFT to fulfil the task of seismic data interpolation. 
	Of course, other tight frame-based transforms can be utilized instead of FFT.
	In the demo files, you can use DCT to replace FFT by letting:


		A = @(x) idct2(x);
		At = @(x) dct2(x);
	

(2)	These demos tell you how to use this MATLAB package to perform sesic interpolation.





#### The main structure of this MATLAB implementation package 

0. Seismic Dataset:
-	shot600x201.mat
-	field512x430.mat

1. Sampling:
-	random_sampling.m
-	gap_sampling.m
-	regular_sampling.m
-	generate_testdata.m
	
2. Thresholding:
-	HardThresh.m 	(Hard Thresholding)
-	SoftThresh.m	(Soft Thresholding)

3. Algorithms;
-	ist_interp.m
-	pocs_interp.m
-	fista_interp.m <----------- next_t.m
-	fpocs_interp.m <----------- next_t.m
		
4. Auxiliary computation:
-	next_t.m	
-	compute_snr.m

5. Display:
-	seismic.m
	Example: pcolor(X),shading interp
		 colormap(seismic(3))

6. Synthetic data producing file taken from SeismicLab 
	(http://www-geo.phys.ualberta.ca/saig/SeismicLab)
-	ricker.m (generate ricker wavelet)
-	hyperbolic_events.m (generate synthetic seismic data including hyperbolic events)