Re: [GRASS-user] i.atcorr results in all nulls for 3 out of 4 bands

2018-06-07 Thread hindo kush
Respected,

i have the same problem for making the same file for Landsat 8, as i have
the MTL text but how to make it as an 6s parameter.
please advice me too.

regards

On Wed, Jun 6, 2018 at 12:31 PM, Rob Deg  wrote:

> Dear Markus and Gwen(all?),
>
> Iam new to grass gis and want to run the atcorr over Planet scene.
> Can you please give me a hint how to enter those values into the parameter
> file (see below)?
>
> Lets say Iam going to run the atcorr over the blue band using the *.csv you
> posted. Here is a partial parameter file:
>
> 13  - geometrical conditions=Rapideye ETM+ (using
> rapideye parameter instead)
> 5 12 09.40 47.455 8.250   - month day hh.ddd longitude latitude ("hh.ddd"
> is
> in decimal hours GMT)
> 2- atmospheric mode=midlatitude summer
> 1- aerosols model=continental
> 15   - visibility [km] (aerosol model
> concentration)
> -0.600   - mean target elevation above sea level [km]
> (here 600m asl)
> -1000- sensor height (here, sensor on board a
> satellite)
> ??? - Sensor band
>
>
> Best,
> Rob
>
>
>
>
>
>
> --
> Sent from: http://osgeo-org.1560.x6.nabble.com/Grass-Users-f3884509.html
> ___
> grass-user mailing list
> grass-user@lists.osgeo.org
> https://lists.osgeo.org/mailman/listinfo/grass-user
___
grass-user mailing list
grass-user@lists.osgeo.org
https://lists.osgeo.org/mailman/listinfo/grass-user

Re: [GRASS-user] [GRASS-dev] GRASS and Anaconda

2018-06-07 Thread Luca Delucchi
Il gio 7 giu 2018, 20:58 Michael Barton  ha scritto:

