[GRASS-user] Landsat 8 Workflow for Grass 7.0.0RC2
Paul Shapley
p.shapley at gmail.com
Wed Feb 18 04:44:00 PST 2015
Hi Grass Users,
I want to use the following Landsat 8 workflow to obtain some pansharpened
images but have a problem when i convert the DN values to 8bit for
pansharpening. Here is my preferred workflow in Grass 7.0.0RC2.
· r.in.gdal (success)
· i.landsat.toar (success)
· r.rescale bands 432 & 8 (convert reflectance 'DN'
values to 0-255, 8bit)
· i.pansharpen bands 432 with 8
· r.rescale (convert 0-255 values back to reflectance
'DN)
· i.colors.enhance (pansharp bands)
· d.rgb
My question...is this a valid workflow? because i'm getting the same two
errors below and the images when displayed look solid grey when viewed with
d.rgb:-
*displaying the pansharpened images error log:-*
(Tue Feb 17 14:57:26 2015)
g.list rast
LC82040242014206LGN00_B1
LC82040242014206LGN00_B10
LC82040242014206LGN00_B11
LC82040242014206LGN00_B2
LC82040242014206LGN00_B3
LC82040242014206LGN00_B4
LC82040242014206LGN00_B5
LC82040242014206LGN00_B6
LC82040242014206LGN00_B7
LC82040242014206LGN00_B8
LC82040242014206LGN00_B9
LC82040242014206LGN00_BQA
LC82040242014206LGN00_refl1
LC82040242014206LGN00_refl10
LC82040242014206LGN00_refl11
LC82040242014206LGN00_refl2
LC82040242014206LGN00_refl3
LC82040242014206LGN00_refl4
LC82040242014206LGN00_refl5
LC82040242014206LGN00_refl6
LC82040242014206LGN00_refl7
LC82040242014206LGN00_refl8
LC82040242014206LGN00_refl9
LC82040242014318LGN00_B1
LC82040242014318LGN00_B10
LC82040242014318LGN00_B11
LC82040242014318LGN00_B2
LC82040242014318LGN00_B3
LC82040242014318LGN00_B4
LC82040242014318LGN00_B5
LC82040242014318LGN00_B6
LC82040242014318LGN00_B7
LC82040242014318LGN00_B8
LC82040242014318LGN00_B9
LC82040242014318LGN00_BQA
LC82040242014318LGN00_refl1
LC82040242014318LGN00_refl10
LC82040242014318LGN00_refl11
LC82040242014318LGN00_refl2
LC82040242014318LGN00_refl2_rescale
LC82040242014318LGN00_refl3
LC82040242014318LGN00_refl3_rescale
LC82040242014318LGN00_refl4
LC82040242014318LGN00_refl4_rescale
LC82040242014318LGN00_refl5
LC82040242014318LGN00_refl6
LC82040242014318LGN00_refl7
LC82040242014318LGN00_refl8
LC82040242014318LGN00_refl8_rescale
LC82040242014318LGN00_refl9
(Tue Feb 17 14:57:27 2015) Command finished (0 sec)
Command 'd.rast
map=LC82040242014318LGN00_pansharp at PERMANENT' failed
Details: Raster map
<LC82040242014318LGN00_pansharp at PERMANENT> not found
*i.colors.enhance error log:-*
(Tue Feb 17 13:06:40 2015)
i.colors.enhance red=LC82040242014318LGN00_B4 at PERMANENT
green=LC82040242014318LGN00_B3 at PERMANENT
blue=LC82040242014318LGN00_refl2 at PERMANENT strength=96
Processing...
Process Process-3:
Traceback (most recent call last):
File
"C:\OSGeo4W\apps\Python27\lib\multiprocessing\process.py",
line 258, in _bootstrap
self.run()
File
"C:\OSGeo4W\apps\Python27\lib\multiprocessing\process.py",
line 114, in run
self._target(*self._args, **self._kwargs)
File "C:\OSGeo4W\apps\grass\grass-7.0.0RC2\scripts\i.color
s.enhance.py", line 103, in get_percentile_mp
result = get_percentile(map, percentiles)
File "C:\OSGeo4W\apps\grass\grass-7.0.0RC2\scripts\i.color
s.enhance.py", line 89, in get_percentile
percentiles = values, quiet = True)
File "C:\OSGeo4W\apps\grass\grass-7.0.0RC2\etc\python\gras
s\script\core.py", line 427, in read_command
return handle_errors(returncode, stdout, args, kwargs)
File "C:\OSGeo4W\apps\grass\grass-7.0.0RC2\etc\python\gras
s\script\core.py", line 310, in handle_errors
returncode=returncode)
CalledModuleError: Module run None ['r.quantile', '--q',
'input=LC82040242014318LGN00_refl2 at PERMANENT',
'percentiles=2,96'] ended with error
Process ended with non-zero return code -1073741571. See
errors in the (error) output.
Process Process-2:
Traceback (most recent call last):
File
"C:\OSGeo4W\apps\Python27\lib\multiprocessing\process.py",
line 258, in _bootstrap
self.run()
File
"C:\OSGeo4W\apps\Python27\lib\multiprocessing\process.py",
line 114, in run
self._target(*self._args, **self._kwargs)
File "C:\OSGeo4W\apps\grass\grass-7.0.0RC2\scripts\i.color
s.enhance.py", line 103, in get_percentile_mp
result = get_percentile(map, percentiles)
File "C:\OSGeo4W\apps\grass\grass-7.0.0RC2\scripts\i.color
s.enhance.py", line 89, in get_percentile
percentiles = values, quiet = True)
File "C:\OSGeo4W\apps\grass\grass-7.0.0RC2\etc\python\gras
s\script\core.py", line 427, in read_command
return handle_errors(returncode, stdout, args, kwargs)
File "C:\OSGeo4W\apps\grass\grass-7.0.0RC2\etc\python\gras
s\script\core.py", line 310, in handle_errors
returncode=returncode)
CalledModuleError: Module run None ['r.quantile', '--q',
'input=LC82040242014318LGN00_B3 at PERMANENT',
'percentiles=2,96'] ended with error
Process ended with non-zero return code -1073741571. See
errors in the (error) output.
Process Process-1:
Traceback (most recent call last):
File
"C:\OSGeo4W\apps\Python27\lib\multiprocessing\process.py",
line 258, in _bootstrap
self.run()
File
"C:\OSGeo4W\apps\Python27\lib\multiprocessing\process.py",
line 114, in run
self._target(*self._args, **self._kwargs)
File "C:\OSGeo4W\apps\grass\grass-7.0.0RC2\scripts\i.color
s.enhance.py", line 103, in get_percentile_mp
result = get_percentile(map, percentiles)
File "C:\OSGeo4W\apps\grass\grass-7.0.0RC2\scripts\i.color
s.enhance.py", line 89, in get_percentile
percentiles = values, quiet = True)
File "C:\OSGeo4W\apps\grass\grass-7.0.0RC2\etc\python\gras
s\script\core.py", line 427, in read_command
return handle_errors(returncode, stdout, args, kwargs)
File "C:\OSGeo4W\apps\grass\grass-7.0.0RC2\etc\python\gras
s\script\core.py", line 310, in handle_errors
returncode=returncode)
CalledModuleError: Module run None ['r.quantile', '--q',
'input=LC82040242014318LGN00_B4 at PERMANENT',
'percentiles=2,96'] ended with error
Process ended with non-zero return code -1073741571. See
errors in the (error) output.
--
*Paul J. Shapley *MSc CGeog (GIS) FRGS
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://lists.osgeo.org/pipermail/grass-user/attachments/20150218/dc8abebd/attachment-0001.html>
More information about the grass-user
mailing list