How to add a left pointing arrow with annotation to a plotted DEM in MATLAB?
I have been trying to add a left pointing arrow with the label ‘Mt Everest’ (both in white) to an extracted portion of a DEM. This is the code I have to get to this point, however the arrow and text does not appear when I run the script despite the rest of the properties working.
% extract_mt_everest_dem
%% Load the Mt Everest DEM data
[DEM, DEM_R] = geotiffread(‘aster_gdem_mosaic_UTM_crop.tif’);
DEM_info = geotiffinfo(‘aster_gdem_mosaic_UTM_crop.tif’);
% Manage NoData cells to make then NaN
DEM(DEM<=0) = NaN;
% Extract the lat and lon values for each cell within the data
% Create an index variable
idx = 1:length(DEM(:));
% Get the indices for each cell in the DEM
[I,J] = ind2sub(size(DEM),idx);
clear idx
% Convert the cell indices into projected coordinates
[x,y] = pix2map(DEM_info.RefMatrix, [I(:) J(:)]);
clear I J
% Reshape the x and y vectors to a grid the same size as the data
x_grid = reshape(x,size(DEM));
y_grid = reshape(y,size(DEM));
% Calculate the lat and lon values from the projected coordinates
[lat,lon] = projinv(DEM_info, x(:), y(:));
clear x y
% Reshape the lat and lon vectors to a grid the same size as the data
lon_grid = reshape(lon,size(DEM));
lat_grid = reshape(lat,size(DEM));
clear lat lon
%% Choose the point of interest (POI): Mt Everest
mt_everest = [27.986065 86.922623];
[mt_everest_X, mt_everest_Y] = projfwd(DEM_info, mt_everest(1), mt_everest(2));
%% Calculate distance from Aoraki Mt Cook
x_diff = x_grid – mt_everest_X;
y_diff = y_grid – mt_everest_Y;
% Find the indices of the bounding box of cells within a specified ditance of 20×20 km
max_dist = 10000; % (m)
row_start = find(abs(y_diff(:,1))<=max_dist, 1);
row_end = find(abs(y_diff(:,1))<=max_dist, 1, ‘last’);
col_start = find(abs(x_diff(1,:))<=max_dist, 1);
col_end = find(abs(x_diff(1,:))<=max_dist, 1, ‘last’);
% Extract the ROI from whole array using indexing
DEM_ROI = DEM(row_start:row_end,col_start:col_end);
lat_grid_ROI = lat_grid(row_start:row_end,col_start:col_end);
lon_grid_ROI = lon_grid(row_start:row_end,col_start:col_end);
%% Make a map to display the subset data
% Set the font size
font_size = 18;
set(0, ‘DefaultAxesFontSize’, font_size);
% Make the figure the size of the full screen
figure(‘units’,’normalized’,’outerposition’,[0 0 1 1]);
% Set the backrgound colour to white
set(gcf,’color’,’w’);
% Make a projected map axis using worldmap
ax = worldmap([nanmin(lat_grid_ROI(:)) nanmax(lat_grid_ROI(:))],[nanmin(lon_grid_ROI(:)) nanmax(lon_grid_ROI(:))]);
% Set axis font size to 18
setm(ax,’FontSize’,font_size);
hold on % to allow for overplotting
pcolorm(lat_grid_ROI, lon_grid_ROI, DEM_ROI);
cmap = landcolor; % get the landcolour colour map
colormap(cmap); % use the landcolour colour map
shadem; % shade the DEM to highlight the relief
hbar = colorbar(‘northoutside’,’fontsize’,font_size); % create a colourbar
set(get(hbar,’Xlabel’),’String’,’Elevation (m)’); % add a colourbar label
% Add a scale bar
scalebar(‘color’,[0 0 0],’FontSize’,font_size,’location’,’southeast’);
% Add a label showing Mt Everest with a left-pointing arrow
b = textm(mt_everest(1), mt_everest(2), 4000,’leftarrow Mt. Everest’,’FontSize’, font_size,’Color’,[1 1 1]);
uistack(b, ‘top’);
close
Anyone know how to get my arrow and label to appear?I have been trying to add a left pointing arrow with the label ‘Mt Everest’ (both in white) to an extracted portion of a DEM. This is the code I have to get to this point, however the arrow and text does not appear when I run the script despite the rest of the properties working.
% extract_mt_everest_dem
%% Load the Mt Everest DEM data
[DEM, DEM_R] = geotiffread(‘aster_gdem_mosaic_UTM_crop.tif’);
DEM_info = geotiffinfo(‘aster_gdem_mosaic_UTM_crop.tif’);
% Manage NoData cells to make then NaN
DEM(DEM<=0) = NaN;
% Extract the lat and lon values for each cell within the data
% Create an index variable
idx = 1:length(DEM(:));
% Get the indices for each cell in the DEM
[I,J] = ind2sub(size(DEM),idx);
clear idx
% Convert the cell indices into projected coordinates
[x,y] = pix2map(DEM_info.RefMatrix, [I(:) J(:)]);
clear I J
% Reshape the x and y vectors to a grid the same size as the data
x_grid = reshape(x,size(DEM));
y_grid = reshape(y,size(DEM));
% Calculate the lat and lon values from the projected coordinates
[lat,lon] = projinv(DEM_info, x(:), y(:));
clear x y
% Reshape the lat and lon vectors to a grid the same size as the data
lon_grid = reshape(lon,size(DEM));
lat_grid = reshape(lat,size(DEM));
clear lat lon
%% Choose the point of interest (POI): Mt Everest
mt_everest = [27.986065 86.922623];
[mt_everest_X, mt_everest_Y] = projfwd(DEM_info, mt_everest(1), mt_everest(2));
%% Calculate distance from Aoraki Mt Cook
x_diff = x_grid – mt_everest_X;
y_diff = y_grid – mt_everest_Y;
% Find the indices of the bounding box of cells within a specified ditance of 20×20 km
max_dist = 10000; % (m)
row_start = find(abs(y_diff(:,1))<=max_dist, 1);
row_end = find(abs(y_diff(:,1))<=max_dist, 1, ‘last’);
col_start = find(abs(x_diff(1,:))<=max_dist, 1);
col_end = find(abs(x_diff(1,:))<=max_dist, 1, ‘last’);
% Extract the ROI from whole array using indexing
DEM_ROI = DEM(row_start:row_end,col_start:col_end);
lat_grid_ROI = lat_grid(row_start:row_end,col_start:col_end);
lon_grid_ROI = lon_grid(row_start:row_end,col_start:col_end);
%% Make a map to display the subset data
% Set the font size
font_size = 18;
set(0, ‘DefaultAxesFontSize’, font_size);
% Make the figure the size of the full screen
figure(‘units’,’normalized’,’outerposition’,[0 0 1 1]);
% Set the backrgound colour to white
set(gcf,’color’,’w’);
% Make a projected map axis using worldmap
ax = worldmap([nanmin(lat_grid_ROI(:)) nanmax(lat_grid_ROI(:))],[nanmin(lon_grid_ROI(:)) nanmax(lon_grid_ROI(:))]);
% Set axis font size to 18
setm(ax,’FontSize’,font_size);
hold on % to allow for overplotting
pcolorm(lat_grid_ROI, lon_grid_ROI, DEM_ROI);
cmap = landcolor; % get the landcolour colour map
colormap(cmap); % use the landcolour colour map
shadem; % shade the DEM to highlight the relief
hbar = colorbar(‘northoutside’,’fontsize’,font_size); % create a colourbar
set(get(hbar,’Xlabel’),’String’,’Elevation (m)’); % add a colourbar label
% Add a scale bar
scalebar(‘color’,[0 0 0],’FontSize’,font_size,’location’,’southeast’);
% Add a label showing Mt Everest with a left-pointing arrow
b = textm(mt_everest(1), mt_everest(2), 4000,’leftarrow Mt. Everest’,’FontSize’, font_size,’Color’,[1 1 1]);
uistack(b, ‘top’);
close
Anyone know how to get my arrow and label to appear? I have been trying to add a left pointing arrow with the label ‘Mt Everest’ (both in white) to an extracted portion of a DEM. This is the code I have to get to this point, however the arrow and text does not appear when I run the script despite the rest of the properties working.
% extract_mt_everest_dem
%% Load the Mt Everest DEM data
[DEM, DEM_R] = geotiffread(‘aster_gdem_mosaic_UTM_crop.tif’);
DEM_info = geotiffinfo(‘aster_gdem_mosaic_UTM_crop.tif’);
% Manage NoData cells to make then NaN
DEM(DEM<=0) = NaN;
% Extract the lat and lon values for each cell within the data
% Create an index variable
idx = 1:length(DEM(:));
% Get the indices for each cell in the DEM
[I,J] = ind2sub(size(DEM),idx);
clear idx
% Convert the cell indices into projected coordinates
[x,y] = pix2map(DEM_info.RefMatrix, [I(:) J(:)]);
clear I J
% Reshape the x and y vectors to a grid the same size as the data
x_grid = reshape(x,size(DEM));
y_grid = reshape(y,size(DEM));
% Calculate the lat and lon values from the projected coordinates
[lat,lon] = projinv(DEM_info, x(:), y(:));
clear x y
% Reshape the lat and lon vectors to a grid the same size as the data
lon_grid = reshape(lon,size(DEM));
lat_grid = reshape(lat,size(DEM));
clear lat lon
%% Choose the point of interest (POI): Mt Everest
mt_everest = [27.986065 86.922623];
[mt_everest_X, mt_everest_Y] = projfwd(DEM_info, mt_everest(1), mt_everest(2));
%% Calculate distance from Aoraki Mt Cook
x_diff = x_grid – mt_everest_X;
y_diff = y_grid – mt_everest_Y;
% Find the indices of the bounding box of cells within a specified ditance of 20×20 km
max_dist = 10000; % (m)
row_start = find(abs(y_diff(:,1))<=max_dist, 1);
row_end = find(abs(y_diff(:,1))<=max_dist, 1, ‘last’);
col_start = find(abs(x_diff(1,:))<=max_dist, 1);
col_end = find(abs(x_diff(1,:))<=max_dist, 1, ‘last’);
% Extract the ROI from whole array using indexing
DEM_ROI = DEM(row_start:row_end,col_start:col_end);
lat_grid_ROI = lat_grid(row_start:row_end,col_start:col_end);
lon_grid_ROI = lon_grid(row_start:row_end,col_start:col_end);
%% Make a map to display the subset data
% Set the font size
font_size = 18;
set(0, ‘DefaultAxesFontSize’, font_size);
% Make the figure the size of the full screen
figure(‘units’,’normalized’,’outerposition’,[0 0 1 1]);
% Set the backrgound colour to white
set(gcf,’color’,’w’);
% Make a projected map axis using worldmap
ax = worldmap([nanmin(lat_grid_ROI(:)) nanmax(lat_grid_ROI(:))],[nanmin(lon_grid_ROI(:)) nanmax(lon_grid_ROI(:))]);
% Set axis font size to 18
setm(ax,’FontSize’,font_size);
hold on % to allow for overplotting
pcolorm(lat_grid_ROI, lon_grid_ROI, DEM_ROI);
cmap = landcolor; % get the landcolour colour map
colormap(cmap); % use the landcolour colour map
shadem; % shade the DEM to highlight the relief
hbar = colorbar(‘northoutside’,’fontsize’,font_size); % create a colourbar
set(get(hbar,’Xlabel’),’String’,’Elevation (m)’); % add a colourbar label
% Add a scale bar
scalebar(‘color’,[0 0 0],’FontSize’,font_size,’location’,’southeast’);
% Add a label showing Mt Everest with a left-pointing arrow
b = textm(mt_everest(1), mt_everest(2), 4000,’leftarrow Mt. Everest’,’FontSize’, font_size,’Color’,[1 1 1]);
uistack(b, ‘top’);
close
Anyone know how to get my arrow and label to appear? uistack, leftarrow, matlab, label MATLAB Answers — New Questions