> Along with helping me set up a new way to compile and package GRASS
> binaries for the Mac, Eric Hutton of the Community Surface Dynamics
> Modeling System (CSDMS – http//csdms.colorado.edu) has created a Conda
> recipe for building GRASS in the Anaconda environment. This is now deployed
> in the csdms stack at GitHub. Several versions of GRASS, for Linux and Mac
> are available. Details and instructions for building GRASS with Anaconda
> can be found at:
>
>
>
> https://github.com/csdms-stack/grass-recipe
>
>
>
> Many thanks to Eric and CSDMS
>
>
>
> Happy GRASSing!
>
> Michael
>
>
>
> __
>
> C. Michael Barton
>
> Director, Center for Social Dynamics & Complexity
>
> Professor of Anthropology, School of Human Evolution & Social Change
>
> Head, Graduate Faculty in Complex Adaptive Systems Science
>
> Arizona State University
>
> Tempe, AZ  85287-2402
>
> USA
>
>
>
> voice:480-965-6262 (SHESC), 480-965-8130/727-9746 (CSDC)
>
> fax:  480-965-7671(SHESC), 480-727-0709 (CSDC)
>
> www:  http://csdc.asu.edu, http://shesc.asu.edu
>
> http://www.public.asu.edu/~cmbarton
>
>
> ___
> grass-dev mailing list
> grass-...@lists.osgeo.org
> https://lists.osgeo.org/mailman/listinfo/grass-dev
___
grass-user mailing list
grass-user@lists.osgeo.org
https://lists.osgeo.org/mailman/listinfo/grass-user

Re: [GRASS-user] [EXTERNAL] [GRASS-dev] GRASS and Anaconda

2018-06-07 Thread Michael Barton
No. These are not available for Anaconda, at least for Mac. I tried. You can 
get lastools, though and use them through the command line, including to change 
a laszip to ASCII for upload to GRASS using r.in.xyz.

Michael

__
C. Michael Barton
Director, Center for Social Dynamics & Complexity
Professor of Anthropology, School of Human Evolution & Social Change
Head, Graduate Faculty in Complex Adaptive Systems Science
Arizona State University
Tempe, AZ  85287-2402
USA

voice:480-965-6262 (SHESC), 480-965-8130/727-9746 (CSDC)
fax:  480-965-7671(SHESC), 480-727-0709 (CSDC)
www:  http://csdc.asu.edu, http://shesc.asu.edu
http://www.public.asu.edu/~cmbarton


From: Doug Newcomb 
Date: Thursday, June 7, 2018 at 12:03 PM
To: Michael Barton 
Cc: GRASS developers list , GRASS users 
, Eric Hutton 
Subject: Re: [EXTERNAL] [GRASS-dev] GRASS and Anaconda

Just glancing over the recipe, it does not look like liblas and laszip are 
included for LiDAR analysis.

Doug

On Thu, Jun 7, 2018 at 2:39 PM, Michael Barton 
mailto:michael.bar...@asu.edu>> wrote:
Along with helping me set up a new way to compile and package GRASS binaries 
for the Mac, Eric Hutton of the Community Surface Dynamics Modeling System 
(CSDMS – 
http//csdms.colorado.edu)
 has created a Conda recipe for building GRASS in the Anaconda environment. 
This is now deployed in the csdms stack at GitHub. Several versions of GRASS, 
for Linux and Mac are available. Details and instructions for building GRASS 
with Anaconda can be found at:

https://github.com/csdms-stack/grass-recipe

Many thanks to Eric and CSDMS

Happy GRASSing!
Michael

__
C. Michael Barton
Director, Center for Social Dynamics & Complexity
Professor of Anthropology, School of Human Evolution & Social Change
Head, Graduate Faculty in Complex Adaptive Systems Science
Arizona State University
Tempe, AZ  85287-2402
USA

voice:480-965-6262 (SHESC), 480-965-8130/727-9746 (CSDC)
fax:  480-965-7671(SHESC), 480-727-0709 (CSDC)
www:  
http://csdc.asu.edu,
 
http://shesc.asu.edu

http://www.public.asu.edu/~cmbarton


___
grass-dev mailing list
grass-...@lists.osgeo.org
https://lists.osgeo.org/mailman/listinfo/grass-dev



--
Doug Newcomb
USFWS
551F Pylon Dr
Raleigh, NC
919-856-4520 ext. 14 doug_newc...@fws.gov
-

NOTE: This email correspondence and any attachments to and from this sender is 
subject to the Freedom of Information Act (FOIA) and may be disclosed to third 
parties.​
___
grass-user mailing list
grass-user@lists.osgeo.org
https://lists.osgeo.org/mailman/listinfo/grass-user

[GRASS-user] GRASS and Anaconda

2018-06-07 Thread Michael Barton
Along with helping me set up a new way to compile and package GRASS binaries 
for the Mac, Eric Hutton of the Community Surface Dynamics Modeling System 
(CSDMS – http//csdms.colorado.edu) has created a Conda recipe for building 
GRASS in the Anaconda environment. This is now deployed in the csdms stack at 
GitHub. Several versions of GRASS, for Linux and Mac are available. Details and 
instructions for building GRASS with Anaconda can be found at:

https://github.com/csdms-stack/grass-recipe

Many thanks to Eric and CSDMS

Happy GRASSing!
Michael

__
C. Michael Barton
Director, Center for Social Dynamics & Complexity
Professor of Anthropology, School of Human Evolution & Social Change
Head, Graduate Faculty in Complex Adaptive Systems Science
Arizona State University
Tempe, AZ  85287-2402
USA

voice:480-965-6262 (SHESC), 480-965-8130/727-9746 (CSDC)
fax:  480-965-7671(SHESC), 480-727-0709 (CSDC)
www:  http://csdc.asu.edu, http://shesc.asu.edu
http://www.public.asu.edu/~cmbarton

___
grass-user mailing list
grass-user@lists.osgeo.org
https://lists.osgeo.org/mailman/listinfo/grass-user

Re: [GRASS-user] r.series weights to zero

2018-06-07 Thread Markus Metz
On Tue, Jun 5, 2018 at 2:51 PM, Frank David  wrote:
>
> Hi all Grass users,
>
> It's seems that r.series weights does not accept zero as weight... Is it
possible ?

Yes, weights must be positive.

> and so, why ?

A zero weight means that all values in this map are ignored. Zero weights
can also cause the sum of weights to be zero, leading to division by zero
and finally to +inf / -inf as a result.

Markus M
>
> I've got the following error message when running r.series with one
weight set to 0
>
> "/usr/local/grass-7.4.0svn/etc/python/grass/script/core.py",
> line 332, in handle_errors
> returncode=returncode)
> grass.exceptions.CalledModuleError: Module run None
>
> Thank you for your comment.
>
> Frank
>
>
> ___
> grass-user mailing list
> grass-user@lists.osgeo.org
> https://lists.osgeo.org/mailman/listinfo/grass-user
___
grass-user mailing list
grass-user@lists.osgeo.org
https://lists.osgeo.org/mailman/listinfo/grass-user

Re: [GRASS-user] i.gensig - multiple signatures into a single file for i.maxlik

2018-06-07 Thread Moritz Lennert


Am 7. Juni 2018 16:53:51 MESZ schrieb Paul Shapley :
>Hi Users,
>
>
>What is the correct format for placing multiple signatures into a
>single ascii file.
>
>Is a '#' required between each signature along with a blank line?
>
>or can i start the new signature with line 2 e.g. #grass
>
>and is there a '#' after the last line to end the file.
>

I'm not in front of a computer right now and cannot check the code, so I don't 
have an answer, but why are you trying to craft this file by hand ? i.gensig is 
the tool dedicated to cresting these files for you.

Moritz


>
>1 #
>2 #water
>3 4186
>4 67.9508 48.7346 37.8915 15.3129 13.8473 12.0855
>5 1.74334
>6 0.439504 2.07267
>7 0.662523 1.63501 4.21189
>8 0.530339 2.40757 5.52857 22.433
>9 0.561184 2.30762 5.18846 20.5364 20.4926
>10 0.393218 1.2184 2.63628 9.61528 9.36025 5.8531411
>
>11 #grass
>
>12 2132
>
>13 etc
>
>*OR*
>
>1 #
>2 #water
>3 4186
>4 67.9508 48.7346 37.8915 15.3129 13.8473 12.0855
>5 1.74334
>6 0.439504 2.07267
>7 0.662523 1.63501 4.21189
>8 0.530339 2.40757 5.52857 22.433
>9 0.561184 2.30762 5.18846 20.5364 20.4926
>10 0.393218 1.2184 2.63628 9.61528 9.36025 5.8531411
>
>11 #
>
>12 #grass
>
>13 etc
>
>Many thanks,
___
grass-user mailing list
grass-user@lists.osgeo.org
https://lists.osgeo.org/mailman/listinfo/grass-user

[GRASS-user] i.gensig - multiple signatures into a single file for i.maxlik

2018-06-07 Thread Paul Shapley
Hi Users,


What is the correct format for placing multiple signatures into a
single ascii file.

Is a '#' required between each signature along with a blank line?

or can i start the new signature with line 2 e.g. #grass

and is there a '#' after the last line to end the file.


1 #
2 #water
3 4186
4 67.9508 48.7346 37.8915 15.3129 13.8473 12.0855
5 1.74334
6 0.439504 2.07267
7 0.662523 1.63501 4.21189
8 0.530339 2.40757 5.52857 22.433
9 0.561184 2.30762 5.18846 20.5364 20.4926
10 0.393218 1.2184 2.63628 9.61528 9.36025 5.8531411

11 #grass

12 2132

13 etc

*OR*

1 #
2 #water
3 4186
4 67.9508 48.7346 37.8915 15.3129 13.8473 12.0855
5 1.74334
6 0.439504 2.07267
7 0.662523 1.63501 4.21189
8 0.530339 2.40757 5.52857 22.433
9 0.561184 2.30762 5.18846 20.5364 20.4926
10 0.393218 1.2184 2.63628 9.61528 9.36025 5.8531411

11 #

12 #grass

13 etc

Many thanks,

-- 
*Paul J. Shapley *MSc CGeog (GIS) FRGS
___
grass-user mailing list
grass-user@lists.osgeo.org
https://lists.osgeo.org/mailman/listinfo/grass-user

Re: [GRASS-user] v.class.mlR Error

2018-06-07 Thread Jamille Haarloo
Hello Moritz,

No worries. Thankful these modules are made available for newbies in RS
like me and also happy these interactions are possible for learning.
Hope to get back soon after some adjustments.

Best,
Jamille

On Thu, Jun 7, 2018 at 10:44 AM, Moritz Lennert <
mlenn...@club.worldonline.be> wrote:

> Thanks
>
> On 07/06/18 15:17, Jamille Haarloo wrote:
>
>> The first 20+ lines of Stats_Training_Ben_test:
>>
>> cat,area,perimeter,compact_circle,compact_square,fd,WV_Benat
>> imofo_1_min,WV_Benatimofo_1_max,WV_Benatimofo_1_range,WV_B
>> enatimofo_1_mean,WV_Benatimofo_1_stddev,WV_Benatimofo_1_
>> variance,WV_Benatimofo_1_coeff_var,WV_Benatimofo_1_sum,
>> WV_Benatimofo_1_first_quart,WV_Benatimofo_1_median,WV_Bena
>> timofo_1_third_quart,WV_Benatimofo_2_min,WV_Benatimofo_2_
>> max,WV_Benatimofo_2_range,WV_Benatimofo_2_mean,WV_Benatimofo
>> _2_stddev,WV_Benatimofo_2_variance,WV_Benatimofo_2_
>> coeff_var,WV_Benatimofo_2_sum,WV_Benatimofo_2_first_quart,
>> WV_Benatimofo_2_median,WV_Benatimofo_2_third_quart,WV_Benati
>> mofo_3_min,WV_Benatimofo_3_max,WV_Benatimofo_3_range,WV_B
>> enatimofo_3_mean,WV_Benatimofo_3_stddev,WV_Benatimofo_3_
>> variance,WV_Benatimofo_3_coeff_var,WV_Benatimofo_3_sum,
>> WV_Benatimofo_3_first_quart,WV_Benatimofo_3_median,WV_Bena
>> timofo_3_third_quart,WV_Benatimofo_4_min,WV_Benatimofo_4_
>> max,WV_Benatimofo_4_range,WV_Benatimofo_4_mean,WV_Benatimofo
>> _4_stddev,WV_Benatimofo_4_variance,WV_Benatimofo_4_
>> coeff_var,WV_Benatimofo_4_sum,WV_Benatimofo_4_first_quart,
>> WV_Benatimofo_4_median,WV_Benatimofo_4_third_quart
>> 1144,3832.00,1256.00,5.723635,0.197144,1.729624,13,7
>> 6,63,46.4097077244259,9.98454911351384,99.69122100017,21.
>> 5139237092391,177842,40,47,53,40,138,98,90.2687891440501,15.
>> 2500825418009,232.565017531741,16.8940812061464,345910,81,
>> 92,100,15,61,46,40.8582985386221,7.82663897784868,61.2562776
>> 895802,19.1555675536767,156569,36,42,47,28,124,96,68.
>> 4253653444676,13.5774536655369,184.347248039801,19.
>> 8427200164517,262206,59,68,77
>> 1145,12092.00,2282.00,5.854120,0.192750,1.645226,13,
>> 94,81,51.386288455177,10.5294376761475,110.869057775874,20.
>> 4907534532914,621363,45,52,59,21,220,199,114.230731061859,
>> 23.3590328249442,545.644414516822,20.4489917973953,1381278,
>> 101,114,128,7,76,69,46.4219318557724,8.42747122371732
>> ,71.0222712265835,18.1540726264915,561334,42,48,52,17,198,
>> 181,97.2732385047966,22.492313569247,505.904169697333,
>> 23.1228176577445,1176228,84,97,110
>>
>> [...]
>
> -
>> All the lines of the output of v.db.select Training_Ben2@haarlooj_Ben_Tes
>> t:
>>
>> cat|id|Type|code
>> 1|4|B29|18
>> 2|5|B31|19
>> 3|3|B28|17
>>
>
>
> Again a lack of clear documentation on my side: both the training and the
> segment info should contain the same attributes, with only additional one
> column ('code' in your case) present in the training data.
>
> It should be possible to do this differently, i.e. provide the module with
> the features of all segments, and only the id/cat of each training segment
> with the relevant class and have the module merge the two, but this is not
> implemented, yet.
>
> I also just notice that you have the word 'Training' in both names.
>
> The segment_file/segment_map contains the info (cat + all feature
> variables) of all segments you wish to classify, either in the form of a
> csv file or in the form of a vector map with the info in the attribute
> table.
>
> The training_file/training_map contains the info (cat + all feature
> variables + class) of the training data. Often this is an extract of the
> former, but not necessarily.
>
> All columns in the training file have to be present in the segment file,
> except for the class column (your 'code').
>
> Sorry for the lack of docs. This module has mostly been used internally
> here and so we are not always aware of the unclear and missing parts.
> Having your feedback has been very useful !
>
> Moritz
>
>
___
grass-user mailing list
grass-user@lists.osgeo.org
https://lists.osgeo.org/mailman/listinfo/grass-user

Re: [GRASS-user] v.class.mlR Error

2018-06-07 Thread Moritz Lennert

Thanks

On 07/06/18 15:17, Jamille Haarloo wrote:

The first 20+ lines of Stats_Training_Ben_test:

cat,area,perimeter,compact_circle,compact_square,fd,WV_Benatimofo_1_min,WV_Benatimofo_1_max,WV_Benatimofo_1_range,WV_Benatimofo_1_mean,WV_Benatimofo_1_stddev,WV_Benatimofo_1_variance,WV_Benatimofo_1_coeff_var,WV_Benatimofo_1_sum,WV_Benatimofo_1_first_quart,WV_Benatimofo_1_median,WV_Benatimofo_1_third_quart,WV_Benatimofo_2_min,WV_Benatimofo_2_max,WV_Benatimofo_2_range,WV_Benatimofo_2_mean,WV_Benatimofo_2_stddev,WV_Benatimofo_2_variance,WV_Benatimofo_2_coeff_var,WV_Benatimofo_2_sum,WV_Benatimofo_2_first_quart,WV_Benatimofo_2_median,WV_Benatimofo_2_third_quart,WV_Benatimofo_3_min,WV_Benatimofo_3_max,WV_Benatimofo_3_range,WV_Benatimofo_3_mean,WV_Benatimofo_3_stddev,WV_Benatimofo_3_variance,WV_Benatimofo_3_coeff_var,WV_Benatimofo_3_sum,WV_Benatimofo_3_first_quart,WV_Benatimofo_3_median,WV_Benatimofo_3_third_quart,WV_Benatimofo_4_min,WV_Benatimofo_4_max,WV_Benatimofo_4_range,WV_Benatimofo_4_mean,WV_Benatimofo_4_stddev,WV_Benatimofo_4_variance,WV_Benatimofo_4_coeff_var,WV_Benatimofo_4_sum,WV_Benatimofo_4_first_quart,WV_Benatimofo_4_median,WV_Benatimofo_4_third_quart
1144,3832.00,1256.00,5.723635,0.197144,1.729624,13,76,63,46.4097077244259,9.98454911351384,99.69122100017,21.5139237092391,177842,40,47,53,40,138,98,90.2687891440501,15.2500825418009,232.565017531741,16.8940812061464,345910,81,92,100,15,61,46,40.8582985386221,7.82663897784868,61.2562776895802,19.1555675536767,156569,36,42,47,28,124,96,68.4253653444676,13.5774536655369,184.347248039801,19.8427200164517,262206,59,68,77
1145,12092.00,2282.00,5.854120,0.192750,1.645226,13,94,81,51.386288455177,10.5294376761475,110.869057775874,20.4907534532914,621363,45,52,59,21,220,199,114.230731061859,23.3590328249442,545.644414516822,20.4489917973953,1381278,101,114,128,7,76,69,46.4219318557724,8.42747122371732,71.0222712265835,18.1540726264915,561334,42,48,52,17,198,181,97.2732385047966,22.492313569247,505.904169697333,23.1228176577445,1176228,84,97,110


[...]


-
All the lines of the output of v.db.select Training_Ben2@haarlooj_Ben_Test:

cat|id|Type|code
1|4|B29|18
2|5|B31|19
3|3|B28|17



Again a lack of clear documentation on my side: both the training and 
the segment info should contain the same attributes, with only 
additional one column ('code' in your case) present in the training data.


It should be possible to do this differently, i.e. provide the module 
with the features of all segments, and only the id/cat of each training 
segment with the relevant class and have the module merge the two, but 
this is not implemented, yet.


I also just notice that you have the word 'Training' in both names.

The segment_file/segment_map contains the info (cat + all feature 
variables) of all segments you wish to classify, either in the form of a 
csv file or in the form of a vector map with the info in the attribute 
table.


The training_file/training_map contains the info (cat + all feature 
variables + class) of the training data. Often this is an extract of the 
former, but not necessarily.


All columns in the training file have to be present in the segment file, 
except for the class column (your 'code').


Sorry for the lack of docs. This module has mostly been used internally 
here and so we are not always aware of the unclear and missing parts. 
Having your feedback has been very useful !


Moritz

___
grass-user mailing list
grass-user@lists.osgeo.org
https://lists.osgeo.org/mailman/listinfo/grass-user

Re: [GRASS-user] v.class.mlR Error

2018-06-07 Thread Jamille Haarloo
The first 20+ lines of Stats_Training_Ben_test:

cat,area,perimeter,compact_circle,compact_square,fd,WV_Benatimofo_1_min,WV_Benatimofo_1_max,WV_Benatimofo_1_range,WV_Benatimofo_1_mean,WV_Benatimofo_1_stddev,WV_Benatimofo_1_variance,WV_Benatimofo_1_coeff_var,WV_Benatimofo_1_sum,WV_Benatimofo_1_first_quart,WV_Benatimofo_1_median,WV_Benatimofo_1_third_quart,WV_Benatimofo_2_min,WV_Benatimofo_2_max,WV_Benatimofo_2_range,WV_Benatimofo_2_mean,WV_Benatimofo_2_stddev,WV_Benatimofo_2_variance,WV_Benatimofo_2_coeff_var,WV_Benatimofo_2_sum,WV_Benatimofo_2_first_quart,WV_Benatimofo_2_median,WV_Benatimofo_2_third_quart,WV_Benatimofo_3_min,WV_Benatimofo_3_max,WV_Benatimofo_3_range,WV_Benatimofo_3_mean,WV_Benatimofo_3_stddev,WV_Benatimofo_3_variance,WV_Benatimofo_3_coeff_var,WV_Benatimofo_3_sum,WV_Benatimofo_3_first_quart,WV_Benatimofo_3_median,WV_Benatimofo_3_third_quart,WV_Benatimofo_4_min,WV_Benatimofo_4_max,WV_Benatimofo_4_range,WV_Benatimofo_4_mean,WV_Benatimofo_4_stddev,WV_Benatimofo_4_variance,WV_Benatimofo_4_coeff_var,WV_Benatimofo_4_sum,WV_Benatimofo_4_first_quart,WV_Benatimofo_4_median,WV_Benatimofo_4_third_quart
1144,3832.00,1256.00,5.723635,0.197144,1.729624,13,76,63,46.4097077244259,9.98454911351384,99.69122100017,21.5139237092391,177842,40,47,53,40,138,98,90.2687891440501,15.2500825418009,232.565017531741,16.8940812061464,345910,81,92,100,15,61,46,40.8582985386221,7.82663897784868,61.2562776895802,19.1555675536767,156569,36,42,47,28,124,96,68.4253653444676,13.5774536655369,184.347248039801,19.8427200164517,262206,59,68,77
1145,12092.00,2282.00,5.854120,0.192750,1.645226,13,94,81,51.386288455177,10.5294376761475,110.869057775874,20.4907534532914,621363,45,52,59,21,220,199,114.230731061859,23.3590328249442,545.644414516822,20.4489917973953,1381278,101,114,128,7,76,69,46.4219318557724,8.42747122371732,71.0222712265835,18.1540726264915,561334,42,48,52,17,198,181,97.2732385047966,22.492313569247,505.904169697333,23.1228176577445,1176228,84,97,110
344,56623.00,5622.00,6.664836,0.169303,1.577907,7,117,110,48.9859950903343,12.2148119091359,149.201629975569,24.9353144436706,2773734,41,49,57,13,195,182,98.0623951397842,23.0181102420013,529.833399112924,23.4729227337246,5552587,83,98,113,5,102,97,42.592692015612,9.50012967113416,90.2524637683636,22.3046002061855,2411726,37,43,49,10,199,189,86.872489977571,23.6622500676856,559.902078265688,27.2379093471303,4918981,71,86,102
345,3711.00,932.00,4.315842,0.261451,1.663776,9,88,79,41.4944758825115,16.4423790499195,270.351828821232,39.6254650775079,153986,28,40,55,15,149,134,69.1158717326866,23.6500249395838,559.323679642938,34.2179362665828,256489,52,70,85,6,78,72,36.4807329560765,14.859354754329,220.800423714999,40.7320619687655,135380,24,36,49,7,110,103,49.5238480194018,19.2285018694077,369.735284141817,38.8267524403085,183783,34,49,63
346,3537.00,822.00,3.898965,0.289405,1.642813,23,122,99,68.4042974271982,11.6900148615688,136.6564474637,17.0895912994507,241946,62,69,76,83,244,161,149.022052586938,23.6287627271847,558.318428017593,15.8558832850594,527091,131,148,166,27,97,70,69.2623692394685,7.79821357477528,60.8121349578094,11.2589471893658,244981,65,70,74,49,255,206,126.945716709075,32.1492735900629,1033.57579236872,25.3252133459061,449007,102,124,148
347,62555.00,7296.00,8.229034,0.137122,1.610873,7,194,187,58.1883142834306,12.9492036231701,167.681874474322,22.2539590339318,3639970,50,59,67,13,254,241,123.435153065303,27.877934167845,777.1792134667,22.5850849418045,7721486,107,125,141,4,180,176,51.6065862041404,11.6910103595584,136.679723227302,22.6541052595734,3228250,45,52,59,6,255,249,107.161601790424,34.1020311883702,1162.94853117258,31.8229950081032,6703494,85,107,129
340,3211.00,1110.00,5.525835,0.204201,1.736889,5,92,87,39.2130177514793,12.7088506576442,161.514885038302,32.4097745758543,125913,30,40,49,19,183,164,81.8589224540642,24.2890122701093,589.956117057518,29.6717957455881,262849,65,83,100,5,68,63,33.9841170974774,10.7587310002139,115.750292734965,31.65811537594,109123,26,35,42,17,189,172,74.7689193397695,22.2286365349259,494.11228220184,29.7297817478318,240083,61,75,88
341,58033.00,7022.00,8.222775,0.137226,1.614913,7,109,102,43.4994916685334,11.4706463726522,131.575728206438,26.369610155582,2524406,36,44,52,17,188,171,89.8843416676718,22.9160732139348,525.146411546421,25.4950670926223,5216258,74,90,105,6,87,81,39.6913652576982,9.36784216838539,87.7564668917794,23.6017131372635,2303409,34,41,46,13,199,186,67.3895025244258,19.4984400094688,380.189162802852,28.933942645446,3910815,54,66,79

Re: [GRASS-user] v.class.mlR Error

2018-06-07 Thread Moritz Lennert

On 07/06/18 13:45, Jamille Haarloo wrote:
I see. That was one of the big problems then. I knew it used R-resources 
but wasn't sure the software itself was required. Must have overlooked 
it in the documentation and other user- topics as I was reading about 
the module.


You are right that it is never explicitely stated that it is the user's 
responsibility to install R. I just committed a patch to the man page 
which explains this in the Dependencies section.




Now I have R-3.5, but still get an error:

v.class.mlR --overwrite training_map=Training_Ben2@haarlooj_Ben_Test 
segments_file=C:\Users\haarlooj\Documents\CELOS\Stats_Training_Ben_test 
raster_segments_map=best_myregion1_at_haarlooj_Ben_Test_rank1@haarlooj_Ben_Test 
train_class_column=code output_class_column=vote output_prob_column=prob 
classifiers=svmRadial,rf,C5.0 folds=5 partitions=10 tunelength=10 
weighting_modes=smv,qbwwv weighting_metric=accuracy 
classification_results=C:\Users\haarlooj\Documents\CELOS\v.class.mIRR_optional_output\Ben_test_Classifier-results 
accuracy_file=C:\Users\haarlooj\Documents\CELOS\v.class.mIRR_optional_output\Ben_test_Classifier-accuracy 
model_details=C:\Users\haarlooj\Documents\CELOS\v.class.mIRR_optional_output\Ben_test_Classifier-module-runs 
bw_plot_file=C:\Users\haarlooj\Documents\CELOS\v.class.mIRR_optional_output\Ben_test_Classifier-performance 
r_script_file=C:\Users\haarlooj\Documents\CELOS\v.class.mIRR_optional_output\Ben_test_R_script 
processes=3


[...]


Error in eval(predvars, data, env) : object 'id' not found


This seems to be the relevant error, but I'm not sure I understand where 
it comes from.


Could you provide us with the first 20 lines of

C:\Users\haarlooj\Documents\CELOS\Stats_Training_Ben_test

and the first 20 lines of the output of

v.db.select Training_Ben2@haarlooj_Ben_Test

?


Calls: data.frame ... predict.train -> model.frame -> 
model.frame.default -> eval -> eval

Execution halted
ERROR: There was an error in the execution of the R script.
Please check the R output.
Error in atexit._run_exitfuncs:
Traceback (most recent call last):
   File "C:\Program Files\GRASS GIS
7.4.0\Python27\lib\atexit.py", line 24, in _run_exitfuncs
     func(*targs, **kargs)
   File "C:\Users\haarlooj\AppData\Roaming\GRASS7\addons/scri
pts/v.class.mlR.py ", line 236, in cleanup
     gscript.try_remove(model_output_csv)
NameError: global name 'model_output_csv' is not defined
Error in sys.exitfunc:
Traceback (most recent call last):
   File "C:\Program Files\GRASS GIS
7.4.0\Python27\lib\atexit.py", line 24, in _run_exitfuncs
     func(*targs, **kargs)
   File "C:\Users\haarlooj\AppData\Roaming\GRASS7\addons/scri
pts/v.class.mlR.py ", line 236, in cleanup
     gscript.try_remove(model_output_csv)
NameError: global name 'model_output_csv' is not defined



This last error should not be happening anymore since yesterday. Could 
you reinstall using g.extension ?


Moritz
___
grass-user mailing list
grass-user@lists.osgeo.org
https://lists.osgeo.org/mailman/listinfo/grass-user

Re: [GRASS-user] v.class.mlR Error

2018-06-07 Thread Jamille Haarloo
I see. That was one of the big problems then. I knew it used R-resources
but wasn't sure the software itself was required. Must have overlooked it
in the documentation and other user- topics as I was reading about the
module.

Now I have R-3.5, but still get an error:

v.class.mlR --overwrite training_map=Training_Ben2@haarlooj_Ben_Test
segments_file=C:\Users\haarlooj\Documents\CELOS\Stats_Training_Ben_test
raster_segments_map=best_myregion1_at_haarlooj_Ben_Test_rank1@haarlooj_Ben_Test
train_class_column=code output_class_column=vote output_prob_column=prob
classifiers=svmRadial,rf,C5.0 folds=5 partitions=10 tunelength=10
weighting_modes=smv,qbwwv weighting_metric=accuracy
classification_results=C:\Users\haarlooj\Documents\CELOS\v.class.mIRR_optional_output\Ben_test_Classifier-results
accuracy_file=C:\Users\haarlooj\Documents\CELOS\v.class.mIRR_optional_output\Ben_test_Classifier-accuracy
model_details=C:\Users\haarlooj\Documents\CELOS\v.class.mIRR_optional_output\Ben_test_Classifier-module-runs
bw_plot_file=C:\Users\haarlooj\Documents\CELOS\v.class.mIRR_optional_output\Ben_test_Classifier-performance
r_script_file=C:\Users\haarlooj\Documents\CELOS\v.class.mIRR_optional_output\Ben_test_R_script
processes=3
Running R now. Following output is R output.
During startup - Warning messages:
1: Setting LC_CTYPE=en_US.cp1252 failed
2: Setting LC_COLLATE=en_US.cp1252 failed
3: Setting LC_TIME=en_US.cp1252 failed
4: Setting LC_MONETARY=en_US.cp1252 failed
Loading required package: caret
Loading required package: lattice
Loading required package: ggplot2
Loading required package: foreach
Loading required package: iterators
Loading required package: parallel
During startup - Warning messages:
1: Setting LC_CTYPE=en_US.cp1252 failed
2: Setting LC_COLLATE=en_US.cp1252 failed
3: Setting LC_TIME=en_US.cp1252 failed
4: Setting LC_MONETARY=en_US.cp1252 failed
During startup - Warning messages:
1: Setting LC_CTYPE=en_US.cp1252 failed
2: Setting LC_COLLATE=en_US.cp1252 failed
3: Setting LC_TIME=en_US.cp1252 failed
4: Setting LC_MONETARY=en_US.cp1252 failed
During startup - Warning messages:
1: Setting LC_CTYPE=en_US.cp1252 failed
2: Setting LC_COLLATE=en_US.cp1252 failed
3: Setting LC_TIME=en_US.cp1252 failed
4: Setting LC_MONETARY=en_US.cp1252 failed
Warning message:
In nominalTrainWorkflow(x = x, y = y, wts = weights, info = trainInfo,  :
  There were missing values in resampled performance measures.
Error in eval(predvars, data, env) : object 'id' not found
Calls: data.frame ... predict.train -> model.frame -> model.frame.default
-> eval -> eval
Execution halted
ERROR: There was an error in the execution of the R script.
Please check the R output.
Error in atexit._run_exitfuncs:
Traceback (most recent call last):
  File "C:\Program Files\GRASS GIS
7.4.0\Python27\lib\atexit.py", line 24, in _run_exitfuncs
func(*targs, **kargs)
  File "C:\Users\haarlooj\AppData\Roaming\GRASS7\addons/scri
pts/v.class.mlR.py", line 236, in cleanup
gscript.try_remove(model_output_csv)
NameError: global name 'model_output_csv' is not defined
Error in sys.exitfunc:
Traceback (most recent call last):
  File "C:\Program Files\GRASS GIS
7.4.0\Python27\lib\atexit.py", line 24, in _run_exitfuncs
func(*targs, **kargs)
  File "C:\Users\haarlooj\AppData\Roaming\GRASS7\addons/scri
pts/v.class.mlR.py", line 236, in cleanup
gscript.try_remove(model_output_csv)
NameError: global name 'model_output_csv' is not defined

On Wed, Jun 6, 2018 at 4:05 AM, Moritz Lennert  wrote:

> On 04/06/18 15:42, Jamille Haarloo wrote:
>
>> It did not work out from another station, but I tried today from the same
>> computer.
>> I first reinstalled the module.
>> This is the command output:
>>
>> Running R now. Following output is R output.
>> Traceback (most recent call last):
>>File "C:\Users\haarlooj\AppData\Roaming\GRASS7\addons/scri
>> pts/v.class.mlR.py ", line 632, in 
>>  main()
>>File "C:\Users\haarlooj\AppData\Roaming\GRASS7\addons/scri
>> pts/v.class.mlR.py ", line 569, in main
>>  subprocess.check_call(['Rscript', r_commands],
>> stderr=subprocess.STDOUT, )
>>File "C:\Program Files\GRASS GIS
>> 7.4.0\Python27\lib\subprocess.py", line 537, in check_call
>>  retcode = call(*popenargs, **kwargs)
>>File "C:\Program Files\GRASS GIS
>> 7.4.0\Python27\lib\subprocess.py", line 524, in call
>>  return Popen(*popenargs, **kwargs).wait()
>>File "C:\Program Files\GRASS GIS
>> 7.4.0\Python27\lib\subprocess.py", line 711, in __init__
>>  errread, errwrite)
>>File "C:\Program Files\GRASS GIS
>> 7.4.0\Python27\lib\subprocess.py", line 948, in
>> _execute_child
>>  startupinfo)
>> WindowsError: [Error 2] The system cannot find the file
>> specified
>>
>
> I still have the feeling that the system cannot find Rscript.
>
> Can you open the GRASS terminal window and type 'Rscript' to see if it can
> be found ?
>
> Just to be sure: you do have R installed ?
>
> Error in 

Re: [GRASS-user] i.maxlik yellow map

2018-06-07 Thread Moritz Lennert

On 07/06/18 11:00, Paul Shapley wrote:

Hi Users,

I am re-posting this to see if anyone else may have some suggestions. My 
thanks to Veronica who suggested checking 'g.region' but this appears to 
be fine. I've also tried the same on Grass 7.2.2/7.4 and nightly but all 
produce the same result.


signature file content:-

#
#Bare_peat
2667621
60.4424 75.4845 87.2587
716.802
831.829 1086.7
869.829 1105.97 1179.92
0 0 0

--

Original message:-

I have created a signature file in the g.gui.iclass (supervised 
classification module) and it looks correct in 'preview' and when i 
inspected the 'sig' file. When the signature is applied using 'i.maxlik' 
the result is a solid yellow colour as a layer. What could have gone wrong?


Do you only have one class ? Then it is normal that all pixels will be 
in the same class. You could ask i.maxlik to output the reject map and 
then define a threshold in that map as a mask to show those areas whose 
probability of being in the class is sufficiently high.


Moritz
___
grass-user mailing list
grass-user@lists.osgeo.org
https://lists.osgeo.org/mailman/listinfo/grass-user

[GRASS-user] i.maxlik yellow map

2018-06-07 Thread Paul Shapley
Hi Users,

I am re-posting this to see if anyone else may have some suggestions. My
thanks to Veronica who suggested checking 'g.region' but this appears to be
fine. I've also tried the same on Grass 7.2.2/7.4 and nightly but all
produce the same result.

signature file content:-

#
#Bare_peat
2667621
60.4424 75.4845 87.2587
716.802
831.829 1086.7
869.829 1105.97 1179.92
0 0 0


--

Original message:-

I have created a signature file in the g.gui.iclass (supervised
classification module) and it looks correct in 'preview' and when i
inspected the 'sig' file. When the signature is applied using 'i.maxlik'
the result is a solid yellow colour as a layer. What could have gone wrong?


-- 
*Paul J. Shapley *MSc CGeog (GIS) FRGS
___
grass-user mailing list
grass-user@lists.osgeo.org
https://lists.osgeo.org/mailman/listinfo/grass-user

Re: [GRASS-user] Using GRASS on a ply file

2018-06-07 Thread marion-brunet

The 
complete traceback is : 

Traceback (most recent call last):
 File "C:\PROGRA~1\QGIS2~1.18\apps\grass\grass-7.2.2\gui\wx
python\gui_core\prompt.py", line 416, in OnKeyPressed

self._runCmd(self.GetCurLine()[0].strip())
 File "C:\PROGRA~1\QGIS2~1.18\apps\grass\grass-7.2.2\gui\wx
python\gui_core\prompt.py", line 123, in _runCmd

self.promptRunCmd.emit(cmd=cmd)
 File "C:\PROGRA~1\QGIS2~1.18\apps\grass\grass-7.2.2\etc\py
thon\grass\pydispatch\signal.py", line 229, in emit

dispatcher.send(signal=self, *args, **kwargs)
 File "C:\PROGRA~1\QGIS2~1.18\apps\grass\grass-7.2.2\etc\py
thon\grass\pydispatch\dispatcher.py", line 349, in send

**named
 File "C:\PROGRA~1\QGIS2~1.18\apps\grass\grass-7.2.2\etc\py
thon\grass\pydispatch\robustapply.py", line 60, in
robustApply

return receiver(*arguments, **named)
 File "C:\PROGRA~1\QGIS2~1.18\apps\grass\grass-7.2.2\gui\wx
python\gui_core\goutput.py", line 114, in lambda

self._gconsole.RunCmd(command=cmd))
 File "C:\PROGRA~1\QGIS2~1.18\apps\grass\grass-7.2.2\gui\wx
