Crear cuenta
scatter(Lon,Lat,30,Temp,'Filled')
pcolor(Lon,Lat,Temp)
a = magic(4);
image(a);
map = hsv(16);
colormap(map)
colorbar
%figure(3)
subplot(1,3,3);
mat = vec2mat(d_est,71)%se hace para que queden matrices del mismo tamaño
z=mat';
contourf(xx,yy,z); %ya solucioné el problema del tamaño de las matrices
axis image
title('Kriging Variable')
%figure(3)
subplot(1,3,3);
mat = vec2mat(d_est,71)%se hace para que queden matrices del mismo tamaño
z=mat';
pcolor(xx,yy,z); %ya solucioné el problema del tamaño de las matrices
shading flat
hold on
plot(pos_known(:,1),pos_known(:,2),'.k')
axis image
title('Kriging Variable')
view(2)
colorbar
%http://sourceforge.net/projects/mgstat/files/mGstat/
%%mGstat : A Geostatistical Matlab toolbox
pos_known=[1 1;1 6;5 1;5 6];
val_known=[24 23 22 20]';
x=0:.1:6;
y=0:.1:7;
[xx,yy]=meshgrid(x,y);
pos_est=[xx(:) yy(:)];
V=deformat_variogram('1 Sph(1)');
[d_est,d_var]=gstat_krig(pos_known,val_known,pos_est,V);
figure(1)
subplot(1,3,1);
scatter(pos_est(:,1),pos_est(:,2),10,d_est)
axis image;
title('Kriging mean')
subplot(1,3,2);
%figure(2)
scatter(pos_est(:,1),pos_est(:,2),10,d_var)
axis image;
title('Kriging variance')
%figure(3)
subplot(1,3,3);
mat = vec2mat(d_est,length(x));%se hace para que queden matrices del mismo tamaño
z=mat;
contour3(xx,yy,z,600)
colormap(jet(9))
colorbar
view(0,0)
hold off>> grafica_de_temperatura
Trying to run GSTAT on gstat.cmd