Re: [GRASS-user] 'Growing' a raster on the edges only (Paulo van Breugel)

Thanks Paulo, your suggestion of combining r.neigbors and r.patch worked. From: Paulo van Breugel > 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 r.neighbors function with

[GRASS-user] 'Growing' a raster on the edges only

Hi, I'm working in GRASS 6.4.2, WIndows 7. I have a raster file that, for various reasons, has shrunk to smaller than its original size and it now now longer crosses the boundaries of a vector layer to which I need to assign values from the raster. So I was wondering how I could 'grow' the raste

Re: [GRASS-user] Rainfall interpolation for a time series

his shouldn't take very 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: > &

[GRASS-user] Rainfall interpolation for a time series

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 gr

Re: [GRASS-user] DEM resolution from r.surf.contour

, Markus Neteler : > On Wed, May 4, 2011 at 10:48 AM, Hanlie Pretorius > 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 raste

[GRASS-user] DEM resolution from r.surf.contour

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, t

Re: [GRASS-user] v.to.rast gives empty raster layer

other 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 : > On 04/05/11 08:42, Hanlie Pretorius wrote: > > Hi Moritz, > > >

Re: [GRASS-user] v.to.rast gives empty raster layer

itz Lennert : > 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, >

Re: [GRASS-user] v.to.rast gives empty raster layer

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 : > g.region vect=contours_20m > > Also 20x20m cells wouldn't be too large?!? > > Maris. > > > 2011/5/2, Hanlie Preto

[GRASS-user] v.to.rast gives empty raster layer

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 inpu

Re: [GRASS-user] One chamfered corner with v.buffer on an extents rectangle

Ok, thanks. The other three corners were rounded, but I eliminated the rounding with the -s option. Hanlie 2011/4/20, Markus Metz : > On Wed, Apr 20, 2011 at 2:21 PM, Hanlie Pretorius > wrote: >> Hi, >> >> I'm working in GRASS 6.4.0 and using v.buffer applied to

[GRASS-user] One chamfered corner with v.buffer on an extents rectangle

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 distance=500

[GRASS-user] Installing r.inund.fluv on Windows

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

[GRASS-user] Assign attributes of start and enpoints to connecting lines

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 in

[GRASS-user] Re: Assign attributes of start and enpoints to connecting lines

Hi Shane, Thanks for the suggestion :-) I think you're starting where I want to end up, that is, with an attribute link between the two tables. I want to create this link based on the spatial arrangement of the lines and the points. So, the initial 'joining' of points to line needs to be based on

Re: [GRASS-user] Assign attributes of start and enpoints to connecting lines

mn=end_node 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 : > On 09/12/2010 16:38, Hanlie Pretorius wrote: >&

[GRASS-user] Assign attributes of start and enpoints to connecting lines

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'

[GRASS-user] (no subject)

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 co

[GRASS-user] Euclidean distance between lines

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 statis

[GRASS-user] Conversion from shapefile to GRASS and back

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 not recognised by GRASS and no parameters found Projection of input dataset and current location a

Re: [GRASS-user] Area weighting for vectors

2010/11/3, Micha Silver : > 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

[GRASS-user] Area weighting for vectors

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 i

Re: [GRASS-user] r.reclass with floating point result

2010/10/19, Glynn Clements : > > 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

[GRASS-user] r.reclass with floating point result

Hi, I'm trying to reclassify a raster with these categories: |#|description |- |1|Sandy Loam |2|Sandy Clay Loam |*|no data To this: |#|description |---

Re: [GRASS-user] Converting raster to polygon creates too many polygons

2010/9/29, Markus Metz : > > 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 exactly one meter, but also > dif

[GRASS-user] Converting raster to polygon creates too many polygons

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 com

Re: [GRASS-user] Script for r.proj and g.list issues

2010/9/13, Sylvain Maillard : > 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>> >> -- >>

[GRASS-user] Script for r.proj and g.list issues

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 (

Re: [GRASS-user] Adding options to grass.run_command

Thanks :-) 2010/9/3, Martin Landa : > 2010/9/3 Hanlie Pretorius : >> 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 do