python\core\gconsole.py", line 482, in RunCmd

task = GUI(show=None).ParseCommand(command)
 File "C:\PROGRA~1\QGIS2~1.18\apps\grass\grass-7.2.2\gui\wx
python\gui_core\forms.py", line 2816, in ParseCommand

raise gcmd.GException(e.value)
AttributeError
:
'exceptions.UnicodeDecodeError' object has no attribute
'value'


I am using GRASS 7.2.2, and I am on Windows 10‌

De : "Moritz 
Lennert"
A : marion-bru...@sfr.fr,"Grass-user"
Envoyé: jeudi 7 juin 2018 09:56
Objet : Re: [GRASS-user] Using GRASS on a ply file

On 07/06/18 09:10, Moritz Lennert wrote:


 Hi Marion,

 Please keep the discussion on the list.


 Am 7. Juni 2018 06:55:33 MESZ schrieb marion-bru...@sfr.fr:
 Thank you very much for your answer!

 However, I have another ERROR when I use the command , even when
 I remove the white spaces. My command is:

 v.in.ply input=path\file.ply output=test

 and I get:



 AttributeError

 :

 'exceptions.UnicodeDecodeError' object has no attribute

 'value'



 I really don't know what it means..


 I have the feeling that you provide us with only an extract of the 
