%SGYH dessine des dispositifs d'enregistrements 
% voir : http://jmmeost.free.fr/Sismique/isis.html
% fig1 = http://jmmeost.free.fr/Sismique/sgyh1.gif  

clear all

ds = 2; xs = 0:ds:23; ns = length(xs);
dg = 2; xg = [-24:dg:-2 2:dg:24]; ng = length(xg); for is = 1:ns, xgs(is,:) = xs(is)+xg; end  
hmax = xs(1)-xg(1);
dy = min([ds dg])/2; xy1 = (xs(1)+xgs(1,1))/2; xy2 = (xs(ns)+xgs(ns,ng))/2; xy = xy1:dy:xy2; ny = length(xy);
for is = 1:ns, xys(is,:) = (xs(is)+xgs(is,:))/2; end  

xx = xgs(1,1):dg:xgs(ns,ng); nx = length(xx);
[XS XG] = meshgrid(xs,xx);
HG = XG-XS; for is = 1:ns, HG([1:(xgs(is,1)-xgs(1,1))/dg (xgs(is,ng)-xgs(1,1))/dg+2:nx],is) = NaN; end

[YS XY] = meshgrid(xs,xy);
HY = 2*(XY-YS); for is = 1:ns, HY([1:(xys(is,1)-xys(1,1))/dy (xys(is,ng)-xys(1,1))/dy+2:ny],is) = NaN; end

