Micha Silver | 1 Apr 07:55 2010
Picon

Re: Re: grass-user Digest, Vol 47, Issue 52

On 04/01/2010 01:11 AM, Nathaniel Iwuchukwu wrote:
Thanks Micha,
I ran r.thin and got just one line. is there a way i can improve it?

Good, so we're making progress...
How many polygon areas did you create in the training map? Does it cover various portions of the roads?
On Wed, Mar 31, 2010 at 1:54 PM, Micha Silver <micha <at> arava.co.il> wrote:
Nathaniel Iwuchukwu wrote:

Also, when i convert from vector to raster, i don't see anything on the new raster map. its all blank. what should i do. thanks
I guess that would explain why you're getting "Training map has no classes".
Can you display the vector which covers sections of the road? Do you see filled polygons?
If not maybe you made boundaries with no centroids? In GRASS, area features must be composed of both a boundary *and* a centroid. In the digitizer there's a special button for adding a centroid. You need one inside each polygon boundary.


On Wed, Mar 31, 2010 at 11:20 AM, Nathaniel Iwuchukwu <ikechinedu <at> gmail.com <mailto:ikechinedu <at> gmail.com>> wrote:

   Thanks Micha,
   I will definately post it once i scale through it.
   i am getting an error when i run the i.gensigset on the training data.
   this is what i typed and the result:
   GRASS 6.4.0RC5 (sis):~/Desktop > i.gensigset trainingmap=traingis
   group=sislabgp subgroup=sislabsub signaturefile=sislab.sig
   Finding training classes...
    100%
   ERROR: Training map has no classes
   do you have any clue why i don't have a class?
        Thanks Man
   2010/3/31 Micha Silver <micha <at> arava.co.il <mailto:micha <at> arava.co.il>>


       On 03/31/2010 01:19 AM, Nathaniel Iwuchukwu wrote:
       Thanks Micha,
       this looks good but am still working on it. would you know
       the RGB for a road? from the link you
       No, I can't tell what the RGB values for a road would be, and
       it will change from image to image, and from road to road.
       sent me, i will need it to capture the pixel values of road
       polygons that i am digitizing.
       You *don't* need to find the actual pixel values for the road
       polygons. You just digitize strips that cover the road, and
       feed these areas (after converting to a raster) as the
       training map. The gensigset module creates the collection of
       pixels values that match "roads"

       thanks Micha. you really gave me a lead.
       Good. Post back to the list when you get it worked out.


                Nathaniel

       On Tue, Mar 30, 2010 at 4:26 PM, Micha Silver
       <micha <at> arava.co.il <mailto:micha <at> arava.co.il>> wrote:

           Nathaniel Iwuchukwu wrote:

               sorry my typo. i mean classification.

           I didn't intend to insult. Sometimes typos come out
           funny, and I can't resist pointing it out...

               but am looking for procedure on road extraction from
               satellite imagery with high resolution

           OK, you'll probably want to use "supervised"
           classification. This means giving the classification
           algorithm a "training" layer. So the process involves
           1-  First create a vector layer using the grass digitizer
           with small polygons covering some sections of the roads
           in the satellite image. Convert the vector to a raster
           (v.to.rast). This will be the training map.
           2- Now group the satellite bands into an imagery group
           (i.group)
           3- Run i.gensigset with the training raster and the
           imagery group.
           4- Now use i.smap with the signature file and imagery
           group as input, and a classification map as output.
           5- If the results look good, run r.thin on this raster,
           and then r.to.vect to see what you get.

           Here's a more detailed explanation that you might find
           helpful:
           http://www.surfaces.co.il/?p=3

           Regards,
           Micha



       --         Micha Silver
       http://www.surfaces.co.il/
       Arava Development Co. +972-52-3665918




This mail was received via Mail-SeCure System.


--
Micha Silver

Arava Development Co. +972-52-3665918


This mail was received via Mail-SeCure System.


-- Micha Silver http://www.surfaces.co.il/ Arava Development Co. +972-52-3665918
_______________________________________________
grass-user mailing list
grass-user <at> lists.osgeo.org
http://lists.osgeo.org/mailman/listinfo/grass-user
Szabó Balázs | 1 Apr 09:32 2010
Picon

Installing GIPE

Hi list!

Is it possible install gipe to grass gis 7.0 latest svn?

When I run install_yanns_grass.sh (as $sudo sh install_yanns_grass.sh)
it return with lot of error (
$main.c: In function ‘main’:
$main.c:164: warning: assignment makes pointer from integer without a
cast 
etc etc....
)
However the script copys everithing on the right dir. 

