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