Ignore:
Timestamp:
Jun 19, 2006, 10:30:57 AM (19 years ago)
Author:
duncan
Message:

added function to pmesh.mesh for adding vertices by passing in a geospatial instance.

Location:
inundation/geospatial_data
Files:
2 edited

Legend:

Unmodified
Added
Removed
  • inundation/geospatial_data/geospatial_data.py

    r3154 r3177  
    175175        # represent the new geo-ref
    176176        if  self.geo_reference is not None:
    177             #print "self.geo_reference",self.geo_reference
    178             #print "geo_reference",geo_reference
    179177            #FIXME: Maybe put out a warning here...
    180178            self.data_points = self.get_data_points \
    181179                               (geo_reference=geo_reference)
    182             #self.data_points = geo_reference.change_points_geo_ref \
    183             #                   (self.data_points,
    184             #                    self.geo_reference)
    185180           
    186181        self.geo_reference = geo_reference
     
    205200        returned coordinates are relative to the internal georeference's
    206201        xll and yll corners.
     202
     203        If a geo_reference is passed the points are returned relative
     204        to that geo_reference.
    207205
    208206        Default: absolute is True.
  • inundation/geospatial_data/test_geospatial_data.py

    r3154 r3177  
    921921        '''
    922922
    923         #FIXME (Ole): I changed the georeference arbitrarily, but the test still passes.
    924        
    925923        # create files
    926924        att_dict1 = {}
     
    934932        att_dict2['elevation'] = array([1.0, 15.0, 1.4])
    935933        att_dict2['brightness'] = array([14.0, 1.0, -12.4])
    936         geo_reference2 = Geo_reference(56, 1.0, 2.0) #FIXME (Ole): I changed this, why does it still pass
    937         #OK - it fails now, after the fix revealed by test_create_from_pts_file_with_geo')
     934        geo_reference2 = Geo_reference(56, 1.0, 2.0)
    938935
    939936        G1 = Geospatial_data(pointlist1, att_dict1, geo_reference1)
Note: See TracChangeset for help on using the changeset viewer.