Aprende a hacer Batimetrias en MATLAB
mayo 24, 2020Los perfiles realizados dentro de este estudio así como las batimetrías fueron tratadas y grafícadas a través de la plataforma de MATLAB, por lo cual aquí te enseñamos como procesar los datos de alguna costa de tu interés.
Lo primero que necesitas son los datos batimetricos de tu zona, es decir: las longitudes, latitudes y profundidades de cada una de tus mediciones recopiladas en un archivo de texto (para esta demostración ocupamos los datos obtenidos del canal de la Bocana en sus diferentes recorridos)
Haz click en el video para ver la explicación
Comandos que puedes utilizar para hacer tus propias batimetrias (Ejemplo del recorrido en zigzag del canal de la Bocana)
%grafica del recorrido en zigzag
batzig=dlmread('BatimetriazigZagcorregidos.txt');
xzig=batzig(:,2)
yzig=batzig(:,1)
zzig=batzig(:,3)
figure(4)
plot(xzig,yzig,'pk')
xlabel('Longitud')
ylabel('Latitud')
title('Recorrido en Zigzag')
set(gca,'FontName','Comic Sans MS','FontSize',14,'LineWidth',2.5);
%%%grafica en tres dimensiones
xxzig=min(xzig):0.00005:max(xzig);
yyzig=min(yzig):0.00005:max(yzig);
[xxzig2,yyzig2]=meshgrid(xxzig,yyzig)
zzzig2=griddata(xzig,yzig,zzig,xxzig2,yyzig2);
figure (5)
pcolor(xxzig2,yyzig2,zzzig2)
shading interp;
c = colorbar
set(c,'location','southoutside','LineWidth',2.5)
colormap('HSV')
caxis([-7 0])
xlabel('Longitud')
ylabel('Latitud')
zlabel('Profundidad')
title('Recorrido en Zigzag')
set(gca,'FontName','Comic Sans MS','FontSize',14,'LineWidth',2.5);
%%%%grafico surf
figure (6)
surf(xxzig2,yyzig2,zzzig2)
shading interp;
c = colorbar
set(c,'location','southoutside','LineWidth',2.5)
colormap('HSV')
caxis([-7 0])
xlabel('Longitud')
ylabel('Latitud')
zlabel('Profundidad')
title('Recorrido en Zigzag')
set(gca,'FontName','Comic Sans MS','FontSize',14,'LineWidth',2.5);
%%
Comandos utilizados para hace la batimetria del canal de la Bocana
batrec=dlmread('Batimetriarectacorregidos.txt');
%recorrido linea recta
figure (1)
x=batrec(:,2);
y=batrec(:,1);
z=batrec(:,3);
plot(x,y,'pk')
xlabel('Longitud')
ylabel('Latitud')
title('Recorrido en linea recta')
set(gca,'FontName','Comic Sans MS','FontSize',14,'LineWidth',2.5);
%tres dimendiones recorrido en linea recta
xx=min(x):0.00005:max(x);
yy=min(y):0.00005:max(y);
[xx2,yy2]=meshgrid(xx,yy)
zz2=griddata(x,y,z,xx2,yy2);
figure (2)
surf(xx2,yy2,zz2)
shading interp;
c = colorbar
set(c,'location','southoutside','LineWidth',2.5)
colormap('HSV')
caxis([-7 0])
xlabel('Longitud')
ylabel('Latitud')
zlabel('Profundidad (m)')
title('Recorrido en linea recta')
set(gca,'FontName','Comic Sans MS','FontSize',14,'LineWidth',2.5);
%plano
figure (7)
pcolor(xx2,yy2,zz2)
shading interp;
c = colorbar
set(c,'location','southoutside','LineWidth',2.5')
colormap('HSV')
caxis([-7 0])
xlabel('Longitud')
ylabel('Latitud')
zlabel('Profundidad (m)')
title('Recorrido en linea recta')
set(gca,'FontName','Comic Sans MS','FontSize',14,'LineWidth',2.5);
%grafica del recorrido en zigzag
batzig=dlmread('BatimetriazigZagcorregidos.txt');
xzig=batzig(:,2)
yzig=batzig(:,1)
zzig=batzig(:,3)
figure(4)
plot(xzig,yzig,'pk')
xlabel('Longitud')
ylabel('Latitud')
title('Recorrido en Zigzag')
set(gca,'FontName','Comic Sans MS','FontSize',14,'LineWidth',2.5);
%%%grafica en tres dimensiones
xxzig=min(xzig):0.00005:max(xzig);
yyzig=min(yzig):0.00005:max(yzig);
[xxzig2,yyzig2]=meshgrid(xxzig,yyzig)
zzzig2=griddata(xzig,yzig,zzig,xxzig2,yyzig2);
figure (5)
pcolor(xxzig2,yyzig2,zzzig2)
shading interp;
c = colorbar
set(c,'location','southoutside','LineWidth',2.5)
colormap('HSV')
caxis([-7 0])
xlabel('Longitud')
ylabel('Latitud')
zlabel('Profundidad')
title('Recorrido en Zigzag')
set(gca,'FontName','Comic Sans MS','FontSize',14,'LineWidth',2.5);
%%%%grafico surf
figure (6)
surf(xxzig2,yyzig2,zzzig2)
shading interp;
c = colorbar
set(c,'location','southoutside','LineWidth',2.5)
colormap('HSV')
caxis([-7 0])
xlabel('Longitud')
ylabel('Latitud')
zlabel('Profundidad')
title('Recorrido en Zigzag')
set(gca,'FontName','Comic Sans MS','FontSize',14,'LineWidth',2.5);
%%
%%
figure (3)
plot(x,y,'pk')
xlabel('Longitud')
ylabel('Latitud')
title('Recorridos en zigzag y linea recta')
set(gca,'FontName','Comic Sans MS','FontSize',14,'LineWidth',2.5);hold on
plot(xzig,yzig,'or')
xlabel('Longitud')
ylabel('Latitud')
set(gca,'FontName','Comic Sans MS','FontSize',14,'LineWidth',2.5);
%%%%figure(8)
surf(xx2,yy2,zz2);hold on
surf(xxzig2,yyzig2,zzzig2)
shading interp;
c = colorbar
set(c,'location','southoutside','LineWidth',2.5)
colormap('HSV')
caxis([-7 0])
xlabel('Longitud')
ylabel('Latitud')
zlabel('Profundidad (m)')
title('Recorrido en linea recta y zigzag')
set(gca,'FontName','Comic Sans MS','FontSize',14,'LineWidth',2.5);
%%
figure(9)
pcolor(xxzig2,yyzig2,zzzig2);hold on
pcolor(xx2,yy2,zz2)
shading interp;
c = colorbar
set(c,'location','southoutside','LineWidth',2.5)
colormap('HSV')
caxis([-7 0])
xlabel('Longitud')
ylabel('Latitud')
zlabel('Profundidad')
title('Ambos recorridos')
set(gca,'FontName','Comic Sans MS','FontSize',14,'LineWidth',2.5);
0 comentarios