Re: [GRASS-user] Reclass: Atributte a string

[Franz Schiller]

Hello Glynn

Thanks for the reply. About using r.categories, can I apply it to a raster
map and then export it as a "normal" raster?
Thanks
Franz

On Fri, Nov 12, 2010 at 6:18 AM, Glynn Clements wrote:

>
> Franz Schiller wrote:
>
> > I'm doing a reclass of rasters based on a rule file. In this file I'm
> doing
> > a correspondence between values (1,2... 100) to other integer values.
> > Is it possible to do this but instead to done this to Numbers do this to
> > Strings just like
> > 1 = "Agriculture"
> > 2,3,4= "Forest"
>
> You can associate category labels with each category using
> r.categories. You can't "reclass" to a string.
>
> --
> Glynn Clements 
>
___
grass-user mailing list
grass-user@lists.osgeo.org
http://lists.osgeo.org/mailman/listinfo/grass-user

[GRASS-user] Re: Exporting Vector data to SHP (v.out.ogr)- takes a lot of time

[Luis Lisboa]

The problem is that, I left it this night processing and, when I
returned it has freezed. So, v.out.ogr is not being able to export
vector files produced by r.to.vect... As I said, maybe it's related
with some option in r.to.vect that is messing this up. Or something
related with tables.

Can anyone give me some tips on this?

Thanks

Luis


>...that´ll be a HUGE vector file!
>So it´s normal, that it takes some minutes.


On Thu, Nov 11, 2010 at 3:08 PM, Luis Lisboa wrote:

> Greetings
>
> I have produced a vector map (with r.to.vector) with the following
> metadata:
>
> v.info test01
>
>  
> ++
>  | Layer:   teste01  |
>  | Mapset:  Regional
>|
>  | Location:Country
>  |
>  | Database:E:\v3
>   |
>  | Title:
>   |
>  | Map scale:   1:1
>   |
>  | Map format:  native
>|
>  | Name of creator: orbit  |
>  | Organization:
>|
>  | Source date: Tue Nov 09 11:11:40 2010
>|
>
>  
> ||
>  |   Type of Map:  vector (level: 2)
>|
>  |
>|
>  |   Number of points:   0   Number of areas:  213680
>   |
>  |   Number of lines:0   Number of islands:73649
>|
>  |   Number of boundaries:   348770  Number of faces:  0
>|
>  |   Number of centroids:213666  Number of kernels:0
>|
>  |
>|
>  |   Map is 3D:  No
>   |
>  |   Number of dblinks:  1
>|
>  |
>|
>  | Projection: Universal Transverse Mercator (zone 0)
>   |
>  |   N:  8832707.83027649S:  7028073.62772369
>   |
>  |   E:  1347458.42001915W:200792.5777576
>   |
>  |
>|
>  |   Digitization threshold: 0
>|
>  |   Comments:
>|
>  |
>|
>
>  
> ++
>
>
> Now I'm trying to export it to SHAPEFILE but, after 10 minutes the
> processing bar is not even in 40%. What might be happening? I mean, is it
> some error with the vector or when I convert to vector (from raster)?
> Because, taking more than 10 minutes to export to vector it's not normal (I
> guess) and taking nearly 7 min to convert to vector also...
>
> Thanks for the reply
> Luis L
>
___
grass-user mailing list
grass-user@lists.osgeo.org
http://lists.osgeo.org/mailman/listinfo/grass-user

Re: [GRASS-user] Re: Exporting Vector data to SHP (v.out.ogr)- takes a lot of time

[Achim Kisseler]

Am 12.11.2010 09:40, schrieb Luis Lisboa:

The problem is that, I left it this night processing and, when I returned it 
has freezed. So, v.out.ogr is not being able to export vector files produced by 
r.to.vect... As I said, maybe it's related with some option in r.to.vect that 
is messing this up. Or something related with tables.

Can anyone give me some tips on this?


