<div dir="ltr">Please unsubscribe me. Thank you!</div><div class="gmail_extra"><br><div class="gmail_quote">On Tue, Jan 20, 2015 at 2:00 PM, <span dir="ltr"><<a href="mailto:grass-user-request@lists.osgeo.org" target="_blank">grass-user-request@lists.osgeo.org</a>></span> wrote:<br><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">Send grass-user mailing list submissions to<br>
<a href="mailto:grass-user@lists.osgeo.org">grass-user@lists.osgeo.org</a><br>
<br>
To subscribe or unsubscribe via the World Wide Web, visit<br>
<a href="http://lists.osgeo.org/mailman/listinfo/grass-user" target="_blank">http://lists.osgeo.org/mailman/listinfo/grass-user</a><br>
or, via email, send a message with subject or body 'help' to<br>
<a href="mailto:grass-user-request@lists.osgeo.org">grass-user-request@lists.osgeo.org</a><br>
<br>
You can reach the person managing the list at<br>
<a href="mailto:grass-user-owner@lists.osgeo.org">grass-user-owner@lists.osgeo.org</a><br>
<br>
When replying, please edit your Subject line so it is more specific<br>
than "Re: Contents of grass-user digest..."<br>
<br>
<br>
Today's Topics:<br>
<br>
1. GRASS/gdal grib2 decoding issue for new NCEP MRMS data<br>
(Thomas Adams)<br>
2. Re: Large Grass Raster Notes (Markus Neteler)<br>
3. Re: Large Grass Raster Notes (Vaclav Petras)<br>
<br>
<br>
----------------------------------------------------------------------<br>
<br>
Message: 1<br>
Date: Mon, 19 Jan 2015 21:57:56 -0700<br>
From: Thomas Adams <<a href="mailto:tea3rd@gmail.com">tea3rd@gmail.com</a>><br>
To: "<a href="mailto:grass-user@lists.osgeo.org">grass-user@lists.osgeo.org</a>" <<a href="mailto:grass-user@lists.osgeo.org">grass-user@lists.osgeo.org</a>><br>
Subject: [GRASS-user] GRASS/gdal grib2 decoding issue for new NCEP<br>
MRMS data<br>
Message-ID:<br>
<<a href="mailto:CAGxgkWjUL%2BLMHcV9EB013KGUAyhLnvcT2b-xMNGiSDw-GAHVFw@mail.gmail.com">CAGxgkWjUL+LMHcV9EB013KGUAyhLnvcT2b-xMNGiSDw-GAHVFw@mail.gmail.com</a>><br>
Content-Type: text/plain; charset="utf-8"<br>
<br>
All:<br>
<br>
I'm not on the <a href="mailto:gdal-dev@lists.osgeo.org">gdal-dev@lists.osgeo.org</a> mailing list, so I can not report<br>
this here. NOAA/NCEP is producing new hourly CONUS NEXRAD radar-based<br>
meteorological data -- my interest is precipitation. The data is available<br>
now via a UCAR LDM process, described in the 1st link. The problem is<br>
(using the downloadable sample data), gdal recognizes the data files as<br>
grib2, but it will not decode the data (I have not previously had any<br>
problems). The Google docs file states that there is an issue and points to<br>
an updated wgrib2 that can handle the grib2 decoding (which it does) -- the<br>
second link.<br>
<br>
wgrib2 v0.2.0.1 7/2014 Wesley Ebisuzaki, Reinoud Bokhorst, Jaakko Hyv?tti,<br>
Dusan Jovic, Kristian Nilssen, Karl Pfeiffer, Pablo Romero, Manfred<br>
Schwarb, Arlindo da Silva, Niklas Sondell, Sergey Varlamov<br>
<br>
<a href="https://docs.google.com/document/d/1LeVcn_taIXZgzZb5JgWqaVr0xVs7GmA6RpHcb8ZGiwk/edit" target="_blank">https://docs.google.com/document/d/1LeVcn_taIXZgzZb5JgWqaVr0xVs7GmA6RpHcb8ZGiwk/edit</a><br>
<br>
<a href="http://www.cpc.noaa.gov/products/wesley/wgrib2/" target="_blank">http://www.cpc.noaa.gov/products/wesley/wgrib2/</a><br>
<br>
When I used wgrib2 to produce a binary file, I could then import the<br>
example data file into GRASS:<br>
<br>
<a href="http://r.info" target="_blank">r.info</a> map=MRMS_GaugeCorr_QPE_01H_00.00_20141208_120000@teaiii<br>
<br>
+----------------------------------------------------------------------------+<br>
| Map: MRMS_GaugeCorr_QPE_01H_00.00_ Date: Mon Jan 19 20:36:05<br>
2015 |<br>
| Mapset: teaiii Login of Creator:<br>
teaiii |<br>
| Location:<br>
world_location |<br>
| DataBase:<br>
/home/teaiii/grass/data |<br>
| Title: ( MRMS_GaugeCorr_QPE_01H_00.00_20141208-120000<br>
) |<br>
| Timestamp:<br>
none |<br>
|----------------------------------------------------------------------------|<br>
|<br>
|<br>
| Type of Map: raster Number of Categories:<br>
0 |<br>
| Data Type:<br>
FCELL |<br>
| Rows:<br>
3500 |<br>
| Columns:<br>
7000 |<br>
| Total Cells:<br>
24500000 |<br>
| Projection:<br>
Latitude-Longitude |<br>
| N: 55N S: 20N Res:<br>
0:00:36 |<br>
| E: 60W W: 130W Res:<br>
0:00:36 |<br>
| Range of data: min = 0 max =<br>
210.5 |<br>
|<br>
|<br>
| Data<br>
Description: |<br>
| generated by<br>
r.in.bin |<br>
|<br>
|<br>
|<br>
Comments: |<br>
| r.in.bin -f<br>
input="/home/teaiii/Downloads/MRMS_GRIB2_SAMPLE-precip/t\ |<br>
| est.bin" output="MRMS_GaugeCorr_QPE_01H_00.00_20141208-120000"<br>
bytes\ |<br>
| =4 header=8 bands=1 order="native" north=55 south=20 east=-60<br>
west=-\ |<br>
| 130 rows=3500 cols=7000 anull=-3<br>
flip="v" |<br>
|<br>
|<br>
+----------------------------------------------------------------------------+<br>
<br>
So, I thought the community would be interested in (1) knowing about the<br>
data and (2) the GRASS/gdal decoding issue.<br>
<br>
Cheers!<br>
Tom<br>
-------------- next part --------------<br>
An HTML attachment was scrubbed...<br>
URL: <<a href="http://lists.osgeo.org/pipermail/grass-user/attachments/20150119/2ce8984c/attachment-0001.html" target="_blank">http://lists.osgeo.org/pipermail/grass-user/attachments/20150119/2ce8984c/attachment-0001.html</a>><br>
<br>
------------------------------<br>
<br>
Message: 2<br>
Date: Tue, 20 Jan 2015 17:35:48 +0100<br>
From: Markus Neteler <<a href="mailto:neteler@osgeo.org">neteler@osgeo.org</a>><br>
To: Jeshua Lacock <<a href="mailto:jeshua@3dtopo.com">jeshua@3dtopo.com</a>><br>
Cc: GRASS user list <<a href="mailto:grass-user@lists.osgeo.org">grass-user@lists.osgeo.org</a>><br>
Subject: Re: [GRASS-user] Large Grass Raster Notes<br>
Message-ID:<br>
<CALFmHhuKm1=<a href="mailto:xmcvDkOVa2vGZnwVC_ZAx7PLh7GxSqebBNVX8cA@mail.gmail.com">xmcvDkOVa2vGZnwVC_ZAx7PLh7GxSqebBNVX8cA@mail.gmail.com</a>><br>
Content-Type: text/plain; charset=UTF-8<br>
<br>
On Sun, Jan 18, 2015 at 1:00 AM, Jeshua Lacock <<a href="mailto:jeshua@3dtopo.com">jeshua@3dtopo.com</a>> wrote:<br>
><br>
> Greetings,<br>
><br>
> In the event my experience working with rather large GRASS rasters may be useful, I thought I would share it. The following notes were compiled running the GRASS 7.0.0beta4-141230 package on Mac OS X 10.10.1. The system is a 6-core Xeon with 32GB of RAM, running in full 64-bit mode.<br>
><br>
> 1. Using an r.external virtual mosaic as the input proved impossible trying to run r.resample on it on trying to produce a nearly teracell raster (I ended up breaking sub-tiles). I closely followed the instructions on the Wiki.<br>
><br>
> 2. On smaller sub-tiles of said image, r.resample was insanely slow even with the external imagery.<br>
<br>
You may check/benchmark the other resampling modules, some of them may<br>
be more modern (ideally faster):<br>
<br>
<a href="http://grass.osgeo.org/grass70/manuals/keywords.html#resample" target="_blank">http://grass.osgeo.org/grass70/manuals/keywords.html#resample</a><br>
<br>
In case you try, please share the outcome.<br>
<br>
cheers,<br>
Markus<br>
<br>
<br>
------------------------------<br>
<br>
Message: 3<br>
Date: Tue, 20 Jan 2015 14:04:05 -0500<br>
From: Vaclav Petras <<a href="mailto:wenzeslaus@gmail.com">wenzeslaus@gmail.com</a>><br>
To: Jeshua Lacock <<a href="mailto:jeshua@3dtopo.com">jeshua@3dtopo.com</a>><br>
Cc: GRASS user list <<a href="mailto:grass-user@lists.osgeo.org">grass-user@lists.osgeo.org</a>><br>
Subject: Re: [GRASS-user] Large Grass Raster Notes<br>
Message-ID:<br>
<<a href="mailto:CABo5uVsfKnMC3O%2B-g%2BkdK-yj4-29yhBLzvryTFeQrOL9VgPuaw@mail.gmail.com">CABo5uVsfKnMC3O+-g+kdK-yj4-29yhBLzvryTFeQrOL9VgPuaw@mail.gmail.com</a>><br>
Content-Type: text/plain; charset="utf-8"<br>
<br>
On Sun, Jan 18, 2015 at 1:03 PM, Jeshua Lacock <<a href="mailto:jeshua@3dtopo.com">jeshua@3dtopo.com</a>> wrote:<br>
<br>
> Actually, I misspoke, the wildcard will not work with long lists either as<br>
> it will crash with some kind argument too long message as well.<br>
><br>
<br>
This is strange, the wild card is interpreted inside the module, so the<br>
argument list shouldn't be too long. If you don't quite the pattern the<br>
shell can do the expansion but this is not what you want.<br>
<br>
Correct:<br>
<br>
g.remove ... pattern="landsat_*"<br>
<br>
Incorrect:<br>
<br>
g.remove ... pattern=landsat_*<br>
<br>
The later will transform the command line to:<br>
<br>
g.remove ... pattern=landsat_01 landsat_02 landsat_03 landsat_04 ...<br>
<br>
in case you are in directory with files or directories named landsat_01,<br>
landsat_02, landsat_03, ...<br>
<br>
What is the message you get? And what is the command return code (do echo<br>
$? to find out)?<br>
-------------- next part --------------<br>
An HTML attachment was scrubbed...<br>
URL: <<a href="http://lists.osgeo.org/pipermail/grass-user/attachments/20150120/6678f2d5/attachment-0001.html" target="_blank">http://lists.osgeo.org/pipermail/grass-user/attachments/20150120/6678f2d5/attachment-0001.html</a>><br>
<br>
------------------------------<br>
<br>
_______________________________________________<br>
grass-user mailing list<br>
<a href="mailto:grass-user@lists.osgeo.org">grass-user@lists.osgeo.org</a><br>
<a href="http://lists.osgeo.org/mailman/listinfo/grass-user" target="_blank">http://lists.osgeo.org/mailman/listinfo/grass-user</a><br>
<br>
End of grass-user Digest, Vol 105, Issue 36<br>
*******************************************<br>
</blockquote></div><br><br clear="all"><div><br></div>-- <br><div class="gmail_signature"><div dir="ltr">Desire can be of greater import than utility; experience more important than closure.</div></div>
</div>