[GRASS-user] Adding options to grass.run_command

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: - grass.run_command("

[GRASS-user] Raster resolution problems

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 grid

[GRASS-user] Compiling GRASS with debugging symbols

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 includ

Re: [GRASS-user] Python scripts: 'command not found'

2010/8/19, 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 /usr/lib/grass64/scrip

Re: [GRASS-user] Python scripts: 'command not found'

2010/8/19, Nikos Alexandris : > 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 directory with no error messages, b

[GRASS-user] Python scripts: 'command not found'

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

[GRASS-user] r.sim.water buffer overflow

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 dxin=li

Re: [GRASS-user] Importing multiple files with r.in.xyz

r txt files 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.n

Re: [GRASS-user] Importing multiple files with r.in.xyz

2010/8/11, Micha Silver : > 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} .txt`; r.in.xyz

Re: [GRASS-user] Importing multiple files with r.in.xyz

2010/8/11, Vincent Bain : > 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 while in

Re: [GRASS-user] Importing multiple files with r.in.xyz

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 alread

Re: [GRASS-user] Re: Importing multiple files with r.in.xyz

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 : > Hanlie > > You might find this useful: > http://www.surfaces.co.il/?p=578 > > It's in qgis, but

[GRASS-user] Re: Importing multiple files with r.in.xyz

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 : >>Hi, >> >>I'm working with WinGRASS 6.4RC6 on Win XP > [...] > > please try also th

Re: [GRASS-user] Clipping rasters

2010/8/10, Nikos Alexandris : > 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: >> - > &g

Re: [GRASS-user] Importing multiple files with r.in.xyz

n.xyz -s -g input=$line output=$line > tmpRegion > myregion= ""`head -n 1 tmpRegion`" > g.region $myregion > r.in.xyz --overwrite input=$line output=$line fs=, > done > > and your list_of_files.txt is > raster1 > raster2 > . > > Hope that helps

[GRASS-user] Clipping rasters

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

[GRASS-user] Importing multiple files with r.in.xyz

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 occurenci

[GRASS-user] Re: Filling Null values in a MODIS image

2010/8/4, Hamish : > 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 >>

[GRASS-user] Filling Null values in a MODIS image

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: - r.fillnul

[GRASS-user] Location of r.sim.water code

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 __

[GRASS-user] Issues when reprojecting from a GCS to a PCS

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:

[GRASS-user] Re: Importing Dems with r.in.xyz

2010/6/7, Markus Metz : >> 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 previously, I assume that th

Re: [GRASS-user] Re: Importing Dems with r.in.xyz

2010/6/7, Hamish : > 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 th

Re: [GRASS-user] Re: Importing Dems with r.in.xyz

2010/6/5, Hamish : > 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 th

[GRASS-user] Re: Importing Dems with r.in.xyz

2010/6/4, Micha Silver : > 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

[GRASS-user] Importing Dems with r.in.xyz

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 enlarg

[GRASS-user] Problem importing netCDF

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: - Dr

[GRASS-user] Projecting a GCS onto a computer screen

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 grass-user@lists.osgeo

[GRASS-user] Re: Using an external hard drive as data source

Thanks for all the replies. I reformatted the external drive as FAT and let Ububtu mount it itself, without editing fstab. Then I started GRASS without the sudo option and by selecting the data folder using the /media/volume_name/grassdata path. 2010/5/24, Glynn Clements : > > Hanlie Pre

[GRASS-user] Using an external hard drive as data source

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 I

[GRASS-user] r.sim.water depth and error ranges are the same

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 r.info)

[GRASS-user] Datum in v.in.ogr

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 +x

Re: [GRASS-user] dems from coordinate lists

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 raste

[GRASS-user] Re: dems from coordinate lists

egion is defined as pixel edge, you need to > expand the region 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, M

[GRASS-user] Re: dems from coordinate lists

resent cell centers so you need to expand your > region 1/2 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: >

[GRASS-user] dems from coordinate lists

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. v.in.asci

Re: [GRASS-user] d.* command problems

Hi Glynn, This worked, thanks. 2010/5/11, Glynn Clements : > > 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

[GRASS-user] Re: maps France and Germany

Try http://www.gadm.org. > Date: Wed, 12 May 2010 06:23:34 +0200 > From: Jos? Miguel Barrios > Subject: [GRASS-user] maps France and Germany > To: grass-user@lists.osgeo.org > Message-ID: > > Content-Type: text/plain; charset="iso-8859-1" > > Hi list; > > Can someone indicate me where can

[GRASS-user] d.* command problems

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

[GRASS-user] r.contour generates no lines

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 |

[GRASS-user] Problem after import using r.in.ascii

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 d

[GRASS-user] Re: r.sim.water crashes (repost)

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 : > Hi, > > This is a repost since I got no reaction the previous time. I'm not > sure if it i

[GRASS-user] r.sim.water crashes (repost)

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 r.sim.wate

[GRASS-user] r.sim.water crashes

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 repor

[GRASS-user] Language in GRASS help files

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 t