v.rast.stats: add where option

[veroandreo]

Would you agree to add where option in G7:v.rast.stats? I find it useful in some occasions where user is interested in querying only some features meeting a certain category. I attach the diff

[mmetz]

Replying to veroandreo:

Would you agree to add where option in G7:v.rast.stats? I find it useful in some occasions where user is interested in querying only some features meeting a certain category. I attach the diff

A where option would make sense. You need to add this where option to the call to G7:v.to.rast in G7:v.rast.stats.

[veroandreo]

Replying to mmetz:

Replying to veroandreo:

Would you agree to add where option in G7:v.rast.stats? I find it useful in some occasions where user is interested in querying only some features meeting a certain category. I attach the diff

A where option would make sense. You need to add this where option to the call to G7:v.to.rast in G7:v.rast.stats.

Right! See new patch attached. I tested with:

v.rast.stats map=urbanarea@PERMANENT where="NAME == 'Raleigh'" raster=aspect@PERMANENT column_prefix=aspect method=number

and got expected results and no errors :)

[veroandreo]

new patch adding where option in v.rast.stats

[veroandreo]

MarkusM, do you think I can commit it now?

[mmetz]

Replying to veroandreo:

MarkusM, do you think I can commit it now?

Have you tested your patch also if the new where option is not given? If this succeeds you can commit.

[veroandreo]

I tested... it does not work. It keeps extracting only where="NAME == 'Raleigh'" even if I remove the column and do not set any where... I added this:

if where == "" or where == " " or where == "\n":
		where = None

but it does not work... There seems to be something I need to clean/remove from previous run, but do not know where nor how. Would you help me, please?

[mmetz]

Replying to veroandreo:

I tested... it does not work. It keeps extracting only where="NAME == 'Raleigh'" even if I remove the column and do not set any where... I added this:

if where == "" or where == " " or where == "\n":
		where = None

but it does not work... There seems to be something I need to clean/remove from previous run, but do not know where nor how. Would you help me, please?

See attached patch.

Test commands for the NC sample dataset:

g.region rast=elev_state_500m
g.copy vect=railroads,myrailroads

v.rast.stats map=myrailroads raster=elev_state_500m column_prefix=ele method=average -d where="cat = 8802"
v.db.select map=myrailroads

v.db.update map=myrailroads col=ele_average val=NULL
v.rast.stats map=myrailroads raster=elev_state_500m column_prefix=ele method=average -d 
v.db.select map=myrailroads

v.db.update map=myrailroads col=ele_average val=NULL
v.rast.stats map=myrailroads raster=elev_state_500m column_prefix=ele method=average
v.db.select map=myrailroads

v.db.update map=myrailroads col=ele_average val=NULL
v.rast.stats map=myrailroads raster=elev_state_500m column_prefix=ele method=average -d where="cat = 8802"
v.db.select map=myrailroads

[mmetz]

patch to add where option o v.rast.stats, working with and without where option used

[veroandreo]

Thanks so much, MarkusM! Committed in r73666. Milestone was set to 7.8, but is it worth a backport to next 7.6 release?

[mmetz]

Replying to veroandreo:

Thanks so much, MarkusM! Committed in r73666. Milestone was set to 7.8, but is it worth a backport to next 7.6 release?

I think it's ok if you backport to 7.6

[veroandreo]

Milestone changed to 76 and change backported. Closing.

Thanks again, markusM :)