Changeset 6517 for branches/numpy/anuga/geospatial_data/geospatial_data.py
- Timestamp:
- Mar 16, 2009, 11:06:22 AM (14 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
branches/numpy/anuga/geospatial_data/geospatial_data.py
r6441 r6517 1806 1806 points = G_small.get_data_points() 1807 1807 1808 #FIXME: Remove points outside boundary polygon 1809 # print 'new point',len(points) 1810 # 1811 # new_points=[] 1812 # new_points=array([],typecode=Float) 1813 # new_points=resize(new_points,(len(points),2)) 1814 # print "BOUNDARY", boundary_poly 1815 # for i,point in enumerate(points): 1816 # if is_inside_polygon(point,boundary_poly, verbose=True): 1817 # new_points[i] = point 1818 # print"WOW",i,new_points[i] 1819 # points = new_points 1820 1808 1821 if verbose: print "Number of points in sample to compare: ", len(points) 1809 1822
Note: See TracChangeset
for help on using the changeset viewer.