Наложите Бостонские дороги на ортофотографии. Обратите внимание на то, что 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)

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

roads = readgeotable("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);

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

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

Импортируйте и отобразите многоугольник, представляющий водоем в середине ортофотографии.

Импортируйте файл форм, содержащий гидрографические данные для Согласия как геопространственная таблица. Создайте подтаблицу, которая содержит 14-ю строку таблицы, которая является строкой, содержащей водоем. Проверьте, что подтаблица содержит данные о многоугольнике путем запроса Shape переменная таблицы.

hydro = readgeotable("concord_hydro_area.shp");
pond = hydro(14,:);
pond.Shape

ans = 
  mappolyshape with properties:

              NumRegions: 1
                NumHoles: 3
                Geometry: "polygon"
    CoordinateSystemType: "planar"
            ProjectedCRS: [1x1 projcrs]

Отобразите многоугольник по ортофотографии.

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');