Стандартное отклонение для географических точек
[latdev,londev] = stdm(lat,lon)
[latdev,londev] = stdm(lat,lon,ellipsoid)
[latdev,londev] = stdm(lat,lon,units)
[latdev,londev] = stdm(lat,lon) возвращает векторы строк отклонений географических стандартов широты и долготы для точек данных, указанных столбцами lat и lon.
[latdev,londev] = stdm(lat,lon,ellipsoid) определяет форму Земли, которая будет использоваться ellipsoid, который может быть referenceSphere, referenceEllipsoid, или oblateSpheroid объект или вектор формы [semimajor_axis eccentricity]. Эллипсоид по умолчанию является единичной сферой. Выходные измерения выражаются в единицах расстояния ellipsoid вектор.
[latdev,londev] = stdm(lat,lon, указывает угловые единицы данных. При пропуске единиц измерения units)'degrees' предполагается. Выходные измерения относятся к этим единицам измерения (как расстояние длины дуги).
Если используется один выходной аргумент, то geodevs = [latdev longdev]. Это особенно полезно, если оригинал lat и lon входами являются векторы столбцов.
Определение отклонений географических данных по широте и долготе сложнее простых отклонений суммы квадратов от средних значений данных. Для широтного отклонения выполняется прямой расчет углового стандартного отклонения по среднему географическому значению, рассчитанному по meanm. Для долгот аналогичный расчет выполняется на основе отклонения данных, а не углового отклонения. См. раздел Географическая статистика местоположений точек на сфере в Руководстве пользователя Mapping Toolbox™.
Создание списков широт и долгот с помощью worldcities набор данных и получение стандартного отклонения расстояния для группы (сравните с примером для stdist):
cities = shaperead('worldcities.shp', 'UseGeoCoords', true);
Paris = strcmp('Paris',{cities(:).Name});
London = strcmp('London',{cities(:).Name});
Rome = strcmp('Rome',{cities(:).Name});
Madrid = strcmp('Madrid',{cities(:).Name});
Berlin = strcmp('Berlin',{cities(:).Name});
Athens = strcmp('Athens',{cities(:).Name});
lat = [cities(Paris).Lat cities(London).Lat...
cities(Rome).Lat cities(Madrid).Lat...
cities(Berlin).Lat cities(Athens).Lat]
lon = [cities(Paris).Lon cities(London).Lon...
cities(Rome).Lon cities(Madrid).Lon...
cities(Berlin).Lon cities(Athens).Lon]
[latstd,lonstd]=stdm(lat,lon)
lat =
48.8708 51.5188 41.9260 40.4312 52.4257 38.0164
lon =
2.4131 -0.1300 12.4951 -3.6788 13.0802 23.5183
latstd =
2.7640
lonstd =
68.7772