Небольшие круги от центра и периметра
[lat,lon] = scircle2(lat1,lon1,lat2,lon2)
[lat,lon] = scircle2(lat1,lon1,lat2,lon2,ellipsoid)
[lat,lon] = scircle2(lat1,lon1,lat2,lon2,units)
[lat,lon] = scircle2(lat1,lon1,lat2,lon2,ellipsoid,units)
[lat,lon] = scircle2(lat1,lon1,lat2,lon2,ellipsoid,units,npts)
[lat,lon] = scircle2(track,...)
mat = scircle2(...)
mat = [lat lon]
[lat,lon] = scircle2(lat1,lon1,lat2,lon2) вычисляет маленькие круги (на сфере) с центрами в точках lat1,lon1 и точки на кругах в lat2,lon2. Входы могут быть скалярными или столбцевыми векторами.
[lat,lon] = scircle2(lat1,lon1,lat2,lon2,ellipsoid) вычисляет малую окружность на эллипсоиде, заданный вход ellipsoid, а не путем принятия сферы. ellipsoid является referenceSphere, referenceEllipsoid, или oblateSpheroid объект или вектор формы [semimajor_axis eccentricity]. По умолчанию это сфера единичного радиуса.
[lat,lon] = scircle2(lat1,lon1,lat2,lon2,units) и
[lat,lon] = scircle2(lat1,lon1,lat2,lon2,ellipsoid,units) являются допустимыми формами вызова, которые используют входные units для определения угловых модулей входов и выходов. Если вы опускаете units, 'degrees' принято.
[lat,lon] = scircle2(lat1,lon1,lat2,lon2,ellipsoid,units,npts) использует скалярный входной npts для определения числа точек на одну вычисленную дорожку. Значение по умолчанию npts равен 100.
[lat,lon] = scircle2(track,...) использует track, чтобы задать большой круг или радиус линии ручьи. Если track' = 'gc', затем вычисляются небольшие круги. Если track = 'rh'затем вычисляют круги с радиусами постоянного расстояния линии револьвера. Если вы опускаете дорожку, 'gc' принято.
mat = scircle2(...) возвращает один выходной аргумент, где mat = [lat lon]. Это полезно, если вычислена одна окружность.
Несколько кругов могут быть заданы из одной центральной точки путем обеспечения скалярной lat1,lon1 входы и векторы-столбцы для точек на окружности, lat2,lon2.
Постройте график локуса всех точек на том же расстоянии от Нью-Дели, что и Катманду.
axesm('mercator','MapLatlimit',[0 40],'MapLonLimit',[60 110]); load coastlines % For reference plotm(coastlat,coastlon,'k'); % New Delhi lat1 = 29; lon1 = 77.5; % Kathmandu lat2 = 27.6; lon2 = 85.5; % Plot the cities plotm([lat1 lat2],[lon1 lon2],'b*') [latc,lonc] = scircle2(lat1,lon1,lat2,lon2); plotm(latc,lonc,'b')
