img_sar_azimuth_compress

 

	int no_samples;
	int i,j;
	float taz;
	float range;
	float arg;
	float a,b,c;
	float rx;
	image_t		img;
	int width,height;
	int order,pts;
	cap_fft_cpx x,y;
	float fmax,fmin;

 

	if( (obufs = NO_OUTPUT_BUFFERS()) < 1 ) {
		fprintf(stderr,"img_sar_azimuth_compress: no output buffers\n");
		CsInfo("img_sar_create: no output buffers");
		return(2);
	}
    dt=(1.0/fs)*0.000001;
	t0=2.0*rp*1000.0/C;
    maxRangeIndex=(int)tau*0.000001/dt;
	dtaz=(1.0/prf);
    tp=tpi*(1.0/prf);
    done=0;
    t=0;
    k=0;
    lamda=C/(fc*1000000.0);
// 	fprintf(stderr,"img_sar_azimuth_compress: width=%d height=%d\n",maxRangeIndex,maxAzimuthIndex); 
    for(j=0; j
    

 


for (no_samples = MIN_AVAIL(); no_samples > 0; --no_samples) {
	IT_IN(0);
	img=x(0);
	height=img.height;
	width=img.width;
	/*
	 * round height to next power of 2
	 * Note we are dealing with Azimuth
	 * Also height is already a power of 2 (should be)
	 */
	order = (int) (log((float)height)/log(2.0)+0.5);
	pts = 1 << order;
	if (pts < height ) {
        	pts = pts*2;
        	order += 1;
	}
	rangeFFTLength=width;
	azimuthFFTLength=pts;
	fprintf(stderr,"img_sar_azimuth_compress: rangeFFTLength=%d height=%d azimuthFFTLength=%d\n",width,height,azimuthFFTLength); 

	ref_P=(cap_fft_cpx *)calloc(azimuthFFTLength,sizeof(cap_fft_cpx));
	azimuth_P=(cap_fft_cpx *)calloc(azimuthFFTLength,sizeof(cap_fft_cpx));
	refFFT_P=(cap_fft_cpx *)calloc(azimuthFFTLength,sizeof(cap_fft_cpx));
	azimuthFFT_P=(cap_fft_cpx *)calloc(azimuthFFTLength,sizeof(cap_fft_cpx));

	if(ref_P == NULL || azimuth_P== NULL || azimuthFFT_P== NULL || refFFT_P== NULL) {
		fprintf(stderr,"img_sar_azimuth_compress: could not allocate space \n");
		CsInfo("img_sar_azimuth_compress: could not allocate space");
		return(7);
	}

	mat_PP = (float**)calloc(azimuthFFTLength,sizeof(float*));
	if(mat_PP == NULL) {
		fprintf(stderr,"img_sar_azimuth_compress: could not allocate space \n");
		CsInfo("img_sar_azimuth_compress: could not allocate space");
		return(5);
	}
    for(i=0; ia ) fmin=a;
            
			mat_PP[j][i]=a;
	     }
#endif 

#if 000
         for(j=0; ja ) fmin=a;
            
			mat_PP[j][i]=a;
	     }
#endif 
	
 
	}
	fprintf(stderr,"img_sar_azimuth_compress: fmin=%f fmax=%f \n",fmin,fmax); 

    for(j=0; j
    

 

if(azimuth_P) free(azimuth_P);
if(azimuthFFT_P) free(azimuthFFT_P);
if(refFFT_P) free(refFFT_P);
if(ref_P) free(ref_P);

/*  Capsim (r) Text Mode Kernel (TMK) Star Library (Blocks)
    Copyright (C) 1989-2017  Silicon DSP Corporation

    This library is free software; you can redistribute it and/or
    modify it under the terms of the GNU Lesser General Public
    License as published by the Free Software Foundation; either
    version 2.1 of the License, or (at your option) any later version.

    This library is distributed in the hope that it will be useful,
    but WITHOUT ANY WARRANTY; without even the implied warranty of
    MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the GNU
    Lesser General Public License for more details.

    You should have received a copy of the GNU Lesser General Public
    License along with this library; if not, write to the Free Software
    Foundation, Inc., 59 Temple Place, Suite 330, Boston, MA  02111-1307  USA

    http://www.silicondsp.com
    Silicon DSP  Corporation
    Las Vegas, Nevada
*/

 

/************************************************************************
 *
 * Azimuth compression
 * SAR Image Formation
 * 
 * (c) 1990-2002 XCAD Corporation
 *									
 *		

img_sar_azimuth_compress


Azimuth compression SAR Image Formation


Sasan Ardalan
 
 ************************************************************************/