Re: [GRASS-user] r.fillnulls error
On Mon, Mar 1, 2010 at 8:13 AM, maven apache apachemav...@gmail.com wrote: HI: When I use the r.fillnulls command to resample a map: r.fillnulls input=aa output=ab tension=40 smooth=0.1 It show me that: taking too long to find points for interpolation--please change the region to area where your points are. Continuing calculations... What does it mean? It is a friendly warning to tell you that you *may* be out of area of interest. But in the case of r.fillnulls it is likely that you have large NULL areas, so you can ignore the message. Be sure to validate the result of course (as usual). Markus ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user
Re: [GRASS-user] r.fillnulls error
maven apache wrote: When I use the r.fillnulls command to resample a map: r.fillnulls input=aa output=ab tension=40 smooth=0.1 It show me that: taking too long to find points for interpolation--please change the region to area where your points are. Continuing calculations... What does it mean? It means you should read the fabulous NOTES section of the r.fillnulls man page. ;-) (the warning is harmless in this context, ignore it) Hamish ps::devs- use r.to.vect -zb in the r.fillnulls script? v.surf.rst doesn't need them and it could save a lot of memory etc if lidar data with many holes is thrown at it. ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user
Re: [GRASS-user] r.fillnulls error
2010/3/1 Hamish hamis...@yahoo.com maven apache wrote: When I use the r.fillnulls command to resample a map: r.fillnulls input=aa output=ab tension=40 smooth=0.1 It show me that: taking too long to find points for interpolation--please change the region to area where your points are. Continuing calculations... What does it mean? It means you should read the fabulous NOTES section of the r.fillnulls man page. ;-) (the warning is harmless in this context, ignore it) Hamish ps::devs- use r.to.vect -zb in the r.fillnulls script? v.surf.rst doesn't need them and it could save a lot of memory etc if lidar data with many holes is thrown at it. Oh, thanks, I see it. However what perplexing me these days is to call grass in command line model in windows platform,forward a message: HI: I want to run grass in a command line model, however when I write a bat,it works just for 6.3 version , and some command can not work for example I can run g.version ,rather than r.fillnulls. The grass variable bat is : --- @echo off rem ***Environment variables*** rem Set GRASS Installation Directory Variable set GRASSDIR=F:\anzhuang\GRASS rem set GISBASE=D:\GrassDataBase set GISRC=C:\Documents and Settings\Administrator\. grassrc6 rem Set Path to utilities (libraries and bynaries) used by GRASS set PATH=%GRASSDIR%\msys\bin;%PATH% set PATH=%GRASSDIR%\extrabin;%GRASSDIR%\extralib;%GRASSDIR%\extrabin;%PATH% set PATH=%GRASSDIR%\tcl-tk\bin;%GRASSDIR%\sqlite\bin;%GRASSDIR%\gpsbabel;%PATH% set PATH=%PATH%;%GRASSDIR%\bin;%GRASSDIR%\lib;%GRASSDIR%\scripts;%GRASSDIR%\etc;%GRASSDIR%\sqlite\lib;%GRASSDIR%\sqlite\bin; rem Path to the proj files (notably the epsg projection list) set GRASS_PROJSHARE=%GRASSDIR%\proj -- And I write a work bat : - r.in.gdal input=D:\data\hdf5\new\OMI-Aura_L3-OMTO3e_2009m0729_v003-2009m0731t020127.tif output=aa location=al g.gisenv set=LOCATION_NAME=al g.region rast=aa r.null map=aa setnull=-1.267e+30--1.268e+30 r.fillnulls input=aa output=ab tension=40 smooth=0.1 r.out.gdal input=ab format=GTiff output=D:/OMI-Aura_L3-OMTO3e_2009m0729_v003-2009m0731t020127.tif - It does not work,. I wonder why? In fact what I am trying to do is to call grass in my web application(like the web processing service),now I want to handle the hdf5 data,the following is my idea: 1) use gdalinfo to get the information of the he5,(the gdal model in grass 6.3 does not support this command, the 6.4 support) 2) change the subdataset to geotiff, use gdal_translate 3) import the tiff to grass and create a new location according the tif 4) log in to the new location use g.gisenv 5) set some sepical cell data to null by using r.null 6) Interpolation the null data using r.fillnulls. 7) export the map to geotiff. So anyone can give some advise? ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user
Re: [GRASS-user] Is there a grassflyer.lyx file available?
On Sun, 2010-02-28 at 03:54 +0100, Νίκος Αλεξανδρής wrote: Is there anybody using LyX who already has a grassflyer.lyx or similar? I have prepared a .lyx version out of grassflyer.tex (en). Should I put it in the wiki? Nikos ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user
Re: [GRASS-user] Is there a grassflyer.lyx file available?
2010/3/1 Νίκος Αλεξανδρής nikos.alexand...@felis.uni-freiburg.de: On Sun, 2010-02-28 at 03:54 +0100, Νίκος Αλεξανδρής wrote: Is there anybody using LyX who already has a grassflyer.lyx or similar? I have prepared a .lyx version out of grassflyer.tex (en). Should I put it in the wiki? The best place would be where the flyer is (SVN, right?). Markus ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user
Re: [GRASS-user] Is there a grassflyer.lyx file available?
nikos: I have prepared a .lyx version out of grassflyer.tex (en). Should I put it in the wiki? Markus: The best place would be where the flyer is (SVN, right?). Sure! Should I send it over to you Markus? Or somebody else? Nikos ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user
[GRASS-user] Export/define region with Map Scale
Greetings At GIS Map Display window it's possible to define map scale in the bottom bar. I would like to know is: 1- Is it possible to export and plot maps (rasters/vectorials) with this map scale definition? 2- Is it possible to set region using this map scale definition in order to have by default this setting? Thanks Best Regfards Franz ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user
[GRASS-user] Move points to nearest feature in a different layer
Dear GRASS users, I'm having some problems trying to find a way to snap/move all points from one layer to the nearest node in another layer. I found that with v.distance, I can find distances to all nodes in the desired layer (the information is in the output window, although I haven't figured out yet how to manipulate this). Also, v.edit has a move and a snap option, although I believe the module expects some input like coordinates, not the name of another layer. So, I'm guessing that the solution will probably involve working with these two funcions, but I haven't figured out how Perhaps I'm just looking in the wrong place and there is a function in GRASS that performs this snap to another layer in one go. I'm sorry if this is kind of a beginner question... Thanks for any help. António ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user
Re: [GRASS-user] Move points to nearest feature in a different layer
I've done it using v.distance. It can output the x and y coords of the nearest feature (check upload option). So then I just created a new point vector with the coords given by v.distance Cheers Daniel 2010/3/1 António M. Rodrigues amrodrig...@fcsh.unl.pt: Dear GRASS users, I'm having some problems trying to find a way to snap/move all points from one layer to the nearest node in another layer. I found that with v.distance, I can find distances to all nodes in the desired layer (the information is in the output window, although I haven't figured out yet how to manipulate this). Also, v.edit has a move and a snap option, although I believe the module expects some input like coordinates, not the name of another layer. So, I'm guessing that the solution will probably involve working with these two funcions, but I haven't figured out how Perhaps I'm just looking in the wrong place and there is a function in GRASS that performs this snap to another layer in one go. I'm sorry if this is kind of a beginner question... Thanks for any help. António ___ 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] Move points to nearest feature in a different layer
Yes, that works. thanks. On Mar 1, 2010, at 11:59 AM, Daniel Victoria wrote: I've done it using v.distance. It can output the x and y coords of the nearest feature (check upload option). So then I just created a new point vector with the coords given by v.distance Cheers Daniel 2010/3/1 António M. Rodrigues amrodrig...@fcsh.unl.pt: Dear GRASS users, I'm having some problems trying to find a way to snap/move all points from one layer to the nearest node in another layer. I found that with v.distance, I can find distances to all nodes in the desired layer (the information is in the output window, although I haven't figured out yet how to manipulate this). Also, v.edit has a move and a snap option, although I believe the module expects some input like coordinates, not the name of another layer. So, I'm guessing that the solution will probably involve working with these two funcions, but I haven't figured out how Perhaps I'm just looking in the wrong place and there is a function in GRASS that performs this snap to another layer in one go. I'm sorry if this is kind of a beginner question... Thanks for any help. António ___ 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] Export/define region with Map Scale
Franz Schiller wrote: At GIS Map Display window it's possible to define map scale in the bottom bar. note that it is probably wrong. in order for it to be right the exact distance across your monitor would have to be known as Scale is the ratio of distance on map to distance on ground. If the distance on the screen (in ground units) isn't known, you can't do that calculation. Typically the system (be it X11 or MS Windows) has a rough idea of the dpi, but it's probably +/- 30% and not very reliable. Some software like the GIMP have a calibration tool where you can make it a bit better (until you change resolution settings or use a different monitor). but most software just give you a very rough approximate and hope you don't notice/care. IMO that scale in the map display window is misleading and should be disabled. (just because lots of other gis software pretend to support this doesn't mean we should too; borrowing the monitor calibration tool idea from GIMP is another option) I would like to know is:1- Is it possible to export and plot maps (rasters/vectorials) with this map scale definition? yes. the ps.map module knows its output plots will be on paper, and it knows how big that piece of paper is. therefore it can correctly calculate and use the scale for PostScript (and so PDF format too) hardcopies. 2- Is it possible to set region using this map scale definition in order to have by default this setting? not directly, but you could write a script wrapping g.region and d.info which after calibration could keep you somewhat close to your desired scale. That's probably more trouble than it's worth though. Hamish ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user
Re: [GRASS-user] Is there a grassflyer.lyx file available?
2010/3/1 Νίκος Αλεξανδρής nikos.alexand...@felis.uni-freiburg.de: nikos: I have prepared a .lyx version out of grassflyer.tex (en). Should I put it in the wiki? Markus: The best place would be where the flyer is (SVN, right?). Sure! Should I send it over to you Markus? Or somebody else? You can send it to me. Markus ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user
Re: [GRASS-user] POW() in SQLite?
On Mon, 1 Mar 2010, Hamish wrote: I am trying to use POW(x,n) with SQLite but it fails: Hamish, Mathematical functions such as POW() are not part of SQL. When you look at what appear to be mathematical functions, e.g., ABS(). MIN(), MAX(), they are actually string manipulations. SQLite stores most values as text strings so this is not surprising. I'd manipulate those values first, then pass them to SQLite. Rich ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user
Re: [GRASS-user] r.fillnulls error
2010/3/1 maven apache apachemav...@gmail.com 2010/3/1 Hamish hamis...@yahoo.com maven apache wrote: When I use the r.fillnulls command to resample a map: r.fillnulls input=aa output=ab tension=40 smooth=0.1 It show me that: taking too long to find points for interpolation--please change the region to area where your points are. Continuing calculations... What does it mean? It means you should read the fabulous NOTES section of the r.fillnulls man page. ;-) (the warning is harmless in this context, ignore it) Hamish ps::devs- use r.to.vect -zb in the r.fillnulls script? v.surf.rst doesn't need them and it could save a lot of memory etc if lidar data with many holes is thrown at it. Oh, thanks, I see it. However what perplexing me these days is to call grass in command line model in windows platform,forward a message: HI: I want to run grass in a command line model, however when I write a bat,it works just for 6.3 version , and some command can not work for example I can run g.version ,rather than r.fillnulls. The grass variable bat is : --- @echo off rem ***Environment variables*** rem Set GRASS Installation Directory Variable set GRASSDIR=F:\anzhuang\GRASS rem set GISBASE=D:\GrassDataBase set GISRC=C:\Documents and Settings\Administrator\. grassrc6 rem Set Path to utilities (libraries and bynaries) used by GRASS set PATH=%GRASSDIR%\msys\bin;%PATH% set PATH=%GRASSDIR%\extrabin;%GRASSDIR%\extralib;%GRASSDIR%\extrabin;%PATH% set PATH=%GRASSDIR%\tcl-tk\bin;%GRASSDIR%\sqlite\bin;%GRASSDIR%\gpsbabel;%PATH% set PATH=%PATH%;%GRASSDIR%\bin;%GRASSDIR%\lib;%GRASSDIR%\scripts;%GRASSDIR%\etc;%GRASSDIR%\sqlite\lib;%GRASSDIR%\sqlite\bin; rem Path to the proj files (notably the epsg projection list) set GRASS_PROJSHARE=%GRASSDIR%\proj -- And I write a work bat : - r.in.gdal input=D:\data\hdf5\new\OMI-Aura_L3-OMTO3e_2009m0729_v003-2009m0731t020127.tif output=aa location=al g.gisenv set=LOCATION_NAME=al g.region rast=aa r.null map=aa setnull=-1.267e+30--1.268e+30 r.fillnulls input=aa output=ab tension=40 smooth=0.1 r.out.gdal input=ab format=GTiff output=D:/OMI-Aura_L3-OMTO3e_2009m0729_v003-2009m0731t020127.tif - It does not work,. I wonder why? In fact what I am trying to do is to call grass in my web application(like the web processing service),now I want to handle the hdf5 data,the following is my idea: 1) use gdalinfo to get the information of the he5,(the gdal model in grass 6.3 does not support this command, the 6.4 support) 2) change the subdataset to geotiff, use gdal_translate 3) import the tiff to grass and create a new location according the tif 4) log in to the new location use g.gisenv 5) set some sepical cell data to null by using r.null 6) Interpolation the null data using r.fillnulls. 7) export the map to geotiff. So anyone can give some advise? I am in hurry, anyone can help me? ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user
[GRASS-user] problem with spgrass6
Hi list, I am using spgrass6 in a grass6.4-R2.9.2 on Ubuntu 9.10 without problems / //x=readRAST6(c(raster.a,raster.b))/ but when i run the same line in grass65svn it gives an error with some rasters(DCELL)... (and other are imported without problems(DCELL too)) /CPL ERROR 4: `/media/COMP/grassdb//hardtke/Bisigato/.tmp/rs7/ndvi_sum' does not exist in the file system, and is not recognised as a supported dataset name. Error en .local(.Object, ...) : `/media/COMP/grassdb//hardtke/Bisigato/.tmp/rs7/ndvi_sum' does not exist in the file system, and is not recognised as a supported dataset name. /Any idea? Is it a grass65 bug? or my fault? Leo Hardtke ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user
[GRASS-user] Terrain Curve Values
Just to confirm: the curve output from r.slope.aspect is positive when the curve is convex and negative when the curve is concave. Correct? There's nothing on the man page about this. Perhaps a sentence could be added. Rich ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user
Re: [GRASS-user] POW() in SQLite?
Hamish wrote: I am trying to use POW(x,n) with SQLite but it fails: G65 v.db.select mapname column=POW(length,3) DBMI-SQLite driver error: Error in sqlite3_prepare():SELECT POW(length,3) FROM mapname no such function: POW Nope; not in SQLite: http://www.sqlite.org/lang_corefunc.html ERROR: Unable to open select cursor Similar problem if I try x^n: ... unrecognized token: ^ SQLite doesn't have a power operator: http://www.sqlite.org/lang_expr.html#binaryops this random webhit suggests that it should be possible: http://osdir.com/ml/sqlite-users/2009-03/msg00248.html That says: * Mathematical SQL Functions - The following mathematical SQL functions are available in addition to the SQLite default: ACOS(), ASIN(), ATAN(), ATAN(), ATAN2(), CEIL(), CEILING(), COS(), COT(), DEGREES(), EXP(), FLOOR(), LN(), LOG(), LOG(), LOG2(), LOG10(), MOD(), PI(), POW(), RADIANS(), SIGN(), SIN(), SQRT(), TAN(), TRUNCATE(). Note: in addition to the SQLite default. IOW, those functions are added by SQLiteSpy, and are not part of SQLite. -- Glynn Clements gl...@gclements.plus.com ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user
[GRASS-user] Features Extraction
Hello I'm ab.derdar I'm new to the Grass-user list, although i've used GRASS GIS for past 4 months i'm interrested in : - road extraction from satellite images - remote sensing - geophysical application of GIS well i've read as much as i can from this list archive about my topics of interest but seems like apart from zero-crossing edge detection (i.zc) there is not much feature detection modules in GRASS i have made a handful of scripts (i'm not good yet with C to make a module for them) that revolves about features extraction and i would be happy to share them : - some Gaussians (and derivatives) convolution masks generators ( to be used with mfilter) - edge extractors scipts (Sobel,Canny ...) - some template based feature exctrations (Duda Operator ..) i'd hope to post those scripts but i don't know how or where to post them if someone is interesed in this topic pls contact me asap i'm warning you i'm that i'm not an expert in Geomatic i'm just an average IT profesionnal so don't expect an accurate implementations of said operators and algorithms ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user
Re: [GRASS-user] problem with spgrass6
Hi Leo, these type of questions are best sent to the grass-st...@list.osgeo.org [1] list (cc'd here), which discusses GRASS-GIS, R spgrass6 issues. What does sessionInfo() give? Best Regards, Daniel. [1] http://www.osgeo.org/mailman/listinfo/grass-stats Leo Hardtke wrote: Hi list, I am using spgrass6 in a grass6.4-R2.9.2 on Ubuntu 9.10 without problems / //x=readRAST6(c(raster.a,raster.b))/ but when i run the same line in grass65svn it gives an error with some rasters(DCELL)... (and other are imported without problems(DCELL too)) /CPL ERROR 4: `/media/COMP/grassdb//hardtke/Bisigato/.tmp/rs7/ndvi_sum' does not exist in the file system, and is not recognised as a supported dataset name. Error en .local(.Object, ...) : `/media/COMP/grassdb//hardtke/Bisigato/.tmp/rs7/ndvi_sum' does not exist in the file system, and is not recognised as a supported dataset name. /Any idea? Is it a grass65 bug? or my fault? Leo Hardtke ___ 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] Display of Terrain Perameter Maps
The output maps from r.slope.aspect. r.resamp.rst, and even r.topidx the output maps have values for each cell rather than trends. This is very difficult to interpret by non-technical viewers. I've seen maps of slopes, profile and planar curvature, and topographic wetness indices that are more smoothly colored. If one wants to determine the curvature of slopes adjacent to a specific cell at its base it would be easier to see this if the values do not change with each cell. When we stand at the foot of a hill we tend to not focus on each detail of the terrain but the overall slope, profile curvature, and planar curvature. It's this type of overview that I'm trying to display. Is there a way to get smoother output maps that are visually more informative to non-specialists? Rich ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user
[GRASS-user] R-GDAL incompatible libraries
Greetings I have just installed R, spgrass6, rgdal into my Ubuntun installation. While I was running an example from GRASS book I got this error precip30n - readVECT6(precip_30ynormals, ignore.stderr=TRUE) ERROR: Incompatible library version for module. You need to rebuild GRASS or untangle multiple installations. Where can I find an identification of the correct libraries versions. Thank you Kim ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user
Re: [GRASS-user] R-GDAL incompatible libraries
Kim Besson: I have just installed R, spgrass6, rgdal into my Ubuntun installation. While I was running an example from GRASS book I got this error precip30n - readVECT6(precip_30ynormals, ignore.stderr=TRUE) ERROR: Incompatible library version for module. You need to rebuild GRASS or untangle multiple installations. Where can I find an identification of the correct libraries versions. Thank you Kim Hi Kim! Did you install grass from source? Did you install R, spgrass6 from source? - If you have installed grass from source, then try to rebuild grass (after having installed R and the rest). - If you have grass from the ubuntu repositories, the you should probably stick with R and the rest from the repositories as well. If the latter is the case and it does not work, then there is a problem. Tell us a bit more... Nikos ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user
[GRASS-user] r.param.scale Question
I've run r.param.scale multiple times to extract the various terrain parameters from both 10m and 25m elevation maps. I set the region (and resolution) to that of the input elevation map. All maps but the tangential (planar) curvature map are properly rendered. This one comes out a solid color. I don't see what I'm doing incorrectly on this one map. Suggestions appreciated. Rich ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user
RE: [GRASS-user] r.param.scale Question
I've run r.param.scale multiple times to extract the various terrain parameters from both 10m and 25m elevation maps. I set the region (and resolution) to that of the input elevation map. All maps but the tangential (planar) curvature map are properly rendered. This one comes out a solid color. I don't see what I'm doing incorrectly on this one map. Suggestions appreciated. Rich It's possible there are outliers in your raster map that are skewing the color table; what is the result of r.info -r? Check to see the range of the raster makes sense. -- Eric ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user
[GRASS-user] Raster attributes
Is it possible to link a raster image to a set of attributes stored in a database using GRASS? I am working with the Harmonized World Soil Database which consist of an image and a database that holds a set of soil attributes. These attributes can be linked/joined to the image since the pixel values are IDs that appear in the database. I would like to output a series of images, each representing a different attribute such as soil unit, texture, soil depth. Is this possible in GRASS? All the best, Ned ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user
RE: [GRASS-user] r.param.scale Question
Eric wrote: It's possible there are outliers in your raster map that are skewing the color table; whenever that is the case 'r.colors -e' can help. Hamish ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user
RE: [GRASS-user] r.param.scale Question
On Mon, 1 Mar 2010, Patton, Eric wrote: It's possible there are outliers in your raster map that are skewing the color table; what is the result of r.info -r? Check to see the range of the raster makes sense. Eric, That's not the situation here: no outliers. I successfully ran this at 5 and 10m resolutions, and once at 25m resolution. The region was incorrect in this last case and after properly setting the region I could no longer produce the multicolored map I did before. I'm not using that map since I _must_ get this report done and found a different module (r.resamp.rst) that also generates tangential curves. Thanks, Rich ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user
[GRASS-user] ps.map issue
This is for my information only now as I don't need an answer right away. I have a ps.map script that uses a raster map, two boundary vline maps, and a vpoints map. There are also clauses for text, colortable, and vlegend. It make for a pretty map. However, if I add one more vline map (the stream channels), the colortable does not display and ghostscript complains about a problem. I cannot figure out why. I changed the sequence in which elements were processed and drawn, but this made no difference. Has anyone seen anything like this before? Rich ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user
Re: [GRASS-user] problem with spgrass6
Leo Hardtke wrote: I am using spgrass6 in a grass6.4-R2.9.2 on Ubuntu 9.10 without problems / //x=readRAST6(c(raster.a,raster.b))/ but when i run the same line in grass65svn it gives an error with some rasters(DCELL)... (and other are imported without problems(DCELL too)) /CPL ERROR 4: `/media/COMP/grassdb//hardtke/Bisigato/.tmp/rs7/ndvi_sum' does not exist in the file system, and is not recognised as a supported dataset name. the // in the middle of the path makes me suspect that there is a $VARIABLE there which is supposed to be set but is empty for some reason. ? (does that file exist?) Hamish ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user
Re: [GRASS-user] Raster attributes
Ned Horning wrote: Is it possible to link a raster image to a set of attributes stored in a database using GRASS? categorical (CELL) raster maps have very basic support via the r.category module; the category value of the cell can hold the key column ID. for something more advanced see: http://grass.osgeo.org/wiki/GRASS_SoC_Ideas_2008#Time_series_management_in_GRASS http://grass.osgeo.org/wiki/Time_series_in_GRASS http://trac.osgeo.org/grass/browser/grass-addons/raster/r.rast4d I am working with the Harmonized World Soil Database which consist of an image and a database that holds a set of soil attributes. These attributes can be linked/joined to the image since the pixel values are IDs that appear in the database. I would like to output a series of images, each representing a different attribute such as soil unit, texture, soil depth. Is this possible in GRASS? I think it would be quite possible to write a nice script that does this with r.reclass + r.category. It would be a nice addition to wiki addons :) I think it could be done without the complication of an extra SQL layer, but whatever works. probably the easiest idea though is to simply run 'r.to.vect -v feature=area' and then v.db.connect. You will be limited by the number of cells though, probably if the raster is bigger than 1600x1600 you might run into memory problems. (depends on how many repeated cell values/ clumps are in the map) good luck, Hamish ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user
Re: [GRASS-user] Terrain Curve Values
On Mon, Mar 1, 2010 at 4:19 PM, Rich Shepard rshep...@appl-ecosys.com wrote: Just to confirm: the curve output from r.slope.aspect is positive when the curve is convex and negative when the curve is concave. Correct? Here the answer from Helena (not on this list): On Tue, Mar 2, 2010 at 1:41 AM, Helena Mitasova hmit...@unity.ncsu.edu wrote: yes - convex is positive and concave is negative. It is the same as in v.surf.rst. There's nothing on the man page about this. Perhaps a sentence could be added. I agree. Markus ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user