Re: [GRASS-dev] [GRASS-SVN] r73769 - grass/trunk/locale/po

[Markus Neteler]

Hi Huidae,

On Wed, Dec 5, 2018 at 1:39 PM  wrote:
>
> Author: hcho
> Date: 2018-12-05 04:39:10 -0800 (Wed, 05 Dec 2018)
> New Revision: 73769
>
> Modified:
>grass/trunk/locale/po/grasslibs_ko.po
>grass/trunk/locale/po/grassmods_ko.po
>grass/trunk/locale/po/grasswxpy_ko.po
> Log:
> i18n: Update Korean POs
>

just a question: did you update this using the transifex_merge script?

Otherwise the work may get lost until we have "tx push" working.

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

[GRASS-dev] Option 'band' in r.external

[Nikos Alexandris]

Dear list,

given a VRT built after
```
gdalbuildvrt -separate landsat8.vrt 
LC08_L1TP_184033_20181028_20181028_01_RT_B1[10].TIF
```

which then is structured like
```
gdalinfo landsat8.vrt |grep -e 'Band\|TIF'
  LC08_L1TP_184033_20181028_20181028_01_RT_B10.TIF
  LC08_L1TP_184033_20181028_20181028_01_RT_B11.TIF
Band 1 Block=128x128 Type=UInt16, ColorInterp=Undefined
Band 2 Block=128x128 Type=UInt16, ColorInterp=Undefined
```

the following works:
```
r.external input=landsat8.vrt -e out=B --o
WARNING: Unknown compression method , using default ZLIB
WARNING: Imagery group  already exists and will be overwritten.
Reading band 1 of 2...
r.external complete. Link to raster map  created.
Reading band 2 of 2...
r.external complete. Link to raster map  created.
Default region for this location updated
Region for the current mapset updated
Imagery group  created
```

Both bands are imported.  The same does `r.in.gdal` with or without the
'band' option, i.e. it imports both bands 1 and 2:
```
r.in.gdal input=landsat8.vrt -e band=1,2 out=B --o
WARNING: Unknown compression method , using default ZLIB
Importing 2 raster bands...
WARNING: Raster map  already exists and will be overwritten
Importing raster map ...
100%
WARNING: Raster map  already exists and will be overwritten
Importing raster map ...
100%
Default region for this location updated
Region for the current mapset updated
```

Why does the 'band' option to `r.external` not behave as expected? See:
```
r.external input=landsat8.vrt -e band=1,2 out=B --o
WARNING: Unknown compression method , using default ZLIB
Reading band 1 of 2...
r.external complete. Link to raster map  created.
Default region for this location updated
Region for the current mapset updated
```

It links only to band 1 of 2. Is it a bug?

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

[GRASS-dev] [GRASS GIS] #3702: r.in.wms: examples return errors

[GRASS GIS]

#3702: r.in.wms: examples return errors
-+-
 Reporter:  lucadelu |  Owner:  grass-dev@…
 Type:  defect   | Status:  new
 Priority:  normal   |  Milestone:  7.6.0
Component:  Raster   |Version:  unspecified
 Keywords:  r.in.wms |CPU:  Unspecified
 Platform:  Unspecified  |
