I am looking to calculate the geographic area of an intersection of a GeoDisk with a geographic entity such as a county.
I was hoping for a simple construct like:
state = Entity["Country", "PuertoRico"];
disk = GeoDisk[{18, -67}, Quantity[50, "Kilometers"]];
Print[ GeoGraphics[{disk, state}, GeoScaleBar -> True] ];
Print[" The overlap area is: ",
GeoArea[RegionIntersection[ Polygon[ state ], disk] ] ]
Unfortunately, the RegionIntersection function doesn't handle GeoGraphics structures, and I can find no such thing as a "GeoRegionIntersection" native function. Is there such a function hidden, but in trial, perhaps?
Alternatively, is there a reasonably straightforward way to convert the GeoDisk to a Polygon and do the calculation with something like:
Area[RegionIntersection[ Polygon[state], function[disk] ] ]?
I hope to apply this to several hundred such disks and am also interested in doing this for about a hundred counties for which I will take a weighted sum of population per county.
Thank you in advance, Steve