Стандартное отклонение для географических точек
[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