[GRASS-user] i.histo.match error

GRASS GRASS user
1

Hello all,

Been trying to run i.histo.match but keep getting this error:

Traceback (most recent call last):
  File "i.histo.match.py", line 223, in <module>
    sys.exit(main())
  File "i.histo.match.py", line 89, in main
    curs.execute(query_create)
sqlite3.OperationalError: unknown database tnormalized

All help is appreciated...
(grass7-win7 64b)
luis

2

Luís Mário Ribeiro wrote:

Hello all,

Hi!

Been trying to run i.histo.match but keep getting this error:

Traceback (most recent call last):
  File "i.histo.match.py", line 223, in <module>
    sys.exit(main())
  File "i.histo.match.py", line 89, in main
    curs.execute(query_create)
sqlite3.OperationalError: unknown database tnormalized

..

(grass7-win7 64b)

Would you mind posting the exact command you are trying to execute?

Thanks, Nikos

3

Hi Nikos, all

Been using histo.match thru the GUI. The command that appears in the output
window is (in this example just 3 images):

i.histo.match
input=normalized.LT52030312007052MPS00_tr_reflec3_int@histomatch,normalized.
LT52030312007052MPS00_tr_reflec4_int@histomatch,normalized.LT52030312007180M
PS00_tr_reflec3_int@histomatch max=1000

best regards
luis

-----Mensagem original-----
De: grass-user-bounces@lists.osgeo.org
[mailto:grass-user-bounces@lists.osgeo.org] Em nome de Nikos Alexandris
Enviada: terça-feira, 5 de Março de 2013 17:09
Para: grass-user@lists.osgeo.org
Assunto: Re: [GRASS-user] i.histo.match error

Luís Mário Ribeiro wrote:

Hello all,

Hi!

Been trying to run i.histo.match but keep getting this error:

Traceback (most recent call last):
  File "i.histo.match.py", line 223, in <module>
    sys.exit(main())
  File "i.histo.match.py", line 89, in main
    curs.execute(query_create)
sqlite3.OperationalError: unknown database tnormalized

..

(grass7-win7 64b)

Would you mind posting the exact command you are trying to execute?

Thanks, Nikos
_______________________________________________
grass-user mailing list
grass-user@lists.osgeo.org
http://lists.osgeo.org/mailman/listinfo/grass-user

4

Luís Mário Ribeiro wrote:
..

Been using histo.match thru the GUI. The command that appears in the output
window is (in this example just 3 images):

i.histo.match
input=normalized.LT52030312007052MPS00_tr_reflec3_int@histomatch,normalized.
LT52030312007052MPS00_tr_reflec4_int@histomatch,normalized.LT52030312007180
M PS00_tr_reflec3_int@histomatch max=1000

..

