void RunRelateTest(String wkt1, String wkt2, IBoundaryNodeRule bnRule, String expectedIM)
{
IGeometry g1 = rdr.Read(wkt1);
IGeometry g2 = rdr.Read(wkt2);
IntersectionMatrix im = RelateOp.Relate(g1, g2, bnRule);
String imStr = im.ToString();
Console.WriteLine(imStr);
Assert.IsTrue(im.Matches(expectedIM));
}