[GRASS-user] Patching Maps

GRASS GRASS user
,
1

Dear Nick and Hamish,

Thanks, when I used the query raster over the individual maps for both black and white most gave coordinates plus a 0 or 1, some however gave an * for value. I suspect that even though all the tiles were from the same source and processed in the same way some of the .tif files have ended up as greyscale images.

I have now run ImageMagick convert -monochrome and will reload the tiles and check.

If I still have problems I will come back to the list.

Kind regards

Tim
Tim Southern
17, Park Close,
Sonning Common,
Oxfordshire
RG4 9RY

On 19 May 2014, at 20:00, grass-user-request@lists.osgeo.org wrote:

Send grass-user mailing list submissions to
  grass-user@lists.osgeo.org

To subscribe or unsubscribe via the World Wide Web, visit
  http://lists.osgeo.org/mailman/listinfo/grass-user
or, via email, send a message with subject or body 'help' to
  grass-user-request@lists.osgeo.org

You can reach the person managing the list at
  grass-user-owner@lists.osgeo.org

When replying, please edit your Subject line so it is more specific
than "Re: Contents of grass-user digest..."

Today's Topics:

  1. Re: patching maps together (Nick Ves)
  2. Re: patching maps together (Hamish)
  3. Re: GRASS and Blender (G. Allegri)
  4. GRASS 7.0beta2 (Jon Eiriksson)

----------------------------------------------------------------------

Message: 1
Date: Mon, 19 May 2014 12:39:23 +0300
From: Nick Ves <vesnikos@gmail.com>
To: Tim Southern <tim.southern@talktalk.net>
Cc: Grass Users <grass-user@lists.osgeo.org>
Subject: Re: [GRASS-user] patching maps together
Message-ID:
  <CADCcPvto6TUC4m_5UemZg-mc6psgvc6BRupYX8TchVcVvG_AhQ@mail.gmail.com>
Content-Type: text/plain; charset=UTF-8

Hey Tim,

If your maps don't overlap with each other you could use r.series in your case :

$ g.region rast=$(g.mlist rast sep=",")

$ r.series input=$(g.mlist rast sep=",") out=myMap method=sum

N

On Mon, May 19, 2014 at 11:31 AM, Tim Southern
<tim.southern@talktalk.net> wrote:

Dear All,

I am trying to patch together a number of 1 km by 1 km tiles of a map whose values are either 0 or 1 (black on white background).

If I use r.patch to put 2 tiles together is appears to work fine but if however I try and patch more, even if it is two tiles of already patched tiles I get one set when patched as black background and white features and the other part as white background and black features. It appears that the second in the case of already patched data or third and subsequent in case of more than 2 tiles is inverted with 1 becoming 0 and 0 becoming 1.

I have tried reducing region to the smallest. Maps are at 0.2m resolution. I am hoping to get a 5 tile by 5 tile block

I have tried this across two different platforms (Mac OS and OpenSuse ) on the latest downloadable binary versions of GRASS.

Am I doing something wrong?

Kind regards

Tim
Tim Southern

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

------------------------------

Message: 2
Date: Tue, 20 May 2014 00:05:36 +1200
From: Hamish <hamish.webmail@gmail.com>
To: Tim Southern <tim.southern@talktalk.net>
Cc: grass-user@lists.osgeo.org
Subject: Re: [GRASS-user] patching maps together
Message-ID: <20140520000536.25baa063189ea1fd178f4252@gmail.com>
Content-Type: text/plain; charset=US-ASCII

Tim wrote:

I am trying to patch together a number of 1 km by 1 km tiles of a map
whose values are either 0 or 1 (black on white background).

If I use r.patch to put 2 tiles together is appears to work fine but
if however I try and patch more, even if it is two tiles of already
patched tiles I get one set when patched as black background and
white features and the other part as white background and black
features. It appears that the second in the case of already patched
data or third and subsequent in case of more than 2 tiles is inverted
with 1 becoming 0 and 0 becoming 1.

I have tried reducing region to the smallest. Maps are at 0.2m
resolution. I am hoping to get a 5 tile by 5 tile block

I have tried this across two different platforms (Mac OS and OpenSuse
) on the latest downloadable binary versions of GRASS.

Am I doing something wrong?

