exponenta event banner

геопики

Создание набора синтетических данных для сферы

Описание

пример

Z = geopeaks(lat,lon) оценивает «пикоподобную» функцию в определенных широтах и долготах на поверхности сферы, возвращая синтетический набор данных Z. Функция непрерывна и плавна во всех точках, включая полюса. Напоминает MATLAB peaks функция, geopeaks мягко волнообразно между значениями -10 и 8, где около полудюжины местных экстремалов.

пример

Z = geopeaks(R) оценивает geopeaks функция в центрах ячеек или выборочных точках проводки, определенных объектом географической растровой ссылки, R.

пример

Z = geopeaks(___,spheroid) вычисляет функцию на определенном сфероиде. Выбор сфероида имеет очень мало значения. Этот вариант существует главным образом для поддержки формального тестирования. Если не указать spheroid и GeographicCRS имущество R не пуст, то geopeaks использует сфероид, содержащийся в Spheroid имущества geocrs объект в GeographicCRS имущество R.

Примеры

свернуть все

Определите значения широты и долготы по меридиану, включающему Париж, Франция.

lon = dms2degrees([2 21 3]);
lat = -90:0.5:90;

Создайте набор данных, указав сфероид.

z = geopeaks(lat, lon, wgs84Ellipsoid);

Просмотрите профиль.

figure
plot(lat,z)
xlabel('Latitude')

Figure contains an axes. The axes contains an object of type line.

Создайте объект ссылки на растр для сетки проводок 181 на 361.

latlim = [-90 90];
lonlim = [-180 180];
sampleSpacing = 1;
R = georefpostings(latlim,lonlim,sampleSpacing,sampleSpacing)
R = 
  GeographicPostingsReference with properties:

              LatitudeLimits: [-90 90]
             LongitudeLimits: [-180 180]
                  RasterSize: [181 361]
        RasterInterpretation: 'postings'
            ColumnsStartFrom: 'south'
               RowsStartFrom: 'west'
     SampleSpacingInLatitude: 1
    SampleSpacingInLongitude: 1
      RasterExtentInLatitude: 180
     RasterExtentInLongitude: 360
            XIntrinsicLimits: [1 361]
            YIntrinsicLimits: [1 181]
        CoordinateSystemType: 'geographic'
               GeographicCRS: []
                   AngleUnit: 'degree'


Создайте набор данных.

Z = geopeaks(R);

Отображение результирующего набора данных на карте.

figure
worldmap world
geoshow(Z,R,'DisplayType','surface','CData',Z,'ZData',zeros(size(Z)))
load coastlines
geoshow(coastlat,coastlon,'Color','k')
colorbar

Входные аргументы

свернуть все

Геодезическая широта одной или нескольких точек, заданная как скалярное значение, вектор или матрица. Значения должны быть в градусах.

lat входной аргумент должен соответствовать lon входной аргумент в размере, если одно из значений не является скалярным (в этом случае оно будет расширяться по размеру в соответствии с другим), или lat является вектором столбца и lon является вектором строки (они расширяются, образуя сетку широты-долготы).

Пример: lat = -90:0.5:90

Типы данных: single | double

Геодезическая долгота одной или нескольких точек, заданная как скалярное значение, вектор или матрица. Значения должны быть в градусах.

lon входной аргумент должен соответствовать lat входной аргумент в размере, если одно из значений не является скалярным (в этом случае оно будет расширяться по размеру в соответствии с другим), или lon является вектором столбца и lat является вектором строки (они расширяются, образуя сетку широты-долготы).

Пример: lon = -180:0.5:180

Типы данных: single | double

Географический растр, указанный как GeographicCellsReference или GeographicPostingsReference объект. Географический растр сохраняет широту и долготу точек.

Пример: R = georefcells([0 80], [-140 60], 0.25, 0.25)

Ссылочный сфероид, указанный как referenceEllipsoid, oblateSpheroid, или referenceSphere объект.

Пример: spheroid = referenceEllipsoid('GRS80')

Выходные аргументы

свернуть все

Синтетический набор данных, возвращаемый как скалярное значение, вектор или матрица класса single или double, в зависимости от класса входных данных. Функция оценивается в каждом элементе lat и lon (после расширения, как указано выше), или в каждом центре ячейки или точке проводки, определенной R.

См. также

Представлен в R2015b