How can I clip features (e.g. Rivers) to a smaller region than that map region?

Please, post this in the PyGMT category.