plot3ddem#
plot3ddem(x, y, z, plottype, cmapcolors)
Description#
Plot x (longitude), y (latitude) and z (elevation) data into a defined mesh
Inputs#
- x
- x (longitude) data extracted from xyz fileSet x=[] if it is not availableIt may be 1d or 2d array
- y
- y (latitude) data extracted from xyz fileSet y=[] if it is not availableIt may be 1d or 2d array
- z
z (elevation) data extracted from xyz file
- plottype=’imagesc’;
- Plot type‘imagesc’: 2 dimensional plot using imagesc or imshow‘pcolor’: 2 dimensional plot using pcolor‘contour’: 2 dimensional contour plot, number of contour=ncolor‘surface’: 3 dimensional surface plot
- cmapcolors=’topocmap’;
- Colormap for z dataTopographic (Water/Land) colormaps:‘topocmap’: colormap developed by Arash Karimpour‘topocmaprelief’: colormap developed by Arash Karimpour‘topocmapocean’: colormap developed by Arash Karimpour‘topocmapearth’: colormap developed by Arash Karimpour‘blueocean’: colormap developed by Arash Karimpour‘blueoceansea’: colormap developed by Arash Karimpour‘greenearth’: colormap developed by Arash Karimpour‘greenearthland’: colormap developed by Arash Karimpour‘blgrtopocmap’: colormap developed by Arash Karimpour‘blrdtopocmap’: colormap developed by Arash Karimpour‘grayearth’: colormap developed by Arash Karimpour‘etopo1’: ETOPO1 colormap, https://www.ngdc.noaa.gov/mgg/global/global.html‘gmtglobe’: GMT_globe colormap, https://www.giss.nasa.gov/tools/panoply/colorbars/‘gmtrelief’: GMT_relief colormap, https://www.giss.nasa.gov/tools/panoply/colorbars/‘aendekerk’: Colormap from Florian Aendekerk, http://www.mathworks.com/matlabcentral/fileexchange/63590-landseacolormap-m-Any other available color map such as ‘cool’, ‘winter’, etc can be usedColormap can be defined by user as [n*3] array in RGB color format between 0 and 255
Outputs#
Examples#
[x,y]=meshgrid(linspace(-10,10,50),linspace(-10,10,50));
r=sqrt(x.^2+y.^2)+1e-10; %Add 1e-10 to prevent divide by 0
z=sin(r)./r;
plot3ddem(x,y,z,'pcolor','topocmap')
[x,y]=meshgrid(linspace(-10,10,21),linspace(-10,10,21));
z=(sin(x)+sin(y))./(x+y+1e-10); %Add 1e-10 to prevent divide by 0
plot3ddem(x,y,z,'pcolor','topocmap')
x(:,1)=10.*rand(1000,1);
y(:,1)=10.*rand(1000,1);
z=x.^2+y.^2;
plot3ddem(x,y,z,'pcolor','topocmap')
References#
Geospatial data
Colormap
Color