RE: [GRASS-user] DXF to UTM32

[Pablo Carreira]

> From: nete...@osgeo.org
> Date: Sun, 22 Nov 2009 12:58:48 +0100

> Pablo,
> 
> you notes are very interesting! Could you add them into a related Wiki page:
> http://grass.osgeo.org/wiki/Import_DXF


I have added them. If someone could review it I would thank. 


Regards.

Pablo.





  
_
Novo windowslive.com.br. Descubra como juntar a galera com os produtos Windows 
Live.
http://www.windowslive.com.br/?ocid=WindowsLive09_MSN_Hotmail_Tagline_out09___
grass-user mailing list
grass-user@lists.osgeo.org
http://lists.osgeo.org/mailman/listinfo/grass-user

Re: [GRASS-user] g.region Options Yield Different Values

[Rich Shepard]

On Sun, 22 Nov 2009, Glynn Clements wrote:


It doesn't make any sense to specify lat/lon boundaries for projected
data, particularly for a non-cylindrical projection such as LCC.


Glynn,

  That's true. I've obviously erred when creating new locations and
importing .shp and .e00 files.

  I've used the wxPython GUI to define new locations and the mapset. When I
specify choosing the projection from the list and select LCC I'm asked to
enter the relevant projection parameters. I then specify NAD83 as the datum
and the Washington-Oregon adjustment. Finally, I'm asked if I want to enter
the default region's extent and resolution. This is where I enter the
lat/long coordinates I have.

  Perhaps it's this last step that I have wrong. Should I decline to enter
the extents and resolution and leave those to the information embedded in
the imput data?

  As I've written before it's been a very long time since I last used GRASS
and quite likely that I've mis-used the new GUI.

Thanks,

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

[GRASS-user] Piping Commands

[Pablo Carreira]

Hi, 

I want to do the following:
v.extract some lines from a vector and them
v.buffer those lines

How do I pipe these commands?





Pablo Torres Carreira


  
_
Novo windowslive.com.br. Descubra como juntar a galera com os produtos Windows 
Live.
http://www.windowslive.com.br/?ocid=WindowsLive09_MSN_Hotmail_Tagline_out09___
grass-user mailing list
grass-user@lists.osgeo.org
http://lists.osgeo.org/mailman/listinfo/grass-user

Re: [GRASS-user] Piping Commands

[Achim Kisseler]

Easiest way:

don't pipe!

Just make a new vectormap from v.extract and buffer that.

(I dont think piping is possible here and v.buffer cannot be restricted
on a selection)

Achim



Pablo Carreira schrieb:
> 
> Hi,
> 
> I want to do the following:
> v.extract some lines from a vector and them
> v.buffer those lines
> 
> How do I pipe these commands?
> 
> 
> 
> 
> 
> Pablo Torres Carreira
> 
> 
> 
> Chegou o Windows 7: Incrivelmente simples! Clique e conheça.
> 
> 
> 
> 
> 
> ___
> 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] Project Area from Larger Region

[Rich Shepard]

  I would appreciate points to create a defined subregion and apply one
raster layer and a few vector layers to that subregion. I don't see an
algorithm to do this in either the GRASS book or on the web site.

  When I tried to define the region using coordinates identified when
displaying one vector layer the values were rejected. In brief what I want
to do is:

  1) Identify the appropriate bounds that encompass the drainage basin in
which the project is located.

  2) Clip the larger regional DEM raster layer (probably using r.mask) to
copy the project area into the new location.

  2) Clip the vector layers (v.extract?) so they, too, can be copied into
the new projection location.

  3) Assign the relevant attribute SQLite data tables (db.connect) to this
new region's vector layers.

  After this I can start designing the model.

  Pointers to what I should read to accomplish this will be appreciated.

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

RE: [GRASS-user] Piping Commands

[Pablo Carreira]

In general, is it possible to connect the "out=" to the "in=" of vector 
commands in GRASS

Thanks.


Pablo Torres Carreira




