[GRASS-user] Exporting Grass raster to netCDF via r.out.gdal gives '0' z-values

[Eric Patton]

Hi guys, I am having trouble exporting any Grass raster to netCDF format and 
preserving the z-values.

Ex:

~> r.info GEBCO2014_Atlantic_Canada_100m -g | grep datatype
datatype=DCELL

~> g.region -p

projection: 99 (unnamed)
zone:   0
datum:  wgs84
ellipsoid:  wgs84
north:  838800
south:  -313000
west:   -4400
east:   1225500
nsres:  100
ewres:  100
rows:   11518
cols:   12299
cells:  141659882

~> r.out.gdal in=GEBCO2014_Atlantic_Canada_100m format=netCDF 
out=GEBCO2014_Grand_Banks_100m.nc type=Float64 --v --o

Exporting to GDAL data type: Float64
Checking GDAL data type and nodata value...
100%
Using GDAL data type 
Exporting raster data to netCDF format...
100%
r.out.gdal complete. File  created.

Now, checking the output .nc file with GMT's grdinfo:

~> grdinfo GEBCO2014_Grand_Banks_100m.nc
GEBCO2014_Grand_Banks_100m.nc: Title: GDAL Band Number 1
GEBCO2014_Grand_Banks_100m.nc: Command: Wed Dec 04 15:25:18 2019: GDAL Create( 
GEBCO2014_Grand_Banks.nc, ... )
GEBCO2014_Grand_Banks_100m.nc: Remark:
GEBCO2014_Grand_Banks_100m.nc: Pixel node registration used [Cartesian grid]
GEBCO2014_Grand_Banks_100m.nc: Grid file format: nd = GMT netCDF format (64-bit 
float), CF-1.7
GEBCO2014_Grand_Banks_100m.nc: x_min: -4400 x_max: 1225500 x_inc: 100 name:  
n_columns: 12299
GEBCO2014_Grand_Banks_100m.nc: y_min: -313000 y_max: 838800 y_inc: 100 name:  
n_rows: 11518
GEBCO2014_Grand_Banks_100m.nc: z_min: 0 z_max: 0 name: GDAL Band Number 1
GEBCO2014_Grand_Banks_100m.nc: scale_factor: 1 add_offset: 0
GEBCO2014_Grand_Banks_100m.nc: format: classic

Zeros for the z-range. gdalinfo also fails to reveal any z-range info. I'm not 
sure what I'm doing wrong here.

System Info:
GRASS version: 7.8.1 (installed from UbuntuGIS-unstable PPA)
Code revision: exported
Build date: 2019-12-03
Build platform: x86_64-pc-linux-gnu
GDAL: 3.0.2
PROJ: 6.2.1
GEOS: 3.8.0
SQLite: 3.22.0
Python: 3.6.9
wxPython: 4.0.1
Platform: Linux-5.3.0-24-generic-x86_64-with-LinuxMint-19.2-tina

Any pointers most appreciated!

~ Eric.

--
Eric Patton
GPG ID: 
rsa4096/0x[B72AE021](http://keys.gnupg.net/pks/lookup?op=get&search=0x24CC056BB72AE021)
Fingerprint: 1B0F 7AFD 9B5B 0E59 DCBC 36B4 24CC 056B B72A E021___
grass-user mailing list
grass-user@lists.osgeo.org
https://lists.osgeo.org/mailman/listinfo/grass-user

Re: [GRASS-user] GRASS 7.8.1 MAC binaries

[Markus Neteler]

Hi,

On Tue, Nov 26, 2019 at 5:57 PM cesaregal  wrote:
>
> I see in the Osgeo downloading page that it is available for downloading
> Grass GIS 7.6.x (old stable) for MAC OS X.
>
> I don't see it available for downloading Grass GIS 7.8.1 for MAC OS X (with
> PROJ 6 and GDAL 3 support).
> Can someone tell me if this version is due to be released?

Seems that the efforts are currently stuck here:

Ticket URL: 

More ideas are welcome.

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

Re: [GRASS-user] SQL: generating numeric class numbers from class text labels?

[Markus Neteler]

Thanks for your answers.
In fact I need it in Python...

Best
Markus