message. Please provide the complete traceback and the exact command you 
used.

 If you launch the command via the GUI you can find the history of commands 
in the console window (up arrow to see previous commands) or you can use the 
'copy' button in the module GUI to copy the command line just before pushing 
run.

Sorry, on my phone I didn't see correctly that you actually provided the
command.

But the traceback is still necessary.

Also which version of GRASS GIS are you using ? On which OS ?

Moritz



___
grass-user mailing list
grass-user@lists.osgeo.org
https://lists.osgeo.org/mailman/listinfo/grass-user

Re: [GRASS-user] Using GRASS on a ply file

2018-06-07 Thread Moritz Lennert

On 07/06/18 09:10, Moritz Lennert wrote:



Hi Marion,

Please keep the discussion on the list.


Am 7. Juni 2018 06:55:33 MESZ schrieb marion-bru...@sfr.fr:

Thank you very much for your answer!
However, I have another ERROR when I use the command , even when
I remove the white spaces. My command is:
v.in.ply input=path\file.ply output=test
and I get:

AttributeError
:
'exceptions.UnicodeDecodeError' object has no attribute
'value'

I really don't know what it means..



I have the feeling that you provide us with only an extract of the message. 
Please provide the complete traceback and the exact command you used.

If you launch the command via the GUI you can find the history of commands in 
the console window (up arrow to see previous commands) or you can use the 
'copy' button in the module GUI to copy the command line just before pushing 
run.