I guess it freezes because it uses so much memory.

Please check if v.out.ogr works with a small vector map. E.g. shrink 
region before r.to.vect or choose only some features from the vector map 
to export.


?

Achim


Thanks

Luis



...that´ll be a HUGE vector file!
So it´s normal, that it takes some minutes.



On Thu, Nov 11, 2010 at 3:08 PM, Luis Lisboa mailto:luislisboa1...@gmail.com>> wrote:

Greetings

I have produced a vector map (with r.to.vector) with the following
metadata:
v.info  test01
  
++
  | Layer:   teste01
  |
  | Mapset:  Regional
|
  | Location:Country
|
  | Database:E:\v3
   |
  | Title:
   |
  | Map scale:   1:1
   |
  | Map format:  native
|
  | Name of creator: orbit
|
  | Organization:
|
  | Source date: Tue Nov 09 11:11:40 2010
|
  
||
  |   Type of Map:  vector (level: 2)
|
  |
|
  |   Number of points:   0   Number of areas:
  213680 |
  |   Number of lines:0   Number of islands:
  73649  |
  |   Number of boundaries:   348770  Number of faces:
  0  |
  |   Number of centroids:213666  Number of kernels:
  0  |
  |
|
  |   Map is 3D:  No
   |
  |   Number of dblinks:  1
|
  |
|
  | Projection: Universal Transverse Mercator (zone 0)
   |
  |   N:  8832707.83027649S:  7028073.62772369
   |
  |   E:  1347458.42001915W:200792.5777576
   |
  |
|
  |   Digitization threshold: 0
|
  |   Comments:
|
  |
|
  
++


Now I'm trying to export it to SHAPEFILE but, after 10 minutes the
processing bar is not even in 40%. What might be happening? I mean,
is it some error with the vector or when I convert to vector (from
raster)? Because, taking more than 10 minutes to export to vector
it's not normal (I guess) and taking nearly 7 min to convert to
vector also...

Thanks for the reply
Luis L




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

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

[GRASS-user] Re: Need some help in converting a contour map to DEM map...

[nikmor]

Hello sir

I am very thankful to all who helped me, I have got the correct output after
setting the resolution and running v.surf.rst. I am attaching the nviz view
of the DEM after exagerating it.Now I am looking forward to drape an aerial
image of the same area on the DEM that I have got. In order to do this I
loaded the DEM map and then draped the imported aerial map (imported using
r.in.gdal, but that single jpeg image(satimage.jpg) got imported into 3
different maps namely satimage.red, satimage.green, satimage.blue, therefore
is just draped only one of it on the DEM). The output view that I get seems
okay but its too dark (black and white-- please see the attachment), though
the original aerial image was colourful which I had imported. Why is it
happening so?? Can any one help me with this.

regards
http://osgeo-org.1803224.n2.nabble.com/file/n5731657/iit_grass_images.tar.gz
iit_grass_images.tar.gz 
Nikhil
-- 
View this message in context: 
http://osgeo-org.1803224.n2.nabble.com/Need-some-help-in-converting-a-contour-map-to-DEM-map-tp5720519p5731657.html
Sent from the Grass - Users mailing list archive at Nabble.com.
___
grass-user mailing list
grass-user@lists.osgeo.org
http://lists.osgeo.org/mailman/listinfo/grass-user

Re: [GRASS-user] Re: Need some help in converting a contour map to DEM map...

[Achim Kisseler]

of the DEM after exagerating it.Now I am looking forward to drape an aerial
image of the same area on the DEM that I have got. In order to do this I
loaded the DEM map and then draped the imported aerial map (imported using
r.in.gdal, but that single jpeg image(satimage.jpg) got imported into 3
different maps namely satimage.red, satimage.green, satimage.blue, therefore
is just draped only one of it on the DEM). The output view that I get seems
okay but its too dark (black and white-- please see the attachment), though
the original aerial image was colourful which I had imported. Why is it
happening so?? Can any one help me with this.


