Наложите дороги Бостона на ортофото. Обратите внимание, что mapshow рисует новый слой в осях, а не заменяет его содержимое.

Отобразите изображение.

figure
mapshow boston.tif
axis image off manual

Преобразуйте дороги Бостона в модули измерения футов и наложения на ортофото.

S = shaperead('boston_roads.shp');
surveyFeetPerMeter = unitsratio('sf','meter');
x = surveyFeetPerMeter * [S.X];
y = surveyFeetPerMeter * [S.Y];
mapshow(x,y)

Считайте векторные данные и отобразите их с помощью пунктирной линии.

roads = shaperead('boston_roads.shp');
figure
mapshow(roads,'LineStyle',':');

Создайте спецификацию символов, чтобы различать различные типы дорог. Например, можно скрыть очень незначительные дороги (CLASS = 6), отключив их видимость и сделав основные дороги (CLASS = 1-4) более видимыми, увеличив их ширины линии. Эта спецификация символа также использует цвет для различения типов дорог.

roadspec = makesymbolspec('Line',...
                          {'ADMIN_TYPE',0,'Color','cyan'}, ...
                          {'ADMIN_TYPE',3,'Color','red'},...
                          {'CLASS',6,'Visible','off'},...
                          {'CLASS',[1 4],'LineWidth',2});

Отобразите векторные данные с помощью спецификации символов.

figure
mapshow('boston_roads.shp','SymbolSpec',roadspec);

Создайте спецификацию символа и укажите цвет по умолчанию, используемый для линий. Как видно из предыдущего примера, значение по умолчанию является синим. Этот пример устанавливает значение по умолчанию черное.

roadspec = makesymbolspec('Line',...
                              {'Default', 'Color', 'black'}, ...
                              {'ADMIN_TYPE',0,'Color','c'}, ...
                              {'ADMIN_TYPE',3,'Color','r'},...
                              {'CLASS',6,'Visible','off'},...
                              {'CLASS',[1 4],'LineWidth',2});

Отобразите векторные данные, используя спецификацию символа. Обратите внимание, что основные дороги, отображенные синим цветом в предыдущем примере, теперь являются черными.

figure
mapshow('boston_roads.shp','SymbolSpec',roadspec);

Создайте спецификацию символов, задав различные свойства.

roadspec = makesymbolspec('Line',...
                              {'Default', 'Color', 'yellow'}, ...
                              {'ADMIN_TYPE',0,'Color','c'}, ...
                              {'ADMIN_TYPE',3,'Color','r'},...
                              {'CLASS',6,'Visible','off'},...
                              {'CLASS',[1 4],'LineWidth',2});

Отобразите векторные данные, задав цвет в командной строке.

figure
mapshow('boston_roads.shp', 'Color', 'black', 'SymbolSpec', roadspec);

Импортируйте ортофото Concord, MA, вместе с объектом ссылки камер карты и палитрой. Отобразите ортофото с помощью mapshow функция.

[ortho,R,cmap] = readgeoraster('concord_ortho_w.tif');
mapshow(ortho,cmap,R) 

Импортируйте файл shapefile, содержащий координаты береговой линии пруда в середине ортофото. Проверьте, что данные представляют многоугольник, запросив их Geometry поле.

pond = shaperead('concord_hydro_area.shp','RecordNumbers',14);
pond.Geometry

ans = 
'Polygon'

Отобразите многоугольник на ортофото. Внешний контур очерчивает водоем, а внутренние контуры - острова.

mapshow(pond,'FaceColor','c')

Ознакомьтесь с данными о местности SDTS для горы Вашингтон. Получите такую информацию, как индикаторы недостающих данных, используя georasterinfo. Замените отсутствующие данные на NaN значения с использованием standardizeMissing.

[Z,R] = readgeoraster('sdts/9129CATD.DDF','OutputType','double');
info = georasterinfo('sdts/9129CATD.DDF');
Z = standardizeMissing(Z,info.MissingDataIndicator);

Просмотрите данные местности как mesh. Применить палитру, подходящую для данных местности, используя demcmap.

figure
mapshow(Z,R,'DisplayType','mesh');
demcmap(Z)

Просмотрите данные местности как поверхность.

figure
mapshow(Z,R,'DisplayType','surface');
demcmap(Z)

Просмотрите данные местности как 3-D поверхность.

view(3);
axis normal

Ознакомьтесь с файлами данных о местности для горы Вашингтон и горы Дартмут. Чтобы построить график данных как поверхности, используя mapshowРастр должен иметь тип single или double. Укажите тип данных для растра, используя 'OutputType' Пара "имя-значение".

[ZWash,RWash] = readgeoraster('MtWashington-ft.grd','OutputType','double');
[ZDart,RDart] = readgeoraster('MountDartmouth-ft.grd','OutputType','double');

Найдите отсутствующие данные с помощью georasterinfo функция. Функция возвращает объект с MissingDataIndicator свойство, которое указывает, какое значение представляет отсутствующие данные. Замените отсутствующие данные на NaN значения с использованием standardizeMissing функция.

infoWash = georasterinfo('MtWashington-ft.grd');
ZWash = standardizeMissing(ZWash,infoWash.MissingDataIndicator);

infoDart = georasterinfo('MountDartmouth-ft.grd');
ZDart = standardizeMissing(ZDart,infoDart.MissingDataIndicator);

Отображение данных местности под контурными линиями и метками путем определения 'ZData' Пара "имя-значение" как матрица нулей. Применить палитру, подходящую для данных местности, используя demcmap.

hold on
mapshow(ZWash,RWash,'DisplayType','surface','ZData',zeros(RWash.RasterSize))
mapshow(ZDart,RDart,'DisplayType','surface','ZData',zeros(RDart.RasterSize))
axis equal
demcmap(ZWash)

Наложение черных контурных линий и меток.

mapshow(ZWash,RWash,'DisplayType','contour',...
    'LineColor','k','ShowText','on');
mapshow(ZDart,RDart,'DisplayType','contour',...
    'LineColor','k','ShowText','on');