This might be a pretty a pretty simple question, but still
in the process of learning about GeoTools – so rather
than agonizing, I thought I’d ask…

I have a polygon shapefile that
imports nicely into geotools as a FeatureResults
containing MultiPolygons (they are counties, so they
are irregular polygons). I also
have an ArrayList of x and y co-ordinates as double
values.

What is the best way of seeing which polygons contain the
set of points?

I know that MultiPolygon has < contains(Geometry) >, and that Polygon has < contains(double
x, double y) >. What would be
the best way for packaging an array of x y double values into a Geometry, or
alternatively is there a way to somehow make use of Polygon’s < contains(double
x, double y) > method?

Thanks,

J