Micha Silver  schrieb am Mi., 4. Dez. 2019, 18:57:

> How about doing this in R? The labels will be read into R as factors, and
> the factor levels can easily be extracted as numbers.
>
>
> Something like this:
>
>
> micha@tp480:~$ v.info -c stations
> Displaying column types/names for database connection of layer <1>:
> INTEGER|cat
> INTEGER|station_num
> TEXT|station_he
> TEXT|station_en
> TEXT|type
> INTEGER|x_coord
> INTEGER|y_coord
> DOUBLE PRECISION|long
> DOUBLE PRECISION|lat
> INTEGER|elev
> TEXT|date_open
> DOUBLE PRECISION|dist
> DOUBLE PRECISION|azim
>
>
> micha@tp480:~$ R
>
>
> R version 3.5.2 (2018-12-20) -- "Eggshell Igloo"
> Copyright (C) 2018 The R Foundation for Statistical Computing
> Platform: x86_64-pc-linux-gnu (64-bit)
> .
>
> > library(rgrass7)
> Loading required package: XML
> GRASS GIS interface loaded with GRASS version: GRASS 7.6.0 (2019)
> and location: ITM
> > use_sf()
> > stations = readVECT("stations")
> WARNING: Vector map  is 3D. Use format specific layer creation
>  options (parameter 'lco') to export   (default).
> Exporting 94 features...
>  100%
> .
>
> > stations['new_station_num'] = as.numeric(stations$station_en)
> > stations$new_station_num
>  [1] 71 26  6 55 54 63  7  8 31 30 46 84 92 38 32 88 27 12 67 62 47 33 53
> 76 89
> [26]  2 86 11 40 65 64 45 13 85 60 59  1 74 73 22 19 15 39 50 56 14 44 23
> 36 83
> [51] 41 42 43 18 17 75 16 82 81 37 48 28 87  3 66 10 34 91 61 93 94 72  5
> 4 68
> [76] 78 77  9 29 51 58 57 49 52 24 25 80 79 35 70 69 90 21 20
>
> > writeVECT(SDF=stations, vname="new_stations")
>
>
> Best regards, Micha
>
>
> On 04/12/2019 19:11, Markus Neteler wrote:
>
> Hi,
>
> I have a landuse map with text labels (forest, street, ...). Forr.learn.ml I 
> need to have them as numeric classes.
> It is not important for me which number is assigned but I search for
> an automated solution, i.e. SQL statement unless there is a different
> way.
>
> So:
>
> cat|label|label_int
> 1|forest|1
> 2|forest|1
> 3|street|2
> 4|forest|1
> 5|street|2
> 6|urban|3
> ...
>
> I guess I have done that already some years ago but I can't remember
> the trick :-)
>
> thanks for a hint,
> Markus
> ___
> grass-user mailing 
> listgrass-user@lists.osgeo.orghttps://lists.osgeo.org/mailman/listinfo/grass-user
>
> --
> Micha Silver
> Ben Gurion Univ.
> Sde Boker, Remote Sensing Lab
> cell: +972-523-665918https://orcid.org/-0002-1128-1325
>
>
___
grass-user mailing list
grass-user@lists.osgeo.org
https://lists.osgeo.org/mailman/listinfo/grass-user

Re: [GRASS-user] SQL: generating numeric class numbers from class text labels?

[Micha Silver]

How about doing this in R? The labels will be read into R as factors, 
and the factor levels can easily be extracted as numbers.


Something like this:


micha@tp480:~$ v.info -c stations
Displaying column types/names for database connection of layer <1>:
INTEGER|cat
INTEGER|station_num
TEXT|station_he
TEXT|station_en
TEXT|type
INTEGER|x_coord
INTEGER|y_coord
DOUBLE PRECISION|long
DOUBLE PRECISION|lat
INTEGER|elev
TEXT|date_open
DOUBLE PRECISION|dist
DOUBLE PRECISION|azim


micha@tp480:~$ R


R version 3.5.2 (2018-12-20) -- "Eggshell Igloo"
Copyright (C) 2018 The R Foundation for Statistical Computing
Platform: x86_64-pc-linux-gnu (64-bit)
.

