linedecimation.c

超声波成像算法

/*           1儔僀儞偺憸傪梌偊傜傟偨僒儞僾儖悢偵娫堷偒   */
/*  丂丂丂丂丂丂丂丂嶰師尦僨乕僞偺攝楍偵奿擺偡傞         */

/*                 LineDecimation                        */
#include <math.h>
#include <stdlib.h>
#include "Reconst.h"

void LineDecimation(int ix,int iy,int nx, int ny, int nz,int ns,int ib, 
					 COMPLEX16 *s,double *v)
/*	堷悢丗
/*		ix,iy		丗3師尦僨乕僞偺x,y傪梌偊傞僀儞僨僢僋僗
/*		int nx,ny,nz丗x,y,z曽岦偺攝楍偺師尦
/*					ix*ny*nz + iy*nz + iz偱傾僋僙僗偡傞
/*		int ns		丗儔僀儞僨乕僞偺攝楍偺悺朄
/*		COMPLEX16 s[i]丗暋慺嶰師尦僨乕僞傪梌偊傞攝楍
/*		double v[i]		丗嶰師尦僨乕僞傪奿擺偡傞攝楍
/************************************************************/

{
	/*寁嶼傪崅懍偵幚峴偡傞偨傔丄梌偊傜傟偨1儔僀儞僨乕僞傪嬫娫
	丂偵暘妱偟丄嬫娫枅偺嵟戝傪戙昞抣偲偡傞丅*/
	int i,istep,iz,ix_bias,iy_bias ;
	double a,b,c,d, max, space, z, max1=0.0;
	//ib=ix*ny*ns+iy*ns;
	//space=(double)ns/(double)nz;//嬫娫偺暆
	space=(double)ns/(double)nz;
	ix_bias=ix*ny*nz;
	iy_bias=iy*nz;
	z = 0.0;
	i = 0;
	for(iz=0;iz<nz;iz++){//墱峴偒偑256側偺偱丄幚嵺偺夋慺抣256屄傪憸嵞惗偝傟偨夋慺2048屄偺拞偐傜慖傇
		z += space;
		//istep=((int)z)%ns;//嬫娫偺忋尷
		istep=((int)z)%ns;
		/*尰嵼偺index偐傜istep傑偱偺嬫娫偺嵟戝抣傪媮傔傞*/
		max=0.0;
		for(;i<=istep;i++){//嬫娫撪偺夋慺偺偆偪嵟戝偺傕偺傪慖傇
			a = s[ib+i].real*s[ib+i].real+s[ib+i].imag*s[ib+i].imag;//愨懳抣傪偲偭偰夋慺抣傪媮傔傞
			//b = s[ib+i+1].real*s[ib+i+1].real+s[ib+i+1].imag*s[ib+i+1].imag;
			//c = s[ib+i+2].real*s[ib+i+2].real+s[ib+i+2].imag*s[ib+i+2].imag;
			if(max<a) max =a;
			//max1+=a;
			//d=(a+b+c)/3;
			//if(max<d) max =d;


		}
		//max=max1/space;//?
		v[ix_bias+iy_bias+iz]=max;//嬫娫撪偺嵟戝夋慺抣傪3師尦夋憸攝楍偵奿擺
	}
	return;
}
/**********************************

   (墘嶼偱媮傔傜傟偨壖憐價乕儉1杮暘偺夋慺悢)        2048
------------------------------------------------ = ------ = 8
丂丂(3師尦夋憸傪偮偔傞偺偵昁梫側墱峴偒夋慺悢乯    丂 256    

  埲忋偐傜丄侾嬫娫偼8屄偺夋慺偱峔惉偝傟傞丅偙偺側偐偐傜夋慺偑嵟戝偱偁傞傕偺傪慖傃
  幚嵺偺3師尦夋憸偱昁梫側夋慺偵側傞丅
  偙偺張棟傪慡壖憐價乕儉杮悢暘偩偗孞傝曉偡丅

*************************************/