GRASS GIS (Geographic Resources Analysis Support System) is a free, open source geographical information system (GIS) capable of handling raster, topological vector, image processing, and graphic data.

learn more… | top users | synonyms

1
vote
1answer
330 views

mdenoise/r.denoise: How to

I want to use mdenoise / r.denoise but i cant get it working. I have Xubuntu 12.04 and Grass 6.4.1 GRASS tells me after installing the plugin: g.extension extension=r.denoise ...
4
votes
4answers
774 views

How to automatically plot intermediate contours with QGIS, GRASS GIS or SAGA GIS?

I'm mapping a steep mountain range. My main contour interval is 50 m and my intermediate contour interval is 10 m. In steep areas, the interval of 10 m is too small, so I need to leave ...
15
votes
3answers
3k views

Visual Impact Assessment in QGIS or GRASS

Is it possible using QGIS or GRASS to perform a Zone of Visual Impact assessment? E.g. where can a 100m high structure be seen at a viewing height of 1.8m within 30km, taking into account elevation ...
3
votes
2answers
1k views

1'' DHM: Hillshade is too rough - how to “generalize”?

I got this hillshade from my DHM (I used the hillshade-module from qgis). For my purpose there is too much relief detail (roughness). How can I get rid of this in QGIS or GRASS? Thanks! Martin
1
vote
1answer
366 views

r.sim.water G_malloc error

I've been trying to run r.sim.water on a 3 meter LiDAR-generated DEM I have, but keep getting G_malloc errors: r.sim.water --verbose elevin=meterDEM.final@Sediment_Simulation ...
1
vote
2answers
233 views

Can GRASS read Spatialite geometry fields?

I've looked at the question: How to import SpatiaLite data into GRASS? already. The links there talk about two separate tasks: Passing SQL to a Spatialite DB from GRASS and connecting to Spatialite ...
1
vote
2answers
840 views

How to calculate the average value of a raster within polygons?

I am using the GRASS plugin in QGIS. I am trying to calculate the average value of a raster dataset within a set of polygons. I have created a new mapset and have imported both the vector and ...
1
vote
1answer
97 views

mapping number of days a site is flooded from daily inundated areas maps

I have a set of daily inundated area polygons and would like combine them into one map that can give me a map showing the number of days a site remained inundated. I am unable to figure out any ...
1
vote
1answer
285 views

Why can't r.patch find the rasters?

I have elevation rasters which loaded just fine in a GRASS session. They represent adjacent areas, so I want to meld them into one raster. I ran r.patch using the GRASS GUI, and it would only accept ...
0
votes
1answer
357 views

GRASS and QGIS- “ not recognized…format” but it worked before

Prologue: I used my Ubuntu file viewing GUI to move data files around and "clean house" a bit. I started a new GRASS session and tried to add NED elevation rasters which I have successfully added ...
0
votes
1answer
49 views

Does 'r.patch' act only on “loaded” rasters?

Can I run the GRASS command 'r.patch' on raster files which are not "loaded" into the current GRASS session? That is, rasters which are simply in my folder but not in a GRASS location & mapset? ...
0
votes
1answer
136 views

GRASS importing same-name rasters fails

