MATLAB: Is there a command in the Mapping Toolbox 2.7.1 (R2008b) that will add filled contours to a map

Mapping Toolbox

I am using a map created with grid2image(map,refvec);

I want to create a filled contour on my map using the output of

[visZ,visrefvec] = viewshed(Z,refvec,lat1,lon1)

If I use contourm(visZ,visrefvec) I will see the contour lines on the map but the contour lines are not filled.

Also, contourf(visZ) to create a filled contour, does not seem to work.

Best Answer

There is no shipped function which allows this functionality in MATLAB 7.9 (R2008b). As a workaround you can use the CONTOURF function, along with its 'Fill' property.

Related Solutions

MATLAB: What this error means …and what i shd do to sort it out???

It looks like you don't have videoinput function in your installation. You need to have Image Acquisition Toolbox to have that function. You can do

>> ver

to verify what you have.

MATLAB: How to overlay VMAP0 data on a topographic relief surface using the Mapping Toolbox

To achieve this, you need to use interpolation to determine the altitude values to place the graphics objects created from VMAP0 data to sit on a topographic surface. More specifically, you need to use the latitude, longitude, and altitude data from the topographic grid as the basis to calculate the altitude values for the given latitude and longitude values from the VMAP0 data.

% here, Z is an array of altitude values.
[latsize,lonsize] = size(Z);
latlocs = linspace(latlim(1),latlim(2),latsize);
lonlocs = linspace(lonlim(1),lonlim(2),lonsize);
%VMAPlines is a structure of line information, created from VMAP0 data
for i=1:length(VMAPlines)
      str = VMAPlines(i).tag;
      if ~isempty(strmatch('Accurate;',str))
         % since the segments of data are separated by NaNs, 
         % you will need to interpolate each segment separately
          NaNLocs = find(isnan(VMAPlines(i).long));
          for n = 1:length(NaNLocs)-1
              % get the segment's latitude and longitude data
              tmplat = VMAPlines(i).lat(NaNLocs(n)+1:NaNLocs(n+1)-1);
              tmplon = VMAPlines(i).long(NaNLocs(n)+1:NaNLocs(n+1)-1);
             % interpolate the altitude
              tmpalt = interp2(lonlocs,latlocs,Z,tmplon,tmplat);
             % add the segment's altitude info, separated by NaNs
              VMAPlines(i).alt = [VMAPlines(i).alt; NaN; tmpalt(:)]; 
          end
          % display the line
          hTRline(i) = displaym(VMAPlines(i));       
      end
end

You can use this approach to "drape" vector data other than VMAP0 over terrain by substituting the appropriate data structures for VMAPlines. Note that there may be locations where the vectors become hidden for short distances.

Best Answer

Related Solutions

MATLAB: What this error means …and what i shd do to sort it out???

MATLAB: How to overlay VMAP0 data on a topographic relief surface using the Mapping Toolbox

Related Question