This is my install_yanns_grass.sh settings:
GRASSSVNDIR=/usr/local/grass-7.0.svn
ADDONSVNDIR=/home/baa/grass-addons

----------------------------------------------
I tryed with g.extension manager, but this command won't work: 

g.extension extension=r.hydro.CASC2D \
svnurl=https://svn.osgeo.org/grass/grass-addons/gipe/

Someone could help me?

Thanks

Balazs

_______________________________________________
grass-user mailing list
grass-user <at> lists.osgeo.org
http://lists.osgeo.org/mailman/listinfo/grass-user

Achim Kisseler | 1 Apr 11:05 2010
Picon

Re: ERROR mapset <mapname> not found


>> Could it be that the mapset is locked? If yes, where can I find out and
>> how can I change that?
>
> A mapset is locked by creating a ".gislock" file in the mapset
> directory (note the leading "."). But the error which you report
> doesn't sound like a locked mapset.

right, thats not the problem.

>
> It would help if you specified the exact command which you run and the
> exact error message (copied and pasted, not re-typed) which it
> generates.

here I use the gui: Config -> GRASS working enviroment -> mapset access

>
> Are you sure that the directory is actually a mapset, i.e. does it
> contain the WIND and VAR files? Do these files have the correct
> ownership and permissions?

Yes, I'm sure.

 >>>>>>>>EDIT>>>>>>>>

Now I used the command line:

--%<--
g.mapsets addmapset=mapname
--%<--

and it worked!

But now: the gui doesnt work anymore:

--%<--
Traceback (most recent call last):
   File "/usr/local/grass-6.5.svn/etc/wxpython/wxgui.py",
line 257, in OnMapsets

dlg = preferences.MapsetAccess(parent=self, id=wx.ID_ANY)
   File "/usr/local/grass-6.5.svn/etc/wxpython/gui_modules/pr
eferences.py", line 2137, in __init__

self.mapsetlb.CheckItem(self.all_mapsets.index(mset), True)
ValueError
:
list.index(x): x not in list
--%<--

BTW:
g.mapsets -s
shows the tcltc dialog!?

Thats ok for me, because it works, but WHY does this happen?

Thanks everyone very much for the hints,
Achim

_______________________________________________
grass-user mailing list
grass-user <at> lists.osgeo.org
http://lists.osgeo.org/mailman/listinfo/grass-user

Micha Silver | 1 Apr 14:23 2010
Picon

Re: Re: grass-user Digest, Vol 47, Issue 52

On 04/01/2010 02:52 PM, Nathaniel Iwuchukwu wrote:
i created 11 polygons from different portions of the road.
So if I understand correctly, the procedure is working, but you're not getting the results you'd expect?
Maybe you could supply to the grass-users list some URL's of images of what the satellite rasters look like, with your training polygons, and the resulting roads. Then perhaps someone will be able to suggest how to get more complete results. 


 
2010/4/1 Micha Silver <micha <at> arava.co.il>
On 04/01/2010 01:11 AM, Nathaniel Iwuchukwu wrote:
Thanks Micha,
I ran r.thin and got just one line. is there a way i can improve it?

Good, so we're making progress...
How many polygon areas did you create in the training map? Does it cover various portions of the roads?
On Wed, Mar 31, 2010 at 1:54 PM, Micha Silver <micha <at> arava.co.il> wrote:

Nathaniel Iwuchukwu wrote:

Also, when i convert from vector to raster, i don't see anything on the new raster map. its all blank. what should i do. thanks
I guess that would explain why you're getting "Training map has no classes".
Can you display the vector which covers sections of the road? Do you see filled polygons?
If not maybe you made boundaries with no centroids? In GRASS, area features must be composed of both a boundary *and* a centroid. In the digitizer there's a special button for adding a centroid. You need one inside each polygon boundary.


