Mrnip2vrms.c

/* Estimate root mean square velocity (VRMS) from RNIP parameter

 This program is based on multifocus formula that allows to estimate
 VRMS from RNIP: vrms = sqrt(2*RNIP*v0/t0)

 Reference: Diffraction imaging by multifocusing, Berkovitch A.,
 Belfer I., Hassin Y., and Landa E., 2009.

 Programer: Rodolfo A C Neves (Dirack) 21/09/2020
 
 Email: rodolfo_profissional@hotmail.com
 
 License: GPL-3.0 <https://www.gnu.org/licenses/gpl-3.0.txt>.
*/

#include "rnip2vrms_lib.h"
#include <rsf.h>
/*^*/

int main(int argc,char* argv[]){
	
	bool verb; // Verbose parameter
	int nt0; // Number of t0's in stacked section
	float dt0; // Time sampling in stacked section
	float ot0; // Origin of the time axis in stacked section
	int nm0; // Number of CMP's in stacked section
	float dm0; // CMP sampling in stacked section
	float om0; // Origin of the CMP axis in stacked section
	float v0; // Near surface velocity
	float** rnip; // (t0,m0,rnip) section
	float** vrmsSection; // (t0,m0,VRMS) section

	sf_file in,out;
	sf_init(argc,argv);

	in = sf_input("in");
	out = sf_output("out");

	/* Read coordinates of the stacked section file */
	if(!sf_histint(in,"n1",&nt0)) sf_error("No n1= in input");
	if(!sf_histfloat(in,"d1",&dt0)) sf_error("No d1= in input");
	if(!sf_histfloat(in,"o1",&ot0)) sf_error("No o1= in input");
	if(!sf_histint(in,"n2",&nm0)) sf_error("No n2= in input");
	if(!sf_histfloat(in,"d2",&dm0)) sf_error("No d2= in input");
	if(!sf_histfloat(in,"o2",&om0)) sf_error("No o2= in input");

	if(!sf_getbool("verb",&verb)) verb=0;
	/*1: verbose; 0: quiet*/

	if(!sf_getfloat("v0",&v0)) v0=1.5;
	/* Near surface velocity (Km/s) */
	
	if(verb){
		sf_warning("Input file (Rnip Section)");
		sf_warning("nt0=%f dt0=%f ot0=%f",nt0,dt0,ot0);
		sf_warning("nm0=%f dm0=%f om0=%f",nm0,dm0,om0);
		sf_warning("Command line parameters");
		sf_warning("v0=%f",v0);
	}

	/* Load rnip values */
	rnip = sf_floatalloc2(nt0,nm0);
	sf_floatread(rnip[0],nm0*nt0,in);
	vrmsSection = sf_floatalloc2(nt0,nm0);

	/* Calculate VRMS from RNIP */
	vrmsSection = calculateVrmsSectionForRnipVector(rnip,nt0,ot0,dt0,nm0,v0);

	/* Output the VRMS section */
	sf_floatwrite(vrmsSection[0],nt0*nm0,out);
}