[GRASS-user] i.pansharpen error

GRASS GRASS user
1

Hi?

I am trying to perform pansharpening on landsat raster, using pygrass. Below is the snippet

···

Kind Regards,
Joseph Kariuki

Geospatial Engineer | GIS / Web Developer

2

ERROR: Raster map <LE07_L1TP_166060_20161016_20161111_01_T1_B8.TIF.aux> not
      found

just guessing, does

    LE07_L1TP_166060_20161016_20161111_01_T1_B8.TIF.aux

exist? is this the correct raster name?

-----
best regards
Helmut
--
View this message in context: http://osgeo-org.1560.x6.nabble.com/i-pansharpen-error-tp5321277p5321290.html
Sent from the Grass - Users mailing list archive at Nabble.com.

3

* Joseph Kariuki <joehene@gmail.com> [2017-05-22 13:35:44 +0300]:

Hi?

I am trying to perform pansharpening on landsat raster, using pygrass.
Below is the snippet

Not answering your question, but please give i.fusion.hpf a try.

Nikos

*n = 0for x in bluear:
mypansharpen.append(grass.run_command("i.pansharpen", red=redar[n],
green=greenar[n], blue=bluear[n], pan=panar[n],
                 output=bluear[n][0:21]+"pansharpen", method="ihs"))
n = n + 1*

Running the code returns the following error

*ERROR: Raster map <LE07_L1TP_166060_20161016_20161111_01_T1_B8.TIF.aux> not
      found
Traceback (most recent call last):
File "/usr/local/grass-7.2.2svn/scripts/i.pansharpen", line 452, in <module>
   main()
File "/usr/local/grass-7.2.2svn/scripts/i.pansharpen", line 119, in main
   kv = grass.raster_info(map=pan)
File "/usr/local/grass-7.2.2svn/etc/python/grass/script/raster.py",
line 78, in raster_info
   s = read_command('r.info <http://r.info>', flags='gre', map=map)
File "/usr/local/grass-7.2.2svn/etc/python/grass/script/core.py",
line 461, in read_command
   return handle_errors(returncode, stdout, args, kwargs)
File "/usr/local/grass-7.2.2svn/etc/python/grass/script/core.py",
line 329, in handle_errors
   returncode=returncode)
grass.exceptions.CalledModuleError: Module run None ['r.info
<http://r.info>', '-gre',
'map=LE07_L1TP_166060_20161016_20161111_01_T1_B8.TIF.aux'] ended with
error
Process ended with non-zero return code 1. See errors in the (error) output.
Traceback (most recent call last):
File "/home/hempire/PycharmProjects/GRASSproject/landsatProcessing.py",
line 185, in <module>
   output=bluear[n][0:21]+"pansharpen", method="ihs"))
File "/usr/local/grass-7.2.2svn/etc/python/grass/script/core.py",
line 410, in run_command
   return handle_errors(returncode, returncode, args, kwargs)
File "/usr/local/grass-7.2.2svn/etc/python/grass/script/core.py",
line 329, in handle_errors
   returncode=returncode)
grass.exceptions.CalledModuleError: Module run None ['i.pansharpen',
'blue=LE07_L1TP_166060_20030303_20170126_01_T1_B2',
'green=LE07_L1TP_166060_20030303_20170126_01_T1_B3',
'pan=LE07_L1TP_166060_20161016_20161111_01_T1_B8.TIF.aux',
'output=LE07_L1TP_166060_2003pansharpen', 'method=ihs',
'red=LE07_L1TP_166060_20030303_20170126_01_T1_B4'] ended with error
Process ended with non-zero return code 1. See errors in the (error) output.*

Note that I have also tried grass.read_command in place of
grass.run_command but still getting the same error output.

Kind Regards,

*Joseph Kariuki*

*Geospatial Engineer** | **GIS / Web Developer*

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

--
Nikos Alexandris | Remote Sensing & Geomatics
GPG Key Fingerprint 6F9D4506F3CA28380974D31A9053534B693C4FB3