I suspect that the images are paletted, and the color palettes differ
between the images. After confirming that you might run the images
through a bulk color palette -> 24bit RGB filter. (I'm not sure what the
right NetPBM or ImageMagick command is, but I'm sure there is one)

after that, gdalbuildvrt is wonderful for import + patch in one step.

good luck,

--
Hamish <hamish.webmail@gmail.com>

------------------------------

Message: 3
Date: Mon, 19 May 2014 16:21:55 +0200
From: "G. Allegri" <giohappy@gmail.com>
To: Vincent Bain <bain@toraval.fr>
Cc: GRASS user list <grass-user@lists.osgeo.org>
Subject: Re: [GRASS-user] GRASS and Blender
Message-ID:
  <CAB4g1=wenN5SiAUe1nUgOiLnfoXXARKVzjMrx7Rev-YERhsGag@mail.gmail.com>
Content-Type: text/plain; charset="utf-8"

Thanks Vincent for sharing your experience.
Indeed there aren't production ready (open source) tools to manage the
whole process.
Point data are relatively easier to manage (though they generally bring
their complexities), but CityGML like models, geological models, etc. lack
a ready to use toolbox.
I know that 3D brings exponential difficulties then 2D, and probably theres
isn't such a widespread need for it to have enough investments to support
its development in the FOSS context.
A good work is being done by Sourcepole, with PostGIS + (SF)CGAL, but the
road is still very long...

giovanni

2014-05-19 11:43 GMT+02:00 Vincent Bain <bain@toraval.fr>:

Le lundi 19 mai 2014 ? 09:57 +0200, G. Allegri a ?crit :

Do you have experimented robust workflows that can do one or both
things?

Hello Giovanni,
yes in a production context (natural hazrd ingeneering) I do this kind
of back and forth quite often, but with limited functionality: for my
personal needs, it is mainly a question of 3d point clouds export/import
with coordinate system switching; unfortunately I cannot give any advice
concerning the consistency of attribute data. To keep record of an
existing database structure within GRASS, perhaps a first step would be
to export 3d objects with their cat values, so one would retrieve it
when importing back to GRASS.
In blender for example I know you can add custom properties to objects
or classes of objects, but don't know how simple it is to link it with
e.g. .vtk "POINT_DATA" or .ply feature "property". I'll dig into that.

V.

--
Giovanni Allegri
http://about.me/giovanniallegri
Twitter: https://twitter.com/_giohappy_
blog: http://blog.spaziogis.it
GEO+ geomatica in Italia http://bit.ly/GEOplus
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://lists.osgeo.org/pipermail/grass-user/attachments/20140519/d03f471c/attachment-0001.html&gt;

------------------------------

Message: 4
Date: Mon, 19 May 2014 14:24:33 +0000
From: "Jon Eiriksson" <jeir@hi.is>
To: "GRASS user list" <grass-user@lists.osgeo.org>
Subject: [GRASS-user] GRASS 7.0beta2
Message-ID: <851236F3-B068-4A4E-9107-AA4B13B86425@hi.is>
Content-Type: text/plain; format=flowed

Hi,

I have been running GRASS 7.0 SVN (and frameworks) from Michael Barton's
site since Feruary on two Macs, MacOSX 10.7.5 and MacOSX Mavericks,
happily, on a daily basis. Tried installing the GARSS 7.0beta2 package
and frameworks on the 10.7.5 machine, installation succeeded, but I got
this message on startup:
---
| ~ @ nfmac129 (je)
| => '/Applications/GRASS/GRASS-7.0.app/Contents/MacOS/grass.sh'; exit
Rebuilding Addon HTML manual pages index...
Rebuilding Addon menu...
Python 2.7.6 found.
arch: posix_spawnp: /opt/local/bin/pythonw2.7: Bad CPU type in
executable
logout
---

In Terminal, [which python] gives:
---
| ~ @ nfmac129 (je)
| => which python
/Library/Frameworks/Python.framework/Versions/2.7/bin/python

---

Tried reinstalling the February Frameworks and GRASS 7.0 SVN package,
but get the same message.

On the Mavericks machine the installation went OK, GRASS 7.0 started
normally and the wxpython gui seems to work, but a few Terminal commands
I tried give error messages (e.g. [g.list type=vect] is not recognised).

Installation of GRASS-7.0beta1 failed when I tried about a month ago,
decided to wait for the next beta release.

As I do not know what the message means (is there something wrong on my
mac, perhaps?) Any ideas would be welcome, have other users had similar
problems?

Jon

------------------------------

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

End of grass-user Digest, Vol 97, Issue 18
******************************************