Re: [GRASS-user] How to calculates linear regression from three or more raster maps?

2013-02-22 Thread Glynn Clements
=?GBK?B?zfW1wrvU?= wrote: How to calculates linear regression from three or more raster maps £¨X1 = a + bX2 + cX3 +dX4 + ...)? r.series with a method of slope, offset and/or detcoeff. -- Glynn Clements gl...@gclements.plus.com ___ grass-user

Re: [GRASS-user] How to calculates linear regression from three or more raster maps?

2013-02-22 Thread 王德辉
I want to caculate linear regression function of X1 raster map to X2,X3,X4. r.series can not do it. Thanks! dehui 在2013-02-22 16:02:00,王德辉deh...@gig.ac.cn写道: =?GBK?B?zfW1wrvU?= wrote: How to calculates linear regression from three or more raster maps £¨X1 = a + bX2 + cX3 +dX4 + ...)?

Re: [GRASS-user] How to calculates linear regression from three or more raster maps?

2013-02-22 Thread Markus Metz
On Fri, Feb 22, 2013 at 9:44 AM, 王德辉 deh...@gig.ac.cn wrote: I want to caculate linear regression function of X1 raster map to X2,X3,X4. The addon r.regression.multi (GRASS 7 only) does that [0]. Markus M [0] http://grasswiki.osgeo.org/wiki/AddOns/GRASS7/raster#r.regression.multi r.series

Re: [GRASS-user] Problem in query vector maps

2013-02-22 Thread Vaclav Petras
Hi, if you are able to test it with the 6.4.3svn, release candidate or the version 7, please do. There were some changes (some only in 7) which may affect your problem. Vaclav On 20 February 2013 14:57, Marcello Benigno benigno.marce...@gmail.com wrote: Just one more information, I made some

[GRASS-user] r.stream.order and accum map

2013-02-22 Thread Micha Silver
I'm trying to use r.stream.order to get a vector with the full set of stream orders and topological attributes, but it errors with the need for a flow accumulation map as follows: r.stream.order --o -a streams=str dirs=fdir elev=dtm_full vect=stream_order

[GRASS-user] Quick note on EVI2 vs NDVI

2013-02-22 Thread Nikos Alexandris
Hi list! Just wanted to drop a note on EVI2. It looks and performs better than NDVI (tested with large Landsat5TM-based mosaics over Greece). It sure helps visual discrimination of dense vegetation from the rest. Will post screenshots at some point laaater (in the Wiki of course) -- using

Re: [GRASS-user] How to calculates linear regression from three or more raster maps?

2013-02-22 Thread Moritz Lennert
On 22/02/13 09:02, Glynn Clements wrote: =?GBK?B?zfW1wrvU?= wrote: How to calculates linear regression from three or more raster maps £¨X1 = a + bX2 + cX3 +dX4 + ...)? r.series with a method of slope, offset and/or detcoeff. Does r.series allow multiple regression ? How would you get

Re: [GRASS-user] How to calculates linear regression from three or more raster maps?

2013-02-22 Thread Moritz Lennert
On 22/02/13 13:32, Moritz Lennert wrote: On 22/02/13 09:02, Glynn Clements wrote: =?GBK?B?zfW1wrvU?= wrote: How to calculates linear regression from three or more raster maps £¨X1 = a + bX2 + cX3 +dX4 + ...)? r.series with a method of slope, offset and/or detcoeff. Does r.series allow

Re: [GRASS-user] Producing map layouts and map composition

2013-02-22 Thread Moritz Lennert
On 21/02/13 18:42, katrin eggert wrote: Dear Nikos, thank you for your reply. Just two questions: 1- Which one is better to produce composites? I mean, more easier and stable? As Nicos said, you have to try out for yourself. Most people will probably find QGIS easier in terms of interactive

[GRASS-user] Run a GRASS module from a custom module (C language)

2013-02-22 Thread José María Michia Roberts
I want to write a module to automatically do a large number of import/process/export tasks , but I cannot figure out how to run v.in.ogr from C program (I do not understand how to process and export, but maybe I can understand it after getting import). I've this functionality working on a Python

Re: [GRASS-user] Run a GRASS module from a custom module (C language)