4

[please keep it on the ML]

I just did a batch import of all raster bands including panchromatic band
therefor it exists.

ERROR: Raster map <LE07_L1TP_166060_20161016_20161111_01_T1_B8.TIF.aux> not
      found

it looks for => TIF.aux;

are you sure that this aux is a valid raster?

please post here the relevant part of g.list type=raster

and r.info-output for LE07_L1TP_166060_20161016_20161111_01_T1_B8.TIF.aux

-----
best regards
Helmut
--
View this message in context: http://osgeo-org.1560.x6.nabble.com/i-pansharpen-error-tp5321277p5321300.html
Sent from the Grass - Users mailing list archive at Nabble.com.

5

I shall give it a try. thanks

···

On Mon, May 22, 2017 at 4:17 PM, Nikos Alexandris <nik@nikosalexandris.net> wrote:

Hi?

I am trying to perform pansharpening on landsat raster, using pygrass.
Below is the snippet

Not answering your question, but please give i.fusion.hpf a try.

Nikos

n = 0for x in bluear:
mypansharpen.append(grass.run_command(“i.pansharpen”, red=redar[n],
green=greenar[n], blue=bluear[n], pan=panar[n],
output=bluear[n][0:21]+“pansharpen”, method=“ihs”))
n = n + 1

Running the code returns the following error

ERROR: Raster map <LE07_L1TP_166060_20161016_20161111_01_T1_B8.TIF.aux> not
found
Traceback (most recent call last):
File “/usr/local/grass-7.2.2svn/scripts/i.pansharpen”, line 452, in
main()
File “/usr/local/grass-7.2.2svn/scripts/i.pansharpen”, line 119, in main
kv = grass.raster_info(map=pan)
File “/usr/local/grass-7.2.2svn/etc/python/grass/script/raster.py”,
line 78, in raster_info
s = read_command(‘r.info <http://r.info>’, flags=‘gre’, map=map)
File “/usr/local/grass-7.2.2svn/etc/python/grass/script/core.py”,
line 461, in read_command
return handle_errors(returncode, stdout, args, kwargs)
File “/usr/local/grass-7.2.2svn/etc/python/grass/script/core.py”,
line 329, in handle_errors
returncode=returncode)
grass.exceptions.CalledModuleError: Module run None [‘r.info
<http://r.info>’, ‘-gre’,
‘map=LE07_L1TP_166060_20161016_20161111_01_T1_B8.TIF.aux’] ended with
error
Process ended with non-zero return code 1. See errors in the (error) output.
Traceback (most recent call last):
File “/home/hempire/PycharmProjects/GRASSproject/landsatProcessing.py”,
line 185, in
output=bluear[n][0:21]+“pansharpen”, method=“ihs”))
File “/usr/local/grass-7.2.2svn/etc/python/grass/script/core.py”,
line 410, in run_command
return handle_errors(returncode, returncode, args, kwargs)
File “/usr/local/grass-7.2.2svn/etc/python/grass/script/core.py”,
line 329, in handle_errors
returncode=returncode)
grass.exceptions.CalledModuleError: Module run None [‘i.pansharpen’,
‘blue=LE07_L1TP_166060_20030303_20170126_01_T1_B2’,
‘green=LE07_L1TP_166060_20030303_20170126_01_T1_B3’,
‘pan=LE07_L1TP_166060_20161016_20161111_01_T1_B8.TIF.aux’,
‘output=LE07_L1TP_166060_2003pansharpen’, ‘method=ihs’,
‘red=LE07_L1TP_166060_20030303_20170126_01_T1_B4’] ended with error
Process ended with non-zero return code 1. See errors in the (error) output.

Note that I have also tried grass.read_command in place of
grass.run_command but still getting the same error output.

Kind Regards,

Joseph Kariuki

Geospatial Engineer* | *GIS / Web Developer

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

–
Nikos Alexandris | Remote Sensing & Geomatics
GPG Key Fingerprint 6F9D4506F3CA28380974D31A9053534B693C4FB3

