I'm having two problems when exporting the output of r.stream.order
to a shapefile. I have attached the created table (saved in my setup in sqlite) to the 'streams' map thru layer 2. Then v.db.select streams layer=2 shows all the details of strahler order, prev_str, etc. First problem: When I try to do v.out.ogr, the column header named "next_stream" fails because it's 10 characters long, too long for a shapefile dbf. As a result all attributes are missing. Here's a patch to io.c in r.stream.order I put in place to correct this: [micha@SL6 r.stream.order]$ diff -u io.c.orig io.c --- io.c.orig 2012-08-28 12:21:46.020275045 +0300 +++ io.c 2012-08-28 12:22:17.951382882 +0300 @@ -292,7 +292,7 @@ /* table definition */ char *tab_cat_col_name = "cat integer"; char *tab_stream = "stream integer"; - char *tab_next_stream = "next_stream integer"; + char *tab_next_stream = "next_str integer"; char *tab_prev_streams; char *tab_strahler = "strahler integer"; char *tab_horton = "horton integer"; @@ -300,7 +300,7 @@ char *tab_hack = "hack integer"; char *tab_length = "length double precision"; char *tab_cumlength = "cum_length double precision"; - char *tab_stright = "stright double precision"; + char *tab_stright = "straight double precision"; char *tab_fractal = "fractal double precision"; char *tab_distance = "out_dist double precision"; char *tab_topo_dim = "topo_dim integer"; (Also corrects a small typo in the column name "stright") After recompiling, then rerunning the r.stream.order addon, the sqlite table now has all column headers <= 10 characters suitable to shapefile. Problem 2, where I'm stuck. I export the vector map using : v.out.ogr streams layer=2 dsn=streams.shp The process finishes OK, and the shape is created with all attrib columns, but *all* values are NULL (including the cat). Any ideas what I've done wrong? Thanks, Micha |
_______________________________________________ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user