import com.vividsolutions.jts.algorithm.locate.IndexedPointInAreaLocator; //导入依赖的package包/类
public void testGrid()
	{
		// Use fixed PM to try and get at least some points hitting the boundary
		GeometryFactory geomFactory = new GeometryFactory(pmFixed_1);
//		GeometryFactory geomFactory = new GeometryFactory();
		
		PerturbedGridPolygonBuilder gridBuilder = new PerturbedGridPolygonBuilder(geomFactory);
		gridBuilder.setNumLines(20);
		gridBuilder.setLineWidth(10.0);
    gridBuilder.setSeed(1185072199562L);
		Geometry area = gridBuilder.getGeometry();
		
//    PointInAreaLocator pia = new IndexedPointInAreaLocator(area); 
    PointOnGeometryLocator pia = new IndexedPointInAreaLocator(area); 

		PointInAreaStressTester gridTester = new PointInAreaStressTester(geomFactory, area);
		gridTester.setNumPoints(100000);
		gridTester.setPIA(pia);
		
		boolean isCorrect = gridTester.run();
		assertTrue(isCorrect);
	}
 

import com.vividsolutions.jts.algorithm.locate.IndexedPointInAreaLocator; //导入依赖的package包/类
/**
 * Determines the {@link Location} of the given {@link Coordinate}
 * in this geometry.
 *
 * @param p the point to test
 * @return the location of the point in the geometry
 */
public int locate(Coordinate pt) {
    if (this.parentGeom instanceof Polygonal && this.parentGeom.getNumGeometries() > 50) {
        // lazily init point locator
        if (this.areaPtLocator == null) {
            this.areaPtLocator = new IndexedPointInAreaLocator(this.parentGeom);
        }
        return this.areaPtLocator.locate(pt);
    }
    return this.ptLocator.locate(pt, this.parentGeom);
}
 

import com.vividsolutions.jts.algorithm.locate.IndexedPointInAreaLocator; //导入依赖的package包/类
public synchronized PointOnGeometryLocator getPointLocator() {
    if (this.pia == null) {
        this.pia = new IndexedPointInAreaLocator(this.getGeometry());
    }

    return this.pia;
}
 

import com.vividsolutions.jts.algorithm.locate.IndexedPointInAreaLocator; //导入依赖的package包/类
/**
 * Sets a polygonal mask.
 *
 * @param mask
 * @throws IllegalArgumentException if the mask is not polygonal
 */
public void setExtent(Geometry mask) {
    if (!(mask instanceof Polygonal)) {
        throw new IllegalArgumentException("Only polygonal extents are supported");
    }
    this.maskPoly = mask;
    this.setExtent(mask.getEnvelopeInternal());
    this.extentLocator = new IndexedPointInAreaLocator(mask);
}
 

import com.vividsolutions.jts.algorithm.locate.IndexedPointInAreaLocator; //导入依赖的package包/类
/**
 * Sets a polygonal mask.
 *
 * @param mask
 * @throws IllegalArgumentException if the mask is not polygonal
 */
public void setExtent(Geometry mask) {
    if (!(mask instanceof Polygonal))
        throw new IllegalArgumentException("Only polygonal extents are supported");
    this.maskPoly = mask;
    setExtent(mask.getEnvelopeInternal());
    extentLocator = new IndexedPointInAreaLocator(mask);
}
 

import com.vividsolutions.jts.algorithm.locate.IndexedPointInAreaLocator; //导入依赖的package包/类
/**
 * Sets a polygonal mask.
 * 
 * @param mask
 * @ throws IllegalArgumentException if the mask is not polygonal
 */
public void setExtent(Geometry mask)
{
	if (! (mask instanceof Polygonal))
		throw new IllegalArgumentException("Only polygonal extents are supported");
	this.maskPoly = mask;
	setExtent(mask.getEnvelopeInternal());
	extentLocator = new IndexedPointInAreaLocator(mask);
}
 

import com.vividsolutions.jts.algorithm.locate.IndexedPointInAreaLocator; //导入依赖的package包/类
/**
 * Sets a polygonal mask.
 * 
 * @param mask
 * @throws IllegalArgumentException if the mask is not polygonal
 */
public void setExtent(Geometry mask)
{
	if (! (mask instanceof Polygonal))
		throw new IllegalArgumentException("Only polygonal extents are supported");
	this.maskPoly = mask;
	setExtent(mask.getEnvelopeInternal());
	extentLocator = new IndexedPointInAreaLocator(mask);
}
 

import com.vividsolutions.jts.algorithm.locate.IndexedPointInAreaLocator; //导入依赖的package包/类
/**
 * rasterize the mask clipped with the Rectangle scaled back to full size with an offset onto a BufferedImage
 */
public int[] rasterMaskJTS(Rectangle rect, int offsetX, int offsetY, double scalingFactor) {

	// create the buffered image of the size of the Rectangle
    int[] mask = new int[rect.width* rect.height];
    GeometryFactory gf = new GeometryFactory();

    // define the clipping region in full scale
    Coordinate[] coords = new Coordinate[]{
        new Coordinate((int) (((double) rect.getMinX() / scalingFactor)), (int) (((double) rect.getMinY() / scalingFactor))),
        new Coordinate((int) (((double) rect.getMaxX() / scalingFactor)), (int) (((double) rect.getMinY() / scalingFactor))),
        new Coordinate((int) (((double) rect.getMaxX() / scalingFactor)), (int) (((double) rect.getMaxY() / scalingFactor))),
        new Coordinate((int) (((double) rect.getMinX() / scalingFactor)), (int) (((double) rect.getMaxY() / scalingFactor))),
        new Coordinate((int) (((double) rect.getMinX() / scalingFactor)), (int) (((double) rect.getMinY() / scalingFactor))),};

    Polygon geom = gf.createPolygon(gf.createLinearRing(coords));
    PreparedPolygon ppol=new PreparedPolygon(geom);

    int numPix=rect.width*rect.height;
    for (Geometry p : maskGeometries) {
        if (ppol.intersects(p)) {
        	Geometry pg=p.intersection(geom).buffer(0);
        	IndexedPointInAreaLocator locator=new IndexedPointInAreaLocator(pg);


        	int x=0;
        	int y=0;

        	for(int ii=0;ii<numPix;ii++){
        		if(ii%rect.width==0){
        			x=0;
        			y++;
        		}
    			//Point point=gf.createPoint(new Coordinate(rect.x+x,rect.y+y));
    			//PreparedPoint ppoint=new PreparedPoint(point);
    			//if(ppoint.within(pg)){
        		int loc=locator.locate(new Coordinate(rect.x+x,rect.y+y));
        		if(loc==Location.INTERIOR||loc==Location.BOUNDARY)
         		try{
         			mask[x]=1;
         		}catch(Exception e){
         			logger.warn(e.getMessage()+"  x:"+x+"  y:"+y);
         		}
    			}
         }
     //}
    }
    return mask;

}