Hi Nikhil,

take a look at
http://grass.fbk.eu/grass64/manuals/html64_user/imageryintro.html
it might give you a good introduction.
especially this part:

Visualizing (true) color composites
To quickly combine the first three channels to a near natural color 
image, the GRASS command d.rgb can be used or the graphical GIS manager 
(gis.m). It assigns each channel to a color which is then mixed while 
displayed. With a bit more work of tuning the grey scales of the 
channels, nearly perfect colors can be achieved. Channel histograms can 
be shown with d.histogram.


Regards,
Achim

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

Re: [GRASS-user] Re: Exporting Vector data to SHP (v.out.ogr)- takes a lot of time

[Luis Lisboa]

Hello Achim... Not really. It didn't work
Let me start from the beggining:
1- My orioginal raster file has the following:
|   Type of Map:  raster   Number of Categories: 11
 |
 |   Data Type:CELL
  |
 |   Rows: 5381
  |
 |   Columns:  4481
  |
 |   Total Cells:  24112261
  |
 |Projection: UTM (zone -24)
 |
 |N: 9659623.5017993S: 8045266.26306668   Res: 300.01063719
  |
 |E: 1077874.87767815W: -266394.50364563   Res: 299.993167
 |
 |   Range of data:min = 1  max = 11
 |
It's a landcover map

2- I convert to vectorusing the following expression
r.to.vect -v input=teste01_l...@regional output=teste05 feature=area
and it returns:
Extracting areas...
Building topology for vector map ...
Registering primiti 907000907256 primitives registered
3527832 vertices registered
Building areas...
321179 areas built
85875 isles built
Attaching islands...
Attaching centroids...
Number of nodes: 671952
Number of primitives: 907256
Number of points: 0
Number of lines: 0
Number of boundaries: 586103
Number of centroids: 321153
Number of areas: 321179
Number of isles: 85875
Number of areas without centroid: 26
r.to.vect complete.
My vector is:
 ||
 |   Type of Map:  vector (level: 2)
 |
 |
 |
 |   Number of points:   0   Number of areas:  321179
  |
 |   Number of lines:0   Number of islands:85875
 |
 |   Number of boundaries:   586103  Number of faces:  0
 |
 |   Number of centroids:321153  Number of kernels:0
 |
 |
 |
 |   Map is 3D:  No
  |
 |   Number of dblinks:  1
 |
 |
 |
 | Projection: Universal Transverse Mercator (zone 0)
  |
 |   N:   9659623.5017993S:  8045266.26306668
  |
 |   E:   911978.65632755W:  -266394.50364563
  |
 |
 |
 |   Digitization threshold: 0
 |
 |   Comments:
 |
 |
 |
 ++
this is a vector composed by areas (as far as I wanted and I believe)
3- I reduced the region in order to export just a piece
4- I exported but, when I open it in QGIS I get an empty vector file
5- Then again: v.out.ogr -c input=test...@regional type=area
dsn=c:\teste07.shp
Exporting 321179 areas (may take some time)...
And after 20 minutes it seems freezed

It this too much polygons?

On Fri, Nov 12, 2010 at 9:06 AM, Achim Kisseler  wrote:

