Documentation |
Intersection points for great and small circle pairs
[newlat,newlong] = gcxsc(gclat,gclong,gcaz,sclat,sclong,scrange)
[newlat,newlong] = gcxsc(...,units)
[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'.
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.