On Wed, Mar 31, 2010 at 11:20 AM, Nathaniel Iwuchukwu <ikechinedu <at> gmail.com <mailto:ikechinedu <at> gmail.com>> wrote:

   Thanks Micha,
   I will definately post it once i scale through it.
   i am getting an error when i run the i.gensigset on the training data.
   this is what i typed and the result:
   GRASS 6.4.0RC5 (sis):~/Desktop > i.gensigset trainingmap=traingis
   group=sislabgp subgroup=sislabsub signaturefile=sislab.sig
   Finding training classes...
    100%
   ERROR: Training map has no classes
   do you have any clue why i don't have a class?
        Thanks Man
   2010/3/31 Micha Silver <micha <at> arava.co.il <mailto:micha <at> arava.co.il>>


       On 03/31/2010 01:19 AM, Nathaniel Iwuchukwu wrote:
       Thanks Micha,
       this looks good but am still working on it. would you know
       the RGB for a road? from the link you
       No, I can't tell what the RGB values for a road would be, and
       it will change from image to image, and from road to road.
       sent me, i will need it to capture the pixel values of road
       polygons that i am digitizing.
       You *don't* need to find the actual pixel values for the road
       polygons. You just digitize strips that cover the road, and
       feed these areas (after converting to a raster) as the
       training map. The gensigset module creates the collection of
       pixels values that match "roads"

       thanks Micha. you really gave me a lead.
       Good. Post back to the list when you get it worked out.


                Nathaniel

       On Tue, Mar 30, 2010 at 4:26 PM, Micha Silver
       <micha <at> arava.co.il <mailto:micha <at> arava.co.il>> wrote:

           Nathaniel Iwuchukwu wrote:

               sorry my typo. i mean classification.

           I didn't intend to insult. Sometimes typos come out
           funny, and I can't resist pointing it out...

               but am looking for procedure on road extraction from
               satellite imagery with high resolution

           OK, you'll probably want to use "supervised"
           classification. This means giving the classification
           algorithm a "training" layer. So the process involves
           1-  First create a vector layer using the grass digitizer
           with small polygons covering some sections of the roads
           in the satellite image. Convert the vector to a raster
           (v.to.rast). This will be the training map.
           2- Now group the satellite bands into an imagery group
           (i.group)
           3- Run i.gensigset with the training raster and the
           imagery group.
           4- Now use i.smap with the signature file and imagery
           group as input, and a classification map as output.
           5- If the results look good, run r.thin on this raster,
           and then r.to.vect to see what you get.

           Here's a more detailed explanation that you might find
           helpful:
           http://www.surfaces.co.il/?p=3

           Regards,
           Micha



       --         Micha Silver
       http://www.surfaces.co.il/
       Arava Development Co. +972-52-3665918




This mail was received via Mail-SeCure System.


--
Micha Silver

Arava Development Co. +972-52-3665918


This mail was received via Mail-SeCure System.


-- Micha Silver http://www.surfaces.co.il/ Arava Development Co. +972-52-3665918


This mail was received via Mail-SeCure System.


-- Micha Silver http://www.surfaces.co.il/ Arava Development Co. +972-52-3665918
_______________________________________________
grass-user mailing list
grass-user <at> lists.osgeo.org
http://lists.osgeo.org/mailman/listinfo/grass-user
M S | 1 Apr 14:55 2010
Picon

Re: 3D lines to raster with elevation

(answering self)
v.to.points seems to be a nice solution, as going from 3D line to raster, and capturing Z of line didnt output raster features.

Data went from 3D lines, to 3D points (verticies), to raster with Z.

On Mon, Mar 29, 2010 at 12:15 PM, M S <mseibel <at> gmail.com> wrote:
Should v.to.rast be able to take a 3D line, and convert to raster with applicable Z values brought over?

This command:
v.to.rast input=line3d <at> PERMANENT output=line3d use=z type=line layer=1 rows=4096 --overwrite

outputs this:
Loading data...
Reading features...
Writing raster map...
Converted areas: 0 of 0
Converted points/lines: 0 of 1
v.to.rast complete. 

When I query the 3D line in GRASS, it has Z values:
Map: line3d
Mapset: PERMANENT
Type: Line
Id: 1
Length: 1619.419511
Line height min: 97.439897
Line height max: 98.149869
Layer: 1
Category: 1


Mark

_______________________________________________
grass-user mailing list
grass-user <at> lists.osgeo.org
http://lists.osgeo.org/mailman/listinfo/grass-user
Brian Walters | 1 Apr 17:52 2010
Picon

i.atcorr

Hello users,
I'm interested in using i.atcorr to do atmospheric corrections.  However, the link on the GRASS add-ons page no longer works.  Is i.atcorr still available?  If so, where can I download it for version 6.4?
Thanks for any response,

Brian

_______________________________________________
grass-user mailing list
grass-user <at> lists.osgeo.org
http://lists.osgeo.org/mailman/listinfo/grass-user
Markus Neteler | 1 Apr 18:14 2010

Re: i.atcorr

Hi Brian,

