Mapping Toolbox

gcxsc

Provide intersection coordinates for great circles paired with small circles

Syntax

• [newlat,newlong] = gcxsc(gclat,gclong,gcaz,sclat,sclong,scrange)
  [newlat,newlong] = gcxsc(gclat,gclong,gcaz,...
                           sclat,sclong,scrange,units)
  

Description

For a pairing of a great circle with a small circle, there are four possible intersection conditions: the circles are identical (possible because great circles are a subset of small circles), they do not intersect, they are tangent to each other (the small circle interior to the great circle) and hence they intersect once, or they intersect twice.

Great circle notation consists of a point on the great circle and the azimuth at that point along which the great circle proceeds.

Small circle notation consists of a center point and a radius in units of angular arc length.

[newlat,newlong] = gcxsc(gclat,gclong,gcaz,sclat,sclong,scrange) returns the points of intersection of a great circle in great circle notation followed by a small circle in small circle notation. For multiple pairings, the inputs must be column vectors. The results are two-column matrices with the coordinates of the intersection points. If the circles do not intersect, or are identical, two NaNs are returned and a warning is displayed. If the two circles are tangent, the single intersection point is repeated twice.

[newlat,newlong] = gcxsc(...,units) specifies the standard angle unit string. The default value is 'degrees'.

Examples

Given a great circle passing through (43ºN,0º) and proceeding on an azimuth of 10º, where does it intersect with a small circle centered at (47ºN,3ºE) with an arc length radius of 12º?

• [newlat,newlong] = gcxsc(43,0,10,47,3,12)
  newlat =
     35.5068   58.9143
  newlong =
     -1.6159    5.4039
  

See Also

gc2sc	Convert great circle to small circle notation
gcxgc scxsc rhxrh crossfix polyxpoly	Other intersection functions

gcxgc

geod2aut