[GRASS-user] Patching Maps

Tim Southern tim.southern at talktalk.net
Mon May 19 13:00:33 PDT 2014


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 at lists.osgeo.org wrote:

> Send grass-user mailing list submissions to
> 	grass-user at 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 at lists.osgeo.org
> 
> You can reach the person managing the list at
> 	grass-user-owner at 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 at gmail.com>
> To: Tim Southern <tim.southern at talktalk.net>
> Cc: Grass Users <grass-user at lists.osgeo.org>
> Subject: Re: [GRASS-user] patching maps together
> Message-ID:
> 	<CADCcPvto6TUC4m_5UemZg-mc6psgvc6BRupYX8TchVcVvG_AhQ at 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 at 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 at 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 at gmail.com>
> To: Tim Southern <tim.southern at talktalk.net>
> Cc: grass-user at lists.osgeo.org
> Subject: Re: [GRASS-user] patching maps together
> Message-ID: <20140520000536.25baa063189ea1fd178f4252 at 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 at gmail.com>
> 
> 
> ------------------------------
> 
> Message: 3
> Date: Mon, 19 May 2014 16:21:55 +0200
> From: "G. Allegri" <giohappy at gmail.com>
> To: Vincent Bain <bain at toraval.fr>
> Cc: GRASS user list <grass-user at lists.osgeo.org>
> Subject: Re: [GRASS-user] GRASS and Blender
> Message-ID:
> 	<CAB4g1=wenN5SiAUe1nUgOiLnfoXXARKVzjMrx7Rev-YERhsGag at 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 at 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>
> 
> ------------------------------
> 
> Message: 4
> Date: Mon, 19 May 2014 14:24:33 +0000
> From: "Jon Eiriksson" <jeir at hi.is>
> To: "GRASS user list" <grass-user at lists.osgeo.org>
> Subject: [GRASS-user] GRASS 7.0beta2
> Message-ID: <851236F3-B068-4A4E-9107-AA4B13B86425 at 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 at lists.osgeo.org
> http://lists.osgeo.org/mailman/listinfo/grass-user
> 
> End of grass-user Digest, Vol 97, Issue 18
> ******************************************



More information about the grass-user mailing list