-+-
 Hi, I'm testing r.in.wms examples from the manual page but most of them
 return errors


 {{{
 r.in.wms url="https://tiles.maps.eox.at/wms?; layers=s2cloudless
 output=sentinel2 format=png
 ERROR 1: KEA was compiled against GDAL 2.2, but the current library
 version is 2.3
 Downloading data from WMS server...
 Reprojecting raster...
 ERROR 1: KEA was compiled against GDAL 2.2, but the current library
 version is 2.3
 ERROR 1: Cannot guess driver for
 /home/lucadelu/grassdata/nc_spm_08/test/.tmp/geowork/14177.0
 ERROR: gdalwarp failed
 }}}


 {{{
 r.in.wms url="http://ows.mundialis.de/services/service?; layer=TOPO-WMS
 output=topo format=png
 ERROR 1: KEA was compiled against GDAL 2.2, but the current library
 version is 2.3
 Downloading data from WMS server...
 RuntimeError: module compiled against API version 0xc but this version of
 numpy is 0xa
 Traceback (most recent call last):
   File "/home/lucadelu/compilati/grass_trunk/dist.x86_64-pc-linux-
 gnu/scripts/r.in.wms", line 241, in 
 sys.exit(main())
   File "/home/lucadelu/compilati/grass_trunk/dist.x86_64-pc-linux-
 gnu/scripts/r.in.wms", line 227, in main
 fetched_map = wms.GetMap(options, flags)
   File "/home/lucadelu/compilati/grass_trunk/dist.x86_64-pc-linux-
 gnu/etc/r.in.wms/wms_base.py", line 215, in GetMap
 self.temp_map = self._download()
   File "/home/lucadelu/compilati/grass_trunk/dist.x86_64-pc-linux-
 gnu/etc/r.in.wms/wms_drv.py", line 194, in _download
 tile_dataset = self._pct2rgb(temp_tile, temp_tile_pct2rgb)
   File "/home/lucadelu/compilati/grass_trunk/dist.x86_64-pc-linux-
 gnu/etc/r.in.wms/wms_drv.py", line 285, in _pct2rgb
 src_data = src_band.ReadAsArray(0, iY, src_ds.RasterXSize, 1)
   File "/usr/lib/python2.7/dist-packages/osgeo/gdal.py", line 2598, in
 ReadAsArray
 import gdalnumeric
   File "/usr/lib/python2.7/dist-packages/osgeo/gdalnumeric.py", line 1, in
 
 from osgeo.gdal_array import *
   File "/usr/lib/python2.7/dist-packages/osgeo/gdal_array.py", line 17, in
 
 _gdal_array = swig_import_helper()
   File "/usr/lib/python2.7/dist-packages/osgeo/gdal_array.py", line 16, in
 swig_import_helper
 return importlib.import_module('_gdal_array')
   File "/usr/lib/python2.7/importlib/__init__.py", line 37, in
 import_module
 __import__(name)
 ImportError: No module named _gdal_array
 }}}


 {{{
 r.in.wms url="http://watzmann-geog.urz.uni-heidelberg.de/cached/osm;
 layers=osm_auto:all output=osm format=png
 ERROR 1: KEA was compiled against GDAL 2.2, but the current library
 version is 2.3
 Downloading data from WMS server...
 ERROR: Unable to fetch capabilities from
. Reason:
Proxy Error
 }}}

-- 
Ticket URL: 
GRASS GIS 

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

Re: [GRASS-dev] t.rast.univar Unable to get statistics for raster map

[Veronica Andreo]

Hi,

El mié., 5 dic. 2018 a las 11:04, Moritz Lennert (<
mlenn...@club.worldonline.be>) escribió:

> On 5/12/18 10:47, Margherita Di Leo wrote:
> > Hi,
> >
> > I have a strds and I try to get statistics for an area within a mask. A
> > couple of images have no data within the masked area, and t.rast.univar
> > gives the warning message:  "Unable to get statistics for raster map".
> > In the resulting csv, the missed dates are skipped, and this is sort of
> > unexpected and creates hurdles when for example I compare it with the
> > behavior of t.rast.what, that keeps the date and associates a no-data to
> > it. Wouldn't it be better if also t.rast.univar would yield a no-data
> > rather than skip the date?
>

+1

>
> This can probably somehow be solved within t.rast.univar, but IIUC the
> actual issue comes from the underlying r.univar call which has the same
> behaviour, for example when call with the 'zones' parameter and the '-t'
> flag: data is just absent if there are only null values in a given zone.
> I've been struggling with that, for example for addon modules such as
> i.segment.stats. I don't know how easy it would be to change that within
> r.univar, but it would be nice.
>
> IIUC, the issue comes from stats.c [1]:
>
> 123 if (stats[z].n == 0)
> 124 continue;
>
> Maybe a flag to fill the stats with NULL values, instead of ignoring
> them, would be appropriate ?


+1!!

But I don't know what an good NULL value would be here.
>

IIRC, r.out.xyz and t.rast.out.xyz use the '*' character.
It would be good to create a ticket for this, no? As to not forget
afterwards...

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

Re: [GRASS-dev] [GRASS GIS] #396: add NULL and MASK support to i.smap

[GRASS GIS]

#396: add NULL and MASK support to i.smap
--+-
  Reporter:  dylan|  Owner:  grass-dev@…
  Type:  enhancement  | Status:  new
  Priority:  minor|  Milestone:  7.8.0
 Component:  Raster   |Version:
