00001
00002
00003
00004
00005
00006
00007
00008
00009
00010
00011
00012
00013
00014
00015
00016
00017
00018
00019
00020
00021
00022
00023
00024
00025
00026
00027
00028
00029
00030
00031
00032
00033 #include <math.h>
00034
00035 #include "image.h"
00036 #include "proto2D.h"
00037 #include "radar.h"
00038
00039 int main(int argc, char *argv[]){
00040
00041 imau2 im0;
00042 imafl im2, moyint, moycarre;
00043 pixfl tab2[65536];
00044
00045
00046 read_ima_t rea;
00047 moy2D_t moy;
00048 moycarre2D_t moycar;
00049 kuan_t kuan;
00050 write_ima_t wri;
00051
00052
00053 param par0, *ptp;
00054 int i, j;
00055 double tmpdb;
00056
00057
00058 param_debut(argc, argv, &par0);
00059 ptp = &par0;
00060 ptp = read_ima_lect(&rea, ptp, ">> image initiale (u2) :");
00061 ptp = kuan_lect(&kuan, ptp, "Kuan :");
00062 moy.dimX = kuan.dimX; moycar.dimX = kuan.dimX;
00063 moy.dimY = kuan.dimY; moycar.dimY = kuan.dimY;
00064 ptp = write_ima_lect(&wri, ptp, ">> image resultat (u2) :");
00065 param_fin(argc, argv, &par0);
00066
00067
00068 read_imau2_init(&rea, &im0);
00069
00070 im2.nc = im0.nc;
00071 im2.nr = im0.nr;
00072 alloc_imafl(&im2);
00073 for(j=0; j<65536; j++)
00074 tab2[j] = (pixfl)j*(pixfl)j;
00075 for(j=0; j<im0.nr; j++)
00076 for(i=0; i<im0.nc; i++)
00077 im2.p[j][i] = tab2[im0.p[j][i]];
00078 moy2Dfl_init(&moy, im2, &moyint);
00079 moycarre2Dfl_init(&moycar, im2, &moycarre);
00080 ;
00081 write_ima_init(&wri);
00082
00083
00084 moy2Dfl_calc(&moy, im2, &moyint);
00085 moycarre2Dfl_calc(&moycar, im2, &moycarre);
00086 kuanfl_calc(&kuan, im2, moyint, moycarre, &im2);
00087
00088 for(j=0; j<im0.nr; j++)
00089 for(i=0; i<im0.nc; i++){
00090 tmpdb = sqrt((double)im2.p[j][i]) +0.5;
00091 if( tmpdb < 65536. )
00092 im0.p[j][i] = (pixu2)tmpdb;
00093 else
00094 im0.p[j][i] = 65535;
00095 }
00096
00097 write_imau2_ferm(&wri, im0);
00098 }
00099
00100
00101