maptriml

Линии для обрезки к четырехугольнику долготы широты

Синтаксис

[lat,lon] = maptriml(lat0,lon0,latlim,lonlim)

Описание

[lat,lon] = maptriml(lat0,lon0,latlim,lonlim) возвращается отфильтровал NaN- разграниченные векторные наборы данных карты, от которых были отброшены все точки, лежащие вне желаемых пределов широты и долготы. Эти пределы заданы двухэлементными векторами latlim и lonlim, которые имеют форму [south-limit north-limit] и [west-limit east-limit], соответственно.

Примеры

свернуть все

Для векторных данных весьма обычно расширить вне географической области в настоящее время интереса. Например, вы можете иметь данные о береговой линии для всего мира, но интересуетесь отображением Австралии только. В этом и других ситуациях, вы можете хотеть устранить ненужные данные из рабочей области и из вычислений для того, чтобы сохранить память или ускорять обработку и отображение. Можно обрезать данные о линии путем простого удаления точек вне необходимой области путем усечения линий в системе координат карты или в некоторую другую заданную область. maptriml функция возвращает переменные, содержащие только те точки, которые лежат в заданной области.

Загрузите coastlines Matfile. Этот файл содержит данные для всего мира.

load coastlines

Задайте необходимую область, сосредоточенную на Австралии.

latlim = [-50 0]; 
lonlim = [105 160];

Удалите все данные о линии вне этих пределов, произведя векторы линии, при помощи maptriml функция.

[linelat,linelon] = maptriml(coastlat,coastlon,latlim,lonlim);

Исследуйте переменные, чтобы видеть, как много данных уменьшалось. Отсеченные данные являются только на 10% столь же большими, как исходный набор данных.

whos
  Name             Size            Bytes  Class     Attributes

  coastlat      9865x1             78920  double              
  coastlon      9865x1             78920  double              
  latlim           1x2                16  double              
  linelat        977x1              7816  double              
  linelon        977x1              7816  double              
  lonlim           1x2                16  double              

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

axesm('MapProjection', 'miller', 'Frame', 'on',...
'FlatLimit', latlim, 'FlonLimit', lonlim)
geoshow(linelat,linelon)

Figure contains an axes object. The axes object contains 2 objects of type patch, line.

Смотрите также

|

Представлено до R2006a