Resolution:   |   Keywords:  i.smap
   CPU:  Other|   Platform:  All
--+-
Changes (by marisn):

 * milestone:  6.4.6 => 7.8.0


Comment:

 This still is a case.

 To make things worse – if NULL cells are inside a block, result differs
 significantly from a different sized block without NULL cells. Thus maps
 with NULL values should use as small blocks as possible.

-- 
Ticket URL: 
GRASS GIS 

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

Re: [GRASS-dev] t.rast.univar Unable to get statistics for raster map

[Moritz Lennert]

On 5/12/18 10:47, Margherita Di Leo wrote:

Hi,

I have a strds and I try to get statistics for an area within a mask. A 
couple of images have no data within the masked area, and t.rast.univar 
gives the warning message:  "Unable to get statistics for raster map". 
In the resulting csv, the missed dates are skipped, and this is sort of 
unexpected and creates hurdles when for example I compare it with the 
behavior of t.rast.what, that keeps the date and associates a no-data to 
it. Wouldn't it be better if also t.rast.univar would yield a no-data 
rather than skip the date?


This can probably somehow be solved within t.rast.univar, but IIUC the 
actual issue comes from the underlying r.univar call which has the same 
behaviour, for example when call with the 'zones' parameter and the '-t' 
flag: data is just absent if there are only null values in a given zone. 
I've been struggling with that, for example for addon modules such as 
i.segment.stats. I don't know how easy it would be to change that within 
r.univar, but it would be nice.


IIUC, the issue comes from stats.c [1]:

123 if (stats[z].n == 0)
124 continue;

Maybe a flag to fill the stats with NULL values, instead of ignoring 
them, would be appropriate ? But I don't know what an good NULL value 
would be here.


Moritz

[1] 
https://trac.osgeo.org/grass/browser/grass/trunk/raster/r.univar/stats.c#L123

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

Re: [GRASS-dev] [GRASS GIS] #3701: i.sentinel.prep fails with southern hemisphere Sentinel 2 level 1 data

[GRASS GIS]

#3701: i.sentinel.prep fails with southern hemisphere Sentinel 2 level 1 data
---+-
  Reporter:  geomax22  |  Owner:  grass-dev@…
  Type:  defect| Status:  new
  Priority:  normal|  Milestone:
 Component:  Addons|Version:  unspecified
Resolution:|   Keywords:
   CPU:  x86-64|   Platform:  MSWindows
---+-

Comment (by mlennert):

 Replying to [ticket:3701 geomax22]:
 > Using Sentinel S2A_MSIL1C from USGS in an area in Western Australia (UTM
 Zone 50 South ) fails when trying to create an input file for
 i.sentinel.mask with the error:
 >
 "C:\Users\Max\AppData\Roaming\GRASS7\addons\scripts\i.sentinel.preproc.py
 -a -t
 
input_dir=E:\S2_Abrholos\Sentinel\S2A_MSIL1C_20181128T021341_N0207_R060_T50JKN_20181128T080405.SAFE
 elevation=bathy_LL_sub_utm50@S2 atmospheric_model=Automatic
 aerosol_model=Continental model aod_value=0.07 suffix=depth
 text_file=mask_prep
 > ERROR: Latitude -29.3983025 is out of range
 > (Wed Dec 05 13:13:03 2018) Command finished (7 sec)
 > All works fine with imagery in the northern hemisphere.

 The '>' and '<' were inverted for negative numbers. Should be corrected in
 r73762. That should be available for Windows by tomorrow the latest.

-- 
Ticket URL: 
GRASS GIS 

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

[GRASS-dev] t.rast.univar Unable to get statistics for raster map

[Margherita Di Leo]

Hi,

I have a strds and I try to get statistics for an area within a mask. A
couple of images have no data within the masked area, and t.rast.univar
gives the warning message:  "Unable to get statistics for raster map". In
the resulting csv, the missed dates are skipped, and this is sort of
unexpected and creates hurdles when for example I compare it with the
behavior of t.rast.what, that keeps the date and associates a no-data to
it. Wouldn't it be better if also t.rast.univar would yield a no-data
rather than skip the date?

Thanks
Regards,

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