% dŽfinition de coordonnŽes sous forme de matrice
% contenant des hauteurs. Ces donnŽes ont ŽtŽ prises
% via Internet ˆ l'adresse: http://swisstopo.ch.

%DHM25-MATRIXMODELL                             (c)BUNDESAMT F. LANDESTOPOGRAPHIE
%MODELLNUMMER                                   1111 TEIL
%MODELLNAME                                     TESTDATENSATZ ALBIS            
%MODELLVERSION            06.12.1995   14:43    ZEITPUNKT DER INTERPOLATION
%PROGRAMM                 dhm_int     V 2.19    INTERPOLATIONSPROGRAMM
%--------------------------------------------------------------------------------
%NORD-WEST ECKE     [M]  679000.0  238000.0     ERSTER HOEHENWERT
%SUED-OST ECKE      [M]  683000.0  234000.0     LETZTER HOEHENWERT
%MASCHENWEITE WE/NS [M]      25.0      25.0
%MATRIXDIMENSIONEN WE/NS    161       161       TOTAL     25921 MATRIXPUNKTE%
%HOEHENBEREICH     [DM]    5121      8930       (6 CHARACTER PRO HOEHENWERT)
%--------------------------------------------------------------------------------
%STAND DES BASISMODELLS    1982                 HOEHENKURVEN
%                          1982                 SEEKONTUREN
%

% on charge les donnŽes
load albis

% on rŽarrange la matrice
B=reshape(A',1,77*340);

% on Žlimine les zŽros de remplissage
i=find(B==0);
C=B(1:i(1)-1);

% la matrice est une matrice carrŽe de 161 sur 161
D=reshape(C',161,161);

% on purge
clear A B C

% on dessine
surfl(D)
shading interp
colormap(gray)
axis off
rotate3d on