> Date: Mon, 23 Nov 2009 16:14:55 +0100
> From: a...@jupiter.uni-freiburg.de
> To: pablotcarre...@hotmail.com
> CC: grass-user@lists.osgeo.org
> Subject: Re: [GRASS-user] Piping Commands
> 
> Easiest way:
> 
> don't pipe!
> 
> Just make a new vectormap from v.extract and buffer that.
> 
> (I dont think piping is possible here and v.buffer cannot be restricted
> on a selection)
> 
> Achim
> 
> 
> 
> Pablo Carreira schrieb:
> > 
> > Hi,
> > 
> > I want to do the following:
> > v.extract some lines from a vector and them
> > v.buffer those lines
> > 
> > How do I pipe these commands?
> > 
> > 
> > 
> > 
> > 
> > Pablo Torres Carreira
> > 
> > 
> > 
> > Chegou o Windows 7: Incrivelmente simples! Clique e conheça.
> > 
> > 
> > 
> > 
> > 
> > ___
> > grass-user mailing list
> > grass-user@lists.osgeo.org
> > http://lists.osgeo.org/mailman/listinfo/grass-user
  
_
Novo windowslive.com.br. Descubra como juntar a galera com os produtos Windows 
Live.
http://www.windowslive.com.br/?ocid=WindowsLive09_MSN_Hotmail_Tagline_out09___
grass-user mailing list
grass-user@lists.osgeo.org
http://lists.osgeo.org/mailman/listinfo/grass-user

[GRASS-user] r.prominence

[Bulent Arikan]

Hi,

I have been using r.param.scale and inverting DEM method for identifying
   peaks. I am specifically interested in finding the high spots on a
   landscape; not just the highest single cell. It has been suggested that
   r.prominence may be of help. I realise that this comes in a file that
needs
   to be compiled (.c extension). Before compiling, I will appreciate any
   insight on what it does.

   Thank you
-- 
BÜLENT ARIKAN
School of Human Evolution and Social Change
Arizona State University
Tempe - AZ
85287-2402
___
grass-user mailing list
grass-user@lists.osgeo.org
http://lists.osgeo.org/mailman/listinfo/grass-user

[GRASS-user] r.horizon problems with large DEM

[Joshua Campbell]

I have been trying to run r.horizon on a SRTM DEM for the continental US,
with no success. The operation runs for about 45 minutes then eventually
just stops with no warnings or errors

I'm running GRASS 6.4-rc5 64-bit in Ubuntu 9.10 on a virtual machine
(Virtual Box on Windows XP 64-bit).

The r.horizon function will consume as much RAM as apply to the machine, up
to 6144 MB, plus an additional 4GB of swap memory.

The DEM is in Geographic coordinates - WGS 84. I tried to run r.horizon in
only 1 direction with maximum distance of 5310 meters.

Is there a file size limitation? Or any suggestions for an optimal size?

Thanks,
Josh

Joshua Campbell
Kansas Applied Remote Sensing Program
University of Kansas
Lawrence, KS
___
grass-user mailing list
grass-user@lists.osgeo.org
http://lists.osgeo.org/mailman/listinfo/grass-user

Re: [GRASS-user] r.prominence

[Benjamin Ducke]

It's very simple. It really just averages the values of raster
cells for a neighbourhood window of a size m x n. It then determines
how much the cell in the window center deviates from that average.
Add some normalization and you get a map that shows you which cells
"stick out" most. Use it on a map with elevation values and you
can (tentatively) call that topographic prominence.
The module lets you choose different window sizes and shapes and
whether to use local (within window) or global (over whole map)
normalization.

Ben

- Original Message -
From: "Bulent Arikan" 
To: grass-user-requ...@lists.osgeo.org, grass-user@lists.osgeo.org
Sent: Monday, November 23, 2009 7:46:52 PM GMT +01:00 Amsterdam / Berlin / Bern 
/ Rome / Stockholm / Vienna
Subject: [GRASS-user] r.prominence


Hi, 

I have been using r.param.scale and inverting DEM method for identifying 
peaks. I am specifically interested in finding the high spots on a 
landscape; not just the highest single cell. It has been suggested that 
r.prominence may be of help. I realise that this comes in a file that needs 
to be compiled (.c extension). Before compiling, I will appreciate any 
insight on what it does. 

Thank you 
-- 
BÜLENT ARIKAN 
School of Human Evolution and Social Change 
Arizona State University 
Tempe - AZ 
85287-2402 

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


--
Files attached to this email may be in ISO 26300 format (OASIS Open Document 
Format). If you have difficulty opening them, please visit http://iso26300.info 
for more information.

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

[GRASS-user] Project Area from Larger Region

[Helmut Kudrnovsky]

