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

Чтобы определить, в каких файлах 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)