gmtselect treatment of line segments
|Target version:||Candidate for next minor release|
When a multi-segment file with line-segments is processed, what should happen if a part of a segment is inside the selected region/polygon:
1) Pass the entire segment
2) Pass just the points of the segment that are inside
3) Clip segment at boundary, possibly splitting segments into many segments
4) All points must be inside for segment to pass.
Seems we need a new option to control what is meant by "inside" for lines, and possibly the same for polygons?
#1 Updated by Paul about 4 years ago
- File selectlines.pdf added
- Status changed from New to Resolved
- Assignee set to Paul
- Target version set to Candidate for next minor release
- % Done changed from 0 to 100
Resolution to this old issue: Since records in gmtselect are read one at the time it is not possible to do what was requested as options (1) and (4) above without a complete i/o rewrite. Option (2) remains the default action, shown in the bottom plot in the attached graph. Basically, input data are considered as individual points. Option (3) however can instead be achieved via gmtspatial -T (add -Fl for lines). Since we wish gmtselect to handle huge datasets we will retain the idea of reading record by record. In v 5.2, r13275.