Changeset 6533 for branches/numpy/anuga/geospatial_data

Timestamp:

Mar 17, 2009, 4:02:54 PM (16 years ago)

Author:

rwilson

Message:

Revert back to 6481, prior to auto-merge of trunk and numpy branch.

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)