2013-02-22 Thread Rainer M Krug
-BEGIN PGP SIGNED MESSAGE- Hash: SHA1 On 22/02/13 14:40, Jos← Mar■a Michia Roberts wrote: I want to write a module to automatically do a large number of import/process/export tasks , but I cannot figure out how to run v.in.ogr from C program (I do not understand how to process and

[GRASS-user] Bar in ps.map

2013-02-22 Thread Martin Album Ytre-Eide
Hello all. is it possible to bars in ps.map - like you would do in v.thematic.chart ? Martin ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user

[GRASS-user] Reproject Reclassed Raster Problems

2013-02-22 Thread Stephen Sefick
I am trying to re-project raster maps that have been re-classed. They are not showing up in the new mapset. What should I do? I can provide any information that is necessary to help solve this problem. many thanks, -- Stephen Sefick ** Auburn

Re: [GRASS-user] r.stream.order and accum map

2013-02-22 Thread jentjr
From a comment on line 74 of main.c it looks like the accumulation map is required when calculating the Strahler stream order. in_acc_opt-key = accum; /* required if strahler stream order is calculated for existing stream network */ in_acc_opt-type = TYPE_STRING; in_acc_opt-required = NO;

Re: [GRASS-user] r.stream.order and accum map

2013-02-22 Thread Micha Silver
Thanks, On 02/22/2013 05:17 PM, jentjr wrote: From a comment on line 74 of main.c it looks like the accumulation map is required when calculating the Strahler stream order. Actually you *can* create a raster stream order map based on

[GRASS-user] how to sample a series at one location?

2013-02-22 Thread Michael Barton
Is there tool somewhere, including in the new temporal GIS modules, to sample the value of a raster series at one cell location? I'd like to get text something like this for a cell that I specify with xy coordinates or a cat: map1, value 1 map2, value 2 map3, value 3 map4, value 4 map5, value 5

Re: [GRASS-user] [GRASS-dev] how to sample a series at one location?

2013-02-22 Thread Markus Neteler
On Fri, Feb 22, 2013 at 6:41 PM, Michael Barton michael.bar...@asu.edu wrote: Is there tool somewhere, including in the new temporal GIS modules, to sample the value of a raster series at one cell location? I'd like to get text something like this for a cell that I specify with xy coordinates

[GRASS-user] reclassed raster no longer showing up

2013-02-22 Thread Stephen Sefick
I have reclassed a raster with the rules 1=115 2=115 and it no longer shows up... qgis will display the raster, but not in grass. Is there a color table that I should use? Or some other problem? -- Stephen Sefick ** Auburn University

Re: [GRASS-user] [GRASS-dev] how to sample a series at one location?

2013-02-22 Thread Michael Barton
But I want to do it with a time series of hundreds or thousands of maps. Michael __ C. Michael Barton Director, Center for Social Dynamics Complexity Professor of Anthropology, School of Human Evolution Social Change Arizona State University Tempe, AZ 85287-2402

[GRASS-user] r.report

2013-02-22 Thread ALT SHN
Hi users, Can someone help with this doubt of mine? When I run r.report on a raster, it says that 48% of my cells have no no data - does this refers to NULL cells? (if so, than it's ok, otherwise I must recheck my data). Thanks in advance --

Re: [GRASS-user] [GRASS-dev] how to sample a series at one location?

2013-02-22 Thread Michael Barton
Thanks Doug, I knew I could do it with a script and v.what.rast. What I was hoping was that there is a shortcut already usable in GRASS modules. Looks like the new temporal GIS tools may be able to do it. Michael __ C. Michael Barton Director, Center for Social

Re: [GRASS-user] [GRASS-dev] how to sample a series at one location?

2013-02-22 Thread Markus Neteler
On Fri, Feb 22, 2013 at 7:23 PM, Michael Barton michael.bar...@asu.edu wrote: But I want to do it with a time series of hundreds or thousands of maps. Sure, I feed the output of g.mlist to r.what (using for example the 1440 MODIS LST maps of a year as input). Markus

Re: [GRASS-user] r.report

2013-02-22 Thread Nikos Alexandris
On Friday 22 of February 2013 19:20:56 ALT SHN wrote: Hi users, Hi! Can someone help with this doubt of mine? When I run r.report on a raster, it says that 48% of my cells have no no data - does this refers to NULL cells? (if so, than it's ok, otherwise I must recheck my data). Yes, it

Re: [GRASS-user] [GRASS-dev] how to sample a series at one location?

2013-02-22 Thread Michael Barton
Very cool. Your graph is what I'm trying to do also. What I'm doing is running a agent-based model (ABM) of small-holder agropastoral land-use in Java that is dynamically coupled to a landscape evolution model in Python and GRASS. One of the outputs is a map of net erosion/deposition in each

Re: [GRASS-user] [GRASS-dev] how to sample a series at one location?

2013-02-22 Thread Sören Gebbert
Hi, Am 22.02.2013 21:44 schrieb Newcomb, Doug doug_newc...@fws.gov: t.create-- t.register-- t.vect.observer.strds in GRASS7 ? Looks really nifty, could be useful with the Landsat Cube data sets http://landsat.usgs.gov/documents/Oct27_29_2009_huang_LST_boston.ppt Yes, thats the idea. I

Re: [GRASS-user] Bar in ps.map

2013-02-22 Thread Hamish
Martin Album Ytre-Eide wrote: is it possible to bars in ps.map - like you would do in v.thematic.chart ? can you post an image of what you mean? thanks, Hamish ___ grass-user mailing list grass-user@lists.osgeo.org

Re: [GRASS-user] reclassed raster no longer showing up

2013-02-22 Thread Nikos Alexandris
Stephen Sefick wrote: I have reclassed a raster with the rules 1=115 2=115 and it no longer shows up... qgis will display the raster, but not in grass. That is weird! Usually it is Quantum that doesn't read, for whatsoever might be the reason, some GRASS raster color tables. You

Re: [GRASS-user] Bar in ps.map

2013-02-22 Thread Vaclav Petras
Hi, d.vect.thematic has the parameter psmap but I've never tried. http://grass.osgeo.org/grass64/manuals/d.vect.thematic.html On 22 February 2013 22:19, Hamish hamis...@yahoo.com wrote: Martin Album Ytre-Eide wrote: is it possible to bars in ps.map - like you would do in v.thematic.chart ?

[GRASS-user] GRASS 70 digitizing question

2013-02-22 Thread E . J . Neafsey
Dear grass-user, I am running the latest svn of GRASS 70 downloaded last night (2/21/2013). When I am editing vectory lines using the digitize function in the wxgui, it is hard to see the vertices or points on lines. There used to be a small + sign that would appear over all vertices when a

Re: [GRASS-user] GRASS 70 digitizing question

2013-02-22 Thread Martin Landa
Hi, 2013/2/22 E.J. Neafsey enea...@lumc.edu: I am running the latest svn of GRASS 70 downloaded last night (2/21/2013). When I am editing vectory lines using the digitize function in the wxgui, it is hard to see the vertices or points on lines. vertices are disabled by default. Go to

Re: [GRASS-user] Interpolate a huge time-series with looping

2013-02-22 Thread Richard Chirgwin
Marcello - That's not so hard. In Grass-GIS 6.4 (the version I'm on ATM) the general form would be: for i in `echo select distinct id from rainfall_time_series | db.select -c`; do ## get year / month / day to use in output map name year=`echo select year from rainfall_time_series where

Re: [GRASS-user] How to calculates linear regression from three or more raster maps?

2013-02-22 Thread 王德辉
That's it.Thank you Markus! dehui 在2013-02-22 18:27:06,王德辉deh...@gig.ac.cn写道: On Fri, Feb 22, 2013 at 9:44 AM, 王德辉 deh...@gig.ac.cn wrote: I want to caculate linear regression function of X1 raster map to X2,X3,X4. The addon r.regression.multi (GRASS 7 only) does that [0]. Markus M

[GRASS-user] v.select query and GEOS error?

2013-02-22 Thread charlie
I attempted to create a v.select query using the “touch” option, which requires GEOS,  and kept   getting an error code stating that the only “supported operation is the overlap feature.”.   the package manager states that libgeos 3.3.3 is installed, so I’m perplexed. Is there anything im

Re: [GRASS-user] v.select query and GEOS error?

2013-02-22 Thread Hamish
charlie wrote: I attempted to create a v.select query using the “touch” option, which requires GEOS,  and kept  getting an error code stating that the only “supported operation is the overlap feature.”.   the package manager states that libgeos 3.3.3 is installed, so I’m perplexed. Is there