>   I would appreciate points to create a defined subregion and apply one
>raster layer and a few vector layers to that subregion. I don't see an
>algorithm to do this in either the GRASS book or on the web site.
>
>   When I tried to define the region using coordinates identified when
>displaying one vector layer the values were rejected. In brief what I want
>to do is:
>
>   1) Identify the appropriate bounds that encompass the drainage basin in
>which the project is located.
>
>   2) Clip the larger regional DEM raster layer (probably using r.mask) to
>copy the project area into the new location.
>
>   2) Clip the vector layers (v.extract?) so they, too, can be copied into
>the new projection location.
>
>   3) Assign the relevant attribute SQLite data tables (db.connect) to this
>new region's vector layers.
>
>   After this I can start designing the model.
>
>   Pointers to what I should read to accomplish this will be appreciated.
>
>Rich

hi Rich,

maybe v.in.region could help you? 
(http://grass.osgeo.org/grass64/manuals/html64_user/v.in.region.html)

- zoom to your area of interest of your larger region

-  v.in.region e.g. "v.in.region output=testarea", so  the extent of your area 
of interest is defined by the vector "testarea"

-  "g.region vect=testarea", by this  the region is set to your area of 
interest for analysis. "This region defines the geographic area in which all 
GRASS displays and raster analyses will be done." 
(http://grass.osgeo.org/grass64/manuals/html64_user/g.region.html)

hope this helps
Helmut



__
GRATIS für alle WEB.DE-Nutzer: Die maxdome Movie-FLAT!
Jetzt freischalten unter http://movieflat.web.de



smime.p7s
Description: S/MIME Cryptographic Signature
___
grass-user mailing list
grass-user@lists.osgeo.org
http://lists.osgeo.org/mailman/listinfo/grass-user

[GRASS-user] EPSG:2192 not working

[Richard Chirgwin]

I really ought to stick to Australia where I know how to define locations!

The issue is this. Using EPSG:2192 to create a European location doesn't 
create the projection file correctly.


1. After the EPSG code selection panel, Grass-GIS offers a choice of 
transformation parameters. I have tried several, with the same result.


2. After the transformation selection, the following errors occur:

g.proj returned the following message:
ERROR: region for current mapset is not set run "g.region"

then

Warning: location  at GISDBASE  
is not a directory or does not exist.


(Note: there are no permissions problems, I am running another instance 
of Grass-GIS in a lat-long location to the same GISDBASE)


3. Back at the startup panel, the new location exists but won't launch:

Warning: PERMANENT is not a valid mapset

(no surprises there, since the region wasn't created properly).

Looking in the directory, I find the following:
- the PROJ file has not been created.
- the DEFAULT_WIND file is:
proj   99
zone   0

Being unfamiliar with Europe, I'd welcome either projection parameters 
suitable for France and Germany (ie, what I would put into the text 
location creation interface to get a suitable Lambert location), or a 
fix for the problem using EPSG:2192.


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

[GRASS-user] Further to location creation problem

[Richard Chirgwin]

A couple more things.

1 - the system is Grass 6.40RC5 on Linux.

2 - There's also a problem using the command line location creation 
interface. If I select UTM, it fails immediately after requesting the 
"location description".


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

RE: [GRASS-user] Piping Commands

[Pablo Carreira]

Hi, Javier!

Let me explain:
I have a very large vector of rivers, with many CATs. Here in Brazil we have a 
law that say that, in general, a 30m area aside of rivers should be preserved 
with trees.
Them I look in a high resolution orbital image and where I find a portion river 
without trees I v.extract that portion of river, v.buffer 30m, v.patch that 
buffer to larger common vector of buffers, g.remove the buffer and the 
extracted vector.

It is not a batch process, but I think it could be indeed scripted, and I would 
apreciate some help with that.

Anyway, I still want to discover is if it is possible to pipe grass commands, I 
don't know much of programing but if I'm not wrong piping mean that the output 
of an application becomes the input of another application.

Achim told me not to pipe, but I am a obstinate kind of person (I think most of 
us are)...

ops, I almost forgot, in the company I use windows (becouse I can't use linux 
there), but at home I use Ubuntu linux.

Regards

Pablo Torres Carreira




From: fjgarciapri...@hotmail.com
To: pablotcarre...@hotmail.com
Subject: RE: [GRASS-user] Piping Commands
Date: Mon, 23 Nov 2009 19:28:29 +








Hi, I'm not an expert or anything else, but maybe I can help you.

You are using Windows or Linux?

Sometime ago I have to upload a large number of images (Ortophotos) to grass 
(more than the command aloud me and I use a Linux script to do that performed 
over the basis of processing each group of files and then perform another 
operations on the result.

Is this what are you looking for?

From: pablotcarre...@hotmail.com
To: grass-user@lists.osgeo.org
Date: Mon, 23 Nov 2009 12:56:46 -0200
Subject: [GRASS-user] Piping Commands









Hi, 

I want to do the following:
v.extract some lines from a vector and them
v.buffer those lines

How do I pipe these commands?





Pablo Torres Carreira


  
Chegou o Windows 7: Incrivelmente simples! Clique e conheça.
  
Windows Live:  Friends get your Flickr, Yelp, and Digg updates when they e-mail 
you.  
_
Novo site do Windows Live: Novidades, dicas dos produtos e muito mais. Conheça!
http://www.windowslive.com.br/?ocid=WindowsLive09_MSN_Hotmail_Tagline_out09___
grass-user mailing list
grass-user@lists.osgeo.org
http://lists.osgeo.org/mailman/listinfo/grass-user

Re: [GRASS-user] Piping Commands

[Alex Mandel]

I would consider it best practice to run 2 commands in sequence, storing
the intermediate product. It can save you a lot of time, help
troubleshoot data errors, make it easier to repeat your steps. You also
won't have to repeat step one if something goes wrong in step two.

An example
v.somecommand out=step1
v.anothercommand in=step1 out=step2

These can easily be turned into a script, and the layer names can be
variable derived if you want to loop through several places doing the
same process.

Alex

Pablo Carreira wrote:
> Hi, Javier!
> 
> Let me explain:
> I have a very large vector of rivers, with many CATs. Here in Brazil we have 
> a law that say that, in general, a 30m area aside of rivers should be 
> preserved with trees.
> Them I look in a high resolution orbital image and where I find a portion 
> river without trees I v.extract that portion of river, v.buffer 30m, v.patch 
> that buffer to larger common vector of buffers, g.remove the buffer and the 
> extracted vector.
> 
> It is not a batch process, but I think it could be indeed scripted, and I 
> would apreciate some help with that.
> 
> Anyway, I still want to discover is if it is possible to pipe grass commands, 
> I don't know much of programing but if I'm not wrong piping mean that the 
> output of an application becomes the input of another application.
> 
> Achim told me not to pipe, but I am a obstinate kind of person (I think most 
> of us are)...
> 
> ops, I almost forgot, in the company I use windows (becouse I can't use linux 
> there), but at home I use Ubuntu linux.
> 
> Regards
> 
> Pablo Torres Carreira
> 
> 
> 
> 
> From: fjgarciapri...@hotmail.com
> To: pablotcarre...@hotmail.com
> Subject: RE: [GRASS-user] Piping Commands
> Date: Mon, 23 Nov 2009 19:28:29 +
> 
> 
> 
> 
> 
> 
> 
> 
> Hi, I'm not an expert or anything else, but maybe I can help you.
> 
> You are using Windows or Linux?
> 
> Sometime ago I have to upload a large number of images (Ortophotos) to grass 
> (more than the command aloud me and I use a Linux script to do that performed 
> over the basis of processing each group of files and then perform another 
> operations on the result.
> 
> Is this what are you looking for?
> 
> From: pablotcarre...@hotmail.com
> To: grass-user@lists.osgeo.org
> Date: Mon, 23 Nov 2009 12:56:46 -0200
> Subject: [GRASS-user] Piping Commands
> 
> 
> 
> 
> 
> 
> 
> 
> 
> Hi, 
> 
> I want to do the following:
> v.extract some lines from a vector and them
> v.buffer those lines
> 
> How do I pipe these commands?
> 
> 
> 
> 
> 
> Pablo Torres Carreira
> 
> 
> 
> Chegou o Windows 7: Incrivelmente simples! Clique e conheça.  
>   
> Windows Live:  Friends get your Flickr, Yelp, and Digg updates when they 
> e-mail you.
> _
> Novo site do Windows Live: Novidades, dicas dos produtos e muito mais. 
> Conheça!
> http://www.windowslive.com.br/?ocid=WindowsLive09_MSN_Hotmail_Tagline_out09
> 
>
___
grass-user mailing list
grass-user@lists.osgeo.org
http://lists.osgeo.org/mailman/listinfo/grass-user

RE: [GRASS-user] Piping Commands

[Rich Shepard]

On Mon, 23 Nov 2009, Pablo Carreira wrote:


It is not a batch process, but I think it could be indeed scripted, ...

Achim told me not to pipe, ...


Pablo,

  Rather than fighting upriver, go with the flow.

  Write a script that uses temporary files. First process writes to that
file, and the second process reads from it. When you are finished you can
have the script delete the temporary file(s). The end result is the same as
piping, but it works how complex, multi-module models are built and run
within GRASS.

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

Re: [GRASS-user] Old GRASS GIS compilation

[Glynn Clements]

Firman Hadi wrote:

> I am compiling GRASS 4.3 on Debian Woody.
> My aim is to use r.le.dist function. I have tried to compile it but there
> is one error that I can't solve. The error message is like below:
> 
> /home/firman/grass4.3src/src/libes/vask
>  make -f OBJ/make.rules
> 
> rm -f OBJ/V_call.o
> gcc   -I/home/firman/grass4.3src/src/include   -DUSE_TERMIO -c V_call.c
> V_call.c: In function 'V_call':
> V_call.c:355: warning: comparison between pointer and integer
> V_call.c: At top level:
> V_call.c:430: error: static declaration of 'fmt' follows non-static
> declaration
> V_call.c:171: error: previous implicit declaration of 'fmt' was here
> V_call.c:449: error: static declaration of 'centered' follows non-static
> declaration
> V_call.c:234: error: previous implicit declaration of 'centered' was here
> make: *** [OBJ/V_call.o] Error 1
> GISGEN failure at STEP: src/libes/vask
> 
> By this email, I would like to ask help from anyone who have experience
> installing GRASS 4.3.

You can fix these specific errors by moving the definitions of fmt()
and centered() from the bottom of the file to the top.

But you may encounter many more such errors trying to compile GRASS
4.3 with a modern C compiler. You could try using an older compiler,
but that may not work with the headers from modern libraries.

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

Re: [GRASS-user] Old GRASS GIS compilation

[Dylan Beaudette]

On Monday 23 November 2009, Glynn Clements wrote:
> Firman Hadi wrote:
> > I am compiling GRASS 4.3 on Debian Woody.
> > My aim is to use r.le.dist function. I have tried to compile it but there
> > is one error that I can't solve. The error message is like below:
> >
> > /home/firman/grass4.3src/src/libes/vask
> >  make -f OBJ/make.rules
> >
> > rm -f OBJ/V_call.o
> > gcc   -I/home/firman/grass4.3src/src/include   -DUSE_TERMIO -c V_call.c
> > V_call.c: In function 'V_call':
> > V_call.c:355: warning: comparison between pointer and integer
> > V_call.c: At top level:
> > V_call.c:430: error: static declaration of 'fmt' follows non-static
> > declaration
> > V_call.c:171: error: previous implicit declaration of 'fmt' was here
> > V_call.c:449: error: static declaration of 'centered' follows non-static
> > declaration
> > V_call.c:234: error: previous implicit declaration of 'centered' was here
> > make: *** [OBJ/V_call.o] Error 1
> > GISGEN failure at STEP: src/libes/vask
> >
> > By this email, I would like to ask help from anyone who have experience
> > installing GRASS 4.3.
>
> You can fix these specific errors by moving the definitions of fmt()
> and centered() from the bottom of the file to the top.
>
> But you may encounter many more such errors trying to compile GRASS
> 4.3 with a modern C compiler. You could try using an older compiler,
> but that may not work with the headers from modern libraries.


Hi,
Any reason not to use the current version of GRASS?

Dylan

-- 
Dylan Beaudette
Soil Resource Laboratory
http://casoilresource.lawr.ucdavis.edu/
University of California at Davis
530.754.7341
___
grass-user mailing list
grass-user@lists.osgeo.org
http://lists.osgeo.org/mailman/listinfo/grass-user

Re: [GRASS-user] g.region Options Yield Different Values

[Glynn Clements]

Rich Shepard wrote:

>That's true. I've obviously erred when creating new locations and
> importing .shp and .e00 files.
> 
>I've used the wxPython GUI to define new locations and the mapset. When I
> specify choosing the projection from the list and select LCC I'm asked to
> enter the relevant projection parameters. I then specify NAD83 as the datum
> and the Washington-Oregon adjustment. Finally, I'm asked if I want to enter
> the default region's extent and resolution. This is where I enter the
> lat/long coordinates I have.
> 
>Perhaps it's this last step that I have wrong.

Yes. If you create a projected location, you need to enter projected
coordinates.

> Should I decline to enter
> the extents and resolution and leave those to the information embedded in
> the imput data?

I think that you may need to enter something here, but it doesn't have
to be correct, and you can change it later.

The current region is usually ignored when importing maps. Vector maps
will contain cartographic coordinates. Raster maps will either have
its bounds stored in metadata, or you will get clearly bogus bounds
and have to set the correct bounds manually with r.region afterwards.

If you have imported a vector map or a raster map with the correct
bounds, you can set the current and/or default region from the map.

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

Re: [GRASS-user] i.pca with multiple input bands

[Νίκος Αλεξανδρής]

On Sun, 2009-11-22 at 07:38 -0800, pdu wrote:
> Hi  there
> 
> Is there any way to specify multiple input bands to i.pca? I'm working with
> aviris data which has 224 bands and it'd be a disaster having to put each
> band like
> 
> i.pca input=ba...@permanent,ba...@permanent, ..., band...@permanent
> I tried input=ba...@permanent and it doesn't work.

> Thanks a lot.


I was messing around for quite some time with i.pca but never tried with
so many "bands". Does this really work for you (224 bands)?

Thanks, Nikos

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

RE: [GRASS-user] Piping Commands

[Glynn Clements]

Pablo Carreira wrote:

> In general, is it possible to connect the "out=" to the "in=" of vector 
> commands in GRASS

No; you can't "pipe" GRASS maps.

GRASS maps consist of multiple files, and they often need to be read
and/or written non-sequentially. Neither of these are possible with a
pipe.

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

Re: [GRASS-user] g.region Options Yield Different Values

[Rich Shepard]

On Mon, 23 Nov 2009, Glynn Clements wrote:


Yes. If you create a projected location, you need to enter projected
coordinates.


Glynn,

  And if they're available only as lat/long geographic coordinates ...?


I think that you may need to enter something here, but it doesn't have to
be correct, and you can change it later.


  OK.


The current region is usually ignored when importing maps. Vector maps
will contain cartographic coordinates. Raster maps will either have its
bounds stored in metadata, or you will get clearly bogus bounds and have
to set the correct bounds manually with r.region afterwards.


  I know that the .prj and .e00 files have the projection information and
cartographic coordinates so I use those. It's creating a new, empty location
that seems to be giving me the problems.


If you have imported a vector map or a raster map with the correct bounds,
you can set the current and/or default region from the map.


  True. And then I can identify the values, too.

Much appreciated,

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

[GRASS-user] Data Source: Interbase/Firebird

[Rich Shepard]

  The USGS National Hydrologic Database is a great resource for those of us
who need these data. The state files are huge (by comparison to most with
which I've worked): Oregon is 425M and Nevada is 330M. Unfortunately
(perhaps fortunately) they're not in ESRI or MS format, but Borland's
Interbase/Firebird format.

  I know that GDAL works on raster data and ogr on vector data from ESRI
.shp and MapInfo. I assume the NHD files are vector. Will v.in.ogr work on
them and is there a linux-based tool to convert the attribute data?

Rich

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

RE: [GRASS-user] Piping Commands

[Pablo Carreira]

Ok, no pipelines... 8(

here is how the script looks like:

#!/bin/sh
echo 'Enter CAT number'
read catnumber more
v.extract input=hidro output=temp_hidro list=$x
v.buffer input=temp_hidro output=temp_buffer buffer=30
v.patch input=temp_buffer output=many_buffers --overwrite --quiet
g.remove vect=temp_hidro,temp_buffer
echo 'Done!'

I have just learned shell script ;). I'm going to test it tomorow.

Thanks to everyone for the help.

Pablo Torres Carreira




> From: gl...@gclements.plus.com
> Date: Mon, 23 Nov 2009 22:54:34 +
> To: pablotcarre...@hotmail.com
> CC: grass-user@lists.osgeo.org
> Subject: RE: [GRASS-user] Piping Commands
> 
> 
> Pablo Carreira wrote:
> 
> > In general, is it possible to connect the "out=" to the "in=" of vector 
> > commands in GRASS
> 
> No; you can't "pipe" GRASS maps.
> 
> GRASS maps consist of multiple files, and they often need to be read
> and/or written non-sequentially. Neither of these are possible with a
> pipe.
> 
> -- 
> Glynn Clements 
  
_
Novo site do Windows Live: Novidades, dicas dos produtos e muito mais. Conheça!
http://www.windowslive.com.br/?ocid=WindowsLive09_MSN_Hotmail_Tagline_out09___
grass-user mailing list
grass-user@lists.osgeo.org
http://lists.osgeo.org/mailman/listinfo/grass-user

Re: [GRASS-user] Data Source: Interbase/Firebird

[John Callahan]

Reading the info on http://nhd.usgs.gov/data.html, it looks like the NHD 
is available as ESRI personal geodatabase (mdb files, which is the 
native format of the NHD data model developed by USGS) as well as 
shapefiles (which includes basic geometry and attributes, no 
relationships, etc...)


OGR can read/write shapefiles and read-only ESRI personal personal 
geodatabase (MDB) files, http://www.gdal.org/ogr/ogr_formats.html.


- John

**
John Callahan, Geospatial Application Developer
Delaware Geological Survey, University of Delaware
Email: john.calla...@udel.edu
**




Rich Shepard wrote:
  The USGS National Hydrologic Database is a great resource for those 
of us

who need these data. The state files are huge (by comparison to most with
which I've worked): Oregon is 425M and Nevada is 330M. Unfortunately
(perhaps fortunately) they're not in ESRI or MS format, but Borland's
Interbase/Firebird format.

  I know that GDAL works on raster data and ogr on vector data from ESRI
.shp and MapInfo. I assume the NHD files are vector. Will v.in.ogr 
work on

them and is there a linux-based tool to convert the attribute data?

Rich

___
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

Re: [GRASS-user] Data Source: Interbase/Firebird

[Rich Shepard]

On Mon, 23 Nov 2009, John Callahan wrote:

Reading the info on http://nhd.usgs.gov/data.html, it looks like the NHD is 
available as ESRI personal geodatabase (mdb files, which is the native format 
of the NHD data model developed by USGS) as well as shapefiles (which 
includes basic geometry and attributes, no relationships, etc...)


John,

  I'm embarrassed to write that I did not look carefully enough at the Web
site. I saw that I could order the files I wanted, but neglected to check
format assuming they'd be .e00 like most federal agency data. Mea culpa!

OGR can read/write shapefiles and read-only ESRI personal personal 
geodatabase (MDB) files, http://www.gdal.org/ogr/ogr_formats.html.


  I know. I've been using v.in.org for the vector data from the Oregon state
GIS repository.

  Thank you ... I'll re-order the data. They put the data on one of their
ftp servers. Oregon data (zipped not tarred and gzipped) is 425M; unzipped
it's 1.1G.

Much appreciated,

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

[GRASS-user] does the grass support windows 7

[maven apache]

does the grass support windows 7 system?
___
grass-user mailing list
grass-user@lists.osgeo.org
http://lists.osgeo.org/mailman/listinfo/grass-user

Re: [GRASS-user] does the grass support windows 7

[Paolo Cavallini]

maven apache ha scritto:
> does the grass support windows 7 system?

We are finding that 7 is more or less the same as vista (marketing... ;) ): 
GRASS
used from QGIS plugin crashes at about the same commands as in vista.
All the best.
-- 
Paolo Cavallini: http://www.faunalia.it/pc
___
grass-user mailing list
grass-user@lists.osgeo.org
http://lists.osgeo.org/mailman/listinfo/grass-user

Re: [GRASS-user] g.region Options Yield Different Values

[Glynn Clements]

Rich Shepard wrote:

> > Yes. If you create a projected location, you need to enter projected
> > coordinates.
> 
> Glynn,
> 
>And if they're available only as lat/long geographic coordinates ...?

Create the location anyhow and set the region from a map once you've
imported one.

If for some reason you need suitable bounds before you've imported any
maps you can use use "m.proj -i ..." to convert lat-lon coordinates to
projected coordinates (but projecting the corners of a lat-lon
rectangle may be an underestimate; e.g. for LCC, the southernmost
point will typically be somewhere in the middle of the southern edge,
not one of the corners).

For a more accurate estimate, create a lat-lon location, set the
region, use v.in.region and v.split to create a sampled boundary,
re-project that to the LCC location with v.proj, and use
"g.region vect=..." to set the default region from the projected
boundary.

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