figure(1), clf, colormap('jet')
subplot(2,1,1), grid on, box on, hold on
for is = 1:ns, plot(xs(is),xs(is),'k*'), plot(xgs(is,:),xs(is),'ko','MarkerSize',4),  end, legend('TIRS','RECEPTEURS','Location','NorthWest')
pcolor(xgs(1,1):dg:xgs(ns,ng),xs,HG'), shading('interp')
xlabel('XG (M)'), ylabel('XS (M)'), 
title(['DISPOSITIF MOBILE SYMETRIQUE 24 TRACES : DS = ' num2str(ds) 'M , DG = ' num2str(dg) 'M, DEPORT = (XG - XS) M'])
set(gca,'XLim',[xx(1) xx(nx)],'YLim',[xs(1) xs(ns)],'DataAspectRatio',[1 1 1],'CLim',[-1 1]*hmax), colorbar

subplot(2,1,2), grid on,  box on, hold on
for is = 1:ns,  plot(xs(is),xs(is),'k*'), plot(xys(is,:),xs(is),'k.'), end, legend('TIRS','POINTS-MILIEUX','Location','NorthWest')
pcolor(xy,xs,HY'), shading('interp')
xlabel('XY (M)'), ylabel('XS (M)'), title(['POINTS MILIEUX : DY = ' num2str(dy) 'M, DEPORT = (XG - XS) M'])
set(gca,'XLim',[xx(1) xx(nx)],'YLim',[xs(1) xs(ns)],'DataAspectRatio',[1 1 1],'CLim',[-1 1]*hmax), colorbar

clear all

ds = 2; xs = 0:ds:23; ns = length(xs);
dg = 2; xg = [-24:dg:-2]; ng = length(xg); for is = 1:ns, xgs(is,:) = xs(is)+xg; end  
hmax = xs(1)-xg(1);
dy = min([ds dg])/2; xy1 = (xs(1)+xgs(1,1))/2; xy2 = (xs(ns)+xgs(ns,ng))/2; xy = xy1:dy:xy2; ny = length(xy);
for is = 1:ns, xys(is,:) = (xs(is)+xgs(is,:))/2; end  

xx = xgs(1,1):dg:xgs(ns,ng); nx = length(xx);
[XS XG] = meshgrid(xs,xx);
HG = XG-XS; for is = 1:ns, HG([1:(xgs(is,1)-xgs(1,1))/dg (xgs(is,ng)-xgs(1,1))/dg+2:nx],is) = NaN; end

[YS XY] = meshgrid(xs,xy);
HY = 2*(XY-YS); for is = 1:ns, HY([1:(xys(is,1)-xys(1,1))/dy (xys(is,ng)-xys(1,1))/dy+2:ny],is) = NaN; end

figure(2), clf, colormap('jet')
subplot(2,1,1), grid on, box on, hold on
for is = 1:ns, plot(xs(is),xs(is),'k*'), plot(xgs(is,:),xs(is),'ko','MarkerSize',4), end, legend('TIRS','RECEPTEURS','Location','NorthWest')
pcolor(xgs(1,1):dg:xgs(ns,ng),xs,HG'), shading('interp')
xlabel('XG (M)'), ylabel('XS (M)'), 
title(['DISPOSITIF MOBILE DISSYMETRIQUE 12 TRACES : DS = ' num2str(ds) 'M , DG = ' num2str(dg) 'M, DEPORT = (XG - XS) M'])
set(gca,'XLim',[xx(1) xx(nx)],'YLim',[xs(1) xs(ns)],'DataAspectRatio',[1 1 1],'CLim',[-1 0]*hmax), colorbar

subplot(2,1,2), grid on,  box on, hold on
for is = 1:ns, plot(xs(is),xs(is),'k*'), plot(xys(is,:),xs(is),'k.'), end, legend('TIRS','POINTS-MILIEUX','Location','NorthWest')
pcolor(xy,xs,HY'), shading('interp')
xlabel('XY (M)'), ylabel('XS (M)'), title(['POINTS MILIEUX : DY = ' num2str(dy) 'M, DEPORT = (XG - XS) M'])
set(gca,'XLim',[xx(1) xx(nx)],'YLim',[xs(1) xs(ns)],'DataAspectRatio',[1 1 1],'CLim',[-1 0]*hmax), colorbar

clear all

ds = 2; xs = -10:ds:32; ns = length(xs);
dg = 2; xg = 0:dg:22; ng = length(xg); for is = 1:ns, xgs(is,:) = xg; end  
hmax = xs(ns)-xg(1);
dy = min([ds dg])/2; xy1 = (xs(1)+xgs(1,1))/2; xy2 = (xs(ns)+xgs(ns,ng))/2; xy = xy1:dy:xy2; ny = length(xy);
for is = 1:ns, xys(is,:) = (xs(is)+xgs(is,:))/2; end  

xx = xgs(1,1):dg:xgs(ns,ng); nx = length(xx);
[XS XG] = meshgrid(xs,xx);
HG = XG-XS; for is = 1:ns, HG([1:(xgs(is,1)-xgs(1,1))/dg (xgs(is,ng)-xgs(1,1))/dg+2:nx],is) = NaN; end

[YS XY] = meshgrid(xs,xy);
HY = 2*(XY-YS); for is = 1:ns, HY([1:(xys(is,1)-xys(1,1))/dy (xys(is,ng)-xys(1,1))/dy+2:ny],is) = NaN; end

figure(3), clf, colormap('jet')
subplot(2,1,1), grid on, box on, hold on
for is = 1:ns,  plot(xs(is),xs(is),'k*'), plot(xgs(is,:),xs(is),'ko','MarkerSize',4), end, legend('TIRS','RECEPTEURS','Location','NorthWestOutside')
pcolor(xgs(1,1):dg:xgs(ns,ng),xs,HG'), shading('interp')
xlabel('XG (M)'), ylabel('XS (M)'), 
title(['DISPOSITIF STATIQUE 12 TRACES : DS = ' num2str(ds) 'M , DG = ' num2str(dg) ' M, DEPORT = (XG - XS) M'])
set(gca,'XLim',[xs(1) xs(ns)],'YLim',[xs(1) xs(ns)],'DataAspectRatio',[1 1 1],'CLim',[-1 1]*hmax), colorbar

subplot(2,1,2), grid on,  box on, hold on
for is = 1:ns, plot(xs(is),xs(is),'k*'), plot(xys(is,:),xs(is),'k.'), end, , legend('TIRS','POINTS-MILIEUX','Location','NorthWestOutside')
pcolor(xy,xs,HY'), shading('interp')
xlabel('XY (M)'), ylabel('XS (M)'), title(['POINTS MILIEUX : DY = ' num2str(dy) 'M, DEPORT = (XG - XS) M'])
set(gca,'XLim',[xs(1) xs(ns)],'YLim',[xs(1) xs(ns)],'DataAspectRatio',[1 1 1],'CLim',[-1 1]*hmax), colorbar