Project

General

Profile

« Previous | Next » 

Revision 3abda63e

Added by Jim Regetz over 12 years ago

  • ID 3abda63e99672025aef15bedcaf37a37e397c7bc

now actually refining before smoothing (GRASS version)

View differences:

terrain/research/oregon/arcgis/v2/multiscalesmooth.py
219 219 
    tmp_rast.add('v.smooth')
220 220 
    tmp_rast.remove('v%d' % NUM_LEVELS)
221 221 

  		




  222
  
  
    
    # smooth, refine and combine each layer in turn


  		




  
  222
  
    
    # refine, smooth, and combine each layer in turn


  		




  223
  223
  
    
    for j in reversed(range(NUM_LEVELS)):


  		




  224
  224
  
    

  		




  225
  225
  
    
        i = j + 1


  		




  226
  226
  
    
        gs.message('Refining from %d to %d' % (i, j), flag='i')


  		




  227
  227
  
    

  		




  228
  228
  
    
        refine_region()


  		




  
  229
  
    
        gs.run_command('r.resample', input=smooth, output=smooth,


  		




  
  230
  
    
            overwrite=True, quiet=True)


  		




  
  231
  
    
        gs.run_command('r.resample', input='v.smooth', output='v.smooth',


  		




  
  232
  
    
            overwrite=True, quiet=True)


  		




  229
  233
  
    

  		




  230
  234
  
    
        # create smoothed higher resolution versions of z and v


  		




  231
  235
  
    
        #TODO: is this the same as circle with radius 2 in arc?


  		












Also available in:  Unified diff