> library(rgrass7)
Loading required package: XML
GRASS GIS interface loaded with GRASS version: GRASS 7.6.0 (2019)
and location: ITM
> use_sf()
> stations = readVECT("stations")
WARNING: Vector map  is 3D. Use format specific layer creation
 options (parameter 'lco') to export  stations['new_station_num'] = as.numeric(stations$station_en)
> stations$new_station_num
 [1] 71 26  6 55 54 63  7  8 31 30 46 84 92 38 32 88 27 12 67 62 47 33 
53 76 89
[26]  2 86 11 40 65 64 45 13 85 60 59  1 74 73 22 19 15 39 50 56 14 44 
23 36 83
[51] 41 42 43 18 17 75 16 82 81 37 48 28 87  3 66 10 34 91 61 93 94 72  
5  4 68

[76] 78 77  9 29 51 58 57 49 52 24 25 80 79 35 70 69 90 21 20

> writeVECT(SDF=stations, vname="new_stations")

Best regards, Micha


On 04/12/2019 19:11, Markus Neteler wrote:

Hi,

I have a landuse map with text labels (forest, street, ...). For
r.learn.ml I need to have them as numeric classes.
It is not important for me which number is assigned but I search for
an automated solution, i.e. SQL statement unless there is a different
way.

So:

cat|label|label_int
1|forest|1
2|forest|1
3|street|2
4|forest|1
5|street|2
6|urban|3
...

I guess I have done that already some years ago but I can't remember
the trick :-)

thanks for a hint,
Markus
___
grass-user mailing list
grass-user@lists.osgeo.org
https://lists.osgeo.org/mailman/listinfo/grass-user


--
Micha Silver
Ben Gurion Univ.
Sde Boker, Remote Sensing Lab
cell: +972-523-665918
https://orcid.org/-0002-1128-1325

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

Re: [GRASS-user] SQL: generating numeric class numbers from class text labels?

[Micha Silver]

  
  
How about doing this in R? The labels will be read into R as
  factors, and the factor levels can easily be extracted as numbers.


Something like this:


micha@tp480:~$ v.info -c stations
  Displaying column types/names for database connection of layer
  <1>:
  INTEGER|cat
  INTEGER|station_num
  TEXT|station_he
  TEXT|station_en
  TEXT|type
  INTEGER|x_coord
  INTEGER|y_coord
  DOUBLE PRECISION|long
  DOUBLE PRECISION|lat
  INTEGER|elev
  TEXT|date_open
  DOUBLE PRECISION|dist
  DOUBLE PRECISION|azim


micha@tp480:~$ R



R version 3.5.2 (2018-12-20) -- "Eggshell Igloo"
  Copyright (C) 2018 The R Foundation for Statistical Computing
  Platform: x86_64-pc-linux-gnu (64-bit)
  .
  
  > library(rgrass7)
  Loading required package: XML
  GRASS GIS interface loaded with GRASS version: GRASS 7.6.0 (2019)
  and location: ITM
  > use_sf()
  > stations = readVECT("stations")
  WARNING: Vector map  is 3D. Use format specific
  layer creation
   options (parameter 'lco') to export 
   (default).
  Exporting 94 features...
   100%
  .

> stations['new_station_num'] =
  as.numeric(stations$station_en)
  > stations$new_station_num
   [1] 71 26  6 55 54 63  7  8 31 30 46 84 92 38 32 88 27 12 67 62
  47 33 53 76 89
  [26]  2 86 11 40 65 64 45 13 85 60 59  1 74 73 22 19 15 39 50 56
  14 44 23 36 83
  [51] 41 42 43 18 17 75 16 82 81 37 48 28 87  3 66 10 34 91 61 93
  94 72  5  4 68
  [76] 78 77  9 29 51 58 57 49 52 24 25 80 79 35 70 69 90 21 20
  

> writeVECT(SDF=stations, vname="new_stations")
   

Best regards, Micha


On 04/12/2019 19:11, Markus Neteler
  wrote:


  Hi,

I have a landuse map with text labels (forest, street, ...). For
r.learn.ml I need to have them as numeric classes.
It is not important for me which number is assigned but I search for
an automated solution, i.e. SQL statement unless there is a different
way.

So:

