Dear Stefan, thank you very much for your feedback, the patch and suggestions.
I have integrated some of your suggestions in t.rast.what. There is now stdout and coordinates support in the latest trunk version (r64368). The column layout bug is now fixed and the documentation was updated according to the new features. Everything is now covered with tests. Hence, if the modules does not run as you expect, please make a test run (in a grass session make python test_what.py in the testsuite folder of t.rast.what) and send me the result of the test, so i can investigate the problem. I did not modified the options in t.rast.what as your patch suggested, since the options should run by default in your suggested configuration. The temporal database layout and its handling (distributed, mapset specific temporal databases) should be identical in the grass7.0 release version and grass7 trunk. Best regards Soeren 2015-01-30 8:49 GMT+01:00 Blumentrath, Stefan <[email protected]>: > Forgot to mention: > > t.rast.what does not seem to write to stdout (as the description suggests), > only to file. Output to stdout could be an improvement for people who want to > read the output directly into R for example (for that purpose the " > one_point_per_row_output"-layout is splendid). > > Carrying also the coordinates-option from r.what to t.rast.what would be nice > for cases where one only has a few coordinates to process and where copying > them from a text file or web and pasting them into the module is more > convenient than creating a vector map for a single use case in advance (and > probably deleting the vector map afterwards)... > > Finally, as written in the r.what ticket, being able to make use of the > site_name column also in "vector-map-mode" would be useful, because it is > more demanding (both for user and computer) to join additional information to > the output (or the other way around) using two columns in double precision > compared to e.g. one cat column(or whatever a user might like to have in > order to carachterize the sites), esp. because coordinates may be given with > different level of precision in r.what and v.to.db for example. But this > applies to r.what first of course... > An intermediate solution could be to pipe something like: > v.to.db map=INPUT_VECTOR_MAP option=coor qcolumn=cat separator=comma -p | awk > -v FS=',' -v OFS=',' '{print $2,$3,$1}' | tail -n +2 > to the coordinates option in r.what... That way the cat column would show up > in the r.what output... > > Thanks again for the excellent tool! > > Cheers, > Stefan > > > > -----Original Message----- > From: [email protected] > [mailto:[email protected]] On Behalf Of Blumentrath, Stefan > Sent: 29. januar 2015 07:30 > To: [email protected] > Subject: Re: [GRASS-dev] [GRASS GIS] #2564: r.what: remove 400 maps limit and > add some new output options > > Great news! Thanks. > I shall test immediately. > > -----Original Message----- > From: [email protected] > [mailto:[email protected]] On Behalf Of GRASS GIS > Sent: 29. januar 2015 00:16 > Subject: Re: [GRASS-dev] [GRASS GIS] #2564: r.what: remove 400 maps limit and > add some new output options > > #2564: r.what: remove 400 maps limit and add some new output options > -------------------------+---------------------------------------------- > -------------------------+---- > Reporter: sbl | Owner: grass-dev@… > Type: enhancement | Status: new > Priority: normal | Milestone: 7.1.0 > Component: Raster | Version: svn-trunk > Keywords: r.what | Platform: Unspecified > Cpu: Unspecified | > -------------------------+---------------------------------------------- > -------------------------+---- > > Comment(by huhabla): > > I have just committed a new temporal module called t.rast.what in trunk > (r64349 - r64351) that utilized r.what to sample space time raster datasets > using vector points. It provides three output layouts that transforms the > r.what output into row or column layouts. It can run several r.what > processes in parallel using a maximum of 400 raster map layer in a single > r.what process. > > Please have a look at it. > > Any feedback about its performance, usefulness and handling is highly > welcome. > > Improving r.what to allow more than 400 maps is simple, just edit the fixed > value in the source. But, editing this value may cause open file handler > limit problems with your OS. However, these limits can be adjusted at kernel > level configuration for Linux. > > Implementing a more intelligent file handling solution requires more effort. > > -- > Ticket URL: <http://trac.osgeo.org/grass/ticket/2564#comment:1> > GRASS GIS <http://grass.osgeo.org> > > _______________________________________________ > 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 > _______________________________________________ > 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