I have NED elevation data in 3 rasters. (That's how it downloaded from USGS NED page.) The files are named "61471006", "61471006_2", and "61471006_3". I succeeded in using the GRASS GUI and menu ...
0
votes
0answers
131 views

Making a Bridge Classification tool in GRASS GIS

I'm working on a tool for classifying bridges with GRASS. I have been using examples online but I've come to a standstill as how to actually make this tool in GRASS. I have my raster of my area and ...
0
votes
1answer
204 views

How to specify column names with v.in.ascii

I have managed to import a CSV file into a points vector file using GRASS v.in.ascii v.in.ascii --o input=test1.csv output=test2 format=point fs=, skip=1 x=2 y=3 z=0 cat=0 This is test1.csv ... ...
2
votes
2answers
587 views

r.mapcalculator and python again

I've been trying to write a script creating binary maps: a=raw_input ("enter a threshold") b=float(a) grass.run_command('r.mapcalculator', formula="**'base_map'>=b**", outfile='map'+str(a), ...
8
votes
1answer
1k views

Open software tool or QGIS plug-in to make land cover classification

Is there any tool or plugin for QGIS/GRASS to make an automatic land cover classification from satellite imagery --- based on CLC nomenclature for example? Of course, it will not be 100% sure but it ...
2
votes
1answer
152 views

Why is GRASS not thinning a terraflow raster properly?

I want to convert a raster map encoding directional flows (the output of r.terraflow) to a vector map. I'm using r.thin to process the output. It terminates successfully with 0 pixels removed during ...
5
votes
1answer
396 views

How to divide a polyline by intersecting isochrones

I have two shapefiles, the first representing a path (a polyline with a single category), the second representing isochrons from the same origin (point 0) which intersect the above mentioned path. I ...
6
votes
2answers
1k views

creating a raster of the residuals of a regression between two rasters

I am trying to create a raster of the residuals of a regression between two rasters. i.e. I would like to carry out a regression of one raster against another of the same extent and plot the residual ...
1
vote
3answers
184 views

Looking for large GML files containing vector polygon data

I have few vector data files with 500 to 700 MB GML files. But I am looking for files larger than these preferably 1 GB, 2 GB. It would be helpful if somebody can provide some pointers regarding the ...
2
votes
2answers
5k views

Viewing .las format LIDAR data in QGIS?

Is there an easy path to visualizing this LIDAR data? I have QGIS, and some USGS LIDAR data in (.las) format downloaded from http://lidar.cr.usgs.gov/ . This means I have both the (.las) and ...
1
vote
2answers
395 views

What I'm doing wrong with the command r.out.ascii in GRASS GIS?

I have raster in .bil format which I want to save in GRASS GIS with r.out.ascii. I get this result in .txt file: north: 1N south: 0 east: 1E west: 0 rows: 1 cols: 1 55537 I suppose it is something ...
3
votes
1answer
417 views

Change raster value for an area based on a class of another raster

I would like to change all raster values in raster A that fall into the locations of category/class 1 of raster B. Here, some more details. I did a calculation on a DEM resulting in different ...
5
votes
2answers
178 views

Vectorize cell boundaries in GRASS

I need to vectorize the cell boundaries of a raster map and I have forgotten how to do this. For instance, if I have a raster map of a region that is just map=1 with a cell size of 500m, I would ...
11
votes
4answers
2k views

How to connect lines in a network?

I have a shapefile of lines and another of polygons. When two or more lines intersect a polygon they should be connected (snapped), but sometimes during manual digitization this doesn't occur. I need ...
4
votes
3answers
2k views

How do you open/display and extract NetCDF files info?

I am working with NetCDF files using R & Matlab and I would like to know if the GIS community is using them. I know that ArcGIS can open & display, however most of the times the results are ...
6
votes
3answers
2k views

Dissolve or unsplit lines on common attributes in PostGIS or GRASS?

I have a road centerline shapefile that has roads split at every intersection. However, I would like to join LINESTRINGs that have a common set of attributes into a single LINESTRING. So long as the ...
5
votes
2answers
385 views

Is it possible to use the Goode Homolosine projection in QGIS?

Is it possible to use the Goode Homolosine projection in QGIS? Both QGIS and GRASS seem to have trouble with the interruptions, squashing the higher latitudes together. I tried +proj=goode ...
3
votes
2answers
168 views

How to make the layer 2 use the same vector categories from layer 1?

I'm having some doubts concerning the layer mechanism. Creating a new layer (ie a new connection between a vector and a DB table), I've imported a vector (rivers) through v.in.ogr, then I've cloned ...
1
vote
2answers
283 views

Automatically apply topology rule to only selected portions of a .shp file?

I am looking for some examples/tutorials on how to apply a topology rule to only a selected portion of a shp file using a script in ArcGIS/GRASS. Using ArcGIS seems to always require an interactive ...
8
votes
2answers
149 views

File naming for Enterprise GIS database?

I'm designing the guidelines for an Enterprise GIS database. Part of this database will include customized geoprocessing tools using a variety of software. I'm interested in knowing whether anyone ...
2
votes
1answer
239 views

How to keep (complex) color table when exporting from GRASS

I understand how to work with color tables within GRASS - but is there any way to preserve color tables when opening the raster in a different program? For example, I created a color table in GRASS ...
2
votes
1answer
236 views

Silence user warnings for GRASS-PYTHON program

I've created a program that runs externally from GRASS v6.4.1, using python 2.6. It imports an ESRI shapefile, and exports a seperate KML for every polygon in this vector file. Polygons with ...
1
vote
2answers
2k views

How to set cell size for vector to raster conversion in QGIS GRASS plugin?

I'm converting a vector (loaded in GRASS plugin in QGIS (1.7.0.-wroclaw) and would like to convert this to a raster with a defined cell size. At the same time I need a value from the attribute tabel ...
3
votes
2answers
330 views

run r.stats as a quiet module in GRASS-Python script

Has anyone else used r.stats as a quiet module in a GRASS-python script? As part of my program I'm using: grass.run_command("r.stats", flags='l', input=<file>, output=<*.csv file>, ...
5
votes
2answers
278 views

Can I split up v.in.ogr into multiple v.clean commands?

I'm trying to load a large set of polygons into GRASS with "v.in.ogr" and am running out of memory at the "Break Boundaries" stage. I'm wondering if I can break into steps all of the different things ...
7
votes
1answer
458 views

Why does GRASS eliminate polygons when using v.generalize?

I am facing a dilemma while trying to use v.generalize to simplify polygons which are similar in size and structure to US Counties or States. Often, I lose polygons from the data set after running ...
2
votes
0answers
117 views

How to fix Graphserver warning “osmdb has no edges!”?

I'm really new to GIS. I have some shapefiles data that I'm trying to import in graphserver. I'm doing: gs_osmdb_compile output3.osm output3.osmdb And the result is: WARNING: osmdb has no edges!. ...
3
votes
2answers
246 views

Hand fit curve and then apply to time series of rasters, make a map. How?

Maybe a data management problem, not too sure, but likely a challenge which someone has dealt with before. I have a time series of rasters. I can summarise the information by taking averages of each ...
4
votes
0answers
153 views

What open source solution would work best for auto extracting features from satellite images? [duplicate]

Possible Duplicate: Feature Extraction - Land cover I have been looking for an open source solution which would allow me to define a region (or collection of images) and extract features ...
1
vote
1answer
132 views

Blank Vector Selection Dropdown in QGIS GRASS Plugin

I realise there are some issues which getting vectors into the GRASS module in QGIS, but is the behaviour i am getting expected? I have loaded the vector shapefile via GRASS, closed G and reopened ...
5
votes
2answers
1k views

Quickly create contour polygons from cost points using OSS for web app

This is related to a previous question I've asked: Create filled contour map with GRASS I've created a point layer using OpenStreetMap data and pgRouting's driving_distance() function, where the ...
1
vote
2answers
93 views

v.dissolve (and other scripts) not in GRASS path

I see that I have the v.dissolve module in my grass64/scripts dir, but it is not found in my path from the GRASS command line. I have the GRASS_ADDON_PATH environment variable set: ...
6
votes
1answer
124 views

Would it be possible to construct a stand-alone grass.script python package?

I would like to migrate dependency on ArcPy to grass.script without requiring the full installation of Grass. Would it be technically possible to have a portable directory containing grass.script and ...
3
votes
2answers
726 views

How to combine polygons in a vector with GRASS?

I need to combine polygons of similar geological classes which are stored in a vector file called geology. I was looking for somethiing like "union selected features" or similar. In the GRASS Book ...
1
vote
1answer
232 views

How can GRASS export an .flt file?

I am using GRASS and connecting it to another application. The other application needs an FLT file of the raster map which I should export from GRASS. The other application only reads FLT files. Can ...
1
vote
1answer
105 views

How do I convert my coordinates to the limit of geographic region (outer edge)?

I am trying to use r.in.bin. The parameters of north, south, east, and west needs values: north=float Northern limit of geographic region (outer edge) south=float Southern limit of geographic region ...
5
votes
1answer
1k views

Create filled contour map with GRASS

Background I'm trying to make contour/isoline/isochrone maps to show what areas are within walking distance to schools for my pet project to help with the Safe Routes to School program. I followed ...
5
votes
1answer
242 views

Need to remove small areas, but only if they do not have adjoining boundary

I'm generalizing some lat/lon data in grass64 and as a precursor to the generalization, I am removing small areas from the data using 'v.clean'. The following command is working well: $ v.clean ...
1
vote
2answers
191 views

Can r.mask be called from a batch file?

I am using a shell script that is called by a batch file. The shell script contains the GRASS scripts that I am using for my program. Every time I use the r.mask script, it displays the error r.mask: ...