cat|label|label_int
1|forest|1
2|forest|1
3|street|2
4|forest|1
5|street|2
6|urban|3
...

I guess I have done that already some years ago but I can't remember
the trick :-)

thanks for a hint,
Markus
___
grass-user mailing list
grass-user@lists.osgeo.org
https://lists.osgeo.org/mailman/listinfo/grass-user

-- 
Micha Silver
Ben Gurion Univ.
Sde Boker, Remote Sensing Lab
cell: +972-523-665918
https://orcid.org/-0002-1128-1325
  

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

Re: [GRASS-user] SQL: generating numeric class numbers from class text labels?

[Hernán De Angelis]

The variable to increment should be of course as corrected below:


# get the unique classes

SELECT distinct(label) FROM table;

-> fetch results in an array "txt_label"

# iterate over array and insert new integer labels in table

*int_label = 0 *

foreach txt_label {

    INSERT INTO table.label_int VALUES int_label WHERE label = 
'txt_label'


*int_label++*

}


/H.





On 2019-12-04 18:11, Markus Neteler wrote:

Hi,

I have a landuse map with text labels (forest, street, ...). For
r.learn.ml I need to have them as numeric classes.
It is not important for me which number is assigned but I search for
an automated solution, i.e. SQL statement unless there is a different
way.

So:

cat|label|label_int
1|forest|1
2|forest|1
3|street|2
4|forest|1
5|street|2
6|urban|3
...

I guess I have done that already some years ago but I can't remember
the trick :-)

thanks for a hint,
Markus
___
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
___
grass-user mailing list
grass-user@lists.osgeo.org
https://lists.osgeo.org/mailman/listinfo/grass-user

Re: [GRASS-user] SQL: generating numeric class numbers from class text labels?

[Hernán De Angelis]

Hi

I cannot think of a short and simple SQL "one liner" statement but can 
think of a short script that will do that using SQL statements.


In pseudo-code, translate to your favorite language:


# get the unique classes

SELECT distinct(label) FROM table;

-> fetch results in an array "txt_label"

# iterate over array and insert new integer labels in table

int_label = 0

foreach txt_label {

    INSERT INTO table.label_int VALUES int_label WHERE label = 'txt_label'

    n++

}


This may seem obvious to you, but since you asked :-)

Cheers and good luck


Hernán





On 2019-12-04 18:11, Markus Neteler wrote:

Hi,

I have a landuse map with text labels (forest, street, ...). For
r.learn.ml I need to have them as numeric classes.
It is not important for me which number is assigned but I search for
an automated solution, i.e. SQL statement unless there is a different
way.

So:

cat|label|label_int
1|forest|1
2|forest|1
3|street|2
4|forest|1
5|street|2
6|urban|3
...

I guess I have done that already some years ago but I can't remember
the trick :-)

thanks for a hint,
Markus
___
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

[GRASS-user] SQL: generating numeric class numbers from class text labels?

[Markus Neteler]

Hi,

I have a landuse map with text labels (forest, street, ...). For
r.learn.ml I need to have them as numeric classes.
It is not important for me which number is assigned but I search for
an automated solution, i.e. SQL statement unless there is a different
way.

So:

cat|label|label_int
1|forest|1
2|forest|1
3|street|2
4|forest|1
5|street|2
6|urban|3
...

I guess I have done that already some years ago but I can't remember
the trick :-)

thanks for a hint,
Markus
___
grass-user mailing list
grass-user@lists.osgeo.org
https://lists.osgeo.org/mailman/listinfo/grass-user

Re: [GRASS-user] Problem with v.class.mIR

[Moritz Lennert]

[Guisseppe sent me the files].

The issue is your training file.

As the manual page states:


"The user provides a set of objects (or segments) to be classified, 
including all feature variables describing these object, and a set of 
objects to be used as training data, including the same feature 
variables as those describing the unknown objects, plus one additional 
column indicating the class each training falls into. The training data 
can, but does not have to be, a subset of the set of objects to be 
classified."


IOW, the training data file has to contain the same feature variables as 
the data you wish to classify.


Your data is:

$ head -1 segmentation_map.csv
cat,area,perimeter,compact_circle,fd

$ head -1 training_map.csv
cat,id,cate

