[GRASS-user] r.thin ?

charlie carowan at yahoo.com
Tue Nov 9 11:21:44 EST 2010


Micha and Achim, thank you for the
 replies.  In response to your questions, I am trying to convert a raster map (a scanned & georeferenced map with only XY data) with a shoreline and some wetland islands into a vector map so I can compare gains/losses over time.

A preview of the map is here: http://historicalcharts.noaa.gov/tiled_jpgs_done/zoomifyURLDrivenWebPage.htm?zoomifyImagePath=579-12-1929

The command I used is as follows:
r.thin input=NOAA1929 output=NOAA1929THIN
 iterations=6000 

I have tried r.thin with 1000 iterations, which did successfully thin the lines down to the features borders quite well, however only 50% of the map was completed.

I bumped the iterations up to 6000, it completed in 4200, but it seemed that it thinned the entire map down to a 1px line.

Thanks

--- On Tue, 11/9/10, Micha Silver
 <micha at arava.co.il> wrote:

From: Micha Silver <micha at arava.co.il>
Subject: Re: [GRASS-user] r.thin ?
To: "charlie" <carowan at yahoo.com>
Cc: grass-user at lists.osgeo.org
Date: Tuesday, November 9, 2010, 5:24 AM

charlie wrote:

> Hi.  I am currently using Grass 6.4 on linux and beginning the the process of vectorizing an old NOAA raster map with r.thin. The map is approximately 12,000 pixels by 8,000 pixels wide.
> 
> According to this post, ( http://osgeo-org.1803224.n2.nabble.com/r-thin-td2718260.html#a2718260 )the iteration variable has to be equal or larger than half the number of cells of the wider line so I
 increased the
 iterations to 6,000, and it would
 up finishing in approximately 4,200 passes.
> 
> Previously, I was getting excellent preliminary results, with fewer (1,000) iterations,  distinct features were visible, however at least 50% of the map was not completed. 
> Now with an increase to 6,000 iterations r.thin completes, but the results are a single 1px wide horizontal line.  Im sure I am overlooking something relatively simple; can anyone give me any advice on how to best run r.thin in this particular case?
> 

What is it you are trying to do? The r.thin module is designed to make line features in a raster exactly one pixel wide so that when converting to a line vector, the line will be clean.
If you run r.thin on a raster that represents areas, and give a high enough iteration value, the areas will eventually be reduced to lines.

> 
> Thanks!
> 
> Charlie
> 
> 
> <a
 href="http://www.wetlandresearch.com">Wetlands</a>
> 
> 
> This mail was received via Mail-SeCure System.
> ------------------------------------------------------------------------
> 
> _______________________________________________
> grass-user mailing list
> grass-user at lists.osgeo.org
> http://lists.osgeo.org/mailman/listinfo/grass-user
> 
> This mail was received via Mail-SeCure System.
> 
> 
>   

-------------- next part --------------
An HTML attachment was scrubbed...
URL: http://lists.osgeo.org/pipermail/grass-user/attachments/20101109/9e9ced4b/attachment-0001.html


More information about the grass-user mailing list