Kind Regards,
Joseph Kariuki

Geospatial Engineer | GIS / Web Developer

6

LE07_L1TP_166060_20030303_20170126_01_T1_B1
LE07_L1TP_166060_20030303_20170126_01_T1_B2
LE07_L1TP_166060_20030303_20170126_01_T1_B3
LE07_L1TP_166060_20030303_20170126_01_T1_B4
LE07_L1TP_166060_20030303_20170126_01_T1_B5
LE07_L1TP_166060_20030303_20170126_01_T1_B61
LE07_L1TP_166060_20030303_20170126_01_T1_B62
LE07_L1TP_166060_20030303_20170126_01_T1_B7
LE07_L1TP_166060_20030303_20170126_01_T1_B8
LE07_L1TP_166060_20030303_20170126_01_T1_BQA
LE07_L1TP_166060_2003composite
LE07_L1TP_166060_20161016_20161111_01_T1_B1
LE07_L1TP_166060_20161016_20161111_01_T1_B2
LE07_L1TP_166060_20161016_20161111_01_T1_B3
LE07_L1TP_166060_20161016_20161111_01_T1_B4
LE07_L1TP_166060_20161016_20161111_01_T1_B5
LE07_L1TP_166060_20161016_20161111_01_T1_B61
LE07_L1TP_166060_20161016_20161111_01_T1_B62
LE07_L1TP_166060_20161016_20161111_01_T1_B7
LE07_L1TP_166060_20161016_20161111_01_T1_B8
LE07_L1TP_166060_20161016_20161111_01_T1_BQA
LE07_L1TP_166060_2016composite

···

On Mon, May 22, 2017 at 4:30 PM, Helmut Kudrnovsky <hellik@web.de> wrote:

[please keep it on the ML]

I just did a batch import of all raster bands including panchromatic band
therefor it exists.

ERROR: Raster map <LE07_L1TP_166060_20161016_20161111_01_T1_B8.TIF.aux> not
found

it looks for => TIF.aux;

are you sure that this aux is a valid raster?

please post here the relevant part of g.list type=raster

and r.info-output for LE07_L1TP_166060_20161016_20161111_01_T1_B8.TIF.aux

best regards
Helmut

View this message in context: http://osgeo-org.1560.x6.nabble.com/i-pansharpen-error-tp5321277p5321300.html
Sent from the Grass - Users mailing list archive at Nabble.com.

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

Kind Regards,
Joseph Kariuki

Geospatial Engineer | GIS / Web Developer

7

Joseph Kariuki wrote

LE07_L1TP_166060_20030303_20170126_01_T1_B1
LE07_L1TP_166060_20030303_20170126_01_T1_B2
LE07_L1TP_166060_20030303_20170126_01_T1_B3
LE07_L1TP_166060_20030303_20170126_01_T1_B4
LE07_L1TP_166060_20030303_20170126_01_T1_B5
LE07_L1TP_166060_20030303_20170126_01_T1_B61
LE07_L1TP_166060_20030303_20170126_01_T1_B62
LE07_L1TP_166060_20030303_20170126_01_T1_B7
LE07_L1TP_166060_20030303_20170126_01_T1_B8
LE07_L1TP_166060_20030303_20170126_01_T1_BQA
LE07_L1TP_166060_2003composite
LE07_L1TP_166060_20161016_20161111_01_T1_B1
LE07_L1TP_166060_20161016_20161111_01_T1_B2
LE07_L1TP_166060_20161016_20161111_01_T1_B3
LE07_L1TP_166060_20161016_20161111_01_T1_B4
LE07_L1TP_166060_20161016_20161111_01_T1_B5
LE07_L1TP_166060_20161016_20161111_01_T1_B61
LE07_L1TP_166060_20161016_20161111_01_T1_B62
LE07_L1TP_166060_20161016_20161111_01_T1_B7
LE07_L1TP_166060_20161016_20161111_01_T1_B8
LE07_L1TP_166060_20161016_20161111_01_T1_BQA
LE07_L1TP_166060_2016composite

