Re: [GRASS-user] r.stats.zonal output raster type

2021-01-24 Thread ming han
Hi everyone I just notice that r.statistics is not in grass add on, how to install it? 于2021年1月24日周日 下午2:51写道: > Try with `r.statistics`. Unfortunately, there are a few modules with > similar functionality, yet they work on different raster data types. > > > Nikos > --- > > > > On

Re: [GRASS-user] v.net.allpairs - reduce number of calculations

2021-01-24 Thread Markus Metz
What about v.net.path to get shortest paths from a start node to a number of end nodes? Markus M On Mon, Jan 11, 2021 at 2:32 PM Daniel McInerney < daniel.o.mciner...@gmail.com> wrote: > > Hi List, > > I've been using v.net.allpairs [1] for the calculation of shortest paths > between pairs of

Re: [GRASS-user] Error in compilation of GRASS GIS

2021-01-24 Thread Markus Metz
On Thu, Jan 21, 2021 at 5:22 AM Maris Nartiss wrote: > > Hello Ashutosh. > > Compilation and development issues are better to ask in dev mailing list. > > Also – do not send errors as images. I have no idea what the error > message you got was and I assume I am not the only one. MS Windows > does

Re: [GRASS-user] r.stats.zonal output raster type

2021-01-24 Thread ming han
Got it, many thanks. 于2021年1月24日周日 下午2:51写道: > Try with `r.statistics`. Unfortunately, there are a few modules with > similar functionality, yet they work on different raster data types. > > > Nikos > --- > > > > On 2021-01-24 21:25, ming han wrote: > > Hi everyone > > Is it possible to let

Re: [GRASS-user] r.stats.zonal output raster type

2021-01-24 Thread nik
Try with `r.statistics`. Unfortunately, there are a few modules with similar functionality, yet they work on different raster data types. Nikos --- On 2021-01-24 21:25, ming han wrote: Hi everyone Is it possible to let r.stats.zonal output raster to be CELL? Thanks Ming

[GRASS-user] r.stats.zonal output raster type

2021-01-24 Thread ming han
Hi everyone Is it possible to let r.stats.zonal output raster to be CELL? Thanks Ming ___ grass-user mailing list grass-user@lists.osgeo.org https://lists.osgeo.org/mailman/listinfo/grass-user

[GRASS-user] use grass within anaconda

2021-01-24 Thread ming han
Hi Everyone Is there a way to use grass within anaconda? Thanks Ming ___ grass-user mailing list grass-user@lists.osgeo.org https://lists.osgeo.org/mailman/listinfo/grass-user

Re: [GRASS-user] [GRASS-dev] compare a DCELL and FCELL question

2021-01-24 Thread Markus Metz
A note about the usage of GRASS_EPSILON = 1.0e-15: The range of single precision floating point is about 1.4012984643e−45 to 3.4028234664e38 That means for large numbers, using GRASS_EPSILON will also detect floating point precision limits, not only meaningful differences. For small numbers,

Re: [GRASS-user] [GRASS-dev] compare a DCELL and FCELL question

2021-01-24 Thread Markus Metz
On Sun, Jan 24, 2021 at 5:06 PM ming han wrote: > > Hi Everyone > >Many thanks for your help. Is if(fabs(map_A - map_B) <= 1.0e-15, ... ) approach may increase the runtime compare to '==' way? The formulas are very simple, I don't think that differences in runtime can be reliably measured.

Re: [GRASS-user] [GRASS-dev] compare a DCELL and FCELL question

2021-01-24 Thread ming han
Hi Everyone Many thanks for your help. Is if(fabs(map_A - map_B) <= 1.0e-15, ... ) approach may increase the runtime compare to '==' way? Thanks Ming Markus Metz 于2021年1月24日周日 上午10:57写道: > Trying to answer the original question: with a DCELL map "cat1_acc_riv" > and a FCELL map

Re: [GRASS-user] [GRASS-dev] compare a DCELL and FCELL question

2021-01-24 Thread Markus Metz
Trying to answer the original question: with a DCELL map "cat1_acc_riv" and a FCELL map "cat1_minacc", why is "float(cat1_acc_riv) == float(cat1_minacc)" not equal to "int(cat1_acc_riv) == int(cat1_minacc)" ? int truncates to integer while float converts to single precision floating point. E.g.

Re: [GRASS-user] compare a DCELL and FCELL question

2021-01-24 Thread ming han
Hi Markus and Micha I am just trying to find grids have the same values in these two rasters, I will try the threshold approach. Thanks Ming Markus Neteler 于2021年1月24日周日 上午6:58写道: > Hi Ming, > > On Sun, Jan 24, 2021 at 10:49 AM ming han wrote: > > > > Hi Micha > > > > Many thanks

Re: [GRASS-user] compare a DCELL and FCELL question

2021-01-24 Thread Markus Neteler
Hi Ming, On Sun, Jan 24, 2021 at 10:49 AM ming han wrote: > > Hi Micha > > Many thanks for your reply. > Here is the command I am using: > > if(float(cat1_acc_riv) == float(cat1_minacc), str_r, null()) > > The str_r is a CELL raster. the result is different when I change it

Re: [GRASS-user] compare a DCELL and FCELL question

2021-01-24 Thread Micha Silver
I can see that the maximum values are different. Have a look at the Range values. What are you trying to do? Maybe there's a simpler approach... On Sun, Jan 24, 2021 at 11:49 AM ming han wrote: > > Hi Micha > > Many thanks for your reply. > Here is the command I am using: > >

Re: [GRASS-user] compare a DCELL and FCELL question

2021-01-24 Thread ming han
Hi Micha Many thanks for your reply. Here is the command I am using: if(float(cat1_acc_riv) == float(cat1_minacc), str_r, null()) The str_r is a CELL raster. the result is different when I change it to: if(int(cat1_acc_riv) == int(cat1_minacc), str_r, null())

Re: [GRASS-user] compare a DCELL and FCELL question

2021-01-24 Thread Micha Silver
Is there some reason that you expect the rasters to be the same? Maybe begin by posting the outputs of `r.info` for both rasters, and what command you used to compare them? On Sun, Jan 24, 2021 at 6:13 AM ming han wrote: > > Hi Everyone > > I tried to compare if grids in two rasters are the