Документация

Чтобы определить, в каких файлах DTED вы нуждаетесь, используйте dteds функция, которая возвращает массив ячеек из символьных векторов:

dteds([41.1 43.9],[-71.9 -69.1])
ans = 
    '\DTED\W072\N41.dt0'
    '\DTED\W071\N41.dt0'
    '\DTED\W070\N41.dt0'
    '\DTED\W072\N42.dt0'
    '\DTED\W071\N42.dt0'
    '\DTED\W070\N42.dt0'
    '\DTED\W072\N43.dt0'
    '\DTED\W071\N43.dt0'
    '\DTED\W070\N43.dt0'

Отметьте три важных соображения об использовании файлов DTED:

Чтобы определить файлы USGS DEM, вы нуждаетесь, используете usgsdems функция:

usgsdems([41.1 43.9],[-71.9 -69.1])
ans = 
    'portland-w'
    'portland-e'
    'bath-w'
    'boston-w'
    'boston-e'
    'providence-w'
    'providence-e'
    'chatham-w'

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

Чтобы определить файлы GTOPO30, вы нуждаетесь, используете gtopo30s функция:

gtopo30s([41.1 43.9],[-71.9 -69.1])
ans = 
    'w100n90'

Задайте сферу интересов и определите файл, который будет получен:

latlim = [ 41.20  41.95];
lonlim = [-70.95 -70.10];

Чтобы определить, в каких файлах DTED вы нуждаетесь, используйте dteds функция, которая возвращает массив ячеек из символьных векторов:

dteds(latlim, lonlim)
ans = 
    'dted\w071\n41.dt0'

В этом примере необходим только один файл DTED, таким образом, ответ является односимвольным вектором. Для получения дополнительной информации о dteds функционируйте, смотрите, Идентифицируют Файлы DEM).

Если у вас нет CD-ROM, содержащего этот файл, загрузите его с источника, обозначенного в следующих технических замечаниях:

Найдите картографические данные онлайн

Исходные данные стали сжатым tar или архивом zip, перед использованием которого необходимо расшириться.

Используйте dted функция, чтобы создать сетку ландшафта и вектор ссылки в рабочей области в полном разрешении. Если больше чем один файл с именем DTED n41.dt0 существует на пути, вашей рабочей папкой должен быть /dted/w071 для того, чтобы быть уверенным тот dted находит правильный файл. Если файл не находится на пути, вам предлагают перейти к n41.dt0 файл dted функция:

samplefactor = 1;
[capeterrain, caperef] = dted('n41.dt0', ...
  samplefactor, latlim, lonlim);

Поскольку файлы DTED не содержат батиметрических глубин, вертикальных изменений уменьшения нуля немного, чтобы представить их с синим, когда палитра сбрасывается:

capeterrain(capeterrain == 0) = -1;

Используйте usamap создать пустую карту осей для области, заданной пределами широты и долготы:

figure;
ax = usamap(latlim,lonlim);

Считайте данные для области, заданной пределами широты и долготы от usastatehi файл форм:

capecoast = shaperead('usastatehi',...
    'UseGeoCoords', true,...
    'BoundingBox', [lonlim' latlim']);

Отобразите береговые линии на осях карты, который был создан с usamap:

geoshow(ax, capecoast, 'FaceColor', 'none');

В этой точке выглядит так карта:

Представьте вертикальные изменения и установите палитру соответственно:

meshm(capeterrain, caperef, size(capeterrain), capeterrain);
demcmap(capeterrain)

Получившаяся карта, показанная ниже, является окном на Кейп-Коде и иллюстрирует относительную грубость данных об уровне 0 DTED.

Чтобы следовать за этим осуществлением, необходимо получить необходимые файлы DTED из Интернета, как описано в Картографических данных Находки Онлайн, или от CD-ROM. Это дает к набору директорий, которые содержат следующие файлы:

/dted
    /w070 
        n41.avg 
        n41.dt0 
        n41.max 
        n41.min 
        n43.avg 
        n43.dt0 
        n43.max 
        n43.min 
    /w071 
        n41.avg 
        n41.dt0 
        n41.max 
        n41.min 
        n42.avg 
        n42.dt0 
        n42.max 
        n42.min 
        n43.avg 
        n43.dt0 
        n43.max 
        n43.min 
    /w072 
        n41.avg 
        n41.dt0 
        n41.max 
        n41.min 
        n42.avg 
        n42.dt0 
        n42.max 
        n42.min 
        n43.avg 
        n43.dt0 
        n43.max 
        n43.min 

Измените свою рабочую папку в папку, которая включает папку DTED верхнего уровня (который всегда называют dted).

Используйте dted функция, указывая, что папка в качестве первого аргумента:

latlim = [ 41.1  43.9];
lonlim = [-71.9 -69.1];
samplefactor = 5;
[capetopo,caperef] = dted(pwd, samplefactor, latlim, lonlim);

Демонстрационное факторное значение 5 указывает, что только каждая пятая ячейка данных, и в широте и в долготе, будет считана из исходного файла DTED. Можно выбрать большее значение, чтобы сохранить обработку памяти и скорости и отображение, за счет разрешения и точности. Размер вашего массива вертикального изменения (capetopo) будет обратно пропорционально квадрату демонстрационного фактора.

Когда DTEDs не содержат батиметрических глубин, перекод все нулевые вертикальные изменения к-1, чтобы позволить водным областям быть представленными правильно:

capetopo(capetopo==0)=-1;

Получите пределы широты и долготы сетки вертикального изменения; используйте их, чтобы чертить карту схемы области, чтобы ориентировать средство просмотра:

[latlim,lonlim] = limitm(capetopo,caperef);

figure;
ax = usamap(latlim,lonlim);
capecoast = shaperead('usastatehi',...
    'UseGeoCoords', true,...
    'BoundingBox', [lonlim' latlim']);
geoshow(ax,capecoast,'FaceColor','None');

Карта теперь выглядит так.

Представьте сетку вертикального изменения с meshm, и затем повторно окрасьте карту с demcmap отобразить гипсометрические цвета (оттенки вертикального изменения):

meshm(capetopo, caperef, size(capetopo), capetopo);
demcmap(capetopo)