On Thu, Apr 1, 2010 at 5:52 PM, Brian Walters <bfwalters83 <at> gmail.com> wrote:
> Hello users,
> I'm interested in using i.atcorr to do atmospheric corrections.  However,
> the link on the GRASS add-ons page no longer works.

where did you find the link so that we can correct it?

>  Is i.atcorr still available?  If so, where can I download it for version 6.4?

i.atcorr is part of "main" GRASS now, it is included in 6.4 or later.

Note: if you have self-compiled GRASS, you need to configure before
compilation with --with-cxx since i.atcorr is written in C++.

Hope this helps,
Markus
_______________________________________________
grass-user mailing list
grass-user <at> lists.osgeo.org
http://lists.osgeo.org/mailman/listinfo/grass-user

Bulent Arikan | 1 Apr 19:41 2010
Picon

'v.surf.rst' VS. 'r.resamp.rst'


Hi,

I'd like to thank everyone for suggestions on my previous question about the interpolation of DEMs from 30 to 10 m res. and  problems I had with 'v.surf.rst' module. I am still experimenting interpolation with it and i hope to get a good interpolated DEM finally. Hamish suggested I use 'r.resamp.rst' instead and I did (with the default values). The results look really impressive; it took equally long time for a single DEM (about 90 minutes). I used the resultant DEM to make slope and aspect maps; nothing seems weird. I also used the mapcalc module to check the difference between 30m and 10m DEMs and it is within -1 and 1 m range. However, I used 'r.resamp.rst' on one imagery only.

So, I am writing to ask for suggestions about which module to use for 30m to 10m interpolation. I understand that 'v.surf.rst' and 'r.resamp.rst' use the same approach (i.e., spline tension) but in different ways (which is not surprising since GRASS offers a lot of alternative ways of manipulating data) and given my limited experience with GIS, I will appreciate suggestions from the more experienced users.       

Thanks,

* I am running a slightly old version (late 2009) of GRASS 6.5 on Mac OS 10.5.8

--
BÜLENT ARIKAN
School of Human Evolution and Social Change
Arizona State University
Tempe - AZ
85287-2402
_______________________________________________
grass-user mailing list
grass-user <at> lists.osgeo.org
http://lists.osgeo.org/mailman/listinfo/grass-user
Matthew Mulbrandon | 2 Apr 03:40 2010
Picon

RE: grass pydev or emacs

Hello,

ran this program from grass_python wiki http://grass.osgeo.org/wiki/GRASS_and_Python
import os
import grass.script as grass
env = grass.gisenv()
print env
#List all environment variables
for key, value in os.environ.items(): print key, value
#Print region extent
r = grass.read_command("g.region", flags='p' )
print r

on second line I keep getting error

Traceback (most recent call last):
File "/home/matthew/workspace/dfdsfds/grass.py", line 2, in <module>
import grass.script as grass
File "/home/matthew/workspace/dfdsfds/grass.py", line 2, in <module>
import grass.script as grass
ImportError: No module named script

I checked variables and made make sure grass.script is there. I can not find problem. Any clues would be helpful.

Ubuntu 9.10, grass 6.4rc5, pydev , eclipse  Version: 3.5.2

Thanks,

Matthew


Hotmail: Trusted email with Microsoft’s powerful SPAM protection. Sign up now.
Hotmail has tools for the New Busy. Search, chat and e-mail from your inbox. Learn More.
_______________________________________________
grass-user mailing list
grass-user <at> lists.osgeo.org
http://lists.osgeo.org/mailman/listinfo/grass-user
Martin Landa | 2 Apr 10:47 2010
Picon

Re: RE: grass pydev or emacs

Hi,

2010/4/2 Matthew Mulbrandon <mulbran <at> hotmail.com>:
> Traceback (most recent call last):
>   File "/home/matthew/workspace/dfdsfds/grass.py", line 2, in <module>
>     import grass.script as grass
>   File "/home/matthew/workspace/dfdsfds/grass.py", line 2, in <module>
>     import grass.script as grass
> ImportError: No module named script
>
> I checked variables and made make sure grass.script is there. I can not find
> problem. Any clues would be helpful.
>
> Ubuntu 9.10, grass 6.4rc5, pydev , eclipse  Version: 3.5.2

Try to update to RC6. Martin

--

-- 
Martin Landa <landa.martin gmail.com> * http://gama.fsv.cvut.cz/~landa
_______________________________________________
grass-user mailing list
grass-user <at> lists.osgeo.org
http://lists.osgeo.org/mailman/listinfo/grass-user


Gmane