I guess (not sure though) it is the same error once pointed to me by Luca
Delucchi: the "." in the input raster map names is not allowed :-(.

You need to rename them... Maybe a quick test first to be sure. Shrink your
region and test. Be sure you have beforehand a MASK or else which you use to
set the computational region. Or, export your current working region in a
vector map? (<-- v.in.region).

Hope it works, Nikos

5

[>]Good morning and thanks... Will try changing the "." (which was added by
Grass when performing i.topo.corr).
Best regards, luis

Luís Mário Ribeiro wrote:
..

Been using histo.match thru the GUI. The command that appears in the
output window is (in this example just 3 images):

i.histo.match

input=normalized.LT52030312007052MPS00_tr_reflec3_int@histomatch,normalized.

LT52030312007052MPS00_tr_reflec4_int@histomatch,normalized.LT520303120
07180 M PS00_tr_reflec3_int@histomatch max=1000

..

I guess (not sure though) it is the same error once pointed to me by Luca
Delucchi: the "." in the input raster map names is not allowed :-(.

You need to rename them... Maybe a quick test first to be sure. Shrink
your region and test. Be sure you have beforehand a MASK or else which you
use to set the computational region. Or, export your current working region
in a vector map? (<-- v.in.region).

Hope it works, Nikos

6

Hello list,

Changed the names of 2 files to not have any "." and rerun histo.match. A
different error showed.
I copy below the output window contents. It refers to an error in " r.stats"
(ERROR: Sorry, <fs> is not a valid parameter)
Does anyone know what I'm doing wrong?

_________
(Wed Mar 06 10:12:17 2013)

i.histo.match --overwrite
input=2030312007052_3@histomatch,2030322007052_3@histomatch max=1000
Description:
Generates area statistics for raster map layers.
Keywords:
raster, statistics
Usage:
r.stats [-1AacplgxrnNCi] input=name[,name,...]
[output=name]
   [separator=character] [nv=string] [nsteps=value]
[--overwrite]
   [--verbose] [--quiet]
Flags:
  -1 One cell (range) per line
  -A Print averaged values instead of intervals
  -a Print area totals
  -c Print cell counts
  -p Print APPROXIMATE percents (total percent may not be
100%)
  -l Print category labels
  -g Print grid coordinates (east and north)
  -x Print x and y (column and row)
  -r Print raw indexes of fp ranges (fp maps only)
  -n Suppress reporting of any NULLs
  -N Suppress reporting of NULLs when all values are NULL
  -C Report for cats fp ranges (fp maps only)
  -i Read fp map as integer (use map's quant rules)
--o Allow output files to overwrite existing files
--v Verbose module output
--q Quiet module output
Parameters:
      input Name of input raster map(s)
     output Name for output file (if omitted or "-" output
to stdout)
  separator Field separator
               Special characters: newline, space, comma,
tab
              default: space
         nv String representing no data cell value
              default: *
     nsteps Number of fp subranges to collect stats from
              default: 255
ERROR: Sorry, <fs> is not a valid parameter
Description:
Generates area statistics for raster map layers.
Keywords:
raster, statistics
Usage:
r.stats [-1AacplgxrnNCi] input=name[,name,...]
[output=name]
   [separator=character] [nv=string] [nsteps=value]
[--overwrite]
   [--verbose] [--quiet]
Flags:
  -1 One cell (range) per line
  -A Print averaged values instead of intervals
  -a Print area totals
  -c Print cell counts
  -p Print APPROXIMATE percents (total percent may not be
100%)
  -l Print category labels
  -g Print grid coordinates (east and north)
  -x Print x and y (column and row)
  -r Print raw indexes of fp ranges (fp maps only)
  -n Suppress reporting of any NULLs
  -N Suppress reporting of NULLs when all values are NULL
  -C Report for cats fp ranges (fp maps only)
  -i Read fp map as integer (use map's quant rules)
--o Allow output files to overwrite existing files
--v Verbose module output
--q Quiet module output
Parameters:
      input Name of input raster map(s)
     output Name for output file (if omitted or "-" output
to stdout)
  separator Field separator
               Special characters: newline, space, comma,
tab
              default: space
         nv String representing no data cell value
              default: *
     nsteps Number of fp subranges to collect stats from
              default: 255
ERROR: Sorry, <fs> is not a valid parameter
(Wed Mar 06 10:12:25 2013) Command finished (8 sec)

_____

-----Mensagem original-----
De: grass-user-bounces@lists.osgeo.org
[mailto:grass-user-bounces@lists.osgeo.org] Em nome de Luís Mário Ribeiro
Enviada: quarta-feira, 6 de Março de 2013 09:30
Para: 'Nikos Alexandris'
Cc: grass-user@lists.osgeo.org
Assunto: Re: [GRASS-user] i.histo.match error

[>]Good morning and thanks... Will try changing the "." (which was added by
Grass when performing i.topo.corr).
Best regards, luis

Luís Mário Ribeiro wrote:
..

Been using histo.match thru the GUI. The command that appears in the
output window is (in this example just 3 images):

i.histo.match

input=normalized.LT52030312007052MPS00_tr_reflec3_int@histomatch,normalized.

LT52030312007052MPS00_tr_reflec4_int@histomatch,normalized.LT520303120
07180 M PS00_tr_reflec3_int@histomatch max=1000

..

I guess (not sure though) it is the same error once pointed to me by Luca
Delucchi: the "." in the input raster map names is not allowed :-(.

You need to rename them... Maybe a quick test first to be sure. Shrink
your region and test. Be sure you have beforehand a MASK or else which you
use to set the computational region. Or, export your current working region
in a vector map? (<-- v.in.region).

Hope it works, Nikos

_______________________________________________
grass-user mailing list
grass-user@lists.osgeo.org
http://lists.osgeo.org/mailman/listinfo/grass-user

7

On Wed, Mar 6, 2013 at 11:16 AM, Luís Mário Ribeiro <luis.mario@adai.pt> wrote:

Hello list,

Changed the names of 2 files to not have any "." and rerun histo.match. A
different error showed.
I copy below the output window contents. It refers to an error in " r.stats"
(ERROR: Sorry, <fs> is not a valid parameter)

It seems that you need a newer version of i.histo.match, this was
fixed couple of weeks ago [1].

Regards,
Anna

[1] http://trac.osgeo.org/grass/changeset/54833/grass-addons/grass7/imagery/i.histo.match/i.histo.match.py

Does anyone know what I'm doing wrong?

_________
(Wed Mar 06 10:12:17 2013)

i.histo.match --overwrite
input=2030312007052_3@histomatch,2030322007052_3@histomatch max=1000
Description:
Generates area statistics for raster map layers.
Keywords:
raster, statistics
Usage:
r.stats [-1AacplgxrnNCi] input=name[,name,...]
[output=name]
   [separator=character] [nv=string] [nsteps=value]
[--overwrite]
   [--verbose] [--quiet]
Flags:
  -1 One cell (range) per line
  -A Print averaged values instead of intervals
  -a Print area totals
  -c Print cell counts
  -p Print APPROXIMATE percents (total percent may not be
100%)
  -l Print category labels
  -g Print grid coordinates (east and north)
  -x Print x and y (column and row)
  -r Print raw indexes of fp ranges (fp maps only)
  -n Suppress reporting of any NULLs
  -N Suppress reporting of NULLs when all values are NULL
  -C Report for cats fp ranges (fp maps only)
  -i Read fp map as integer (use map's quant rules)
--o Allow output files to overwrite existing files
--v Verbose module output
--q Quiet module output
Parameters:
      input Name of input raster map(s)
     output Name for output file (if omitted or "-" output
to stdout)
  separator Field separator
               Special characters: newline, space, comma,
tab
              default: space
         nv String representing no data cell value
              default: *
     nsteps Number of fp subranges to collect stats from
              default: 255
ERROR: Sorry, <fs> is not a valid parameter
Description:
Generates area statistics for raster map layers.
Keywords:
raster, statistics
Usage:
r.stats [-1AacplgxrnNCi] input=name[,name,...]
[output=name]
   [separator=character] [nv=string] [nsteps=value]
[--overwrite]
   [--verbose] [--quiet]
Flags:
  -1 One cell (range) per line
  -A Print averaged values instead of intervals
  -a Print area totals
  -c Print cell counts
  -p Print APPROXIMATE percents (total percent may not be
100%)
  -l Print category labels
  -g Print grid coordinates (east and north)
  -x Print x and y (column and row)
  -r Print raw indexes of fp ranges (fp maps only)
  -n Suppress reporting of any NULLs
  -N Suppress reporting of NULLs when all values are NULL
  -C Report for cats fp ranges (fp maps only)
  -i Read fp map as integer (use map's quant rules)
--o Allow output files to overwrite existing files
--v Verbose module output
--q Quiet module output
Parameters:
      input Name of input raster map(s)
     output Name for output file (if omitted or "-" output
to stdout)
  separator Field separator
               Special characters: newline, space, comma,
tab
              default: space
         nv String representing no data cell value
              default: *
     nsteps Number of fp subranges to collect stats from
              default: 255
ERROR: Sorry, <fs> is not a valid parameter
(Wed Mar 06 10:12:25 2013) Command finished (8 sec)

_____

-----Mensagem original-----
De: grass-user-bounces@lists.osgeo.org
[mailto:grass-user-bounces@lists.osgeo.org] Em nome de Luís Mário Ribeiro
Enviada: quarta-feira, 6 de Março de 2013 09:30
Para: 'Nikos Alexandris'
Cc: grass-user@lists.osgeo.org
Assunto: Re: [GRASS-user] i.histo.match error

[>]Good morning and thanks... Will try changing the "." (which was added by
Grass when performing i.topo.corr).
Best regards, luis

Luís Mário Ribeiro wrote:
..

Been using histo.match thru the GUI. The command that appears in the
output window is (in this example just 3 images):

i.histo.match

input=normalized.LT52030312007052MPS00_tr_reflec3_int@histomatch,normalized.

LT52030312007052MPS00_tr_reflec4_int@histomatch,normalized.LT520303120
07180 M PS00_tr_reflec3_int@histomatch max=1000

..

I guess (not sure though) it is the same error once pointed to me by Luca
Delucchi: the "." in the input raster map names is not allowed :-(.

You need to rename them... Maybe a quick test first to be sure. Shrink
your region and test. Be sure you have beforehand a MASK or else which you
use to set the computational region. Or, export your current working region
in a vector map? (<-- v.in.region).

Hope it works, Nikos

_______________________________________________
grass-user mailing list
grass-user@lists.osgeo.org
http://lists.osgeo.org/mailman/listinfo/grass-user

_______________________________________________
grass-user mailing list
grass-user@lists.osgeo.org
http://lists.osgeo.org/mailman/listinfo/grass-user

8

Thank you Ana and Nikos,

The version I have was instaled yesterday again and the error subsisted. I
found the log of the last change
(http://trac.osgeo.org/grass/changeset/54833) and manually edited my
"i.histo.match.py", changing "fs" for "separator".

All is working now but I have no idea how or why grass instaled an old
version :slight_smile:

Best regards
luis

-----Mensagem original-----
De: Anna Kratochvílová [mailto:kratochanna@gmail.com]
Enviada: quarta-feira, 6 de Março de 2013 19:57
Para: Luís Mário Ribeiro
Cc: Nikos Alexandris; GRASS user list
Assunto: Re: [GRASS-user] i.histo.match error

On Wed, Mar 6, 2013 at 11:16 AM, Luís Mário Ribeiro <luis.mario@adai.pt>
wrote:

Hello list,

Changed the names of 2 files to not have any "." and rerun
histo.match. A different error showed.
I copy below the output window contents. It refers to an error in "

r.stats"

(ERROR: Sorry, <fs> is not a valid parameter)

It seems that you need a newer version of i.histo.match, this was fixed
couple of weeks ago [1].

Regards,
Anna

[1]
http://trac.osgeo.org/grass/changeset/54833/grass-addons/grass7/imagery/i.hi
sto.match/i.histo.match.py

Does anyone know what I'm doing wrong?

_________
(Wed Mar 06 10:12:17 2013)

i.histo.match --overwrite
input=2030312007052_3@histomatch,2030322007052_3@histomatch max=1000
Description:
Generates area statistics for raster map layers.
Keywords:
raster, statistics
Usage:
r.stats [-1AacplgxrnNCi] input=name[,name,...] [output=name]
   [separator=character] [nv=string] [nsteps=value] [--overwrite]
   [--verbose] [--quiet]
Flags:
  -1 One cell (range) per line
  -A Print averaged values instead of intervals
  -a Print area totals
  -c Print cell counts
  -p Print APPROXIMATE percents (total percent may not be
100%)
  -l Print category labels
  -g Print grid coordinates (east and north)
  -x Print x and y (column and row)
  -r Print raw indexes of fp ranges (fp maps only)
  -n Suppress reporting of any NULLs
  -N Suppress reporting of NULLs when all values are NULL
  -C Report for cats fp ranges (fp maps only)
  -i Read fp map as integer (use map's quant rules)
--o Allow output files to overwrite existing files
--v Verbose module output
--q Quiet module output
Parameters:
      input Name of input raster map(s)
     output Name for output file (if omitted or "-" output
to stdout)
  separator Field separator
               Special characters: newline, space, comma, tab
              default: space
         nv String representing no data cell value
              default: *
     nsteps Number of fp subranges to collect stats from
              default: 255
ERROR: Sorry, <fs> is not a valid parameter
Description:
Generates area statistics for raster map layers.
Keywords:
raster, statistics
Usage:
r.stats [-1AacplgxrnNCi] input=name[,name,...] [output=name]
   [separator=character] [nv=string] [nsteps=value] [--overwrite]
   [--verbose] [--quiet]
Flags:
  -1 One cell (range) per line
  -A Print averaged values instead of intervals
  -a Print area totals
  -c Print cell counts
  -p Print APPROXIMATE percents (total percent may not be
100%)
  -l Print category labels
  -g Print grid coordinates (east and north)
  -x Print x and y (column and row)
  -r Print raw indexes of fp ranges (fp maps only)
  -n Suppress reporting of any NULLs
  -N Suppress reporting of NULLs when all values are NULL
  -C Report for cats fp ranges (fp maps only)
  -i Read fp map as integer (use map's quant rules)
--o Allow output files to overwrite existing files
--v Verbose module output
--q Quiet module output
Parameters:
      input Name of input raster map(s)
     output Name for output file (if omitted or "-" output
to stdout)
  separator Field separator
               Special characters: newline, space, comma, tab
              default: space
         nv String representing no data cell value
              default: *
     nsteps Number of fp subranges to collect stats from
              default: 255
ERROR: Sorry, <fs> is not a valid parameter (Wed Mar 06 10:12:25 2013)
Command finished (8 sec)

_____

-----Mensagem original-----
De: grass-user-bounces@lists.osgeo.org
[mailto:grass-user-bounces@lists.osgeo.org] Em nome de Luís Mário
Ribeiro
Enviada: quarta-feira, 6 de Março de 2013 09:30
Para: 'Nikos Alexandris'
Cc: grass-user@lists.osgeo.org
Assunto: Re: [GRASS-user] i.histo.match error

[>]Good morning and thanks... Will try changing the "." (which was
added by Grass when performing i.topo.corr).
Best regards, luis

Luís Mário Ribeiro wrote:
..

Been using histo.match thru the GUI. The command that appears in the
output window is (in this example just 3 images):

i.histo.match

input=normalized.LT52030312007052MPS00_tr_reflec3_int@histomatch,normalized.

LT52030312007052MPS00_tr_reflec4_int@histomatch,normalized.LT52030312
0
07180 M PS00_tr_reflec3_int@histomatch max=1000

..

I guess (not sure though) it is the same error once pointed to me by
Luca
Delucchi: the "." in the input raster map names is not allowed :-(.

You need to rename them... Maybe a quick test first to be sure.
Shrink your region and test. Be sure you have beforehand a MASK or
else which you use to set the computational region. Or, export your
current working region in a vector map? (<-- v.in.region).

Hope it works, Nikos

_______________________________________________
grass-user mailing list
grass-user@lists.osgeo.org
http://lists.osgeo.org/mailman/listinfo/grass-user

_______________________________________________
grass-user mailing list
grass-user@lists.osgeo.org
http://lists.osgeo.org/mailman/listinfo/grass-user