Le 29 novembre 2016 22:33:39 GMT+01:00, Markus Metz <markus.metz.gisw...@gmail.com> a écrit : > For >reprojection, something like r.in.xyz could work:
I'm currently going through this process and it raises a series of issues: > >1. convert raster cells to vector points with raster cell value as >vector z >value. I use r.to.vect for this. The -b flag allows +/- fast conversion, but it takes into account all cells. It would be great to have a flag telling it to only convert non-null cells. Or at least a hint to set a mask, which seems to have the desired effect. >2. reproject these vector points to the target location I locally modified v.proj to add a 'Don't build topology' flag. Is there any indication against this ? >3. export the vector points with v.out.ascii format=point >4. import the exported points with r.in.xyz method=sum This works nicely. Thanks ! If I find the time I will wrap this into a r.proj.aggregate (?) script... > >The method options of r.in.xyz could be added to r.proj for spatial >aggregation, but that would mean that most of the code of r.in.xyz >would >need to be copied to r.proj, then adapted. That would be an interesting >enhancement, also considering that gdalwarp offers as resampling >methods >near, bilinear, cubic, cubicspline, lanczos, average, mode, max, min, >med, >Q1, Q3. But it doesn't offer sum.... Moritz
_______________________________________________ grass-dev mailing list grass-dev@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-dev