[GRASS-user] (no subject)
Favor desuscribirme de esta lista, gracias RODOLFO A. DEL PRADO QUIROZ Arquitecto Universidad Nacional de Colombia MSc en Medio Ambiente y Desarrollo Profesional Universitario Alcaldia de Medellín Celular 3007794381 ___ grass-user mailing list grass-user@lists.osgeo.org https://lists.osgeo.org/mailman/listinfo/grass-user
Re: [GRASS-user] (no subject)
Hi, On Thu, Dec 6, 2018 at 9:38 AM SHIVAM RAI wrote: > > Hi, > please don't delay. > please answer the question that i have asked recently you still did not send which operating system you use. In another email to me you stated that it is GRASS GIS 7.4.0. > when executing "d.tile" and "d.what.rast" in console or cmd Please elaborate "in console or cmd". > error: not implemented in the wxGUI. try to add as a layer command. > thanks Markus ___ grass-user mailing list grass-user@lists.osgeo.org https://lists.osgeo.org/mailman/listinfo/grass-user
Re: [GRASS-user] (no subject)
Dear Shivam rai There is no module called "d.tile", So it won't work. What did you try with d.what.rast? Please post the command you tried !! Regards Sajid On Thu, Dec 6, 2018 at 9:38 AM SHIVAM RAI wrote: > Hi, > please don't delay. > please answer the question that i have asked recently > when executing "d.tile" and "d.what.rast" in console or cmd > error: not implemented in the wxGUI. try to add as a layer command. > thanks > ___ > grass-user mailing list > grass-user@lists.osgeo.org > https://lists.osgeo.org/mailman/listinfo/grass-user ___ grass-user mailing list grass-user@lists.osgeo.org https://lists.osgeo.org/mailman/listinfo/grass-user
Re: [GRASS-user] (no subject)
Hi Please Reply Me On Thu, Dec 6, 2018 at 2:07 PM SHIVAM RAI wrote: > Hi, > please don't delay. > please answer the question that i have asked recently > when executing "d.tile" and "d.what.rast" in console or cmd > error: not implemented in the wxGUI. try to add as a layer command. > thanks > ___ grass-user mailing list grass-user@lists.osgeo.org https://lists.osgeo.org/mailman/listinfo/grass-user
[GRASS-user] (no subject)
Hi, please don't delay. please answer the question that i have asked recently when executing "d.tile" and "d.what.rast" in console or cmd error: not implemented in the wxGUI. try to add as a layer command. thanks ___ grass-user mailing list grass-user@lists.osgeo.org https://lists.osgeo.org/mailman/listinfo/grass-user
[GRASS-user] (no subject)
I’ve tried unsuccessfully to add a satellite to the i.atcorr tool. I followed the steps listed on the page, but I cannot get the actual integration of the tool to work. I have all the needed information to add the tool. It’s for Planet’s Dove satellites. Any help is greatly appreciated! (they have a few types of satellites, and for now at least I only need the one that starts with 0f in the excel document. Also I have the band information for each of the 4 bands Planet Blue band B1 (455nm - 415nm) Planet Green band B2 (500nm - 590nm) Planet Red band B3 (590nm - 670nm) Planet NIR band B4 (780nm - 860nm) Thanks for the help, Gwen relative_spectral_response_median_per_ps_sat_group.xlsx Description: relative_spectral_response_median_per_ps_sat_group.xlsx ___ grass-user mailing list grass-user@lists.osgeo.org https://lists.osgeo.org/mailman/listinfo/grass-user
Re: [GRASS-user] (no subject)
Julita Dudek: Hello, Hello Julita, I use i.atcorr modul, but I have one question: - in the "run_i_atcorr(radiance_flag, input_band, input_range, elevation, visibility, parameters, output, output_range)" party, what is the input_range output_range, you take this from MTL.txt or from Radiance min and max values, or you fixed it 1:255 (because of DN)? I think the reason for the input and output range(s), being set by default to [0,255], is the kind of "need" at the time of the module's creation. Most programs and applications were somewhat and somehow "restricted" to 8-bit. In other words, yes, "because of DN"s ranging in [0,255] for various sensors. This is, in my view, something we need to rethink about, given the many different sensors and the need to update programs to be able to deal with more than 8-bit data. Input range If the input data are: - Spectral Radiance values, it greatly varies and depends on the sensor and the band. Since you mention "MTL.txt", which implies you are handling Landsat data, in my attempt to automate the process for Landsat imagery, (see https://github.com/NikosAlexandris/i.landsat.atcorr/blob/master/i.landsat.atcorr.py#L419) the actual minimum and maximum values for each input band are set. I think this is no harm. Though having some nominal input range might be better (?). ( Looking at the definition of Spectral Radiance (which derives from radiance), the units are: W / m^2 / str / nm, that is: Watts per Square meter per Steradian per Nanometer. Obviously, the "dynamic" range of spectral radiance values, depend on a sensor's radiometric characteristics. Worth mention (?), a sensor's capacity may be appreciated via its bitness: how many bits is the sensor capable of (digitising the incident electromagnetic energy)? 6, 8, 10, 11 or even 12? And, practically, for users, what is the bitness of the distributed products (usually 8 or 16-bit)? ) - reflectance values, then it is [0,1]. Output range The output of i.atcorr, is *always* Spectral Reflectance. The option to rescale is, im my understanding, a "helper" to get the data ready for subsequent processing with modules or else that expect values ranging in [0,255]. In this sense, the default [0,255] is "correct" as long as the actual input is equal or less than 8-bit and/or an 8-bit output is wanted _and_ there is less care to loose radiometric detail in case the input data are of more than 8-bits. In general, I keep the output in [0,1], as is the module's intended behaviour (see https://github.com/NikosAlexandris/i.landsat.atcorr/blob/master/i.landsat.atcorr.py#L435). The question is "what is the (end-)user's need?" In case the input is of >8-bit and the subsequent purpose is more than simple visualisations, then it is better to think about how to proceed and keep the output being (corrected) spectral radiance ranging in [0,1]. Rescale accordingly, later on depending on if and what exactly is needed. Or directly (re)set the output to something larger than [0,255] via the rescale option. Directed to @grass-dev I think it is good to update the description of the module's input and output elements as well as rethink the output scale. I suggest at the very least something like: --%<--- svn diff main.cpp Index: main.cpp === --- main.cpp(revision 69474) +++ main.cpp(working copy) @@ -463,6 +463,7 @@ struct Options opts; opts.iimg = G_define_standard_option(G_OPT_R_INPUT); +opts.iimg->description = _("Spectral radiance (or reflectance with -r) input image"); opts.iscl = G_define_option(); opts.iscl->key = "range"; @@ -490,12 +491,13 @@ opts.icnd->description = _("Name of input text file with 6S parameters"); opts.oimg = G_define_standard_option(G_OPT_R_OUTPUT); +opts.oimg->description = _("Spectral reflectance output image"); opts.oscl = G_define_option(); opts.oscl->key = "rescale"; opts.oscl->type = TYPE_INTEGER; opts.oscl->key_desc = "min,max"; -opts.oscl->answer = "0,255"; +opts.oscl->answer = "0,1"; opts.oscl->required = NO; opts.oscl->description = _("Rescale output raster map"); opts.oscl->guisection = _("Output"); --->%-- I suggest to further customize descriptions of "input" and "output" options, i.e. to modify the "Name of..." (which derives from the standard set of options), for most imagery modules. Why not being explicit about what the expected input is? Nikos Thank you very much for your help. Best regards -- Julita Dudek UMR 5600 Environnement Ville Société, CNRS Centre de recherche en géographie et aménagement, LCRE, Université Lyon 3 0649 07 25 75 ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user
[GRASS-user] (no subject)
Hello, I use i.atcorr modul, but I have one question: - in the "run_i_atcorr(radiance_flag, input_band, input_range, elevation, visibility, parameters, output, output_range)" party, what is the input_range output_range, you take this from MTL.txt or from Radiance min and max values, or you fixed it 1:255 (because of DN)? Thank you very much for your help. Best regards -- Julita Dudek UMR 5600 Environnement Ville Société, CNRS Centre de recherche en géographie et aménagement, LCRE, Université Lyon 3 0649 07 25 75 ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user
Re: [GRASS-user] (no subject)
Thanks Anna, it worked! On Fri, Aug 12, 2016 at 11:31 AM, Anna Petrášováwrote: > On Fri, Aug 12, 2016 at 1:05 PM, alassane toure wrote: > > Group, > > I am running a grass-python script but when i replace a variable argument > > within the command the script does not work > > > > The following works... > > grass.run_command('r.in.gdal', > > input='/home/alassane/Data/SanbornCD/tile/newconstr_mask.tif', > > output='construction') > > > > But not this... > > command="'r.in.gdal', input='"+sys.argv[5]+"', output='construction'" > > print 'r.in.gdal command argument... = %s\n ' % command > > grass.run_command(command) > > > > print 'r.in.gdal command argument... = %s\n ' % command > > r.in.gdal command argument... = 'r.in.gdal', > > input='/home/alassane/Data/SanbornCD/tile/newconstr_mask.tif', > > output='construction' > > > > > > The following also works... > > command ="/home/alassane/Programs/gdal/demclass1"+" "+infile1+" > "+infile2+" > > "+sys.argv[3]+" "+outfile > > os.system(command) > > > > Your help is appreciated. > > Perhaps this? > > grass.run_command('r.in.gdal', input=sys.argv[5], output='construction') > > Best, > Anna > > > > > Regards, > > Alassane > > > > ___ > > 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] (no subject)
On Fri, Aug 12, 2016 at 1:05 PM, alassane tourewrote: > Group, > I am running a grass-python script but when i replace a variable argument > within the command the script does not work > > The following works... > grass.run_command('r.in.gdal', > input='/home/alassane/Data/SanbornCD/tile/newconstr_mask.tif', > output='construction') > > But not this... > command="'r.in.gdal', input='"+sys.argv[5]+"', output='construction'" > print 'r.in.gdal command argument... = %s\n ' % command > grass.run_command(command) > > print 'r.in.gdal command argument... = %s\n ' % command > r.in.gdal command argument... = 'r.in.gdal', > input='/home/alassane/Data/SanbornCD/tile/newconstr_mask.tif', > output='construction' > > > The following also works... > command ="/home/alassane/Programs/gdal/demclass1"+" "+infile1+" "+infile2+" > "+sys.argv[3]+" "+outfile > os.system(command) > > Your help is appreciated. Perhaps this? grass.run_command('r.in.gdal', input=sys.argv[5], output='construction') Best, Anna > > Regards, > Alassane > > ___ > 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] (no subject)
Group, I am running a grass-python script but when i replace a variable argument within the command the script does not work The following works... grass.run_command('r.in.gdal', input='/home/alassane/Data/ SanbornCD/tile/newconstr_mask.tif', output='construction') But not this... command="'r.in.gdal', input='"+sys.argv[5]+"', output='construction'" print 'r.in.gdal command argument... = %s\n ' % command grass.run_command(command) print 'r.in.gdal command argument... = %s\n ' % command r.in.gdal command argument... = 'r.in.gdal', input='/home/alassane/Data/ SanbornCD/tile/newconstr_mask.tif', output='construction' The following also works... command ="/home/alassane/Programs/gdal/demclass1"+" "+infile1+" "+infile2+" "+sys.argv[3]+" "+outfile os.system(command) Your help is appreciated. Regards, Alassane ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user
[GRASS-user] (no subject)
Hello, I've tried to use r.surf.idw on a large raster. It works well in 6.4.4. I've moved to 7.1 SVN, where everything is marvellous, at least what I use, except r.surf.idw which creates the same output as the impot was. I've also svn-ed and compiled 7.0.1, where r.surf.idw did the same. Namely: r.surf.idw --overwrite --verbose input=output_temp1 output=output_temp2 npoints=7 when finished, output_temp2 is the same as output_temp1 and no interpolation seems to be done. One more big difference: 6.4.4 uses less than 2G ram and finishes in a few minutes. In 7.0 and 7.1 it eats all 6GB RAM, 4 GB SWAP and runs much-much longer, understandably... Should I stay with 6.4.4? thnaks Robert ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user
[GRASS-user] (no subject)
Hello i want to involve in grass gevelopment.i am a S4 B.Tech CS student. which version is prefered, grass gis from ubuntu repo or grass latest 7.0. How to build from source. ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user
Re: [GRASS-user] (no subject)
How to build from source. have a look in the wiki: http://grasswiki.osgeo.org/wiki/Compile_and_Install - best regards Helmut -- View this message in context: http://osgeo-org.1560.n6.nabble.com/no-subject-tp5033202p5033203.html Sent from the Grass - Users mailing list archive at Nabble.com. ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user
Re: [GRASS-user] (no subject)
On 2/9/13, Afsal Iqubal afsal.iqu...@gmail.com wrote: Hello i want to involve in grass gevelopment.i am a S4 B.Tech CS student. which version is prefered, grass gis from ubuntu repo or grass latest 7.0. How to build from source. Please download the GRASS 7 source code from SVN trunk which is the development version. Your contributions will be appreciated. Markus ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user
[GRASS-user] (no subject)
http://www.ganamossucaso.com/dfgqp.html___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user
[GRASS-user] (no subject)
Hi, I need to obtain this format from shapefile, how do I? ESRI ASCII Grid Format (Copied from the ArcWorkstation 8.3 Help File) The ASCII file must consist of header information containing a set of keywords, followed by cell values in row-major order. The file format is NCOLS xxx NROWS xxx XLLCENTER xxx | XLLCORNER xxx YLLCENTER xxx | YLLCORNER xxx CELLSIZE xxx {NODATA_VALUE xxx} row 1 row 2 ... row n where xxx is a number, and the keyword nodata_value is optional and defaults to -. Row 1 of the data is at the top of the grid, row 2 is just under row 1 and so on. For example: ncols 386 nrows 286 xllcorner -128.66338 yllcorner 13.7502065 cellsize 0.2 NODATA_value - - - -123 -123 -123 - - - - - ... - - -123 -123 -123 - - - - - ... - - -117 -117 -117 -119 -119 -119 -119 -119 - ... ... The nodata_value is the value in the ASCII file to be assigned to those cells whose true value is unknown. Cell values should be delimited by spaces. No carriage returns are necessary at the end of each row in the grid. The number of columns in the header is used to determine when a new row begins. The number of cell values must be equal to the number of rows times the number of columns. I'm working GRASS 6.4.2 (2012) OSX 10.6.8 Many thanks, Jose C. Guerrero -- Dr. José Carlos Guerrero Laboratorio de Desarrollo Sustentable y Gestión Ambiental del Territorio (LDSGAT) Instituto de Ecología y Ciencias Ambientales (IECA) Facultad de Ciencias Universidad de la República Montevideo-URUGUAY ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user
Re: [GRASS-user] (no subject)
José wrote: I need to obtain this format from shapefile, how do I? ESRI ASCII Grid Format r.out.arc but shapefiles are vector data model and Arc ASCII grids are raster data model, so you would have to translate between the two with v.to.rast or v.out.ascii+r.in.xyz, or interpolate with v.surf.* first. Hamish ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user
[GRASS-user] (no subject)
a href=http://dehradun.org.in/wp-content/plugins/extended-comment-options/02efpk.html; http://dehradun.org.in/wp-content/plugins/extended-comment-options/02efpk.html/a___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user
[GRASS-user] (no subject)
Roger Bivand: I don't know how you you started the R session. I hope you started R within a running GRASS session, as you seem to be using an existing location. in http://grass.osgeo.org/wiki/R_statistics#MS_Windows there are some hints (Usage II) how you can start R within a GRASS session in windows with the wingrass6.4.x-installer (http://grass.osgeo.org/grass64/binary/mswindows/native/ or http://wingrass.fsv.cvut.cz/grass64/). the example above in the wiki is with wingrass64x in text-mode. FYI with wingrass70 (http://wingrass.fsv.cvut.cz/grass70/) you can start a R-session within a GRASS session also in wingrass-GUI-mode. best regards Helmut ___ Schon gehört? WEB.DE hat einen genialen Phishing-Filter in die Toolbar eingebaut! http://produkte.web.de/go/toolbar smime.p7s Description: S/MIME Cryptographic Signature ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user
[GRASS-user] (no subject)
Pedro Camilo Alcantara Concepcion ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user
[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 columns that v.rast.stats adds to the vector layer. These polygons are not the smallest or on the edges and I've used this procedure successfully before. Below are the region, v.info, r.info v.rast.stats outputs. Can anyone perhaps see why this happens? g.region -p projection: 99 (Transverse Mercator) zone: 0 datum: ** unknown (default: WGS84) ** ellipsoid: wgs84 north: -3122190 south: -3155760 west: -66510 east: -33390 nsres: 90 ewres: 90 rows: 373 cols: 368 cells: 137264 v.info map=subs22_clean...@chi ++ | Layer: subs22_clean...@chi | Mapset: CHI | Location:c83a_tm_29deg_e | Database:F:\Hanlie\grassdata | Title: | Map scale: 1:1 | Map format: native | Name of creator: hanlie | Organization: | Source date: Wed Nov 17 13:24:57 2010 || | Type of Map: vector (level: 2) | | Number of points: 0 Number of areas: 22 | Number of lines:0 Number of islands:1 | Number of boundaries: 822 Number of faces: 0 | Number of centroids:22 Number of kernels:0 | | Map is 3D: No | Number of dblinks: 1 | | Projection: Transverse Mercator | N: -3122280S: -3155670 | E:-33570W:-66510 | | Digitization threshold: 0 | Comments: | ++ r.info map=c83a_dem_srtm_90m_sl...@permanent ++ | Layer:c83a_dem_srtm_90m_sl...@perma Date: Wed Nov 03 16:21:11 2010 | Mapset: PERMANENT Login of Creator: Administrator | Location: c83a_tm_29deg_e | DataBase: F:\Hanlie\grassdata | Title:percent slope ( c83a_dem_srtm_90m_slope ) | Timestamp: none |--- | | Type of Map: raster Number of Categories: 255 | Data Type:FCELL | Rows: 373 | Columns: 368 | Total Cells: 137264 |Projection: Transverse Mercator |N: -3122190S: -3155760 Res:90 |E: -33390W: -66510 Res:90 | Range of data:min = 0.023858 max = 110.669037 | | Data Source: |raster elevation file c83a_dem_srtm_...@permanent | | | Data Description: |generated by r.slope.aspect | | Comments: |slope map elev = c83a_dem_srtm_...@permanent |zfactor = 1.00 format = percent |min_slp_allowed = 0.00 | ++ v.rast.stats raster=c83a_dem_srtm_90m_sl...@permanent vector=subs22_clean...@chi =area =1 colprefix=slope Statistics calculated from raster map and uploaded to attribute table of vector map . Done. Successfully finished Thanks Hanlie ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user
[GRASS-user] (no subject)
Hi, It's a return of an old bug. It was fixed 8 months a go. [http://trac.osgeo.org/grass/ticket/980] WinGRASS gurus - please change any other places where HTML browser is hardcoded to some .exe file. See ticket #980 for details. Maris. I've fixed ticket 980 in all cases in the nsis-installer, so the available standalone WinGrass-installer should have incorporated all fixes. AFAIK this is an osgeo4w issue. in the grass-starting routine there are following variables set in C:\OSGeo4W\etc\ini\grass set GISBASE=%OSGEO4W_ROOT%\apps\grass\grass-6.4.0svn set WINGISBASE=%OSGEO4W_ROOT%\apps\grass\grass-6.4.0svn set GRASS_SH=%OSGEO4W_ROOT%\apps\msys\bin\sh.exe set GRASS_WISH=%OSGEO4W_ROOT%\bin\wish.exe set GRASS_PYTHON=%OSGEO4W_ROOT%\bin\python.exe set GRASS_PROJSHARE=%OSGEO4W_ROOT%\share\proj set GRASS_HTML_BROWSER=%PROGRAMFILES%\Internet Explorer\iexplore.exe so this should be fixed there. please open a ticket there. best regards Helmut ___ WEB.DE DSL Doppel-Flat ab 19,99 euro;/mtl.! Jetzt auch mit gratis Notebook-Flat! http://produkte.web.de/go/DSL_Doppel_Flatrate/2 smime.p7s Description: S/MIME Cryptographic Signature ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user
[GRASS-user] (no subject)
I used the r.in.gdal to import the raster map, and then i pressed the button add raster layer, to add the raster file. But anything happened..the output says that everything matches but is not possible to add any raster file and even to display it. ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user
Re: [GRASS-user] (no subject)
On Fri, 29 Jan 2010, Matteo Toro wrote: I used the r.in.gdal to import the raster map, and then i pressed the button add raster layer, to add the raster file. But anything happened..the output says that everything matches but is not possible to add any raster file and even to display it. Matteo, I assume that you're using the wxPython UI for this. There are three distinct steps involved: 1) import the raster map; 2) indicate by the checkbox that you want the imported map added to the list of loaded layers; 3) right-click on the layer(s) you want displayed and select the option to zoom and display the selected layers. The map is not automatically displayed upon import. Rich ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user
[GRASS-user] (no subject)
dear all, I have installed GRASS and QGIS on Fedora 11. The GRASS problem is quoted as below. It gets the following message and once it has crashed too. = Starting GRASS ... /usr/lib/python2.6/site-packages/wx-2.8-gtk2-unicode/wx/_core.py:14450: UserWarning: wxPython/wxWidgets release number mismatch warnings.warn(wxPython/wxWidgets release number mismatch) == The Other QGIS issue is that it shows that it has been installed but it wont run from the Applications menu. And, as I try to run it from console by typing qgis it says Segmentation fault. === Thanks! Ambrish Dhaka _ We all see it as it is. But on MSN India, the difference lies in perspective. http://in.msn.com___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user
[GRASS-user] (no subject)
Hamish wrote: Wim de Vries wrote: (v.5 on suse 10.3 as well as v.6.2 on debian lenny). I came to the point of defining a new location, but then the list for projections does not show Lambert conformal or mercator. it's there. lcc, merc, etc. Ok.Thanks a lot. I'll have another try.. Next, I cannot get out of the list. Had to hit CTRL-C to completely drop out of the app. q this is something which to me is rather annoying in the debian package. they/we change the default pager (press any key for next page) to be the 'less' program instead of grass's default which is 'more'. the result is that you have to be used to it and know to press 'q' to get out of the list. with 'more' all you have to do is bash the spacebar and enter enough times and you get out automatically. sorry for the inconvenience. you can change it back to 'more' or whatever by setting the GRASS_PAGER variable, see the grass variables help page. for lenny you should get the latest grass, gdal, and proj packages from debian's www.backports.org. then you can have a pretty recent version to play with without much more than an apt-get call. Hamish ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user
[GRASS-user] (no subject)
Markus, Thanks for the feedback. The advice is very helpful. I've been slowly getting used to GRASS but still have a way to go. If I can stick with it I hope to learn enough to create a week-long GRASS for conservation workshop that we can offer at our research station in southeastern Arizona and in our projects in SE Asia. Here are some clarifying remarks and a few more questions if anyone can help. It will work. Everything takes time but I don't expect problems. GRASS 6.2 has been used for 10GB images, since then more testing and fixing happened. Best is to use the latest 6.4.0RC2 version. I'll see if I can get this version set up so it doesn't conflict with the GRASS version bundled with QGIS. My linux (Ubuntu) skill are also still quite poor but I'm learning (again). 4) Georeference shrub / non-shrub map to reference Landsat images [might do this with ENVI or ERDAS unless GRASS is a good choice] Also unclear to me: above you say that you generate the shrub / non-shrub map in GRASS, then it is already geocoded. Or not? Sorry - I wasn't very clear. The images I'm working with are geocoded but the image-to-image registration is poor. I plan to project all of the images into the same projection (some are in UTM and some use different datums) and then do an image-to-image registration. I wasn't sure if GRASS would be an efficient tool to do the image-to-image registration. I should just try it. 5) Create percent shrub cover map using regression tree algorithm, Landsat imagery, and shrub location data from the high resolution shrub / non-shrub map. I will probably do this using proprietary software unless I can do it easily in GRASS since a method has been established for another project. Of course GRASS-R-extention comes to mind for all kinds of statistics. I'm considering using GRASS for this but would like to know if anyone has used GRASS with R as the regression tree engine. I have a vague recollection of someone doing this several years ago but didn't find anything after a quick search. I suppose it would require getting training data from GRASS to R to generate the regression tree and then port the tree results back to GRASS to apply the rules to the dataset. It would be great if an interface already existed to do this but if not I can probably figure it out. I'm just not sure how long it will take me. All the best, Ned ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user
[GRASS-user] (no subject)
unsubscribe ¡Todo sobre Amor y Sexo! La guía completa para tu vida en Mujer de Hoy. http://mx.mujer.yahoo.com/ ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user
Re: [GRASS-user] (no subject)
On Mon, 2008-12-01 at 13:46 -0800, Avila Jose wrote: unsubscribe ¡Todo sobre Amor y Sexo! La guía completa para tu vida en Mujer de Hoy. http://mx.mujer.yahoo.com/ ___ grass-user mailing list grass-user@lists.osgeo.org [*] http://lists.osgeo.org/mailman/listinfo/grass-user Dear Jose, you can use the link [*] above to unsubscribe yourself :-). Once you visit this page you will find the unsubscribe button on the bottom. Kind regards, Nikos ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user
[GRASS-user] (no subject)
Dear GRASS-users, does someone of you know how to export a GRASS-polygone layer to CAD polygones? I used GRASS6.3.0 to create closed polygones from a 'unclean' dwg file and want to export this GIS-file to go on drawing within a CAD-programm (we're an architecure/urban planner office and my collegues are not used to use GIS). I tried v.out.dxf , but this does explode the polygones to single lines, so it ain't working. Dwg might be a solution but there's no v.out.dwg in my GRASS, converters from dxf to dwg are not free of license and the registration in opendwg also is not free of charge. Thanks for your help, Patrick Schirmer -- GMX startet ShortView.de. Hier findest Du Leute mit Deinen Interessen! Jetzt dabei sein: http://www.shortview.de/[EMAIL PROTECTED] ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user
Re: [GRASS-user] (no subject)
Can your CAD system read shapefiles? Mark ~ iPhone'd On Sep 29, 2008, at 3:33 AM, Patrick Schirmer [EMAIL PROTECTED] wrote: Dear GRASS-users, does someone of you know how to export a GRASS-polygone layer to CAD polygones? I used GRASS6.3.0 to create closed polygones from a 'unclean' dwg file and want to export this GIS-file to go on drawing within a CAD- programm (we're an architecure/urban planner office and my collegues are not used to use GIS). I tried v.out.dxf , but this does explode the polygones to single lines, so it ain't working. Dwg might be a solution but there's no v.out.dwg in my GRASS, converters from dxf to dwg are not free of license and the registration in opendwg also is not free of charge. Thanks for your help, Patrick Schirmer -- GMX startet ShortView.de. Hier findest Du Leute mit Deinen Interessen! Jetzt dabei sein: http://www.shortview.de/[EMAIL PROTECTED] ___ 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] (no subject)
Dear GRASS users, could somebody help me to find some solution for the following task? I have a raster map representing many distinct areas (MAP_1) and a second int raster with specific values (MAP_2). The result map should contain the same areas of MAP_1 but with the max MAP_2 value of each area as IDs. I tried some approaches with r.mapcalc but I don't know how to code the region aspect of MAP_1 into r.mapcalc. MAP_1 MAP_2 ---MAP_MAX - - --- 32 5266723974 95 111222 1429784155 ---777555 111322 4152763923 777955 12 4125718443 75 Thanks a lot for any help! Regards, Christian. This mail was sent through http://webmail.uni-jena.de ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user
RE: [GRASS-user] (no subject)
Hi Dylan, It looks like r.param.scale is the module he's looking for...thanks! Bob -Original Message- From: Dylan Beaudette [mailto:[EMAIL PROTECTED] Sent: Tuesday, June 03, 2008 2:32 PM To: grass-user@lists.osgeo.org Cc: Moskovitz, Bob Subject: Re: [GRASS-user] (no subject) Hi Bob, Have a look at r.param.scale -- this module allows you to specify the window size. Cheers, Dylan On Tuesday 03 June 2008, Moskovitz, Bob wrote: Hello list, A coworker was asking me if I know of a slope curvature algorithm that lets you adjust the kernel size. Does GRASS have such a program? Bob Bob Moskovitz Seismic Hazard Zonation Project California Geological Survey http://www.conservation.ca.gov/cgs/shzp?xml:namespace prefix = o ns = urn:schemas-microsoft-com:office:office / CONFIDENTIALITY NOTICE: This communication is intended only for the use of the individual or entity to which it is addressed. This message contains information from the State of California, California Geological Survey, which may be privileged, confidential and exempt from disclosure under applicable law, including the Electronic Communications Privacy Act. If the reader of this communication is not the intended recipient, you are hereby notified that any dissemination, distribution, or copying of this communication is strictly prohibited. -- Dylan Beaudette Soil Resource Laboratory http://casoilresource.lawr.ucdavis.edu/ University of California at Davis 530.754.7341 ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user
[GRASS-user] (no subject)
Hello list, A coworker was asking me if I know of a slope curvature algorithm that lets you adjust the kernel size. Does GRASS have such a program? Bob Bob Moskovitz Seismic Hazard Zonation Project California Geological Survey http://www.conservation.ca.gov/cgs/shzp?xml:namespace prefix = o ns = urn:schemas-microsoft-com:office:office / CONFIDENTIALITY NOTICE: This communication is intended only for the use of the individual or entity to which it is addressed. This message contains information from the State of California, California Geological Survey, which may be privileged, confidential and exempt from disclosure under applicable law, including the Electronic Communications Privacy Act. If the reader of this communication is not the intended recipient, you are hereby notified that any dissemination, distribution, or copying of this communication is strictly prohibited. ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user
Re: [GRASS-user] (no subject)
Hi Bob, Have a look at r.param.scale -- this module allows you to specify the window size. Cheers, Dylan On Tuesday 03 June 2008, Moskovitz, Bob wrote: Hello list, A coworker was asking me if I know of a slope curvature algorithm that lets you adjust the kernel size. Does GRASS have such a program? Bob Bob Moskovitz Seismic Hazard Zonation Project California Geological Survey http://www.conservation.ca.gov/cgs/shzp?xml:namespace prefix = o ns = urn:schemas-microsoft-com:office:office / CONFIDENTIALITY NOTICE: This communication is intended only for the use of the individual or entity to which it is addressed. This message contains information from the State of California, California Geological Survey, which may be privileged, confidential and exempt from disclosure under applicable law, including the Electronic Communications Privacy Act. If the reader of this communication is not the intended recipient, you are hereby notified that any dissemination, distribution, or copying of this communication is strictly prohibited. -- Dylan Beaudette Soil Resource Laboratory http://casoilresource.lawr.ucdavis.edu/ University of California at Davis 530.754.7341 ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user
[GRASS-user] (no subject)
Hello and thanks for your input Is it possible to draw a Profile based on a digital hight model and a selectet line? Sincerely Lukas -- -- Arbeit: Tel: +41 44 632 55 98 ETH Zürich Institut für terrestrische Ökosysteme Professur Waldökologie Universitätstrasse 22 CH-8092 Zürich privat: Lukas Rohrer c/o Dietschi Hohlstrasse 212 8004 Zürich +41 43 322 03 26 +41 78 624 93 20 [EMAIL PROTECTED] May the forest be with you! ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user
Re: [GRASS-user] (no subject)
Lukas Rohrer wrote: Hello and thanks for your input Is it possible to draw a Profile based on a digital hight model and a selectet line? Sincerely Lukas -- -- Arbeit: Tel: +41 44 632 55 98 ETH Zürich Institut für terrestrische Ökosysteme Professur Waldökologie Universitätstrasse 22 CH-8092 Zürich privat: Lukas Rohrer c/o Dietschi Hohlstrasse 212 8004 Zürich +41 43 322 03 26 +41 78 624 93 20 [EMAIL PROTECTED] May the forest be with you! ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user We made a script that draw a profile dem based tracking a line The script export the profile in .dxf format with distances and level labels and write a text file with the profile metadata There are some errors and bugs because it's in development. Andrea Ricci, Pierpaolo Chiraz gfosservices.it #!/bin/bash # # # # MODULE: r.section for GRASS 6.0 (GRASS Shell Script) # AUTHOR(S):Andrea Ricci, Pierpaolo Chiraz # # PURPOSE: Creates dxf territorial profile and dxf export; dem based # COPYRIGHT: (C) 2007 by the GFOSSERVICES # REQUIREMENTS: # # This program is free software under the GNU General Public # License (=v2). Read the file COPYING that comes with GRASS # for details. # ### #OUT=/tmp/out_$$ #set -m #mkfifo $OUT #%Module #% description: Create ascii_dxf profile from a raster dem. #%End #%option #% key: cella #% type: integer #% answer: 10 #% description: (raster cells resolution) #% required : yes #%end #%option #% key: DEM #% type: string #% gisprompt: old,cell,raster #% description: (input raster map) #% required : yes #%end #%option #% key: OUT2 #% type: string #% answer: ascii #% description: (Nome del file ascii di output) #% required : yes #%end #%option #% key: grass_vect #% type: string #% gisprompt: new,vector,vector #% answer: sezione1 #% description: (Nome del file vettoriale di Grass) #% required : yes #%end #%option #% key: dxf_vect #% type: string #% answer: sezione1 #% description: (Nome del file .dxf) #% required : yes #%end #%option #% key: percorso_dxf #% type: string #% description: (Percorso di salvataggio del file dxf) #% gisprompt: new_file,file,output #% required : yes #%end #%option #% key: fattore #% type: integer #% answer: 1 #% description: (fattore di esagerazione delle quote) #% required : no #%end #%option #% key: griglia #% type: integer #% answer: 100 #% description: (passo della griglia x,y) #% required : no #%end #-Input testuale metadato- #%option #% key: meta #% type: string #% answer: #% description: (file con metadati della sezione) #% gisprompt: new_file,file,input #% required : no #%end if [ -z $GISBASE ] ; then echo You must be in GRASS GIS to run this program. exit 1 fi if [ $1 != @ARGS_PARSED@ ] ; then exec g.parser $0 $@ fi DEM=$GIS_OPT_DEM OUT2=$GIS_OPT_OUT2 grass_vect=$GIS_OPT_GRASS_VECT dxf_vect=$GIS_OPT_DXF_VECT percorso_dxf=$GIS_OPT_PERCORSO_DXF cella=$GIS_OPT_CELLA OUT=$GIS_OPT_OUT meta=$GIS_OPT_META fattore=$GIS_OPT_FATTORE griglia=$GIS_OPT_GRIGLIA if [ -z $cella ] ; then cella=1;echo -n Risoluzione del raster (Dimensioni del lato della cella)= $cella metri fi g.region res=$cella #Punti della sezione (linea) tramite la quale si calcolano i valori di giacitura del piano xmessage ' INSERISCI 2 PUNTI PER INDIVIDUARE LA TRACCIA DI SEZIONE ' d.what.rast -t punti_sezione cat punti_sezione | grep -vpunti_sezione_corretti N_punti_corretti=` cat punti_sezione_corretti | wc -l` z=$(($N_punti_corretti/2)) #-- #ciclo per ottenere la X e la Y dei due punti individuati #-- for ((x=1; x=$z; x++)) do #creazione di un file di estrapolazione dati corretto (cioe' senza linee vuote) b=$(($x*2)) cat punti_sezione_corretti | sed 's/:/ /g' | awk -F '( NR == '$b' ) { print $1, $2 }' XY_PUNTI_FINALE done #-- #ciclo per ottenere la Z dei due punti individuati #-- cat punti_sezione | grep actual | sed 's/,/ /g' | sed 's/:/ /g' z_punti_sezione_corretti NR_Z=` cat z_punti_sezione_corretti | wc -l` for ((x=1; x=$NR_Z; x++)) do cat z_punti_sezione_corretti | awk -F '( NR == '$x' ) { print $3 }' Z_PUNTI_FINALE done #-- #rimozione di alcuni files testuali rm punti_sezione punti_sezione_corretti z_punti_sezione_corretti
Re: [GRASS-user] (no subject)
Lukas Rohrer wrote: Is it possible to draw a Profile based on a digital hight model and a selectet line? sure, see r.profile and r.transect modules. If your line is already a vector line and not just a collection of coordinates, you could use 'v.out.ascii format=standard' first to get the line's coordinates. r.profile will take input from stdin, so with a little awk magic you could pipe v.out.ascii directly into r.profile. If you figure that out, maybe post it back here so we could add it as an example on the r.profile help page. see also the profile tool in the GIS.m GUI menu for interactive profile creation. (you could display the line and trace over it) Hamish Be a better friend, newshound, and know-it-all with Yahoo! Mobile. Try it now. http://mobile.yahoo.com/;_ylt=Ahu06i62sR8HDtDypao8Wcj9tAcJ ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user
[GRASS-user] (no subject)
I am curious if any of you experienced developers has ever thought about using GPU as a CPU (or Co-CPU) in order to boost calculation speed... ? Is it or is not possible? http://www.gpgpu.org/ Cheers, Nikos signature.asc Description: This is a digitally signed message part ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user
Re: [GRASS-user] (no subject)
Nikos Alexandris wrote: I am curious if any of you experienced developers has ever thought about using GPU as a CPU (or Co-CPU) in order to boost calculation speed... ? Is it or is not possible? Most of GRASS' performance issues arise from the fact that code simplicity takes precedence over perfomance. Writing optimised code requires a fair amount of effort, and maintaining optimised code is even worse. Consequently, most of GRASS is written without even basic optimisations, let alone more extreme measures such as using a GPU as a co-processor. -- Glynn Clements [EMAIL PROTECTED] ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user
[GRASS-user] (no subject)
hi, can anyone tell that whether i can add rivers and trees on the road map for finding shortest path___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user