>
>
> Am 12.11.2010 09:40, schrieb Luis Lisboa:
>
>  The problem is that, I left it this night processing and, when I returned
>> it has freezed. So, v.out.ogr is not being able to export vector files
>> produced by r.to.vect... As I said, maybe it's related with some option in
>> r.to.vect that is messing this up. Or something related with tables.
>>
>> Can anyone give me some tips on this?
>>
>
> I guess it freezes because it uses so much memory.
>
> Please check if v.out.ogr works with a small vector map. E.g. shrink region
> before r.to.vect or choose only some features from the vector map to export.
>
> ?
>
> Achim
>
>>
>> Thanks
>>
>> Luis
>>
>>
>>  ...that´ll be a HUGE vector file!
>>> So it´s normal, that it takes some minutes.
>>>
>>
>>
>> On Thu, Nov 11, 2010 at 3:08 PM, Luis Lisboa > > wrote:
>>
>>Greetings
>>
>>I have produced a vector map (with r.to.vector) with the following
>>metadata:
>>v.info  test01
>>
>>
>>  
>> ++
>>  | Layer:   teste01
>>  |
>>  | Mapset:  Regional
>>|
>>  | Location:Country
>>|
>>  | Database:E:\v3
>>   |
>>  | Title:
>>   |
>>  | Map scale:   1:1
>>   |
>>  | Map format:  native
>>|
>>  | Name of creator: orbit
>>|
>>  | Organization:
>>|
>>  | Source date: Tue Nov 09 11:11:40 2010
>>|
>>
>>  
>> ||
>>  |   Type of Map:  vector (level: 2)
>>|
>>  |
>>|
>>  |   Number of points:   0   Number of areas:
>>  213680 |
>>  |   Number of lines:0   Number of islands:
>>  73649  |
>>  |   Number of boundaries:   348770  Number of faces:
>>  0  |
>>  |   Number of centroids:213666  Number of kernels:
>>  0  |
>>  |
>>|
>>  |   Map is 3D:  No
>>   |
>>  |   Number of dblinks:  1
>>|
>>  |
>>|
>>  | Projection: Universal Transverse Mercator (zone 0)
>>   |
>>  |   N:  8832707.83027649S:  7028073.62772369
>>

[GRASS-user] current state of r.li

[Paolo Cavallini]

Hi all.
Sorry for crossposting. As some of you know, the r.li suite of GRASS 
commands allows landscape analyses[0]. Its interface is rather complex, 
and is still in TclTk, not ported to either wxpython or qgis. As such, 
it is now more difficult to use than it should be, and it will become 
unusable when TclTk support will be dropped.
The possible solution (thanks Radim) is to rewrite the interface as a 
qgis python plugin. It should not be a huge work (we provisionally 
estimate 2-3 weeks).
The question is: is there anybody willing to invest either his/her time, 
or some money, to write such a plugin?

We (Faunalia) would be happy to help if necessary.
All the best.
--
http://www.faunalia.it/pc

[0]http://grass.fbk.eu/gdp/html_grass64/r.li.setup.html
___
grass-user mailing list
grass-user@lists.osgeo.org
http://lists.osgeo.org/mailman/listinfo/grass-user

[GRASS-user] Re: Need some help in converting a contour map to DEM map...

[nikmor]

Hello sir 

Thanks for the link. I finally got the coloured aerial image draped on the
DEM map and which can be viewed in nviz. Since the aerial image was getting
imported in an rgb files (3 different files) I first needed to combine them
into one raster map, which could be done using r.composite. which takes in
these 3 files (.red, .green, .blue) and one output map name of users choice.
Thanks once again for your help sir.

Nikhil Morajkar
Junior Software Engineer,
Department of Computer Science and Engineering,
Geo-Spatial Information Science & Engineering (GISE),
Advance Research Lab,
Indian Institute of Technology, Bombay (IITB)
Powai, Mumbai
-- 
View this message in context: 
http://osgeo-org.1803224.n2.nabble.com/Need-some-help-in-converting-a-contour-map-to-DEM-map-tp5720519p5732236.html
Sent from the Grass - Users mailing list archive at Nabble.com.
___
grass-user mailing list
grass-user@lists.osgeo.org
http://lists.osgeo.org/mailman/listinfo/grass-user

[GRASS-user] How to delete deadends?

[Patrick_schirmer]

Dear group,

I have a lot of lines that are linked to another forming various nets. 
Several of those lines are deadends. Now I search for a option to delete 
"dead-ends" within those networks. It would be perfect to delete the 
lines, or to search for the nodes that link to more than one line.
I was searching in v.generalize, v.net, v.to.point, v.to.db but won't 
find the proper approach.