Sorry, on my phone I didn't see correctly that you actually provided the 
command.


But the traceback is still necessary.

Also which version of GRASS GIS are you using ? On which OS ?

Moritz

___
grass-user mailing list
grass-user@lists.osgeo.org
https://lists.osgeo.org/mailman/listinfo/grass-user

Re: [GRASS-user] Using GRASS on a ply file

2018-06-07 Thread Moritz Lennert


Hi Marion,

Please keep the discussion on the list.


Am 7. Juni 2018 06:55:33 MESZ schrieb marion-bru...@sfr.fr:
>Thank you very much for your answer!
>However, I have another ERROR when I use the command , even when
>I remove the white spaces. My command is:
>v.in.ply input=path\file.ply output=test
>and I get:
>
>AttributeError
>:
>'exceptions.UnicodeDecodeError' object has no attribute
>'value'
>
>I really don't know what it means..


I have the feeling that you provide us with only an extract of the message. 
Please provide the complete traceback and the exact command you used. 

If you launch the command via the GUI you can find the history of commands in 
the console window (up arrow to see previous commands) or you can use the 
'copy' button in the module GUI to copy the command line just before pushing 
run.

Moritz




.
>
>
>De :
>"Moritz Lennert"
>A : marion-bru...@sfr.fr,grass-user@lists.osgeo.org
>Envoyé: mercredi 6 juin 2018 13:49
>Objet : Re: [GRASS-user] Using GRASS on a ply file
>
>Hello Marion and welcome to the GRASS community
>!
>
>On 06/06/18 12:03, marion-brunet wrote:
> ‌Hello,
>
> I have two point clouds of an open-pit mine taken a few weeks
>apart. I
> would like to estimate the volume that has been excavated during
>that
> interval, i.e the difference between the two point clouds. For
>that I
> used CloudCompare to get a new point cloud corresponding to the
>distance
> between the first two clouds. Now I need to know the
> Z information of each point of this last cloud. I want to use
>GRASS to
> do it, but I can't find a method that is working...
>
> The first method I was suggested is to use*v.in.ply* and
>then*r.report*,
> but I have several issues:
> - When I use *v.in.ply* I have this error: *ERROR: Wrong number
>of
> properties*. Is it a common error? How can I solve it? I have
>attached a
> part of my .ply file to this email so you can see how it is
>formatted.
>
>[Marion sent me an extract of the data offline]
>
>It seems that every data line in your .ply file has a trailing
>white
>space at the end. This is counted as an extra property and explains
>the
>message. When I erase these trailing white spaces, I don't get this
>error anymore. I don't know the format enough to know whether this is
>an
>issue with v.in.ply not tolerating these white spaces, or an issue
>with
>the program creating the .ply file which shouldn't add a white space
>at
>the end.
>
> - After importing my .ply file, I would like to use something
>like
> *r.report *to get the value of each point. However, I don't know
>how to
> convert my file into a format that this command will work with...
>Could
> you help me?
>
>r.report gives you a summary report about a raster, not a value by
>point.
>
>Generally, if you want to cover an entire zone (your entire
>open-pit
>mine, you will probably want to go the raster way. So, with the
>points
>imported with v.in.ply, you could do something like the following:
>
>g.region vect=points res=1 -ap #set the xy extension of your raster
>to
>the xy extension of the point cloud, and the resolution to 1 (you
>have
>to adjust that depending on the size of your zone)
>v.surf.rst in=points out=mined_surface #interpolate a continous
>height
>surface across the entire area
>r.volume mined_surface #calculate volume
>
>
>
> The second method I read about is *r.in.xyz*, but I don't know how
>it
> works, or what format of input I should use. Can I have more
>information
> about it?
>
>r.in.xyz is another option. It takes XYZ data and aggregates it by
>pixel
>using a user-determined aggregation function. If you have a higher
>point
>density than your target resolution, then you could use just the
>data
>part of your .ply file, i.e. the part which has lines such as
>
>[...]
>97.08156 121.6752 99.27532 0 0 0 12 9 2 255
>97.08911 121.718 99.50104 0 0 0 46 40 24 255
>93.25124 115.0372 102.7768 0 0 0 62 53 20 255
>97.96031 123.8419 103.3994 0 0 0 150 184 211 255
>97.10852 121.6866 99.42679 0 0 0 9 9 0 255
>97.11196 121.7049 99.59032 0 0 0 19 12 4 255
>97.13282 121.688 99.43796 0 0 0 38 28 16 255
>97.17371 121.7129 99.41496 0 0 0 43 31 15 255
>97.1614 121.6344 99.27691 0 0 0 9 6 1 255
>[...]
>
>Knowing that there are 17 lines of other info at the top of the
>extract
>you sent me, I could import the data using:
>
>r.in.xyz in=BROL/BascFin2-extrait.ply skip=17 out=xyz_mean
>method=mean
>sep=space
>
>However, as the extent of the cloud is very large, but there are only
>a
>few points, most of the pixels in the resulting raster are filled
>with
>NULL values as there are no points falling within these pixels.
>
>Moritz
>
>
>
___
grass-user mailing list
grass-user@lists.osgeo.org
https://lists.osgeo.org/mailman/listinfo/grass-user