Richard wrote:
> Being more a scripter than a programmer, I suspect
> g.pnmcomp is beyond me ... I can't really see how to pass a
> map to the module, and can't really find any examples to help.
#spearfish example:
r.shaded.relief map=elevation.dem shade=elev.shade
g.region rast=elev.shade
d.mon
Tom,
OK, you either want to use v.patch or v.overlay, and you want to do it
via a script so you don't have to handle things manually. Well, ok, six
files isn't too many ...
Although it can show some bad habits, v.patch is probably the better
tool here ... eg (in Australia not the US):
v.p
Hamish,
Being more a scripter than a programmer, I suspect g.pnmcomp is beyond
me ... I can't really see how to pass a map to the module, and can't
really find any examples to help.
So here's what I would like (ideally) a script to do:
1. set a variable (eg state)
for i in NSW QLD VIC; do
g
Richard,
Yeah... Merge might not be the right word. I'm importing a lot of files
downloaded from the USGS site, for example, state boundaries. There are 6
separate files, but I want them all to go into the same database table/map
layer. The import vector function always wants to put each imp
Daniel wrote:
> Tried recompiling grass but that did
> not change anything. So I ran ldd
> `which r.what` and it did not show any links to a gdal
> library... This is what I got
>
> GRASS 6.5.svn (br_wgs84):~ > ldd `which r.what`
> linux-gate.so.1 =>
> (0x00f71000)
> libgrass_gis.so =>
>
Micha wrote:
> We're used to perceiving shadows below and highlights above
> what is being lit. The same holds for a hill-shaded terrain.
> When we see shadows on a slope, it seems to be south facing,
> and highlights appear north facing. So if you "place" the
> sun at 130° (SE) then you are lighti
David wrote:
> Using GRASS 6.4.0 on windows vista I type "g.extension
> extension=r.stream.order" in the command line at the bottom of
> the layer manager window. The following error is returned:
>
> "Error: svn client required. Please install Subversion first."
>
> I have looked through many pos
Hi,
I am fairly new to GRASS and have what many of you would consider a trivial
or newb problem. Regardless any help would be appreciated.
Using GRASS 6.4.0 on windows vista I type "g.extension
extension=r.stream.order" in the command line at the bottom of the layer
manager window. The following
>>>Pedro Roma:
>>>I want to develop a script that crops a raster in order to only have the
>>>pixels that are located inside a vector polygon.
>>Marcello Gorini:
>> There is probably a more elegant solution, but this will do:
>> g.region vect=your_vector
>> v.to.rast input=your_vector output
It depends on what you mean by "merge".
v.patch is the "official" tool to merge two vectors, but you could use
v.overlay (option=or).
Richard
On 24/09/10 10:56 PM, Thom DeCarlo wrote:
Richard,
Thanks! That was exactly the information I needed.
Now to figure out how to merge layers...
Thom
Pedro Roma:
> > >I want to develop a script that crops a raster in order to only have the
> > >pixels that are located inside a vector polygon.
Marcello Gorini:
> > There is probably a more elegant solution, but this will do:
> >
> > g.region vect=your_vector
> > v.to.rast input=your_vector ou
Pedro Roma:
> >I want to develop a script that crops a raster in order to only have the
> >pixels that are located inside a vector polygon.
Marcello Gorini:
> There is probably a more elegant solution, but this will do:
> g.region vect=your_vector
> v.to.rast input=your_vector output=your_vec
stephen sefick wrote:
> I have a bunch of floating point grids that I need to bring into
> GRASS. They say that they are 0.3m res where the documentation says
> they are 1m- I will ferrett this out of the data providers. My
> question is I will bring these layers into a 1m resolution mapset.
>
Kim Besson wrote:
> So your suggestion would be
> rass.mapcalc("$L=$lmin+
> (($lmax-$lmin)/($QCALmax-$QCALmin)) ($raster_file-$QCALmin)", L=tempL,
> lmin=lmin, lmax=lmax, QCALmax=QCALmax, QCALmin=QCALmin,
> raster_file=raster_file)
>
> and at the end delete L file?
My suggestion would be to repl
I have a bunch of floating point grids that I need to bring into
GRASS. They say that they are 0.3m res where the documentation says
they are 1m- I will ferrett this out of the data providers. My
question is I will bring these layers into a 1m resolution mapset.
Will they be automatically resamp
I have a bunch of files that are floating grid data. They have the
file extension .flt and there is a header file .hdr. r.in.gdal reads
in the .flt file fine, but errors out with the hdr. The gdalinfo
spits out a res=0.3, but the attached metadata says that the
resolution is 1m. Is there a way
>Pedro Roma:
>
>I want to develop a script that crops a raster in order to only have the
>pixels that are located inside a vector polygon.
There is probably a more elegant solution, but this will do:
g.region vect=your_vector
v.to.rast input=your_vector output=your_vector_rasterized use=val
r.m
> Date: Thu, 23 Sep 2010 11:06:21 +0300
> From: Maris Nartiss
> Subject: Re: [GRASS-user] Re: QGIS GRASS plugin shadedrelief issue
> To: nunosousa84
>
> Sorry for You.
> You do made backup of important data before messing around with QGIS, right?
>
> You can try to launch GRASS and rebuild vect
So your suggestion would be
rass.mapcalc("$L=$lmin+
(($lmax-$lmin)/($QCALmax-$QCALmin)) ($raster_file-$QCALmin)", L=tempL,
lmin=lmin, lmax=lmax, QCALmax=QCALmax, QCALmin=QCALmin,
raster_file=raster_file)
and at the end delete L file?
Thanks
Kim
2010/9/24 Glynn Clements
>
> Kim Besson wrote:
>
>
Greetings
I want to develop a script that crops a raster in order to only have the
pixels that are located inside a vector polygon. Which functionality shall I
use?
Thanks
pedro
___
grass-user mailing list
grass-user@lists.osgeo.org
http://lists.osgeo.or
On Fri, Jul 30, 2010 at 7:20 PM, Daniel Victoria
wrote:
> Hi all,
>
> Reading on the Atmospheric Correction Wiki [1] and the i.atcorr man
> page [2] I noticed some inconsistencies. In the manual, the exemple
> uses the lsat7_2002_40 image from the NC dataset directlly but, on the
> wiki tutorial,
On Fri, Jul 16, 2010 at 4:08 AM, Glynn Clements
wrote:
>
> Jarek Jasiewicz wrote:
>
>> In GRASS65 I create dir in raster or vector directory, program module
>> in I want to create, next make, sudo make install, ready
>>
>> in grass 7 not.
>> after make install: "no rules to make install. Stop". O
Richard,
Thanks! That was exactly the information I needed.
Now to figure out how to merge layers...
Thom
Sent via DROID on Verizon Wireless
-Original message-
From: Richard Chirgwin
To: grass-user@lists.osgeo.org
Sent: 2010 Sep, Thu, 23 07:20:51 GMT+00:00
Subject: Re: [GRASS-user]
On 09/24/2010 12:45 AM, Bryan Keith wrote:
Orange boundaries are duplicates - overlapping lines - that GRASS
topology doesn't allow.
You should be able to get rid of these with:
"v.clean in=... out=... tool=break,bpol,rmdupl"
But I've never had good luck with this.
The alternative that mi
Tried recompiling grass but that did not change anything. So I ran ldd
`which r.what` and it did not show any links to a gdal library... This
is what I got
GRASS 6.5.svn (br_wgs84):~ > ldd `which r.what`
linux-gate.so.1 => (0x00f71000)
libgrass_gis.so => /usr/local/grass-6.5.svn/l
Kim Besson wrote:
> I'm doing a Python script to calibrate a raster image. Since I need to do 3
> different calculations I'm using a tempfile (tempL)
>
> This is the code:
>
> tempL = grass.read_command("g.tempfile", pid = "1")
> grass.mapcalc("$L=$lmin+ (($lmax-$lmin)/($QCALmax-$QCALm
Greetings
I'm doing a Python script to calibrate a raster image. Since I need to do 3
different calculations I'm using a tempfile (tempL)
This is the code:
tempL = grass.read_command("g.tempfile", pid = "1")
grass.mapcalc("$L=$lmin+ (($lmax-$lmin)/($QCALmax-$QCALmin))*
($raster_file-$QC
It worked
Thanks
Helena
On Thu, Sep 23, 2010 at 3:49 PM, Glynn Clements wrote:
>
> Helena Herrera wrote:
>
> > I have sent an email a couple of days ago reporting an error with
> g.mapset
> > on Windows. Now I have also a problem with g.mapset on linux (a different
> > one)
> > When I do:
> > g.m
Thanks for all the explanations. Now I know better. Below [1], I tried
to bring all relevant snippets together. Wouldn't it make sense to add
some of these details to g.region.html, or so?
Hermann
[1] Glynn and Hamish gave these helpful explanations:
> to sum up:
>
> * r.in.gdal differs fr
29 matches
Mail list logo
