Changes between Version 2 and Version 3 of GSoC/2014/ImplementationOf3DRasterFlowLine


Ignore:
Timestamp:
May 30, 2014, 8:03:16 PM (10 years ago)
Author:
annakrat
Comment:

week 2

Legend:

Unmodified
Added
Removed
Modified

  • GSoC/2014/ImplementationOf3DRasterFlowLine

    v2 v3  
    6767
    6868I will probably keep the code in sandbox and once it is usable, I will put it in addons
     69
     70'''Week 2'''
     71
     72I started to implement the prototype in Python:
     73 * trilinear interpolation implemented and tested against SciPy method, unit test is included in the code
     74 * Runge-Kutta integration of 4th order implemented too, but that has to be tested more
     75 * I started `r3.flow` itself but the prototype is not ready yet.
     76The code is in [source:sandbox/annakrat/r3.flow sandbox].
     77Things to consider:
     78 * Soeren suggests using interpolation approach by [http://www.isws.illinois.edu/pubdoc/b/iswsb-65.pdf Prickett, 1981] instead of trilinear interpolation from voxel centers.
     79 * r3.flow should be able to compute the gradient so that the input can be also just one 3D raster map instead of 3 maps representing the vector field. For now I could probably use !NumPy's method `gradient` based on central differences. When porting to C, there is [http://grass.osgeo.org/programming7/gpdelib.html GPDE] library in GRASS.