82 long meno,i,j,k,row,col,n,ii,ncol;
84 double r,x,y,dx,dy,coordx,coordy;
123 r=sqrt(pow(x,(
double)2)+pow(y,(
double)2));
125 var->
co[i][j]=
gamma2(r, int_scale, variance);
141 variogramma2(E, N , Z0, variogr, UV->
U, int_scale, variance);
150 gam->
co[k]=variogr->
co[ii][k];
166 if(Z0->
co[i][j]!=UV->
V->
co[2]){
168 weights->
co[ii][k]=gam->
co[k];
178 weights->
co[ii][1]=variogr->
co[ii][2]/(variogr->
co[ii][1]+variogr->
co[ii][2]);
179 weights->
co[ii][2]=variogr->
co[ii][1]/(variogr->
co[ii][1]+variogr->
co[ii][2]);
189 for(i=1;i<=Z0->
nrh;i++){
190 for(j=1;j<=Z0->
nch;j++){
192 if(Z0->
co[i][j]!=UV->
V->
co[2]){
193 weights->
co[ii][1]=1;
199 for(i=1;i<=ncol;i++){
201 if(weights->
co[i][j]<=0) weights->
co[i][j]=0;
247 double X,Y,rx,ry,rz,h2;
250 for(i=1;i<=Z0->
nrh;i++){
251 Y=U->
co[3]+(Z0->
nrh-i)*U->
co[2]+(U->
co[2]/(
double)2);
252 for(j=1;j<=Z0->
nch;j++){
254 for(k=1;k<=E->
nh;k++){
255 X=U->
co[4]+(j-1)*U->
co[1]+(U->
co[1]/(
double)2);
259 h2=sqrt(pow(rx,(
double)2)+pow(ry,(
double)2)+pow(rz,(
double)2));
260 if(h2!=0)variogr->
co[ii][k]=
gamma2(h2,int_scale,variance);
294 double gamma2(
double r,
double scala_integr,
double varianza)
299 gamma2=varianza*(1.0-exp(-r/scala_integr));