Re: [GRASS-user] Long post - Remote Sensing - Classifications
mitch wrote: GUI: I got stuck at the i.cluster window. How do I set the sampling interval values? I'm not allow to use any kind of character but numbers...how do I differentiate row and column value? http://n2.nabble.com/file/n2369388/sampling.png If I leave 1 value I get this error: ERROR: option sample must be provided in multiples of 2 You provided 1 items: 50 that's a bug, it should let you give two values. --help says: [sample=row_interval,col_interval] I've notice that i.group window has been changed: shouldn't “Name of imagery sub-group” and “Name of the raster map(s) to include in group” be under the Required tab rather than Optional? 'group=name' is the only required option. you can list maps in the group without adding any, so the subgroup and input raster maps are not needed sometimes. also subgroup really is optional. btw (for better or worse) those required/optional groupings happen automatically based on module requirements. I'm not really a fan of that but others think it's a useful feature. shrug (don't get me wrong, I don't want to criticize in a “bad” way, it's because maybe you receive a lot of feedbacks about other parts and this module is not so usedI haven't seen so many Remote Sensing questions lately...that's why). polite and researched feedback is always welcome! once you are comfortable that an issue is really a bug, please file a ticket, one for each issue. if just left in the mailing list it'll more easily be forgotten in the noise. -Supervised- I couldn't figure out how to make it work. ...but nothing more than that. Shouldn't I be prompted to a display where I can draw my region of interest? (I don't know, maybe I'm expecting something that it's not supposed to be). for sure old xmon based interactive monitors are gone for wxPy on Windows. not sure about what to expect from these modules, but supervised says to me some sort of map interaction is required... GUI: Imagery → Classify image → Interactive input for supervised classification (requires Xterm) it doesn't just need an terminal window for the text based menu data entry, in needs an Xmonitor as well, which is unavailable in the native MS Windows build. a blank Monitor display and a new shell appear and nothing else. we need a better sorry, won't work on MS Windows (so far) error message there. Are all of the GUI issues due to this http://grass.osgeo.org/wiki/WxPython-based_GUI_for_GRASS#Imagery_tools (Status: development not started yet) perhaps better stated as redevelopment not started in earnest yet, Some behind-the-scenes code has been prepared, but not much. Volunteers welcome! Hamish ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user
[GRASS-user] problem with export of raster data
Hello I try to export a raster integer map whose values range from - to 5880712. How can i specify that i want integer 4 bytes and not integer 2 bytes (it makes integer 2 bytes)?? Thank you. Stéphanie ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user
[GRASS-user] problem with r.surf.idw
Dear all, I want to interpolate a raster map that contains real values. I KNOW the values are real because I output the raster values with r.out.ascii command and the file shows real values. The problem is when I interpolate the raster with the r.surf.idw function and output the interpolated raster value, the file only contains integer (probably truncated) values. the command I type is : r.surf.idw input=original_raster output=intepolated_raster --o Obviously, typing a d.what.rast on the intepolated raster, I only see integer values (as in the ascii file generated). So, why interpolating produces only integer values?? thank you in advance, -- Alberto Pettazzi MeteoGalicia - Departamento de Climatología y Observación Consellería de Medio Ambiente e Desenvolvemento Sostible Rúa de Roma, 6 15707 Santiago de Compostela. A Coruña Teléfono: +34-881-999646 e-mail: alberto.petta...@meteogalicia.es mailto:alberto.petta...@meteogalicia.es ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user
Re: [GRASS-user] Invert values of raster image
Thanks Moritz, r.recode did exactly what I needed. much appreciated Wes Moritz Lennert mlenn...@club.worldonline.be 02/23/09 1:02 PM On 23/02/09 11:16, Wesley Roberts wrote: Dear list. I have a raster file with values coded in 8 bit grey scale 0-255. I would like to invert the raster values so that 255 is now 0 and 0 is now 255. I have not found an out of the box solution in grass but have tried the r.reclass function. How about r.recode with a rule such as: old_low:old_high:new_low:new_high 0:255:255:0 ? I generated a reclass rule file but am not sure about the form of the file. You need a '=', so 0=255 1=254 etc Moritz -- This message is subject to the CSIR's copyright terms and conditions, e-mail legal notice, and implemented Open Document Format (ODF) standard. The full disclaimer details can be found at http://www.csir.co.za/disclaimer.html. This message has been scanned for viruses and dangerous content by MailScanner, and is believed to be clean. MailScanner thanks Transtec Computers for their support. ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user
Re: [GRASS-user] problem with export of raster data
On Mon, Feb 23, 2009 at 12:15 PM, FAROUX STEPHANIE stephanie.far...@meteo.fr wrote: Hello I try to export a raster integer map whose values range from - to 5880712. How can i specify that i want integer 4 bytes and not integer 2 bytes (it makes integer 2 bytes)?? If you export with r.out.gdal you can choose among type File type options: Byte,Int16,UInt16,Int32,UInt32,Float32,Float64,CInt16, CInt32,CFloat32,CFloat64 This is the range which GRASS support for export to GDAL supported formats (of course it depends if the format is compatible with the selected type, see GDAL manual pages for details). Markus ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user
[GRASS-user] problem in projection of shape files
hi I have a shape file which I have to import in GRASS. Shape file is consists of a shape file, a dbf file and a prj file. The prj file has following information. PROJCS[BTM, GEOGCS[GCS_Everest_1830, DATUM[D_Everest_1830, SPHEROID[Everest_1830,6377299.36,300.8017]], PRIMEM[Greenwich,0.0], UNIT[Degree,0.0174532925199433]], PROJECTION[Transverse_Mercator], PARAMETER[False_Easting,50.0], PARAMETER[False_Northing,-200.0], PARAMETER[Central_Meridian,90.0], PARAMETER[Scale_Factor,0.9996], PARAMETER[Latitude_Of_Origin,0.0], UNIT[Meter,1.0]] I will be very grateful if anyone could help me to find the ESPG code and any other way to re-project the shape file and import into GRASS. Best Md. Zahid Parvez ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user
[GRASS-user] OSGeo GRASS GIS Report 2008
Dear Community, we are supposed to complete our GRASS GIS Report 2008 for OSGeo, please help to expand the draft: http://wiki.osgeo.org/wiki/GRASS_GIS_Report_2008 thanks Markus ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user
Re: [GRASS-user] problem in projection of shape files
I would use ogr2ogr http://www.gdal.org/ogr2ogr.html from GDAL to reproject using a proj4 string. For example to wgs84: ogr2ogr -s_srs +proj=tmerc +lat_0=0 +lon_0=90 +k=0.9996 +x_0=50 +y_0=-200 +a=6377276.345 +b=6356075.41314024 +units=m +no_defs -t_srs wgs84 out.shp in.shp -Jamie On Mon, Feb 23, 2009 at 1:19 PM, Zahid Parvez grassgisbanglad...@gmail.comwrote: hi I have a shape file which I have to import in GRASS. Shape file is consists of a shape file, a dbf file and a prj file. The prj file has following information. PROJCS[BTM, GEOGCS[GCS_Everest_1830, DATUM[D_Everest_1830, SPHEROID[Everest_1830,6377299.36,300.8017]], PRIMEM[Greenwich,0.0], UNIT[Degree,0.0174532925199433]], PROJECTION[Transverse_Mercator], PARAMETER[False_Easting,50.0], PARAMETER[False_Northing,-200.0], PARAMETER[Central_Meridian,90.0], PARAMETER[Scale_Factor,0.9996], PARAMETER[Latitude_Of_Origin,0.0], UNIT[Meter,1.0]] I will be very grateful if anyone could help me to find the ESPG code and any other way to re-project the shape file and import into GRASS. Best Md. Zahid Parvez ___ 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] problem in projection of shape files
On Mon, Feb 23, 2009 at 10:19 PM, Zahid Parvez grassgisbanglad...@gmail.com wrote: hi I have a shape file which I have to import in GRASS. Shape file is consists of a shape file, a dbf file and a prj file. The prj file has following information. PROJCS[BTM, GEOGCS[GCS_Everest_1830, DATUM[D_Everest_1830, SPHEROID[Everest_1830,6377299.36,300.8017]], PRIMEM[Greenwich,0.0], UNIT[Degree,0.0174532925199433]], PROJECTION[Transverse_Mercator], PARAMETER[False_Easting,50.0], PARAMETER[False_Northing,-200.0], PARAMETER[Central_Meridian,90.0], PARAMETER[Scale_Factor,0.9996], PARAMETER[Latitude_Of_Origin,0.0], UNIT[Meter,1.0]] I will be very grateful if anyone could help me to find the ESPG code and the best is to ask in the PROJ mailing list which code it is: http://trac.osgeo.org/proj/ - Mailing List Just to be sure that it is right. any other way to re-project the shape file and import into GRASS. Then ogr2ogr or a proper GRASS location will do it. Markus ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user
[GRASS-user] Grass module to use to create a bounding box
Hi List I have a raster and a vector layer but would want to limit the analysis to the bounding box of the vector layer - how do I accomplish this. Is there a vector module to create a vector bounding box (v.*). Thanks -- Sent from my mobile device Sindile Bidla ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user
Re: [GRASS-user] Grass module to use to create a bounding box
You should be able to set the region to the vector layer and then your set. You plan to do the analysis on the raster, correct? John On 2/23/09, Intengu Technologies sindile.bi...@gmail.com wrote: Hi List I have a raster and a vector layer but would want to limit the analysis to the bounding box of the vector layer - how do I accomplish this. Is there a vector module to create a vector bounding box (v.*). Thanks -- Sent from my mobile device Sindile Bidla ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user -- -- John C. Tull, Ph.D. Conservation Director Nevada Wilderness Project 8550 White Fir Street Reno, NV 89523 775.746.7851 (office) 775.224.2947 (mobile) www.wildnevada.org ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user
[GRASS-user] Re: Grass module to use to create a bounding box
Yes analysis on the raster - how do I do that John as I have a region defined already, which is larger than the extent of the vector. On 2/24/09, John C. Tull john.t...@wildnevada.org wrote: You should be able to set the region to the vector layer and then your set. You plan to do the analysis on the raster, correct? John On 2/23/09, Intengu Technologies sindile.bi...@gmail.com wrote: Hi List I have a raster and a vector layer but would want to limit the analysis to the bounding box of the vector layer - how do I accomplish this. Is there a vector module to create a vector bounding box (v.*). Thanks -- Sent from my mobile device Sindile Bidla ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user -- -- John C. Tull, Ph.D. Conservation Director Nevada Wilderness Project 8550 White Fir Street Reno, NV 89523 775.746.7851 (office) 775.224.2947 (mobile) www.wildnevada.org -- Sent from my mobile device Sindile Bidla ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user
[GRASS-user] v.external does not create vector layer
Hello, I am running GRASS 6.3.0 on Linux. I would like to use a postgis table as a read only vector layer in GRASS. One time I had v.external doing more than it does now. Here is a typical execution: # v.external 'dsn=PG:host=localhost user=postgres dbname=spatial_db' output=gpxpints_pg Data source contains 1 layers: trackpnt This would be the behavior I would expect if output was not defined. There is no vector layer created in GRASS. What is the next debugging step here? John ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user
Re: [GRASS-user] problem in projection of shape files
Zahid Parvez wrote: I have a shape file which I have to import in GRASS. Shape file is consists of a shape file, a dbf file and a prj file. The prj file has following information. I will be very grateful if anyone could help me to find the ESPG code and any other way to re-project the shape file and import into GRASS. on the grass (tcl/tk) start-up screen use the Define new location with ... [Georeferenced file] button, then to reproject run v.proj from the location you wish to reproject into. Knowing the EPSG code is not critical. You can also use the nice EPSG browse/search tool from the Define new location with ... [EPSG codes] window to track it down. (or search through the epsg file itself with your favorite text editor/viewer. Hamish ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user
Re: [GRASS-user] r.info doesn't print out all category values
William Hudspeth wrote: I have used to r.average to calculate mean pixel values by county. This involved converting the county vector layer to a raster, then using the county raster layer as the base in r.average. So far so, good. The r.average routine does what I want, but the resultant output raster map shows the following output with r.info: [r.info? are you sure? what's the exact command line params you used for each command?] ... 207740:207740:0.4507908000 207741:207741:0.563316 207742:207742:0.5043974000 207743:207743:0.432117 ... There are 52 categories, but the r.info output only shows about 45 of them, leaving the rest unreported. If I redirect to file, the same thing. If I look at the cats file under the cats subdirectory, it also fails to show all of the categories. How do I show everything that is generated by r.average??? can you reproduce the problem using the Spearfish or North Carolina sample datasets? then we can try to reproduce the effect. Hamish ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user
[GRASS-user] Problem using i.cluster on GRASS GIS 6.4.0RC3
I get the following error when using i.cluster: ERROR: Unable to create report file (null) This is how I wrote the commands i.group group=lsat1 subgroup=lsat1 in=L71171078_07820001007_B10,L71171078_07820001007_B20,L71171078_07820001007_B30,L71171078_07820001007_B40,L71171078_07820001007_B50,L72171078_07820001007_B70 i.cluster group=lsat1 subgroup=lsat1 sig=sig.cluster classes=15 sep=1.5 Windows XP - GRASS GIS 6.4.0RC3 I have seen the following tickets http://trac.osgeo.org/grass/ticket/174 and http://trac.osgeo.org/grass/ticket/70 which seem to suggest that the problem has been solved / there is a work around. Kindly advise on how I must proceed -- Sindile Bidla ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user
Re: [GRASS-user] problem with r.surf.idw
Alberto Pettazzi wrote: I want to interpolate a raster map that contains real values. I KNOW the values are real because I output the raster values with r.out.ascii command and the file shows real values. The problem is when I interpolate the raster with the r.surf.idw function and output the interpolated raster value, the file only contains integer (probably truncated) values. r.surf.idw and r.surf.idw2 only produce CELL (integer) maps. try using (r.to.vect +) v.surf.idw or v.surf.rst instead. Hamish ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user