Is there a contourc equivalent for contourm?

Hello,

I would like to compute contours on geotiffs without producing a plot. For regular rectangluar data 'contour' will compute the contours and plot them, and 'contourc' function will just compute the contours and return the contour matrix without plotting. It appears the geographic contour function 'contourm' works like 'contour' and produces a plot, I would like to get the contours without generating a plot.

Right now my solution is this:

set(0,'DefaultFigureVisible','off');  
altContour = contourm(geoData,geoRef,'LevelList',geoAltitude,'Visible','off');
close(gcf);
set(0,'DefaultFigureVisible','on');

But this is still incurring overhead penalties for setting up the figure... I would ideally like a function like 'contourmc'. Is there some function, even internal or undocumented, that I can use to achieve this?

Answer 1

[Pavan Sahith]

Hi Mark,

I understand you would like to compute contours on geotiffs without producing a plot.

You can use the ‘contourm’ function in this way to achieve the contour matrix without plotting.

% Example data
[X, Y, Z] = peaks;
% Define the desired contour levels
contourLevels = [-5, -3, -1, 1, 3, 5];

% Create a temporary invisible figure to capture the contour data
hFig = figure('Visible', 'off');
% Use contourm to capture the contour data
[ccontourData] = contourm(Y, X, Z, contourLevels)

ccontourData = 2×504

-5.0000 0.6250 0.5522 0.5000 0.3750 0.2500 0.1250 0.0376 0 -0.1182 -0.1250 -0.1893 -0.1978 -0.1335 -0.1250 0 0.1250 0.2500 0.3750 0.5000 0.5191 0.6240 0.6250 0.6594 0.6424 0.6250 -3.0000 -1.3750 -1.4003 -1.3938 25.0000 -1.4696 -1.3750 -1.3385 -1.2959 -1.2973 -1.3326 -1.3750 -1.3997 -1.5000 -1.5106 -1.6250 -1.7500 -1.8750 -1.8831 -1.9593 -1.9912 -1.9893 -1.9564 -1.8900 -1.8750 -1.7500 -1.7472 -1.6250 -1.5000 -1.4696 7.0000 0.2817 0.2500 0.1250

% Close the temporary figure
close(hFig);
% contourData contains the contour lines and levels without displaying the plot

Please refer to the MathWorks Documentation to know more about

‘contourm’ - https://in.mathworks.com/help/map/ref/contourm.html

Hope it helps .