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