Doas anybody of you have an idea?

Thanks for your feedback,

Patrick Schirmer
___
grass-user mailing list
grass-user@lists.osgeo.org
http://lists.osgeo.org/mailman/listinfo/grass-user

Re: [GRASS-user] How to delete deadends?

[Markus Metz]

On Fri, Nov 12, 2010 at 2:07 PM, Patrick_schirmer  wrote:
> Dear group,
>
> I have a lot of lines that are linked to another forming various nets.
> Several of those lines are deadends. Now I search for a option to delete
> "dead-ends" within those networks. It would be perfect to delete the lines,
> or to search for the nodes that link to more than one line.
> I was searching in v.generalize, v.net, v.to.point, v.to.db but won't find
> the proper approach.

Try to remove dangles with v.clean tool=rmdangle
If threshold is < 0, all dangles will be removed, otherwise only
dangles shorter than threshold.

Markus M
___
grass-user mailing list
grass-user@lists.osgeo.org
http://lists.osgeo.org/mailman/listinfo/grass-user

[GRASS-user] i.smap and NULL (Landsat 7)

[Nick Jachowski]

Dear Grass List,

I am trying to do radiometric classifications on Landsat 7 SLC-Off imagery
(i.e. it has lots of NULL data).  My question is regarding i.smap.  The
documentation says:

"The module i.smap does not support MASKed or NULL cells. Therefore it might
be necessary to create a copy of the classification results using e.g.
r.mapcalc."

Does this mean i.smap will not work correctly with Landsat 7 SLC-Off
imagery?  When I try i.gensigset and i.smap I get results, but I'm not sure
if I can trust them due to the disclaimer.

Thanks for any help or advice,

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

Re: [GRASS-user] How to delete deadends?

[Hamish]

Patrick_schirmer wrote:
> I have a lot of lines that are linked to another
> forming various nets. Several of those lines are
> deadends. Now I search for a option to delete
> "dead-ends" within those networks. It would be
> perfect to delete the lines, or to search for the
> nodes that link to more than one line.
> I was searching in v.generalize, v.net,
> v.to.point, v.to.db but won't find the proper
> approach.

not a v.net specific thing, but in general for
vector maps you can use v.clean to remove (or isolate)
line "dangles" with the 'rmdangle' (or 'chdangle')
tools.


maybe that helps,
Hamish



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

Re: [GRASS-user] Re: Exporting Vector data to SHP (v.out.ogr)- takes a lot of time

[Hamish]

> Luis wrote:
> 1- My orioginal raster file has the following:
...
> Registering primiti 907000
> 907256 primitives registered

900k is a lot, but not crazy.

> 3527832 vertices registered
> Building areas...
> 321179 areas built

300k areas is not small, but it is not huge either.
...
> Exporting 321179 areas (may take some time)...
>
> And after 20 minutes it seems freezed

wait longer... maybe hours or overnight.
export to shapefile with v.out.ogr is known to be
quite slow.

are you sure it is frozen? check in process monitor
("top" on mac or linux, ctrl-alt-del on Windows)
if cpu usage by v.out.ogr is 100%, and how much
memory it is using. if it is frozen cpu will be 0%,
if it is still working cpu will be 100%.


> It this too much polygons?

perhaps. maybe resample raster to a coarser grid
or try running v.build.polylines on the vector map
prior to export.


Hamish



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

Re: [GRASS-user] Reclass: Atributte a string

[Glynn Clements]

Franz Schiller wrote:

> Thanks for the reply. About using r.categories, can I apply it to a raster
> map and then export it as a "normal" raster?

category labels are attached to an existing map; there's nothing to
export.

It might help if you explain what your ultimate goal is.

-- 
Glynn Clements 
___
grass-user mailing list
grass-user@lists.osgeo.org
http://lists.osgeo.org/mailman/listinfo/grass-user