Hi Anna, i have added the function Rast3d_location2coord_double() to the raster3d library (r61282) that expects x, y and z as double values.
While doing this i discovered a really horrible bug in the coordinate transformation for rows. The reason for this bug were wrong coordinate tests in the raster3d library tests. Wrong indexing leaded me to wrong northern coordinate to row conversion. This bug appeared in case that coordinates are located exactly at the north/south border between two voxels. This is hopefully fixed now. It should be backported to grass7.0 as well, after heavy testing. Best regards Soeren 2014-07-19 19:04 GMT+02:00 Anna Petrášová <[email protected]>: > Hi, > > I have a question about Rast3d_location2coord [1]. It is casting the col, > row, depth values to integer like this: > > (int) col > > instead of using > > (int)floor(col) > > as a result, col = 0.1 => col = 0 and col = -0.1 => col = 0, which doesn't > seem correct to me. Any opinions? > > Thanks, Anna > > BTW, similar function Rast_northing_to_row from raster library returns > double. > > > [1] http://grass.osgeo.org/programming7/region_8c_source.html#l00291 > > _______________________________________________ > grass-dev mailing list > [email protected] > http://lists.osgeo.org/mailman/listinfo/grass-dev _______________________________________________ grass-dev mailing list [email protected] http://lists.osgeo.org/mailman/listinfo/grass-dev
