[GRASS-user] Question about generating signature files with i.gensig
Greetings I have one question regarding usage of i.gensig If I have training pixel that is not available in all images of input group, will that pixel be considered? Thanks Jenny ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user
[GRASS-user] i.gensig for pixels without values
Greetings I have a raster with a few classes (1,2,3,4,5) and some pixels exists where I don't have values in my group images (no-values). When generating signature files, are those files considered in the calculation of average, mean and numver of pixels? THanks Jenny ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user
[GRASS-user] Using r.resamp.stats
I have one question regarding null values: I'm using r.resamp.stats with method average. My question is, if I have a Target pixel based on 10 null values and 1 non-null value, the target pixel will be based only on valid pixel? I mean, it only requires at least 1 value to produce an output value? Thanks Jenny ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user
[GRASS-user] Using a r.mapcalc expression in grass.mapcalc
I need to use this expression output=@input (input has labels and I want to export them so i need to do this). How can I do this in grass.mapcalc in Python? Does it accept @? Thanks ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user
[GRASS-user] Error while running v.db.addcol
Greetings I0'm using winGRASS. When I run v.db.addcol in my WinXP I get this error: Default driver / database set to: driver: dbf database: $GISDBASE/$LOCATION_NAME/$MAPSET/dbf/ 0 [main] sh 5240 open_stackdumpfile: Dumping stack trace to sh.exe.stackdump 0 [main] sh 5436 open_stackdumpfile: Dumping stack trace to sh.exe.stackdump c:/grass6/scripts/v.db.addcol: [: : integer expression expected 0 [main] sh 2264 open_stackdumpfile: Dumping stack trace to sh.exe.stackdump 0 [main] sh 3416 open_stackdumpfile: Dumping stack trace to sh.exe.stackdump ERROR: Column cat not found But when I run in Win7 or other Win Machine (XP also) i do not get this error. Any explanation for this weird incident? Thanks Jenny ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user
[GRASS-user] Compiling a module in Win OsGEO4w
Greetings I have developed a function, which also includes a shared lib (libsegmentimage.dll) and I have compiled with success in OSGEO4W in winGRASS. (before that I needed to place libsegmentimage.dll in C:\OSGeo4W\usr\src\grass-6.4.svn_src\dist.i686-pc-mingw32\lib . It has places my binary (i.segment.exe) in c:\C:\OSGeo4W\usr\src\grass-6.4.svn_src\dist.i686-pc-mingw32\bin\ and now I want to place it in C:\OSGeo4W\apps\grass\grass-6.4.1. TO do that I just copied the exe to bin folder but the bin it seems not to be working because I get this error: Traceback (most recent call last): File C:/OSGeo4W/apps/grass/grass-6.4.1svn/etc/wxpython/wx gui.py, line 486, in OnRunCmd self.goutput.RunCmd(cmd, switchPage=False) File C:\OSGeo4W\apps\grass\grass-6.4.1svn\etc\wxpython\gu i_modules\goutput.py, line 368, in RunCmd menuform.GUI().ParseCommand(cmdlist, parentframe=self) File C:\OSGeo4W\apps\grass\grass-6.4.1svn\etc\wxpython\gu i_modules\menuform.py, line 1866, in ParseCommand xml.sax.parseString(getInterfaceDescription(cmd[0]).decode(e nc).split('\n',1)[1].replace('', '?xml version=1.0 encoding=utf-8?\n', 1).encode(utf-8), IndexError : list index out of range By the way I have complied like this: make -C i.segment MODULE_TOPDIR=/C/OSGeo4W/usr/src/grass-6.4.svn Can anyone help me? Thanks Jenny ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user
[GRASS-user] Error while running a compiled wingrass module in WinXP and win/
Greetings After Antonio's reply I tried to run my compiled module in WinGRASS (in winXP) and in Win7 and it runned ok in WinXP but I got this error in Win7: This application has requested the Runtime to terminate it in an unusual way. Please contact the application's support team for more information. And the process was terminated. Why might be wrong? is there any known limitation of win7 when compared with winXP? ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user
[GRASS-user] Retrieve list of values in Python
Greetings I would like to know how can I retrieve the list of available values in a Raster map in a Python Script? Thanks Best regards, Jenny Turner ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user
Re: [GRASS-user] GRASS'c Python core function to use as g.mlist
Hello Ricardo Thanks for your feedback but the problem is that mapset cannot be used as input. Because I got this error: p=grass.mlist_grouped ('rast', pattern=pattern, mapset='PERMANENT') TypeError: mlist_grouped() got an unexpected keyword argument 'mapset'. If I avoid using mapset I get a dictionary variable in which, each entry is the list of rasters for each MAPSET. Am I doing something wrong? Thanks Jenny On Fri, Mar 11, 2011 at 2:29 PM, Ricardo Filipe Soares Garcia da ricardo.garcia.si...@gmail.com wrote: Hi Jenny Long story short: use your command like this (note the asterisk on the 'pattern' variable): # python code p=grass.mlist_grouped ('rast', pattern='output*', mapset='Mapping') I am using the Ipython shell to script grass. It has a cool feature, that if you type the name of a function followed by two question marks, you can see the source code of the function (as long as it is coded in python itself). So, in doing: grass.mlist_grouped?? It is possible to see that the mlist_grouped function is internally calling the 'g.mlist' command with only the 'm' flag. Reading the documentation for the g.mlist command [1], I realize that, without using the -r or -e flags, this command accepts wildcards[2] as the 'pattern' string. The asterisk (*) is the wildcard character for substituting for any zero or more characters. So, in order to get a string that starts with 'output' and after that can have anything, you just type 'output*' [1] - http://grass.fbk.eu/grass64/manuals/html64_user/g.mlist.html [2] - http://en.wikipedia.org/wiki/Wildcard_character#Computing On Fri, Mar 11, 2011 at 12:10 PM, Jenny Turner jennyturner1...@gmail.com wrote: Greetings I'm doing a GRASS Python Script where I need to list a set of raster maps but I'm not using its exact name but a basic regular expression (r flag in g.mlist function). So I want to find a list of rasters with output on its name: p=grass.mlist_grouped ('rast', pattern='output', mapset='Mapping') But I get null entries because no raster map is named output only output01 output02 output03 ... I want to list all raster maps with output on itsd name. How can I do this using grass.mlist_grouped? Thanks Jenny ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user -- ___ ___ __ Ricardo Garcia Silva ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user
[GRASS-user] GRASS'c Python core function to use as g.mlist
Greetings I'm doing a GRASS Python Script where I need to list a set of raster maps but I'm not using its exact name but a basic regular expression (r flag in g.mlist function). So I want to find a list of rasters with output on its name: p=grass.mlist_grouped ('rast', pattern='output', mapset='Mapping') But I get null entries because no raster map is named output only output01 output02 output03 ... I want to list all raster maps with output on itsd name. How can I do this using grass.mlist_grouped? Thanks Jenny ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user
Re: [GRASS-user] using try/except with grass.mapcalc in a Python Script
Hello Glynn What I want to do is run grass.mapcalc like this: try: grass.mapcalc() except: grass.fatal(_(An error occurred ...)) I know that this doesn't work. What I want to catch is not syntax erros but errors such as an raster input map that does not exist and is being used in grass.mapcal should raise an error or an excpetion in order to stop the process. Any sugfgestions? Thanks Jenny On Tue, Mar 1, 2011 at 7:26 PM, Glynn Clements gl...@gclements.plus.comwrote: Jenny Turner wrote: I want to include grass.mapcalc function in try except but I'm not catching the errors. Has anyone used this for grass.mapcalc? What sort of error are you expecting? If you call set_raise_on_error(), a non-zero exit status should raise a ScriptException instead of calling sys.exit() (note that sys.exit() raises a SystemExit exception, which can be caught if desired). Note that r.mapcalc will only return a non-zero exit status for fundamental errors such as syntax errors in the expression or failure to open a map. Calculation errors (e.g. division by zero or domain errors) result in null values. -- 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] using try/except with grass.mapcalc in a Python Script
Greetings I want to include grass.mapcalc function in try except but I'm not catching the errors. Has anyone used this for grass.mapcalc? Thanks Jenny ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user
[GRASS-user] Eliminating values based on a null
Greetings in my base map (called X) I have a few null values and, in Y map I want to eliminate pixels that match with null values in X. How can I do this? besides using mapcalc Thanks Jenny ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user
[GRASS-user] Re: Help/manual on Georectifier 6.4.0
Hi I have sent this message a couple of days agor but I didn't get any reply So, my only question now is (at least, please can anyone answer this: My active mapset shall be which one? the one where I have my reference images (already georrectified) or mapset where I have my not-georreferenced images)? Thanks Jenny On Thu, Nov 18, 2010 at 11:22 PM, Jenny Turner jennyturner1...@gmail.comwrote: Greetings I'm writting a tutorial of Georrectify tool (from 6.4.0 not latest launch) for my students. I was only able to find this: http://gama.fsv.cvut.cz/~landa/grass/swf/georect.html. Is there any other Wiki or document? Thanks Jenny ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user
Re: [GRASS-user] Re: Help/manual on Georectifier 6.4.0
Hi Markus Let me just say that, according to snapshots, you have done a great job for new WxGeorrect :) So, my only question now is (at least, please can anyone answer this: My active mapset shall be which one? the one where I have my reference images (already georrectified) or mapset where I have my not-georreferenced images)? The active location is the target location with a mapset with your reference images (indicated by the very first question of the georectifier where you have to specify the source and not the target location). My question was not related with Lcoation (because both maps (not-georreference and georreferenced) are in the same Location but in different mapsets In my example: Mapset:Local: not-georreference and to be georreferenced (How do you call this?) Mapset: Global: georreferenced (how do you call this?) In the first panel it's asked about Source Location and source Mapset. You mean that here I shall define where is my not-georreferenced image? In this case mapset:Local? Second panel asks for a group of images to be goerrectified. Third asks for an image to be displayed. in this case is a georreferenced or not-georreferenced image? So, in the normal Map diSplay window I will put an already georreferenced image right? in this case from Mapset Global? Thanks and sorry for such a dumb questions Jenny Source and target location can be identical. The manual of the new georectifier [1] may be able to answer further questions since all the functionality of the 6.4.0 georectifier is included in the georectifier in 6.4.1 and above. Markus M [1] http://grass.osgeo.org/grass64/manuals/html64_user/wxGUI.GCP_Manager.html Thanks Jenny On Thu, Nov 18, 2010 at 11:22 PM, Jenny Turner jennyturner1...@gmail.com wrote: Greetings I'm writting a tutorial of Georrectify tool (from 6.4.0 not latest launch) for my students. I was only able to find this:http://gama.fsv.cvut.cz/~landa/grass/swf/georect.html. Is there any other Wiki or document? Thanks Jenny ___ 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] Help/manual on Georectifier 6.4.0
Greetings I'm writting a tutorial of Georrectify tool (from 6.4.0 not latest launch) for my students. I was only able to find this: http://gama.fsv.cvut.cz/~landa/grass/swf/georect.html. Is there any other Wiki or document? Thanks Jenny ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user
[GRASS-user] Python Script is not able to find a file in a different Drive
Greetings I have a python script that needs to get a file (import a raster). I have tried in Ubuntu and also in Windows and it works but, When I have that file in a different hard drive (D: instead of C: where WinGRASS is installed) os.path.isfile(input) is false... Why? is it some limitation of GRASS? Thanjs~ Jenny ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user
[GRASS-user] Precision error in r.mapcalc
Greetings I have a raster (FCELL ) that has a minimum value of 0.1312103122 and If I try to have an image with just those values r.mapcalc minim= if(image==0.1312103122,10,0) I get an all Zero image. Then I opened the minim metadata and i see that in comments I have: if(clustered.1 == 0.13121031, 10, 0) r.mapcalc has eliminated the last 2 digits.How can I avoid this to happen? Thanks Jenny ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user
Re: [GRASS-user] Error in GUI's sub-group definition in i.maxlik and others modules
This is a follow-up of a message that I sent last week. I have compiled last weekly snapshot and the problem maintains. in functions just like i.maxlik I can select subgroup but it's not inserted in the command So all functions that requires sub-groups cannot be runned from graphic-interfaces. As well as I tried to add images to a group/sub-group and I cannot open the dropdown box for image selection in i.group. Shall I add this to the Trac list? On Wed, Jun 30, 2010 at 3:35 PM, Jenny Turner jennyturner1...@gmail.comwrote: Hi It might be something to do with that. But for instance, I have used GRASS6.4.RC5 from Jan/Feb and it was ok regarding this issue. and now it's not. Or am I doing sometyhing wrong? On Tue, Jun 29, 2010 at 12:42 PM, Markus Neteler nete...@osgeo.orgwrote: Hi, could it be related to http://trac.osgeo.org/grass/ticket/780 ? Markus On Tue, Jun 29, 2010 at 10:07 AM, Jenny Turner jennyturner1...@gmail.com wrote: Hi again This error not only occurs in i.maxlik but also in the others modules where users should select sub-group (except georrectify GU). Dones anyone has the same problem? Thanks Jenny On Mon, Jun 28, 2010 at 9:51 AM, Jenny Turner jennyturner1...@gmail.com wrote: Hello Martin I got this error in GRASS6.4.0.RC6, from two weeks ago in LinuxUbuntu9.10 but also in Windows. Jenny On Fri, Jun 25, 2010 at 9:57 PM, Martin Landa landa.mar...@gmail.com wrote: Hi, 2010/6/24 Jenny Turner jennyturner1...@gmail.com: I have been testing i.maxlik and I'm getting a strange error. When I select the group, in group parameter, automatically, sub-group list includes the available sub-groups but, when i select one, it's selected in the frame but it's not inserted in the command to be runned. In the bottom it's still REQUIRED. I can only run i.maxlik by inserting directly the parameters in the commandline-output. Has anyone having the same problem? GRASS version? OS? Martin -- Martin Landa landa.martin gmail.com * http://gama.fsv.cvut.cz/~landa ___ 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] Commandline for creating a vectorial
Hi In my case I'm using winGRASS. So If I do v.edit map=vector_file type=boundary,centroid tool=create It's exacly the same as I select Crete New Vector map (from Vector/Develp Vector Map/) without Creating attribute table? Thanks Jenny On Wed, Jun 30, 2010 at 8:09 PM, Markus Neteler nete...@osgeo.org wrote: On Wed, Jun 30, 2010 at 6:11 PM, Jenny Turner jennyturner1...@gmail.com wrote: Greetings I'm using GRASS6.4.0.RC6 in Windows and I need to create a script that creates vectorials. By using WxGUI menu is quite easy Vector/Develop vector map/Create new Vector map and then name and options. is there any commandline to perform exacly this? You can use v.edit for this: http://grass.osgeo.org/grass64/manuals/html64_user/v.edit.html Markus ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user
Re: [GRASS-user] Error in GUI's sub-group definition in i.maxlik and others modules
Hi It might be something to do with that. But for instance, I have used GRASS6.4.RC5 from Jan/Feb and it was ok regarding this issue. and now it's not. Or am I doing sometyhing wrong? On Tue, Jun 29, 2010 at 12:42 PM, Markus Neteler nete...@osgeo.org wrote: Hi, could it be related to http://trac.osgeo.org/grass/ticket/780 ? Markus On Tue, Jun 29, 2010 at 10:07 AM, Jenny Turner jennyturner1...@gmail.com wrote: Hi again This error not only occurs in i.maxlik but also in the others modules where users should select sub-group (except georrectify GU). Dones anyone has the same problem? Thanks Jenny On Mon, Jun 28, 2010 at 9:51 AM, Jenny Turner jennyturner1...@gmail.com wrote: Hello Martin I got this error in GRASS6.4.0.RC6, from two weeks ago in LinuxUbuntu9.10 but also in Windows. Jenny On Fri, Jun 25, 2010 at 9:57 PM, Martin Landa landa.mar...@gmail.com wrote: Hi, 2010/6/24 Jenny Turner jennyturner1...@gmail.com: I have been testing i.maxlik and I'm getting a strange error. When I select the group, in group parameter, automatically, sub-group list includes the available sub-groups but, when i select one, it's selected in the frame but it's not inserted in the command to be runned. In the bottom it's still REQUIRED. I can only run i.maxlik by inserting directly the parameters in the commandline-output. Has anyone having the same problem? GRASS version? OS? Martin -- Martin Landa landa.martin gmail.com * http://gama.fsv.cvut.cz/~landa ___ 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] Commandline for creating a vectorial
Greetings I'm using GRASS6.4.0.RC6 in Windows and I need to create a script that creates vectorials. By using WxGUI menu is quite easy Vector/Develop vector map/Create new Vector map and then name and options. is there any commandline to perform exacly this? Thanks Jenny ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user
Re: [GRASS-user] Error in GUI's sub-group definition in i.maxlik and others modules
Hello Martin I got this error in GRASS6.4.0.RC6, from two weeks ago in LinuxUbuntu9.10 but also in Windows. Jenny On Fri, Jun 25, 2010 at 9:57 PM, Martin Landa landa.mar...@gmail.comwrote: Hi, 2010/6/24 Jenny Turner jennyturner1...@gmail.com: I have been testing i.maxlik and I'm getting a strange error. When I select the group, in group parameter, automatically, sub-group list includes the available sub-groups but, when i select one, it's selected in the frame but it's not inserted in the command to be runned. In the bottom it's still REQUIRED. I can only run i.maxlik by inserting directly the parameters in the commandline-output. Has anyone having the same problem? GRASS version? OS? Martin -- Martin Landa landa.martin gmail.com * http://gama.fsv.cvut.cz/~landa ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user
[GRASS-user] Error in GUI's sub-group definition in i.maxlik and others modules
Greetings I have been testing i.maxlik and I'm getting a strange error. When I select the group, in group parameter, automatically, sub-group list includes the available sub-groups but, when i select one, it's selected in the frame but it's not inserted in the command to be runned. In the bottom it's still REQUIRED. I can only run i.maxlik by inserting directly the parameters in the commandline-output. Has anyone having the same problem? Thanks Jenny ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user
[GRASS-user] Converting raster to point-vector- Question
Hi there I have a few training areas, in raster format, that I have used for Clasification (i.maxlik). Thuis means that pixels have information like this: 633376.703089|222904.095941||9|Category 9 or 642250.637155|224541.035235||10|Category 10 I i convert using this expression: r.to.vect input=lsa7_trainin...@landsat output=points00 feature=point but with this method my cat number corresponds to an Id instead of category. What I need is to have cat = to value category. How can I change ths? Thanks Jenny ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user
[GRASS-user] Error while compiling r.in.gdal
Greetings I'm doing some testing in r.in.gdal (by changing a few stuff just to test) in LINUX with GRASS6.4.0svnRC6 but I go this error: gcc -L/usr/local/grass-6.4.0svn/lib -Wl,--export-dynamic -Wl,-rpath-link,/usr/local/grass-6.4.0svn/lib-o /usr/local/grass-6.4.0svn/bin/r.in.gdal OBJ.i686-pc-linux-gnu/main.o -lgrass_gproj -lgrass_gis -lgrass_datetime -lz -lproj -L/usr/lib -lgdal1.6.0 -lgrass_gis -lgrass_datetime -lz -lgrass_I -lgrass_gis -lgrass_datetime -lz -L/usr/lib -lgdal1.6.0 -lgrass_gmath -lgrass_gis -lgrass_datetime -lz -lm -lz /usr/bin/ld: cannot find -lgdal1.6.0 collect2: ld returned 1 exit status make: *** [/usr/local/grass-6.4.0svn/bin/r.in.gdal] Error 1 make: Leaving directory `/home/anrr/Desktop/dwe_devep/r.in.testing' Can anyone explain me what might be happening? Should I reinstall GRASS? Thanks jenny T. ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user
[GRASS-user] Re: Error while compiling r.in.gdal
Ups.. my mistake. the error is: make: Entering directory `/home/anrr/Desktop/dwe_devep/r.in.testing' gcc -L/usr/local/grass-6.4.0svn/lib -Wl,--export-dynamic -Wl,-rpath-link,/usr/local/grass-6.4.0svn/lib-o /usr/local/grass-6.4.0svn/bin/r.in.testing OBJ.i686-pc-linux-gnu/main.o -lgrass_gproj -lgrass_gis -lgrass_datetime -lz -lproj -L/usr/lib -lgdal1.6.0 -lgrass_gis -lgrass_datetime -lz -lgrass_I -lgrass_gis -lgrass_datetime -lz -L/usr/lib -lgdal1.6.0 -lgrass_gmath -lgrass_gis -lgrass_datetime -lz -lm -lz /usr/bin/ld: cannot find -lgdal1.6.0 collect2: ld returned 1 exit status make: *** [/usr/local/grass-6.4.0svn/bin/r.in.testing] Error 1 make: Leaving directory `/home/anrr/Desktop/dwe_devep/r.in.testing' On Mon, Jun 21, 2010 at 4:09 PM, Jenny Turner jennyturner1...@gmail.comwrote: Greetings I'm doing some testing in r.in.gdal (by changing a few stuff just to test) in LINUX with GRASS6.4.0svnRC6 but I go this error: gcc -L/usr/local/grass-6.4.0svn/lib -Wl,--export-dynamic -Wl,-rpath-link,/usr/local/grass-6.4.0svn/lib-o /usr/local/grass-6.4.0svn/bin/r.in.gdal OBJ.i686-pc-linux-gnu/main.o -lgrass_gproj -lgrass_gis -lgrass_datetime -lz -lproj -L/usr/lib -lgdal1.6.0 -lgrass_gis -lgrass_datetime -lz -lgrass_I -lgrass_gis -lgrass_datetime -lz -L/usr/lib -lgdal1.6.0 -lgrass_gmath -lgrass_gis -lgrass_datetime -lz -lm -lz /usr/bin/ld: cannot find -lgdal1.6.0 collect2: ld returned 1 exit status make: *** [/usr/local/grass-6.4.0svn/bin/r.in.gdal] Error 1 make: Leaving directory `/home/anrr/Desktop/dwe_devep/r.in.testing' Can anyone explain me what might be happening? Should I reinstall GRASS? Thanks jenny T. ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user
Re: [GRASS-user] Using Tck/tk vector digitizer for training areas
Thanks Micha for your help but I got an error while doing this 1- Create new empty vectorial 2- run v.db.addtable command that you kindly sent me 3- use v.digit as defined previously 4- Inside tcltk's digitizing menu select Digitize new Centroid 5- Leave layer as 1 6- Change Mode to Manual Entry 7- First digitize boundary (by selecting boundary) by selecting points with left-mouse-button and close polygon with right-mouse-button 8- A panel with FORM title appears with a tab saying Attributes saying: New record was created saying CAT : 1 (this cat refers to the land cover category right?); train_lbl value and GRASSRGB; 9- I place the values (water,1) for train_lbl and GRASSRGB 10- a box ask me for ASSUME DATA ENCODING AS a list appears. I should select which one? ASCII? 11- In trhis case I select ASCII and press SUBMIT. In the bottom of that box Record successfully updated appeared. 12- Insert centroid-1 and press submit in Form panel 13- Did the same for a few more polygons of cat 1 14- Did the same for a few more polygons of cat 2 and 3 15- Pressed Save and Exit Then I realized that at my v.digit panel command output I had this warning in there: *Could not set Tcl system encoding to 'ascii' (unknown encoding ascii)* *Unable to read vector map* Questions: 1- Is this the right procedure? 2- What encoding should I use? 3- Regarding GRASSRGB field, what information usually you place in there? Color code or color description? 4- For Supervised classification algorithms in GRASS, Vector type (accessible by using Query raster/vector tool) should be Boundary? Because, for instance, if I press a defined boundary, in some places I get Type: AREA and in others: Boundary (as shown bellow) (...) Map: vector06 Mapset: landsat Type: Boundary Id: 9 Layer: 1 Category: 3 Driver: dbf Database: C:\GRASS6\grassdata/North-Carolina/landsat/dbf/ Table: vector06 Key column: cat cat : 3 train_lbl : urban GRASSRGB : 3 (Thu Jun 17 10:34:04 2010) Command finished (0 sec) (Thu Jun 17 10:34:06 2010) v.what --q -a map=vecto...@landsat east_north=639279.744997,227730.071998 (...) Map: vector06 Mapset: landsat Type: Area Sq Meters: 2158429.741 Hectares: 215.843 Acres: 533.360 Sq Miles: 0.8334 Layer: 1 Category: 3 Driver: dbf Database: C:\GRASS6\grassdata/North-Carolina/landsat/dbf/ Table: vector06 Key column: cat cat : 3 train_lbl : urban GRASSRGB : 3 Thanks for your help and support Jenny On Wed, Jun 16, 2010 at 7:59 PM, Micha Silver mi...@arava.co.il wrote: On 06/16/2010 06:55 PM, Jenny Turner wrote: Hi Micha and GRASS-users/list Let me see if I get this straight 1- use v.digit (e.g. v.digit -n map=vector00 bgcmd=d.rast map=lsat5_1987_40 ) 2- Inside tcltk's digitizing menu select Digitize new Centroid 3- Leave layer as 1? 4- Change Mode to Manual Entry? All the above looks good. 5- And for each Centroid I select a Cat right? Yes, always using the *same* cat for training areas of the same type. (Is this the proceeding?) About this, I have two questions: a- I have to leave Insert New record into table checked right? Yes b- After I left Insert New record into table checked, Layer 1 and Cat 1 and Mode: Manual Entry, when I presss in the image a Error Window appeared stating: Database table for this layer is not defined. Where shall i define the database table for a new vectorial that I digitizing? Ah, v.db.addtable vector00 col=train_lbl varchar(16), GRASSRGB varchar(16) THis should setup an attrib table with two columns, one for a label and one for an RGB color for the training area. c- When I press the image a point is created but how can I define the boundaries of the centroid? First digitize the boundary. Here you need not enter any attributes. The cat values are also not relevant. NExt digitize a centroid inside each boundary. Give each centroid a cat value that matches the type of training area. So all the forests will be i.e cat=1, and all the urban will be i.e cat=10, etc. Thanks Jenny On Wed, Jun 9, 2010 at 7:06 PM, Micha Silver mi...@arava.co.il wrote: polygon to a certain category? This indeed is not obvious. What I've done is change the category mode from Next not used to Manual entry. This is set in the lower part of the v.digit window (in tcltk interface). Now each time you digitize a *centroid* for a training area, set its cat value to some standard value you choose. For example, you might decide: train_areacat forest 1 agriculture2 urban 9 water 10 Now for each training polygon which is a forest, set the cat value of its centroid to 1, agri areas will get value 2, and so on for all the centroids of all the training polygons. (Note that the cat values for the boundary lines are pretty much irrelevant in GRASS's vector model). When you're done digitizing training areas, you might add another column or two to the attrib table for the training areas
Re: [GRASS-user] Using Tck/tk vector digitizer for training areas
Hi Micha and GRASS-users/list Let me see if I get this straight 1- use v.digit (e.g. v.digit -n map=vector00 bgcmd=d.rast map=lsat5_1987_40 ) 2- Inside tcltk's digitizing menu select Digitize new Centroid 3- Leave layer as 1? 4- Change Mode to Manual Entry? 5- And for each Centroid I select a Cat right? (Is this the proceeding?) About this, I have two questions: a- I have to leave Insert New record into table checked right? b- After I left Insert New record into table checked, Layer 1 and Cat 1 and Mode: Manual Entry, when I presss in the image a Error Window appeared stating: Database table for this layer is not defined. Where shall i define the database table for a new vectorial that I digitizing? c- When I press the image a point is created but how can I define the boundaries of the centroid? Thanks Jenny On Wed, Jun 9, 2010 at 7:06 PM, Micha Silver mi...@arava.co.il wrote: polygon to a certain category? This indeed is not obvious. What I've done is change the category mode from Next not used to Manual entry. This is set in the lower part of the v.digit window (in tcltk interface). Now each time you digitize a *centroid* for a training area, set its cat value to some standard value you choose. For example, you might decide: train_areacat forest 1 agriculture2 urban 9 water 10 Now for each training polygon which is a forest, set the cat value of its centroid to 1, agri areas will get value 2, and so on for all the centroids of all the training polygons. (Note that the cat values for the boundary lines are pretty much irrelevant in GRASS's vector model). When you're done digitizing training areas, you might add another column or two to the attrib table for the training areas - v.db.addcol train_areas col=train_label varchar(16), RGBCOLOR varchar(16) then give each centroid cat value a suitable label and an RGB color combo, etc. Makes life easier down the line... This mail was received via Mail-SeCure System. ___ grass-user mailing listgrass-u...@lists.osgeo.orghttp://lists.osgeo.org/mailman/listinfo/grass-user This mail was received via Mail-SeCure System. -- Micha Silver Arava Development Co. +972-52-3665918http://www.surfaces.co.il ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user
[GRASS-user] Using Tck/tk vector digitizer for training areas
Greetings I'm using WinGRASS (from OSGEO4W installer) to do some supervised classification. Since I'm using winGRASS i.class function is not available neither wxGUI's vector digitizer. So I'm using tcl/tk but I'm a bit confused on using it. About the background images I used a discussion from last week regarding this but now.. How can I define a category? I mean, How can atribute a polygon to a certain category? ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user
Re: [GRASS-user] Georrectify and nearest neighbor
Hi Ok, but at least at georectify tool, Neares is used right? as well as copying a map from one mapset to another Regarding the link, there is no explanation regarding the tutorial video. I mean, step by step what I do have to do Thanks On Thu, Jun 3, 2010 at 10:59 PM, Markus Neteler nete...@osgeo.org wrote: On Wed, Jun 2, 2010 at 10:04 AM, Jenny Turner jennyturner1...@gmail.com wrote: My name is Jenny Turner and I've been using georrectify graphical module in WxpyTHon GRASS6.4. I saw a video of the proceeding but I would like to know if there is a description or a tutorial of what to do and how. Please take a look here: http://grass.osgeo.org/wiki/Georeferencing I have understand what I have to do but I need to explain to other persons (non-grass users :) ). Just one more thing: Is it used any resampling to the data or only nearest neighbor? Currently nearest neighbor. There may be a chance to have a choice in future in GRASS 7 to also use other methods. Please note that many resampling algorithms are already present in GRASS, just not yet integrated with the image processing part. Markus ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user
[GRASS-user] Georrectify and nearest neighbor
My name is Jenny Turner and I've been using georrectify graphical module in WxpyTHon GRASS6.4. I saw a video of the proceeding but I would like to know if there is a description or a tutorial of what to do and how. I have understand what I have to do but I need to explain to other persons (non-grass users :) ). Just one more thing: Is it used any resampling to the data or only nearest neighbor? Thanks Jenny ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user
[GRASS-user] Re: GRASS GUI assessment
First of all, let me thank you for all your help on describing and commentying the GRASS GUI's. Regarding Benjamin Ducke: - Wxpython Vector Digitizer is not yet available in Windows right? - Regarding QGIS, if a GRASS function is not available at QGIS I have heard that it's quite easy to put in QGIS (just edit a python file and XML as far as I remember). right? - And comparing QGIS and WXPYTHON: Which one is more stable and usable? Thanks, Jenny ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user
[GRASS-user] GRASS GUI assessment
Greetings I have sent an email a couple of weeks ago regarding a quick compare/assessment regarding GRASS versions. Now, I need to do the same but for GRASS GUI's. As far as I can see I can identify 3 GUIS: - Tcl/tk - WxPython - QGIS is there any comparison (Wikis, documents, ...) between those in order to assess and identify the most adequate one? Thanks Best regards Jenny ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user
Re: [GRASS-user] Importing Landsat images with different coordinates
Micha You are absolutely right. I was running it in a shared folder, and for some reason (not space beucase it's not full), i was getting that error. Now, I tried in a Desktop folder and it worked. Thanks 2010/5/26 Micha Silver mi...@arava.co.il This all looks fine. And you said that the destination tiff was smaller than the original? Can you display the tiff files? Maybe it's some disk problem (in the log you sent last time, it seems there was an error reading the file...)? On 05/25/2010 04:10 PM, Jenny Turner wrote: Sorry Micha My original Landsat image: Driver: GTiff/GeoTIFF Files: L71203033_0332601_B10.TIF L71203033_0332601_B10.aux L71203033_0332601_B10.rrd Size is 8161, 7141 Coordinate System is: PROJCS[WGS 84 / UTM zone 29N, GEOGCS[WGS 84, DATUM[WGS_1984, SPHEROID[WGS 84,6378137,298.257223563, AUTHORITY[EPSG,7030]], AUTHORITY[EPSG,6326]], PRIMEM[Greenwich,0], UNIT[degree,0.0174532925199433], AUTHORITY[EPSG,4326]], PROJECTION[Transverse_Mercator], PARAMETER[latitude_of_origin,0], PARAMETER[central_meridian,-9], PARAMETER[scale_factor,0.9996], PARAMETER[false_easting,50], PARAMETER[false_northing,0], UNIT[metre,1, AUTHORITY[EPSG,9001]], AUTHORITY[EPSG,32629]] Origin = (521999.9941792,4414200.000) Pixel Size = (30.000,-30.000) Metadata: AREA_OR_POINT=Point Image Structure Metadata: INTERLEAVE=BAND Corner Coordinates: Upper Left ( 522000.000, 4414200.000) ( 8d44'33.81W, 39d52'39.24N) Lower Left ( 522000.000, 4199970.000) ( 8d44'58.57W, 37d56'49.39N) Upper Right ( 766830.000, 4414200.000) ( 5d52'54.42W, 39d50'9.31N) Lower Right ( 766830.000, 4199970.000) ( 5d57'53.91W, 37d54'29.39N) Center ( 644415.000, 4307085.000) ( 7d20'4.58W, 38d54'2.83N) Band 1 Block=8161x1 Type=Byte, ColorInterp=Gray Min=0.000 Max=255.000 Minimum=0.000, Maximum=255.000, Mean=60.374, StdDev=43.530 Overviews: 2041x1786, 1021x893, 511x447, 256x224, 128x112, 64x56 Metadata: STATISTICS_MINIMUM=0 STATISTICS_MAXIMUM=255 STATISTICS_MEAN=60.374122874202 STATISTICS_MEDIAN=8.5986839406512e-287 STATISTICS_MODE=1.7801025473311e-306 STATISTICS_STDDEV=43.529818998446 LAYER_TYPE=athematic And the obtained file: Driver: GTiff/GeoTIFF Files: destination.TIF Size is 8229, 7218 Coordinate System is: PROJCS[ETRS89 / Portugal TM06, GEOGCS[ETRS89, DATUM[European_Terrestrial_Reference_System_1989, SPHEROID[GRS 1980,6378137,298.2572221010042, AUTHORITY[EPSG,7019]], AUTHORITY[EPSG,6258]], PRIMEM[Greenwich,0], UNIT[degree,0.0174532925199433], AUTHORITY[EPSG,4258]], PROJECTION[Transverse_Mercator], PARAMETER[latitude_of_origin,39.668258], PARAMETER[central_meridian,-8.1331084], PARAMETER[scale_factor,1], PARAMETER[false_easting,0], PARAMETER[false_northing,0], UNIT[metre,1, AUTHORITY[EPSG,9001]], AUTHORITY[EPSG,3763]] Origin = (-54187.411212536622770,23417.314857349840167) Pixel Size = (30.006036956002887,-30.006036956002887) Metadata: AREA_OR_POINT=Area Image Structure Metadata: INTERLEAVE=BAND Corner Coordinates: Upper Left ( -54187.411, 23417.315) ( 8d45'59.52W, 39d52'38.78N) Lower Left ( -54187.411, -193166.260) ( 8d44'57.96W, 37d55'35.73N) Upper Right ( 192732.267, 23417.315) ( 5d52'51.29W, 39d51'26.27N) Lower Right ( 192732.267, -193166.260) ( 5d56'29.99W, 37d54'28.07N) Center ( 69272.428, -84874.473) ( 7d20'4.40W, 38d54'3.73N) Band 1 Block=8229x1 Type=Byte, ColorInterp=Gray Hope this helps you (and mee :D ) Thank you Jenny 2010/5/25 Micha Silver mi...@arava.co.il On 05/25/2010 11:17 AM, Jenny Turner wrote: Yes Sure :) Attached to this email you can find the log I was thinking of just the results of gdalinfo. This way we can see if the tiff is being recognized, and if its CRS is is available in the tiff header. Here's an example: ~gdalinfo geodata/DSASC/landsat/p174r038_7t20020308_z36_nn10.tif Driver: GTiff/GeoTIFF Files: geodata/DSASC/landsat/p174r038_7t20020308_z36_nn10.tif Size is 8577, 7573 Coordinate System is: PROJCS[WGS 84 / UTM zone 36N, GEOGCS[WGS 84, DATUM[WGS_1984, SPHEROID[WGS 84,6378137,298.257223563, AUTHORITY[EPSG,7030]], AUTHORITY[EPSG,6326]], PRIMEM[Greenwich,0], UNIT[degree,0.0174532925199433], AUTHORITY[EPSG,4326]], PROJECTION[Transverse_Mercator], PARAMETER[latitude_of_origin,0], PARAMETER[central_meridian,33], PARAMETER[scale_factor,0.9996], PARAMETER[false_easting,50], PARAMETER[false_northing,0], UNIT[metre,1, AUTHORITY[EPSG,9001]], AUTHORITY[EPSG,32636]] Origin
Re: [GRASS-user] Importing Landsat images with different coordinates
Sorry Micha My original Landsat image: Driver: GTiff/GeoTIFF Files: L71203033_0332601_B10.TIF L71203033_0332601_B10.aux L71203033_0332601_B10.rrd Size is 8161, 7141 Coordinate System is: PROJCS[WGS 84 / UTM zone 29N, GEOGCS[WGS 84, DATUM[WGS_1984, SPHEROID[WGS 84,6378137,298.257223563, AUTHORITY[EPSG,7030]], AUTHORITY[EPSG,6326]], PRIMEM[Greenwich,0], UNIT[degree,0.0174532925199433], AUTHORITY[EPSG,4326]], PROJECTION[Transverse_Mercator], PARAMETER[latitude_of_origin,0], PARAMETER[central_meridian,-9], PARAMETER[scale_factor,0.9996], PARAMETER[false_easting,50], PARAMETER[false_northing,0], UNIT[metre,1, AUTHORITY[EPSG,9001]], AUTHORITY[EPSG,32629]] Origin = (521999.9941792,4414200.000) Pixel Size = (30.000,-30.000) Metadata: AREA_OR_POINT=Point Image Structure Metadata: INTERLEAVE=BAND Corner Coordinates: Upper Left ( 522000.000, 4414200.000) ( 8d44'33.81W, 39d52'39.24N) Lower Left ( 522000.000, 4199970.000) ( 8d44'58.57W, 37d56'49.39N) Upper Right ( 766830.000, 4414200.000) ( 5d52'54.42W, 39d50'9.31N) Lower Right ( 766830.000, 4199970.000) ( 5d57'53.91W, 37d54'29.39N) Center ( 644415.000, 4307085.000) ( 7d20'4.58W, 38d54'2.83N) Band 1 Block=8161x1 Type=Byte, ColorInterp=Gray Min=0.000 Max=255.000 Minimum=0.000, Maximum=255.000, Mean=60.374, StdDev=43.530 Overviews: 2041x1786, 1021x893, 511x447, 256x224, 128x112, 64x56 Metadata: STATISTICS_MINIMUM=0 STATISTICS_MAXIMUM=255 STATISTICS_MEAN=60.374122874202 STATISTICS_MEDIAN=8.5986839406512e-287 STATISTICS_MODE=1.7801025473311e-306 STATISTICS_STDDEV=43.529818998446 LAYER_TYPE=athematic And the obtained file: Driver: GTiff/GeoTIFF Files: destination.TIF Size is 8229, 7218 Coordinate System is: PROJCS[ETRS89 / Portugal TM06, GEOGCS[ETRS89, DATUM[European_Terrestrial_Reference_System_1989, SPHEROID[GRS 1980,6378137,298.2572221010042, AUTHORITY[EPSG,7019]], AUTHORITY[EPSG,6258]], PRIMEM[Greenwich,0], UNIT[degree,0.0174532925199433], AUTHORITY[EPSG,4258]], PROJECTION[Transverse_Mercator], PARAMETER[latitude_of_origin,39.668258], PARAMETER[central_meridian,-8.1331084], PARAMETER[scale_factor,1], PARAMETER[false_easting,0], PARAMETER[false_northing,0], UNIT[metre,1, AUTHORITY[EPSG,9001]], AUTHORITY[EPSG,3763]] Origin = (-54187.411212536622770,23417.314857349840167) Pixel Size = (30.006036956002887,-30.006036956002887) Metadata: AREA_OR_POINT=Area Image Structure Metadata: INTERLEAVE=BAND Corner Coordinates: Upper Left ( -54187.411, 23417.315) ( 8d45'59.52W, 39d52'38.78N) Lower Left ( -54187.411, -193166.260) ( 8d44'57.96W, 37d55'35.73N) Upper Right ( 192732.267, 23417.315) ( 5d52'51.29W, 39d51'26.27N) Lower Right ( 192732.267, -193166.260) ( 5d56'29.99W, 37d54'28.07N) Center ( 69272.428, -84874.473) ( 7d20'4.40W, 38d54'3.73N) Band 1 Block=8229x1 Type=Byte, ColorInterp=Gray Hope this helps you (and mee :D ) Thank you Jenny 2010/5/25 Micha Silver mi...@arava.co.il On 05/25/2010 11:17 AM, Jenny Turner wrote: Yes Sure :) Attached to this email you can find the log I was thinking of just the results of gdalinfo. This way we can see if the tiff is being recognized, and if its CRS is is available in the tiff header. Here's an example: ~gdalinfo geodata/DSASC/landsat/p174r038_7t20020308_z36_nn10.tif Driver: GTiff/GeoTIFF Files: geodata/DSASC/landsat/p174r038_7t20020308_z36_nn10.tif Size is 8577, 7573 Coordinate System is: PROJCS[WGS 84 / UTM zone 36N, GEOGCS[WGS 84, DATUM[WGS_1984, SPHEROID[WGS 84,6378137,298.257223563, AUTHORITY[EPSG,7030]], AUTHORITY[EPSG,6326]], PRIMEM[Greenwich,0], UNIT[degree,0.0174532925199433], AUTHORITY[EPSG,4326]], PROJECTION[Transverse_Mercator], PARAMETER[latitude_of_origin,0], PARAMETER[central_meridian,33], PARAMETER[scale_factor,0.9996], PARAMETER[false_easting,50], PARAMETER[false_northing,0], UNIT[metre,1, AUTHORITY[EPSG,9001]], AUTHORITY[EPSG,32636]] Origin = (603487.500,3623062.500) Pixel Size = (28.500,-28.500) Metadata: TIFFTAG_XRESOLUTION=72 TIFFTAG_YRESOLUTION=72 TIFFTAG_RESOLUTIONUNIT=2 (pixels/inch) AREA_OR_POINT=Point Image Structure Metadata: INTERLEAVE=BAND Corner Coordinates: Upper Left ( 603487.500, 3623062.500) ( 34d 6'16.43E, 32d44'25.92N) Lower Left ( 603487.500, 3407232.000) ( 34d 4'54.01E, 30d47'37.09N) Upper Right ( 847932.000, 3623062.500) ( 36d42'38.00E, 32d41'25.66N) Lower Right ( 847932.000, 3407232.000) ( 36d38'1.81E, 30d44'49.94N) Center ( 725709.750, 3515147.250
Re: [GRASS-user] Importing Landsat images with different coordinates
Hi Following this topic I have tested this methodology by converting from UTM 29N to a specific local coordinate system (EPSG:3763) (available at PROJ4) gdalwarp -t_srs EPSG:3763 band_landsat.TIF destination.tif With that I got a segmentation error, one TIFFReadDirectory a lot of TIFFWriteEncodedStrip and the destination.tif has a different size. Any tips of what migh be wrong? Thank you Jenny On Tue, May 18, 2010 at 8:42 PM, Micha Silver mi...@arava.co.il wrote: Jenny Turner wrote: Greetings Firman, but imagine that I have rasters in UTM for differente zones (e.g. huge countries like USA, etc). How can I produce an automatic importing tool that uses gdalwarp to project an image to WGS84 and import to GRASS? Using gdalwarp you don't have to specify the source coordinate system, *if* the originals have the correct projection defined in their headers. For Landsat images, this is most likely the case, so you should be able to get away with just setting the -t_srs option of gdalwarp and each tile will be re-projected to the WGS84 geographic CRS. Thanks for your help Best regards, Jenny On Mon, May 17, 2010 at 11:10 PM, Firman Hadi jalmibur...@gmail.commailto: jalmibur...@gmail.com wrote: Hi, As far as I know you must define the source coordinate system. I think it is pretty straightforward when you are using gdalwarp. If there are so many images, maybe you can create a shell script to run the process. Cheers, Firman Hadi Center for Remote Sensing - ITB Jl. Ganesha No. 10, Labtek IX-C, 3rd Floor Bandung - 40132 INDONESIA Phone: +62-22-2530701 Fax: +62-22-2530702 Website : http://crs.itb.ac.id ; www.sigro.org http://www.sigro.org Blog : http://jalmiburung.wordpress.com JepratJepret: http://jalmiburung.fotografer.net This mail was received via Mail-SeCure System. ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user This mail was received via Mail-SeCure System. -- Micha Silver Arava Development Co. +972-52-3665918 http://www.surfaces.co.il ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user
Re: [GRASS-user] Importing Landsat images with different coordinates
Greetings Firman, but imagine that I have rasters in UTM for differente zones (e.g. huge countries like USA, etc). How can I produce an automatic importing tool that uses gdalwarp to project an image to WGS84 and import to GRASS? Thanks for your help Best regards, Jenny On Mon, May 17, 2010 at 11:10 PM, Firman Hadi jalmibur...@gmail.com wrote: Hi, As far as I know you must define the source coordinate system. I think it is pretty straightforward when you are using gdalwarp. If there are so many images, maybe you can create a shell script to run the process. Cheers, Firman Hadi Center for Remote Sensing - ITB Jl. Ganesha No. 10, Labtek IX-C, 3rd Floor Bandung - 40132 INDONESIA Phone: +62-22-2530701 Fax: +62-22-2530702 Website : http://crs.itb.ac.id ; www.sigro.org Blog : http://jalmiburung.wordpress.com JepratJepret: http://jalmiburung.fotografer.net ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user
[GRASS-user] GRASS Script for importing data with different source_coordinates/projection
Greetings This almost a follow-up from a previous message that I sent regarding Landsat-7 image importing. I need to develop a script to import images from multi-SOURCE_COORDINATE_SYSTEM to defined target_COORDINATE_SYSTEM (WGS84). Since my input images come from different UTM zones it won't be possible to import images using a pre-defined gdalwarp command to convert coordinate/projection systems, since source coordinates are more than one. Firman made suggestion of using a GRASS Script to read header and get coordinate system. Has anyone developed something similar to this? Thyank you Best regards, Jenny ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user
Re: [GRASS-user] Importing Landsat images with different coordinates
Greetings Thanks It really worked but now I have one particualr question: I need to have an automatic proceeding but, in this case my input images are not in the same coordinate system (e.g. different UTM Zones). Can I apply this methdology without defining source_coordinate system? Thanks Jenny On Mon, May 10, 2010 at 12:21 PM, Firman Hadi jalmibur...@gmail.com wrote: On 10 May 2010, at 18:00, Jenny Turner wrote: Greetings I need to import Landsat images (r.in.gdal) to my location but, my Landsat images are in a specific Coordinate system (WGS84) that is different from mine. Could anyone explain me (briefly ) the steps to import this images? Thanks Jenny Dear Jenny, If your images are in TIFF, you can convert the coordinate system to the one that you have, with gdalwarp. As an example, if my coordinate system is UTM Zone 48S (32748) and my data is in lat-long (4326), I will use gdalwarp command like this: gdalwarp -t_srs EPSG:32748 -s_srs EPSG:4326 source.tif target.tif Cheers, Firman Hadi Center for Remote Sensing - ITB Jl. Ganesha No. 10, Labtek IX-C, 3rd Floor Bandung - 40132 INDONESIA Phone: +62-22-2530701 Fax: +62-22-2530702 Website : http://crs.itb.ac.id ; www.sigro.org Blog : http://jalmiburung.wordpress.com JepratJepret: http://jalmiburung.fotografer.net ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user
[GRASS-user] Importing Landsat images with different coordinates
Greetings I need to import Landsat images (r.in.gdal) to my location but, my Landsat images are in a specific Coordinate system (WGS84) that is different from mine. Could anyone explain me (briefly ) the steps to import this images? Thanks Jenny ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user
Re: [GRASS-user] Compare GRASS versions
Thanks Markus So the main reason for still supporting GRASS 6.4 is its stability right? For instance, do you plan to launch 6.5 as a end-product? Thanks Jenny 6.4 and 6.5 are rather close, 7 is starting to deviate. What I need to know is if there is any document that compares and assesses these 3 versions in order to characterizes each one. To my knowledge there is none. The Wiki would be the right place to start it. For instance GRASS6.5 is able to have dynamic wxpython scripts (e.g. v.krige) that is an interesting funcionality. Yes - but for stability reasons (in case you need it in production), you may consider 6.4. If functionality is missing and available in later versions, it could be discussed to backport it for 6.4.1 if not too complicated. Otherwise, take the next higher version number. Markus ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user
[GRASS-user] Compare GRASS versions
Greetings I'm preparing a small GIS tool to High School students to do some image processing, image vizualization etc and I'm thinking of using GRASS since is free and pretty much do what I need. But I realized that 6.4, 6.5 and 7 have different characteristics, functionalities. What I need to know is if there is any document that compares and assesses these 3 versions in order to characterizes each one. For instance GRASS6.5 is able to have dynamic wxpython scripts (e.g. v.krige) that is an interesting funcionality. Thanks Jenny ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user
re: [GRASS-user] Gorreferencing using GRASS
Greetings This was the right email at the right time. I was just doing some georreferencing and I have one or two questions regarding Martin video (awesome, by the way): 1- Right in the beggining, you have a raster, named mrizka, displayed in Display window1. I want to know: - How did you put the grid; - What is the purpose of this image? Because it is empty... It's totally blank 2- The group you create is the group of images you want to georrectify right? 3- the last panel, before, georrectification, regards image to display for GCP creatioN. This means what? This image is the one that I want to georrectify or the one that is already rectified? Because when it's opened, it already displays a grid and points/numbers. 4- XY are coordinates from not-georreferenced images and E/N are form the Georreferenced images right? This may sound lke dumb questions but I'm not use to this proceeding Thank you for your help Jenny ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user
[GRASS-user] Compare GRASS versions
Greetings I'm preparing a small GIS tool to High School students to do some image processing, image vizualization etc and I'm thinking of using GRASS since is free and pretty much do what I need. But I realized that 6.4, 6.5 and 7 have different characteristics, functionalities. What I need to know is if there is any document that compares and assesses these 3 versions in order to characterizes each one. For instance GRASS6.5 is able to have dynamic wxpython scripts (e.g. v.krige) that is an interesting funcionality. Thanks Jenny ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user
[GRASS-user] Can't initialize wxpython GRASS6.4
Greetings I have just installed in my old laptop a Ubuntu Karmic version of GRASS (6.4). I know it's not new but I'm getting an ODD error that it's not suppose to happen. I have placed nc_spm_08 data in /home/jenny/data folder and everytime I do sudo grass64 -wxpython When it's open GRASS GIS layer Manager I get an error window stating: Unable to get current geographic extent. Force quiting wxGUI. Please run manually g.region to fix the problem. This is weird because it's suppose to work right? Any ideas? Jenny ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user