Hi Jeff,
I finally had a chance to try this. I can't get it to work but I think
I'm close - for some reason, the way I'm creating the geos polygons
seems to always intersect the boundary polygon. It's hard to think of a
good minimal example for this so I've attached an example that
illustrate
Thank you Jeff. I'll try out this solution.
Gary.
Jeff Whitaker wrote:
> Gary: You might be able to use the _geoslib module to compute the
> intersections of those polygons with the map boundary. I do a similar
> thing with the coastline polygons in the _readboundarydata function.
> The
Gary Ruben wrote:
> I'm plotting a coverage map of a sphere using the Mollweide plot in
> basemap. The attachment is an example that is produced by sending an
> array of polygons (one polygon per row described as four corners, one
> per column) described using polar (theta) and azimuthal (phi) angl
