Changeset 6517 for branches/numpy/anuga/geospatial_data/geospatial_data.py

Timestamp:

Mar 16, 2009, 11:06:22 AM (14 years ago)

Author:

rwilson

Message:

Hand-merged recent changes in main trunk. Still work to be done in shallow_water.

File:

• branches/numpy/anuga/geospatial_data/geospatial_data.py (modified) (1 diff)

branches/numpy/anuga/geospatial_data/geospatial_data.py

@@ -1806 +1806 @@
     points = G_small.get_data_points()
 
+    #FIXME: Remove points outside boundary polygon
+#    print 'new point',len(points)
+#
+#    new_points=[]
+#    new_points=array([],typecode=Float)
+#    new_points=resize(new_points,(len(points),2))
+#    print "BOUNDARY", boundary_poly
+#    for i,point in enumerate(points):
+#        if is_inside_polygon(point,boundary_poly, verbose=True):
+#            new_points[i] = point
+#            print"WOW",i,new_points[i]
+#    points = new_points
+
     if verbose: print "Number of points in sample to compare: ", len(points)