and

*ERROR: Raster map <LE07_L1TP_166060_20161016_20161111_01_T1_B8.TIF.aux>

not

      found

your script is looking for

LE07_L1TP_166060_20161016_20161111_01_T1_B8.TIF.aux

but it isn't in the g.list output.

-----
best regards
Helmut
--
View this message in context: http://osgeo-org.1560.x6.nabble.com/i-pansharpen-error-tp5321277p5321306.html
Sent from the Grass - Users mailing list archive at Nabble.com.

8

what does

r.info map=LE07_L1TP_166060_20161016_20161111_01_T1_B8.TIF.aux

?

-----
best regards
Helmut
--
View this message in context: http://osgeo-org.1560.x6.nabble.com/i-pansharpen-error-tp5321277p5321309.html
Sent from the Grass - Users mailing list archive at Nabble.com.

9

ERROR: Raster map <LE07_L1TP_166060_20161016_20161111_01_T1_B8.TIF.aux> not
found

If so how do I import a .aux file yet it is not among extensions supported by r.in.gdal command

···

On Mon, May 22, 2017 at 4:43 PM, Helmut Kudrnovsky <hellik@web.de> wrote:

what does

r.info map=LE07_L1TP_166060_20161016_20161111_01_T1_B8.TIF.aux

?

best regards
Helmut

View this message in context: http://osgeo-org.1560.x6.nabble.com/i-pansharpen-error-tp5321277p5321309.html

Sent from the Grass - Users mailing list archive at Nabble.com.

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

Kind Regards,
Joseph Kariuki

Geospatial Engineer | GIS / Web Developer

10

Joseph Kariuki wrote

ERROR: Raster map
<LE07_L1TP_166060_20161016_20161111_01_T1_B8.TIF.aux>
not
       found

If so how do I import a .aux file yet it is not among extensions supported
by r.in.gdal command

have a look into the geotiff format:

http://www.gdal.org/frmt_gtiff.html

aux.xml-files are:

"PAM (Persitant Auxiliary metadata) .aux.xml sidecar file"

there is no need to import .aux.xml-files as these files only cover
metadata, not raster data.

-----
best regards
Helmut
--
View this message in context: http://osgeo-org.1560.x6.nabble.com/i-pansharpen-error-tp5321277p5321313.html
Sent from the Grass - Users mailing list archive at Nabble.com.

11

Ok. I am trying to understand why i.pansharpen is looking for that .aux.xml file

···

On Mon, May 22, 2017 at 4:50 PM, Helmut Kudrnovsky <hellik@web.de> wrote:

Joseph Kariuki wrote

ERROR: Raster map
<LE07_L1TP_166060_20161016_20161111_01_T1_B8.TIF.aux>
not
found

If so how do I import a .aux file yet it is not among extensions supported
by r.in.gdal command

have a look into the geotiff format:

http://www.gdal.org/frmt_gtiff.html

aux.xml-files are:

“PAM (Persitant Auxiliary metadata) .aux.xml sidecar file”

there is no need to import .aux.xml-files as these files only cover
metadata, not raster data.

best regards
Helmut

View this message in context: http://osgeo-org.1560.x6.nabble.com/i-pansharpen-error-tp5321277p5321313.html

Sent from the Grass - Users mailing list archive at Nabble.com.

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

Kind Regards,
Joseph Kariuki

Geospatial Engineer | GIS / Web Developer

12

On 22/05/17 16:00, Joseph Kariuki wrote:

Ok. I am trying to understand why i.pansharpen is looking for that
.aux.xml file

It's not i.pansharpen, it's your script:

mypansharpen.append(grass.run_command("i.pansharpen", red=redar[n],
  green=greenar[n], blue=bluear[n], pan=panar[n],
  output=bluear[n][0:21]+"pansharpen", method="ihs"))

It seems to me that the aux file must have gotten into your panar list...

Moritz