Небольшие круги от центра и периметра
[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' = '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')
