Bonsoir à tous!

Donc voilà le petit soucis qui me tracasse.

Je souhaite faire une carte à l'échelle de l'Europe en me servant de la fonction etopo. Afin d'avoir ce que je souhaite en terme de "rendu", j'utilise ensuite geoshow et demcmap pour le graphique (en PJ pour l'exemple)

Là où je bloque, c'est que je souhaite faire 'disparaitre' l'océan, afin de ne conserver que la partie terrestre... et je bloque...

Merci de votre aide!!!

Code : Sélectionner tout - Visualiser dans une fenêtre à part
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
clear all
close all
clc
 
samplefactor=4;
[Z,refvec]=etopo('etopo1_ice_c.flt',samplefactor);
refvec=[60/samplefactor, 72, -12];
 
pas=samplefactor/60;
lat=[-90:pas:90-(pas)]';
long=[-180:pas:180-(pas)]';
 
% Filtre des latitudes
clear f
f= lat(:)>=29 & lat(:)<=76;
Z2=Z(f,:);
 
% Filtre des longitudes
clear f
f=find(long(:)>=-13 & long(:)<=61);
Z2=Z2(:,f);
 
close
figure1 = figure('Color',[1 1 1]);
ag1=worldmap('World');
set(gcf,'Color',[1,1,1])
h = worldmap([30 72],[-12 60]);
geoshow(Z2, refvec, 'DisplayType', 'texturemap');
demcmap(Z2)


J'ai trouvé la solution... Il fait utiliser une extension de la fonction demcmap:

Code : Sélectionner tout - Visualiser dans une fenêtre à part
demcmap(Z2,256,[1 1 1],[])