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