Thanks Paulo, your suggestion of combining r.neigbors and r.patch worked.
From: Paulo van Breugel p.vanbreu...@gmail.com
Subject: Re: [GRASS-user] 'Growing' a raster on the edges only
Message-ID: 516ba41a.4060...@gmail.com
Content-Type: text/plain; charset=UTF-8; format=flowed
Try
Hi,
I have rainfall measurements from 6 gauges that I want to interpolate
to an areal value (a 'surface'), so that I can compare the the
interpolated gauge values to a satellite rainfall estimate that covers
a grid cell of 28kmx28km. Two of the gauges are outside, but close to
the border of the
long to compute for 2 years of
data. You should be able to fairly easily output your desired time series to
a text file.
Best,
Nick
On Mon, Aug 15, 2011 at 3:00 PM, Hanlie Pretorius
hanlie.pretor...@gmail.com wrote:
Hi,
I have rainfall measurements from 6 gauges that I want
, Markus Neteler nete...@osgeo.org:
On Wed, May 4, 2011 at 10:48 AM, Hanlie Pretorius
hanlie.pretor...@gmail.com wrote:
Hi,
I'd like to ask a question related to another question I was asking
recently about converting contour vectors to rasters and creating a
DEM from the raster contours using
:
On 03/05/11 15:47, Micha Silver wrote:
On 05/03/2011 04:34 PM, Moritz Lennert wrote:
On 03/05/11 10:20, Micha Silver wrote:
On 05/03/2011 09:11 AM, Hanlie Pretorius wrote:
Hi Maris,
Thanks, your suggestion worked. With the region's resolution set to
5m, I got some output.
Hi Maris
reason why I though the layer was empty. I
adjusted the colour table an re-rendered the map with no effect.
Today, however, it works and I can see it.
Sorry about the confusion.
2011/5/4, Moritz Lennert mlenn...@club.worldonline.be:
On 04/05/11 08:42, Hanlie Pretorius wrote:
Hi Moritz
Hi,
I'd like to ask a question related to another question I was asking
recently about converting contour vectors to rasters and creating a
DEM from the raster contours using r.surf.contour.
My supplied vector contour layer has 20m intervals between contour
lines. When I convert this to raster,
Hi Maris,
Thanks, your suggestion worked. With the region's resolution set to
5m, I got some output.
Regards
Hanlie
2011/5/2, Maris Nartiss maris@gmail.com:
g.region vect=contours_20m
Also 20x20m cells wouldn't be too large?!?
Maris.
2011/5/2, Hanlie Pretorius hanlie.pretor
Hi,
I'm working in GRASS 6.4.1 on Windows XP.
I want to use r.surf.contour to convert contour lines into a DEM.
However, I first need to convert my vector contour lines into raster
contour lines. My region is set to the contours_20m vector layer, with
a 20m resolution.
v.to.rast --overwrite
Hi,
I'm working in GRASS 6.4.0 and using v.buffer applied to an extents
rectangle. For some reason v.buffer is producing one chamfered corner
in the buffered vector.
The command I used was:
v.buffer -s -c --overwrite input=c83a_extents@PERMANENT
output=c83a_extents_5km_buffer@20m_DEM
Ok, thanks.
The other three corners were rounded, but I eliminated the rounding
with the -s option.
Hanlie
2011/4/20, Markus Metz markus.metz.gisw...@googlemail.com:
On Wed, Apr 20, 2011 at 2:21 PM, Hanlie Pretorius
hanlie.pretor...@gmail.com wrote:
Hi,
I'm working in GRASS 6.4.0 and using
Hi,
I'm working with GRASS 6.4.0 under QGIS in Windows XP. I'm posting
here because I got no responses on the GRASS Windows email list.
I'm trying to install the add-on r.inund.fluv
(https://svn.osgeo.org/grass/grass-addons/raster/r.inund.fluv/).
The readme file seems aimed at Linux users, so I
Hi Richard,
Your suggestion worked in principle, but the whole operation needs
more work. I ended up with the same node names in the start and end
node columns, so I think a script is needed to identify the start and
end nodes based on the direction of the flow of the river. An
elevation column
Hi,
I'm working in GRASS 6.4.0 on Windows XP.
I'm creating a hydrological model and am using GIS to prepare the data.
I've split my river layer into about 100 segments and I've digitised
the start and end point of each river segment on a points layer.
For the points layer, I calculate a 'name'
to_column=Jname
I don't get an error but no updates either.
I've found another (proprietary) program that can do this, but I was
curious to see if it's possible in GRASS.
Regards
Hanlie
2010/12/9, Micha Silver mi...@arava.co.il:
On 09/12/2010 16:38, Hanlie Pretorius wrote:
Hi,
I'm working
Hi,
I'm running GRASS through the QGIS (v1.0.5) interface on Windows XP.
The GRASS folder in the QGIS installation says 'grass-6.4.0svn'.
I have a polygon layer with 22 polygons and a raster slope layer. When
I run v.rast.stats with these two layers, two of the polygons have
NULL values in the
Hi,
I need to calculate the maximum Euclidean distance between lines in a
vector layer.
In ArcGIS I use Spatial Analyst to achieve this; it has a Euclidean
distance tool that takes a line vector layer and produces a raster
output, either as table or as an image. The output contains several
Hi,
I've imported a shapefile with 14 polygons using v.in.ogr:
-
v.in.ogr
dsn=/media/0847147784/data/dwaf/liebenbergsvlei/myne/sa_tm_29_deg_east/C83A/c83a_14_subcatchments.shp
output=test
Datum unknown not recognised by GRASS and no parameters found
Projection of input dataset and current
2010/11/3, Micha Silver mi...@arava.co.il:
On 11/03/2010 04:55 PM, Hanlie Pretorius wrote:
Hi,
I need to add some attributes from a soil polygon layer to a catchment
polygon layer using area weighting. In other words, the contribution
of the soil layer needs to be weighted by the amount
Hi,
I need to add some attributes from a soil polygon layer to a catchment
polygon layer using area weighting. In other words, the contribution
of the soil layer needs to be weighted by the amount with which its
polygons overlap the catchment polygons.
For example, suppose I want to weight soil
2010/10/19, Glynn Clements gl...@gclements.plus.com:
Hanlie Pretorius wrote:
So, I have two problems here: 1) I can't get r.reclass to output an
FCELL raster even though its help file suggest that only the input
file needs to be a CELL raster.
A reclass table maps integers to integers. You
Hi,
I'm trying to reclassify a raster with these categories:
|#|description
|-
|1|Sandy Loam
|2|Sandy Clay Loam
|*|no data
To this:
|#|description
Hi,
I've created subcatchments using r.watershed and I would like to
convert these to polygons. My result at the moment contains too many
polygons - 194. The original raster had only 9 areas, which should all
have been contiguous seeing that they came from r.watershed.
Here is the sequence of
2010/9/29, Markus Metz markus.metz.gisw...@googlemail.com:
In GRASS, you have full control of the computational region,
consequently you must be able to explain (to the client/supervisor)
every aspect of the computational region. Why is the horizontal
resolution not only slightly off from
Hi,
I'm working in GRASS 6.4.0RC6 and Ubuntu 10.04.
I'm trying to write a Python script to project hundreds of files from
one location to another. Since r.proj -l outputs text that will have
to be parsed, I decided to rather follow the method explained at this
webpage
2010/9/13, Sylvain Maillard sylvain.maill...@gmail.com:
rasters in mapset 2000_02_february_trmm:
-
g.list type=rast mapset=2000_02_february_trmm
--
no raster files available in mapset 2000_02_february_trmm
Hi,
I'm writing my first GRASS Python script and am having trouble adding
options to commands that I want to run. For example, the following
works:
-
grass.run_command(r.in.xyz, input=directory+'/'+file,
output=rainraster, fs=',', x=2, y=1)
-
but this doesn't:
-
Thanks :-)
2010/9/3, Martin Landa landa.mar...@gmail.com:
2010/9/3 Hanlie Pretorius hanlie.pretor...@gmail.com:
grass.run_command(r.in.xyz, input=directory+'/'+file,
output=rainraster, fs=',', x=2, y=1)
better `input = os.path.join(directory, file)`
but this doesn't
Hi,
I'm using GRASS6.4SVN on Windows XP.
I want to import TRMM precipitation data into a GRASS projected
(Transverse Mercator) location. I have previously successfully import
this data to a GCS (WGS84) location using the procedure at
http://grass.osgeo.org/wiki/Import_XYZ.
My problem is the
Hi
I would like to compile GRASS with extra debugging information and
apparently one needs to enable 'debugging symbols'.
Which of the configuration options would add debugging symbols to the
compilation? I can't see any in the list produced by
-
./configure --help
-
that obviously
2010/8/19, razmjoo...@faunalia.co.uk razmjoo...@faunalia.co.uk:
In Ubuntu, you can right click on the file and tick the box to make it
executable
Right Click Properties Permissions
Allow executing files as program
Cheers
Sab
Thanks, tried this and it workedfor the script in the
Hi,
I'm working in GRASS 6.4RC6 on Ubuntu 10.04 and I'm trying to learn
how to script GRASS using Python.
I have taken the example at
http://grass.osgeo.org/programming6/pythonlib.html that checks if a
vector is 3D and saved it to my home folder. When I try to run it
though, GRASS can't find the
2010/8/19, Nikos Alexandris nikos.alexand...@felis.uni-freiburg.de:
Better not to use sudo (as a general advice). The script needs to be made
executable for your user-name, i.e.:
sudo chown hanlie:hanlie v.check3d
Ran 'sudo chown hanlie:hanlie v.check3d' from the
/usr/lib/grass64/scripts
Hi,
I'm using GRASS64 RC6 on Linux 10.04 (64-bit) and I'm trying to run
r.sim.water on a catchment about 5000 km2 in size.
I started with a very fine resolution in my region (00:00:00.18) and
got the following:
-
GRASS 6.4.0RC6 (world_wgs84):~ r.sim.water
elevin=lieb_dem_25m_clipped
Thanks for the reply. I'll try that if I don't manage with the Linux version.
However, previously when I worked with the SVN version, I was advised
to rather use the stable version.
2010/8/10, Helmut Kudrnovsky hel...@web.de:
Hi,
I'm working with WinGRASS 6.4RC6 on Win XP
[...]
please try
Had a look at the article. Interesting process.
However, the TRMM rainfall data is already spatially distributed,
albeit at a coarse resolution (0.25 degrees).
2010/8/11, razmjoo...@faunalia.co.uk razmjoo...@faunalia.co.uk:
Hanlie
You might find this useful:
http://www.surfaces.co.il/?p=578
My apologies for my ignorance, but how do I get this to run in Linux?
Am using Ubuntu 10.04 and GRASS 6.4RC6.
I entered the following into the GRASS text window and am testing it
with a list of two files (I removed the region settings because the
files all have the same region, which I have
2010/8/11, Vincent Bain b...@toraval.fr:
Hi, short intrusion in this current thread, without knowing exactly what
it is related to : I only mean to suggest you to close the while loop...
This prompt means the shell expects the end of your statement ; here
you should type done to close the
2010/8/11, Micha Silver mi...@arava.co.il:
Samber's method will surely work, but you might more simply try as follows:
# start grass in a location which matches the txt file data
# change to the directory where your txt files are, then do
for infile in *.txt; do outrast=`basename ${infile}
are, then do
for infile in *.txt; do outrast=`basename ${infile} .txt`; r.in.xyz
in=${infile} out=${outrast} fs=,; done
Hanlie Pretorius wrote:
This works nicely, except that my input files have names like:
3B42.000201.0.6.nc.lieb.txt
3B42.000201.3.6.nc.lieb.txt
and I want the output rasters
Hi,
I'm working with WinGRASS 6.4RC6 on Win XP.
I'm following the steps at http://grass.osgeo.org/wiki/Import_XYZ to
import hundreds files containing TRMM rainfall data.
When I try to follow the instruction:
-
Just amend above procedure to use wildcards. Change in above example
all
Hi,
I thought I created a clipped DEM by rasterising an irregular vector
boundary and then using a mask to create the clipped DEM from a
rectangular DEM.
For example:
-
# convert vector border to raster
v.to.rast input=border output=border use=cat type=area layer=1
value=1 rows=4096
is
raster1
raster2
.
Hope that helps
Saber
On Tue, 2010-08-10 at 14:59 +0200, Hanlie Pretorius wrote:
Hi,
I'm working with WinGRASS 6.4RC6 on Win XP.
I'm following the steps at http://grass.osgeo.org/wiki/Import_XYZ to
import hundreds files containing TRMM rainfall data.
When I
2010/8/10, Nikos Alexandris nikos.alexand...@felis.uni-freiburg.de:
Hanlie Pretorius wrote:
I thought I created a clipped DEM by rasterising an irregular vector
boundary and then using a mask to create the clipped DEM from a
rectangular DEM.
For example:
-
#first, check and set your
Hi,
I have a MODIS land cover image that contains NULL values according to r.univar:
-
r.univar map=modis_land_cover_type1_2001_prim...@permanent
total null and non-null cells: 41884
total null cells: 17145
-
So, I tried to fill it with r.fillnulls and I get the following:
-
2010/8/4, Hamish hamis...@yahoo.com:
Hanlie wrote:
| Rows: 11957
| Columns: 6277
...
I have tried to see the NULLS by displaying everything in
white, except for NULLS displayed in red, but I didn't see
anything.
Does anyone know how I can find out if there are actually
Hi,
I work with GRASS 6.4RC6 and would like to look at the code for the
r.sim.water module. I have been looking at
http://trac.osgeo.org/grass/browser/grass/branches/releasebranch_6_4/raster/,
but I can't find it there. Can someone please tell me where I the code
is kept?
Thanks
Hanlie
Hi Everyone,
I'm working in GRASS 6.4RC6 on Win XP and am trying to reproject TRMM
rainfall data from its native GCS (WGS84) to a local PCS. I'm using
the PCS because my other important data, specifically DEMs, are
provided in it. TRMM data come in 0.25°x0.25° squares.
My workfow is as follows:
2010/6/5, Hamish hamis...@yahoo.com:
Hanlie wrote:
At this point, g.region reports 1146474 cells in the region, while I
have 1146370 lines of coordinates in my file.
...
So it looks like there are about 100 coordinates missing from the ASCII
ASCII file.
0.01% ..
Maybe holes in the
2010/6/7, Hamish hamis...@yahoo.com:
Hanlie Pretorius wrote:
I don't think it's this bug because this bug discards only one line of
data. I don't get any data in because the number of coordinate pairs
in the file is less than the number of cells in the defined region.
Markus Metz:
Weird
2010/6/7, Markus Metz markus.metz.gisw...@googlemail.com:
g.region -p
north: -49312.5
south: -74587.5
west: -3015862.5
east: -2987512.5
nsres: 25
ewres: 25
rows: 1011
cols: 1134
cells: 1146474
From the region settings you posted
2010/6/4, Micha Silver mi...@arava.co.il:
On 06/04/2010 04:20 PM, Hanlie Pretorius wrote:
Hi,
I've been following the procedure at
http://grass.osgeo.org/wiki/Import_XYZ to import DEMs into GRASS.
It has worked for two DEMS, but I have a problem with the third one at
the step where one
Hi,
I was just wondering how the GRASS code displays coordinates in a
geographic coordinate system on a computer screen. Do they get
projected in some way? Or is there another trick to it?
Thanks
Hanlie
___
grass-user mailing list
Hi,
I'm using Win XP and GRASS 6.4RC6 to import a TRMM netCDF file using
the command:
-
r.in.gdal -o input =c:/data/3B42.000202.0.6.nc output=test band=precipitation
-
and am getting the error:
-
G_set_window(): Illegal latitude for North
-
gdalinfo for the file outputs:
-
Hi,
I've been following the procedure at
http://grass.osgeo.org/wiki/Import_XYZ to import DEMs into GRASS.
It has worked for two DEMS, but I have a problem with the third one at
the step where one verifies that the number of rows in the ASCII file
corresponds to the number of cells in the
...@gclements.plus.com:
Hanlie Pretorius wrote:
I'm using GRASS on two machines in different locations (one Win XP,
the other Ubuntu 10.04), so I would like to keep my grassdata folder
on an external drive. However, when I try to run grass in Ubuntu using
the external hard drive (NTFS file system), I get
Hi,
I ran r.sim.water on a DEM using the process in the manual for
generating the rainfall, manning and infiltration maps. The height in
the DEM vary from 1429.79 m to 1740.20 m. I opted to output the
rainfall depth and the error maps.
The ranges of both these output maps (as reported by
Hi,
I'm using GRASS 6.4 on two machines in different locations (one Win XP,
the other Ubuntu 10.04), so I would like to keep my grassdata folder
on an external drive. However, when I try to run grass in Ubuntu using
the external hard drive (NTFS file system), I get the following error
as soon as
Hi,
I'm using QGIS to export shapefiles from WGS84 GCS to a local
projected coordinate system so that I can import the projected
shapefile into GRASS.
The PROJ4 definition that I use to export the shapefile (and that I
used to create the GRASS location) is:
+proj=tmerc +lat_0=0 +lon_0=29 +k=1
1/2 a pixel in each direction to ensure the points are in
the center of the pixel. At 25m spacing, your region should be set as:
g.region n=1012.5 s=487.5 w=487.5 e=1012.5 res=25
- Jamie
On Thu, May 13, 2010 at 9:02 AM, Hanlie Pretorius
hanlie.pretor...@gmail.com wrote:
Hi Jamie,
I
Hi Micha,
I tried your suggestion after setting the region to 20m instead of the
raster DEM's 25m.:
v.surf.rst input=dem_2628cc_...@c83 layer=0
elev=dem_2628cc_rst_elev tension=40. segmax=40 npmin=120
dmin=9.998022 dmax=49.990111 zm ult=1.0
This worked, but the differences between the raster
Hi,
I've obtained DEMs in text files with columns X, Y and Z at 25m
spacing. The first three entries in the text file are:
-
X,Y,Z
99550,2.9883e+06,1473.47
99550,2.98828e+06,1473.57
99550,2.98825e+06,1473.63
-
To me it seems the easiest way to import these is the following:
1.
pixel in all four directions. Also look at using r.in.xyz, it
works directly on this type of data.
Jamie
On Thu, May 13, 2010 at 7:26 AM, Hanlie Pretorius
hanlie.pretor...@gmail.com wrote:
Hi,
I've obtained DEMs in text files with columns X, Y and Z at 25m
spacing. The first three entries
Try http://www.gadm.org.
Date: Wed, 12 May 2010 06:23:34 +0200
From: Jos? Miguel Barrios jmbarri...@gmail.com
Subject: [GRASS-user] maps France and Germany
To: grass-user@lists.osgeo.org
Message-ID:
aanlktimmcisrcoamohjfuhsbdi2l7wdqejzfsr2ra...@mail.gmail.com
Content-Type:
Hi Glynn,
This worked, thanks.
2010/5/11, Glynn Clements gl...@gclements.plus.com:
Hanlie Pretorius wrote:
I just pasted the code you sent at the end of the d.correlate script
That won't work. The code is a patch, indicating changes to be made
to the file. I've attached a complete
Hi,
I'm running GRASS 6.4.0RC6 on Windows XP SP1.
Working my way through the book Open Source GIS a GRASS GIS Approach
3rd ed. by Markus Neteler and Helena Mitasova, I have noticed problems
with some d.* commands.
For example, when I type d.histogram into the GUI layer manager
command line, I
Hi,
I tried to create contours from a DEM with the following information:
| Type of Map: raster Number of Categories: 255
| Data Type:FCELL
| Rows: 9856
| Columns: 4844
| Total Cells: 47742464
|Projection: Transverse Mercator
|
Greetings Herbivores,
This is the first time that I've created a location and mapset fo my
own and I'm trying to import TRMM precipitation data into GRASS using
r.in.ascii.
I have an existing CSV file that has the headers lat, long and precip.
The precipitation is a floating point number to two
Hi,
This is a repost since I got no reaction the previous time. I'm not
sure if it is because people can't help. Is it ok to contact the
authors at the bottom of the r.sim.water page if no one on this list
can help me?
I'm working on Windows XP SP3 and Grass 6.4.0RC6. When I try to run
I suppose this must be related to ticket 617
(http://trac.osgeo.org/grass/ticket/617). Guess I'll have to switch to
Linux if I want to work with it now.
2010/4/9, Hanlie Pretorius hanlie.pretor...@gmail.com:
Hi,
This is a repost since I got no reaction the previous time. I'm not
sure
Hi,
I'm working on Windows XP SP3 and Grass 6.4.0RC6. When I try to run
r.sim.water according to the example on page 163 in the book Open
Source GIS a GRASS GIS approach (3rd Ed.), I get one of those Windows
errors that want to 'tell Microsoft about this problem'.
The details of the error report
Hi,
I've just installed GRASS 6.4.0svn for Windows (Running on Windows XP
SP3) and I find a peculiar thing happening in the help files. I was
looking at the help page for the TOPMODEL module and at the top of the
page the text displays in a Slavic language. Further down on the page
some English
72 matches
Mail list logo