i.e. you training data only includes the id and the class, but not the 
variables describing the objects.


The reason this is organised in this manner is that you might have 
training data which is independent of the data you wish to classify, 
i.e. the ids of the training data might not be linked to the ids in the 
object data.


Another issue is that in your command, you do not provide a good name 
for the column that contains the actual class label in the training data:


v.class.mlR segments_map=gbbb@PERMANENT training_map=tainning@PERMANENT 
train_class_column=fd [...]


v.class.mlR segments_map=gbbb training_map=tainning 
train_class_column=class [...]


Your training data does not contain a column 'fd' nor a column 'class'. 
I suppose that in your data it is the column 'cate'.


Moritz

On 4/12/19 10:38, Moritz Lennert wrote:

Sorry for ignoring you up to now, just to much on the platter right now.

Would it be feasible for you to make the data available, privately if
necessary ? The easiest (and lightest) would be just the csv output of
v.db.select on the training_map and the segments_maps.

Moritz

On 2/12/19 16:31, Giuseppe Cillis wrote:

No one can help me?
GC

Il giorno mer 20 nov 2019 alle ore 14:06 Giuseppe Cillis
mailto:giucil...@gmail.com>> ha scritto:

 Thanks for the answer._
 _
 _
 _
 _v.class.mlR segments_map=gbbb@PERMANENT
 training_map=tainning@PERMANENT train_class_column=fd
 output_class_column=vote output_prob_column=prob folds=5
 partitions=10 tunelength=10 weighting_metric=accuracy_

 _and the same error with this command:
 _
 _v.class.mlR segments_map=gbbb training_map=tainning
 train_class_column=class weighting_mode=smv,swv,qbwwv -i_

 The vector used was elaborated with i.segment.stats module.
 An extract is here:
 catareaperimeter   fd
 1  958 208 1.555034
 2  24  22  1.945217
 3  160 80  1.726846
 4  2036242 1.440904
 5  25  26  2.024344
 6  222 96  1.68966
 7  843510121.530879


 /cat /type Integer (9-0)
 /fd /type double/real (25-9)
 For the training area (t/ainning/), I use a simple vector with
 numerical categories (named/"cate"/) and there are 11 features (just
 for a trying). The type is in integer64 (10)

 Il giorno mer 20 nov 2019 alle ore 12:18 Moritz Lennert
 mailto:mlenn...@club.worldonline.be>>
 ha scritto:

 On 20/11/19 11:55, Giuseppe Cillis wrote:
  > Hi,
  > I'm trying to use this module for classification of an old
 aerial photos.
  > After a segmentation (i.segment and i.segment.stats), I would
 to use a
  > machine learning approach for the real classification.
  > I tried with v.class.mIR which use also R.
  > But there is an error and I don't know how to solve it:
  > */Durante l'avvio - Warning messages:
  > 1: Setting LC_CTYPE=it_IT.cp1252 failed
  > 2: Setting LC_COLLATE=it_IT.cp1252 failed
  > 3: Setting LC_TIME=it_IT.cp1252 failed
  > 4: Setting LC_MONETARY=it_IT.cp1252 failed
  > Carico il pacchetto richiesto: caret
  > Carico il pacchetto richiesto: lattice
  > Carico il pacchetto richiesto: ggplot2
  > Error in `$<-.data.frame`(`*tmp*`, fd, value = integer(0)) :
  >    replacement has 0 rows, data has 11
  > Calls: $<- -> $<-.data.frame/*
  > */
  > /*
  > A part is in italian, I'm sorry.
  > I'm a beginner in grass and R.

 R messages are often not very explicit, unfortunately.

 Could you provide us with the exact command line you use (If you
 use the
 GUI, you can get the command line by clicking the 'Copy' button
 once
 you've filled out all the parameters, or in the history of the
 'Console'
 in the Layer manager.) and an extract of the attribute data you
 feed
 into the module ?

 Moritz


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

Re: [GRASS-user] Datum not recognized by Grass

[Zoltan Szecsei]

Hi Helmut,
Thanks for your comments.

I installed everything with OSGeo4W64, and QGIS get the EPSG:2932 but 
Grass not.


Perhaps I have a PATH or some other setting problem?
Perhaps let me know what search paths Grass uses for proj4, and what 
proj files and locations I should scan for.


Regards,
Zoltan

On 2019/12/04 01:30, Helmut Kudrnovsky wrote:

ignore the warning and use GRASS with PROJ6, granted that authority name

(e.g. EPSG) and authority code (e.g. 2932) are known for both CRS's in case
of reprojection

this issue is already by GRASS with PROJ6, see


GRASS version: 7.8.1
Code revision: c865432c9
Build date: 2019-11-10
Build platform: x86_64-w64-mingw32
GDAL: 3.0.2
PROJ: 6.2.1  <=
GEOS: 3.8.0
SQLite: 3.29.0
Python: 3.7.0
wxPython: 4.0.7
Platform: Windows-10-10.0.18362-SP0 (OSGeo4W)


and the output from the underlying PROJ: 6.2.1


C:\>projinfo EPSG:2932 -o PROJ,WKT2_2018
PROJ.4 string:
+proj=tmerc +lat_0=24.45 +lon_0=51.21667 +k=0.9 +x_0=20
+y_0=30 +ellps=intl
+towgs84=-119.4248,-303.65872,-11.00061,1.164298,0.174458,1.096259,3.657065
+units=m +no_defs +type=crs

WKT2_2018 string:
PROJCRS["QND95 / Qatar National Grid",
 BASEGEOGCRS["QND95",
 DATUM["Qatar National Datum 1995",
 ELLIPSOID["International 1924",6378388,297,
 LENGTHUNIT["metre",1]]],
 PRIMEM["Greenwich",0,
 ANGLEUNIT["degree",0.0174532925199433]],
 ID["EPSG",4614]],
 CONVERSION["Qatar National Grid",
 METHOD["Transverse Mercator",
 ID["EPSG",9807]],
 PARAMETER["Latitude of natural origin",24.45,
 ANGLEUNIT["degree",0.0174532925199433],
 ID["EPSG",8801]],
 PARAMETER["Longitude of natural origin",51.21667,
 ANGLEUNIT["degree",0.0174532925199433],
 ID["EPSG",8802]],
 PARAMETER["Scale factor at natural origin",0.9,
 SCALEUNIT["unity",1],
 ID["EPSG",8805]],
 PARAMETER["False easting",20,
 LENGTHUNIT["metre",1],
 ID["EPSG",8806]],
 PARAMETER["False northing",30,
 LENGTHUNIT["metre",1],
 ID["EPSG",8807]]],
 CS[Cartesian,2],
 AXIS["(E)",east,
 ORDER[1],
 LENGTHUNIT["metre",1]],
 AXIS["(N)",north,
 ORDER[2],
 LENGTHUNIT["metre",1]],
 USAGE[
 SCOPE["unknown"],
 AREA["Qatar - onshore"],
 BBOX[24.55,50.69,26.2,51.68]],
 ID["EPSG",2932]]



The problem is that GRASS still assumes datum transformation from X to

WGS84, whereas with PROJ6,

WGS84 is no longer required as pivot datum. Datum transformations from

datum X to datum Y can

sometimes (often) be done without going through WGS84. The requirement is

to have EPSG or other

authority codes and to use PROJ 6.

this warning is just from a simple v.in.ogr of a EPSG:2932-based shapefile
into a EPSG:2932-created GRASS location.

so, there may be some more adaptions needed for actions like v.in.ogr and
co?



-
best regards
Helmut
--
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


--

=
Zoltan Szecsei GPrGISc 0031
Geograph (Pty) Ltd.
GIS and Photogrammetric Services

P.O. Box 7, Muizenberg 7950, South Africa.

Mobile: +27-83-6004028 (WhatsApp only)
Qatar:  +974 5083 2722 www.geograph.co.za
=

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

Re: [GRASS-user] Datum not recognized by Grass

[Zoltan Szecsei]

Hi Markus,
Thanks for the interest in my problem.

Herewith g.version output:

Kind regards,
Zoltan

(Wed Dec  4 10:23:18 2019)
D:\sandbox\\tests\PERMANENT\.tmp/unknown\1608.0.py
version=7.8.1
date=2019
revision=c865432c9
build_date=2019-11-10
build_platform=x86_64-w64-mingw32
build_off_t_size=8
libgis_revision=0
libgis_date="?"
proj=6.2.1
gdal=3.0.2
geos=3.8.0
sqlite=3.29.0
(Wed Dec  4 10:23:18 2019) Command finished (0 sec)





On 2019/12/03 20:06, Markus Neteler wrote:

On Mon, Dec 2, 2019 at 9:49 AM Zoltan Szecsei  wrote:

Hi,
I'm using EPSG 2932 and QGIS etc all are OK with it.

I have installed everything using OSGeo4W64, so how come Grass does not
use the same projections database (as QGis etc)?

Please post the output of

g.version -rge

Best
Markus


--

=
Zoltan Szecsei GPrGISc 0031
Geograph (Pty) Ltd.
GIS and Photogrammetric Services

P.O. Box 7, Muizenberg 7950, South Africa.

Mobile: +27-83-6004028 (WhatsApp only)
Qatar:  +974 5083 2722 www.geograph.co.za
=

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

Re: [GRASS-user] Problem with v.class.mIR

[Moritz Lennert]

Sorry for ignoring you up to now, just to much on the platter right now.

Would it be feasible for you to make the data available, privately if 
necessary ? The easiest (and lightest) would be just the csv output of 
v.db.select on the training_map and the segments_maps.


Moritz

On 2/12/19 16:31, Giuseppe Cillis wrote:

No one can help me?
GC

Il giorno mer 20 nov 2019 alle ore 14:06 Giuseppe Cillis 
mailto:giucil...@gmail.com>> ha scritto:


Thanks for the answer._
_
_
_
_v.class.mlR segments_map=gbbb@PERMANENT
training_map=tainning@PERMANENT train_class_column=fd
output_class_column=vote output_prob_column=prob folds=5
partitions=10 tunelength=10 weighting_metric=accuracy_

_and the same error with this command:
_
_v.class.mlR segments_map=gbbb training_map=tainning
train_class_column=class weighting_mode=smv,swv,qbwwv -i_

The vector used was elaborated with i.segment.stats module.
An extract is here:
cat areaperimeter   fd
1   958 208 1.555034
2   24  22  1.945217
3   160 80  1.726846
4   2036242 1.440904
5   25  26  2.024344
6   222 96  1.68966
7   843510121.530879


/cat /type Integer (9-0)
/fd /type double/real (25-9)
For the training area (t/ainning/), I use a simple vector with
numerical categories (named/"cate"/) and there are 11 features (just
for a trying). The type is in integer64 (10)

Il giorno mer 20 nov 2019 alle ore 12:18 Moritz Lennert
mailto:mlenn...@club.worldonline.be>>
ha scritto:

On 20/11/19 11:55, Giuseppe Cillis wrote:
 > Hi,
 > I'm trying to use this module for classification of an old
aerial photos.
 > After a segmentation (i.segment and i.segment.stats), I would
to use a
 > machine learning approach for the real classification.
 > I tried with v.class.mIR which use also R.
 > But there is an error and I don't know how to solve it:
 > */Durante l'avvio - Warning messages:
 > 1: Setting LC_CTYPE=it_IT.cp1252 failed
 > 2: Setting LC_COLLATE=it_IT.cp1252 failed
 > 3: Setting LC_TIME=it_IT.cp1252 failed
 > 4: Setting LC_MONETARY=it_IT.cp1252 failed
 > Carico il pacchetto richiesto: caret
 > Carico il pacchetto richiesto: lattice
 > Carico il pacchetto richiesto: ggplot2
 > Error in `$<-.data.frame`(`*tmp*`, fd, value = integer(0)) :
 >    replacement has 0 rows, data has 11
 > Calls: $<- -> $<-.data.frame/*
 > */
 > /*
 > A part is in italian, I'm sorry.
 > I'm a beginner in grass and R.

R messages are often not very explicit, unfortunately.

Could you provide us with the exact command line you use (If you
use the
GUI, you can get the command line by clicking the 'Copy' button
once
you've filled out all the parameters, or in the history of the
'Console'
in the Layer manager.) and an extract of the attribute data you
feed
into the module ?

Moritz


___
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