). Because I generated the lines from high raster-values
(river-network), I could choose the cell with the highest value as a
direction for moving node and verts with 'v.edit'.
Thanks for patience, especially Hamish,
Achim
achim schrieb:
Hi Hamish,
thanks for your suggestions. It is very useful
, that solves it:
Ok, I would understand that a matlab struct is expected (the same as we can
get from r.out.mat).
But saving of the struct is not supported in format Matlab 4 (-v4):
Greets,
Achim
___
grass-user mailing list
grass-user@lists.osgeo.org
Hi,
maybe its the database syntax: try 'not cat=16'
Achim
Tim Holland schrieb:
Hello,
I am trying to add labels to a map, where most of the labels get printed
beside the shape (in this case, it is a point vector), while one label
(point #16) gets printed above the shape to prevent label
how can I tell the nodes at
line ends to snap to area boundaries.
Thank a lot in advance,
Achim
inline: select_thresh.png___
grass-user mailing list
grass-user@lists.osgeo.org
http://lists.osgeo.org/mailman/listinfo/grass-user
map=yourtablesmap ##automatically columns=cat INTEGER
will be ran and each category of the map gets an entry; default is
layer=1, so categories for layer one will be made; the tables name will
be yourtablesmap; for layer=2 yourtablesmap_2...
Greets,
Achim
PS:
little script to renew categories
a missing option.
Achim
___
grass-user mailing list
grass-user@lists.osgeo.org
http://lists.osgeo.org/mailman/listinfo/grass-user
: Having
raster-values in a column value is not clean, due to the fact, that I
have eg. to remove small trapped cells, calc center of gravity and more.
Thanks,
Achim
___
grass-user mailing list
grass-user@lists.osgeo.org
http://lists.osgeo.org/mailman/listinfo
at the ends of polylines
Cheers,
Achim
Giovanni Manghi schrieb:
Hi all,
a friend is asking me if there is a straight way in GRASS to make a
buffer like the one represented in this arcgis screenshot
http://evora.sigaberto.org/downloads/buffer.jpg
the buffer would be calculated along
...
Achim
Giovanni Manghi schrieb:
On Fri, 2009-10-09 at 19:32 +0200, Achim Kisseler wrote:
It seems, that there are line-segments and a buffer is made for each. Do
you mean that?
yes, I meant that
cheers
-- Giovanni --
___
grass-user
, because I used
three different kinds of objects coded as positive, negative and
negative multiplied with 10.
Achim
Hamish schrieb:
Achim wrote:
from a raster map I get vector polygons with r.to.vect.
As I use the -v flag, I get one category-number for each
value. This is what I want.
My
- update old table
Is there another, better and faster way? Shouldn't be measuring of
distances a common tool? Do I just don't find it? Do I understand
v.distance wrong (it just returns the cat of the nearest feature)?
Achim
___
grass-user mailing list
grass
table (with stored category-numbers)
-adding column length
-v.to.db option=length
-update original table
Achim
Achim Kisseler schrieb:
Hi grass-users,
I am going to measure the distance from one centroid to another for many
polygons in a layer in order to update the table with the distance
/html64_user/v.to.rast.html)
Try to make unique numbers from your string-column and convert them with
column=... and labelcolumn...
Achim
Dylan Beaudette schrieb:
Hi,
v.to.rast is letting me know that it does not accept text columns in an
attribute table. Has anyone figured out a simple approach
from the threshold of r.watershed.
Greets,
Achim
M S schrieb:
Upon first glance, r.basins.fill seems like a great tool to generate
subbasins within a watershed. Although I must not be understanding
something obvious and essential.
If one has to digitize in the internal ridges in a given
The man page seems to imply that internal divided need digitized?
Sorry, again I didn't get the point. Can you explain a bit more
detailed, please?
Achim
___
grass-user mailing list
grass-user@lists.osgeo.org
http://lists.osgeo.org/mailman/listinfo
(in an script)
Thanks,
Achim
___
grass-user mailing list
grass-user@lists.osgeo.org
http://lists.osgeo.org/mailman/listinfo/grass-user
Thanks for response, Markus!
I just did that a minute ago.
But I found out that I don't have to handle that exception in my case: I
do a r.patch further on, which works even if one raster-file is missing.
So there is no need of an null-map anymore.
Thanks again,
Achim
Markus Metz schrieb
Hi,
I guess the broken areas in the grass topology have dangles inside.
(zoom a broken area)
use v.clean:
-remove small areas and angles
-delete dangles
Hope it helps,
Achim
Νίκος Αλεξανδρής schrieb:
Felix Schalck wrote:
While still working on my huge map of Europe, I've noticed many
Hi,
is there a way in GRASS to build an accumulation map from a directions
map without a dem map?
Best regards,
Achim
___
grass-user mailing list
grass-user@lists.osgeo.org
http://lists.osgeo.org/mailman/listinfo/grass-user
I could find a tool like this
(inside GRASS it does not seem to be possible), I would be very grateful.
Best regards,
Achim
___
grass-user mailing list
grass-user@lists.osgeo.org
http://lists.osgeo.org/mailman/listinfo/grass-user
Hi,
why dont you import one map from one location to another? Some
information will be lost, that would be the same like r.resamp.stats.
http://grass.itc.it/grass64/manuals/html64_user/r.proj.html
Achim
Patrice Dumas schrieb:
Hello,
I have 2 locations, one with a lambert projection
Thanks, Jaroslaw for response,
there shoud be no problem if you have oryginal DEM
I know, but I do NOT have the DEM. Thats the problem.
Achim
Jarosław Jasiewicz schrieb:
Achim Kisseler pisze:
Hello again,
I am still searching a method to generate a flow accumulation map from a
flow
Hi Silvia,
thank you very much. I will use h.tca from JGRASS to get the desired.
Best,
Achim
Silvia Franceschi schrieb:
Hi Achim,
if you are trying to use an open source GIS for hydrological
applications consider to take a look at JGrass. There are a lot of
geomorphological modules.
Best
Dear Markus,
I follow Silvias suggestion and use h.tca in JGRASS.
BTW: Isn't is possible to integrate the horton-tools into grass?
Thanks,
Achim
Markus Neteler schrieb:
2009/11/19 Achim Kisseler a...@jupiter.uni-freiburg.de:
Thanks, Jaroslaw for response,
there shoud be no problem if you
Once I was in the same situation. Lambert projection is good for some
analyzes, because of the equal size of the grid cells.
To get the size of the grid cells, you could convert a N-S row of the
rastermap into vector polygons and get the size.
Not a very direct way, but it sould work.
Achim
:
r.mapcalc ns_col=if(col()==1 , row() , null())
Achim
___
grass-user mailing list
grass-user@lists.osgeo.org
http://lists.osgeo.org/mailman/listinfo/grass-user
Easiest way:
don't pipe!
Just make a new vectormap from v.extract and buffer that.
(I dont think piping is possible here and v.buffer cannot be restricted
on a selection)
Achim
Pablo Carreira schrieb:
Hi,
I want to do the following:
v.extract some lines from a vector and them
Moin grass-user,
I'd like to get the result of r.coin into a table-format to import it in
a database. Before I start fighting with sed and awk:
Do someone know a faster way?
Thanks in advance,
Achim
PS:
Also r.report or r.stat would work, but r.coin most looks like a table
it looks like
Hi all,
I wrote an awk-sed-script (for my specific situation) to do it...
Afterwards I recognized, that it is now general enough. So I used a
combination of r.mapcalc preprocessing and r.stats and combined the
result within the database.
All the best,
Achim
Achim Kisseler schrieb:
Moin grass
something wrong?
Thanks,
Achim
PS:
already did v.clean, r.mask -r (- should not be have an effect)
___
grass-user mailing list
grass-user@lists.osgeo.org
http://lists.osgeo.org/mailman/listinfo/grass-user
Hi Jitu,
a way would be to convert the lines to raster (v.to.rast) and calculate
the statistics (r.stats).
As every line has a unique id, the raster-lines will have.
Hope it helps,
Achim
Jitendra Kumar schrieb:
Hi,
v.what.rast allows to extract the values from the raster at the vector
Yes, this is quite tricky...
Have you seen the literature on eg.
http://www.gdf-hannover.de/media.php?id=0lg=en
?
Try the spearfish dataset to get an idea, how to begin.
Achim
Ross Benisch schrieb:
Hi, I have downloaded and unzip all the data for the 6.4 version of
grass for windows vista
I think v.dissolve corrupted it.
I solved it by making a raster map and then a vector map again.
Thanks, Markus and Pablo,
Achim
Achim Kisseler schrieb:
I made the grid with:
-%--
v.mkgrid --overwrite map=$1_grid_all grid=$north_south,$east_west
position=coor coor=$southwest_corner box
Hi all,
is it possible in grass to calculate the center of gravity?
Given two raster maps:
1) map where areas are defined (can be vector map, too)
2) raster map with gravity values
I cannot find something like this in grass.
Any suggestion welcome,
Achim
Is it possible to run r.sum for many areas AT ONCE?
r.sum does it for a whole region (or masked area)
like: one raster-map with areas, the other with values to sum?
Many thanks,
Achim
Achim Kisseler schrieb:
Hi all,
is it possible in grass to calculate the center of gravity?
Given two
[$3]+=$2*$4; d[$3]+=$4}
END{for (i in a){print i,a[i]/b[i],c[i]/d[i]}}
' temp_stats temp_stats2
-%--
temp_stats2 looks like
cat lon lat
14 49.123 -11.123
15 50.123 -10.123
...
(or whatever)
Cheers,
Achim
___
grass-user mailing list
grass-user
the same size.
Otherwise a lat-depending weighting could help... I will try it soon.
Thanks for your comment,
Achim
___
grass-user mailing list
grass-user@lists.osgeo.org
http://lists.osgeo.org/mailman/listinfo/grass-user
Easy:
just run them to get to know them or read the
-%--
DESCRIPTION
r.mapcalculator Provides a GUI frontend to r.mapcalc. Up to 5 maps can
be combined using simple expressions. Expert mode runs r.mapcalc for
more complex expressions.
-%--
Achim
Gilbert Ferrara schrieb:
Greetings all
I
Hi Will,
for lines you can use
v.to.db -p map=... option=start
-p to print your result
option=start to get start points
option=end to get end points
the result looks like
cat|x|y|z
...
without -p you can upload the results to the table with:
columns=x,y(,z)
Hope it helps,
Achim
William
.
Achim
William Temperley schrieb:
Hi Achim
Thanks - but what I really want is a list of unique nodes (i.e. where
three boundaries meet, there would be one node, with an id) and assign
this id to the relevant boundaries. Perhaps that isn't possible. It
might be easier to do this in PostGIS
that I'm new at GRASS :)
...keep on. Its worth!
Regards,
Achim
___
grass-user mailing list
grass-user@lists.osgeo.org
http://lists.osgeo.org/mailman/listinfo/grass-user
Hi,
I guess
doy = day of year
lat = latitude
While doy in (1-365)
And lat in (0-90)?
Should be. But why name??
Achim
Tim Michelsen schrieb:
Hello,
what do the inputs for i.sunshine in GRASS7 mean?
i.sunhours doy=name lat=name
Thanks in advance,
Timmie
Ah, so maybe you can have different doy on every raster-cell?
Generating a map full of 45-integers would give you valentines day hours
of sunshine everywhere on the map..., guess.
Achim
Tim Michelsen schrieb:
I guess
doy = day of year
I thought the same.
But the doc is strange:
doy=name
Did you do:
g.copy rast=boundary_raster,MASK
or alternatively:
r.mask input=boundary_raster -o
(-o for overwrite)
?
Tim Michelsen schrieb:
Hello,
I while ago I asked how to cut a small map from a larger one using the
region settings as boundary for the small map:
And more:
-%--
...The category values in the new raster output map layer will be the
same as those in the original, except that the resolution and extent of
the new raster output map layer will match those of the current
geographic region settings (see g.region)
-%--
Achim
Tim
is:
syntax error, unexpected NAME, expecting INTEGER or '-'
why does dir_x/dir_y does not return an integer here? Using int(dir_x)
also does not work.
Until now I use a r.mapcalc command with many if, checking the
direction of every pixel every time. But this is quite slow.
Achim
I see,
so variables are functions in this case.
I've been searching for a while, but I cannot find any tool like grasss
r.mapcalc that is a bit more flexible, eg. that can handle recursive tasks.
Thanks a lot,
Achim
Marco Lechner - FOSSGIS e.V. schrieb:
Hi Achim,
the neighborhood modifier
Hi,
can you post the hole command you used?
Achim
ciccp...@libero.it schrieb:
Hi everybody,
I am trying to import a raster file that is in a location SRTM, into another
location LatLong-WGS84 using r.in.gdal
I am inside the LatLong-WGS84 location, and once i imported the file i press
Hi,
labeling vectors does not work at my grass 6.5 (kubuntu64bit).
All settings in d.vect are accepted, but I see no labels.
Building a label map with v.label and loading it, I get an error: map
not found.
Do I miss a point?
Thanks for any help,
Achim
,
Achim
___
grass-user mailing list
grass-user@lists.osgeo.org
http://lists.osgeo.org/mailman/listinfo/grass-user
for hints in advance,
Achim
___
grass-user mailing list
grass-user@lists.osgeo.org
http://lists.osgeo.org/mailman/listinfo/grass-user
:
text labels for line features)
--%--
Error
Map not found.
--%--
right, that is a bug in the wxGUI. filed in the trac'er as #1014
Good to know...what a pitty..
Hamish
Achim
___
grass-user mailing list
grass-user@lists.osgeo.org
http
Hi,
that mean, as I understand, that you have to have a paleo DEM.
You have to ask the question:
On what data basis do I want to reconstruct paleo landscape?
Geomorphological changes cannot be reconstructed just with a present day
DEM.
Achim
On 17.03.2010 14:21, Marc Stonestreet wrote
.
Achim
On 17.03.2010 22:12, Maris Nartiss wrote:
One can't turn back time, still I also have been thinking about
playing with erosion models in backwards mode. If it's possible to
calculate current erosion rate and eroded material ammount, it should
be possible to add erroded material to DEM
By the way: it depends on the detail!
Pre glacial reconstruction of watersheds from isostatical
land-level-changes (ger. Ausgleichsbewegungen) were done for noth
america. I think it was last year.
Regards,
Achim
On 17.03.2010 22:52, Achim Kisseler wrote:
That erosion processes are highly
r.mapcalc result=if( isnull(map),0,map)
On 18.03.2010 09:46, Gary Nobles wrote:
Hi I have a raster
values are:
null
1
3
24
27
etc...
I need null = 0
As i want to add rasters together
so I want to do null=0
I'm sure this is simple
Thank you
Gary
=if( isnull(Bones_AVE),0,Bones_AVE)
Achim
___
grass-user mailing list
grass-user@lists.osgeo.org
http://lists.osgeo.org/mailman/listinfo/grass-user
.
Sometimes I try out QGIS, where importing often works out of the box,
but doesnt create the grass stable topology.
Achim
___
grass-user mailing list
grass-user@lists.osgeo.org
http://lists.osgeo.org/mailman/listinfo/grass-user
some more information about the data you are using?
And: can you post the commands you are used? (You can find the command
line equivalent for the gui in the tab Command Output)
Regards,
Achim
___
grass-user mailing list
grass-user@lists.osgeo.org
http
it be?
Thanks in advance for ideas,
Achim
___
grass-user mailing list
grass-user@lists.osgeo.org
http://lists.osgeo.org/mailman/listinfo/grass-user
that?
The folder has the same permissions as the others.
??
Achim
On 30.03.2010 10:37, Achim Kisseler wrote:
Hi all,
I have a problem with setting the mapset access for one mapset. All
other mapsets can be accessed without problems and from different
mapsets. I get:
--%--
ERROR mapset mapname not found
not the meters-miles thing. Exporting
and importing helps, but that's not a real solution.
Achim
___
grass-user mailing list
grass-user@lists.osgeo.org
http://lists.osgeo.org/mailman/listinfo/grass-user
Sometimes GUI output is not like commandline commands have to be. Maybe
its just the missing .
?
v.out.ogr -c input=poly...@mapset type=area dsn=PG:host=localhost
dbname=postgis user=users password=users format=PostgreSQL
v.out.ogr -c input=poly...@mapset type=area
not the meters-miles thing. Exporting
and importing helps, but that's not a real solution.
Achim
___
grass-user mailing list
grass-user@lists.osgeo.org
http://lists.osgeo.org/mailman/listinfo/grass-user
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@lists.osgeo.org
http://lists.osgeo.org/mailman
,
Stream burning: Thats what I thought about.
Further I suggest to weight the DEM with the contour information.
A guess how this could be done:
v.to.rast contourmap
r.fillnulls (to interpolate between contour lines raster map)?
r.mapcalc dem*contourrast
then r.watershed
Best regards,
Achim
PS: Check
Furthermore, in version 6.5 and up r.watershed becomes veryvery fast!
Am 21.05.2010 14:05, schrieb M S:
I believe that MFD in r.watershed is in versions 6.5 and up.
Mark
On Fri, May 21, 2010 at 7:27 AM, Kristian Foersterk.foers...@tu-bs.de wrote:
Hi,
I am using GRASS 6.4.0RC5+39438
Hi,
in v.split:
I wonder what units are be used as length in a lat-lon environment?
Are these degrees?
Why isn't there an option for m, km,..?
Thanks in advance,
Achim
___
grass-user mailing list
grass-user@lists.osgeo.org
http://lists.osgeo.org
the layer settings from the shell again ;-) )
Achim
___
grass-user mailing list
grass-user@lists.osgeo.org
http://lists.osgeo.org/mailman/listinfo/grass-user
Am 24.05.2010 23:35, schrieb Markus Neteler:
Please open an enhancement report in trac for this.
I'll do.
Thanks for response, Markus!
Achim
___
grass-user mailing list
grass-user@lists.osgeo.org
http://lists.osgeo.org/mailman/listinfo/grass-user
a table update and use 2- for identifying the right items.
I hope you like it,
Achim
Am 10.06.2010 19:40, schrieb Dave Roberts:
Friends,
I have an vector area map with numerous polygons with centroids. When I do
v.what.rast vector=A raster=B column=C
It appears to run, but doesn't actually enter
Hi,
try v.db.reconnect.all or v.db.connect for single vector layers.
Achim
Am 25.10.2010 18:54, schrieb Adam Dershowitz, Ph.D., P.E.:
This seems like it must be really easy and I am just missing something
I have a project but I accidentally created the sqlite database file in the
wrong
Hi Chethan,
working on a 64bit system (kubuntu, kernel 2.6.32), GRASS 6.4 seems
stable here.
On image processing: it'll be a bit faster, but it only makes sense, if
you really work a lot with (large) rasters.
Regards,
Achim
Am 30.10.2010 19:45, schrieb Chethan S:
Greetings everyone
Hi Sharon,
that sould be the radius, but:
due to changing distances on latitude, the circle is not a true circle,
as you can see when you reproject the map.
The buffer distance is in degree.
Hope, everything is correct,
regards,
Achim
Am 03.11.2010 06:44, schrieb Sharon M:
Hi,
I'm
Hi,
there are no areas, because there are no centroids. You can add
centroids via v.centroids.
You always have to choose a use-parameter, otherwise the module does
not know, which numbers it has to generate.
Achim
Am 05.11.2010 14:02, schrieb Luisa Peña:
| Number of points: 0
=your_column) or
use=cat to use the categogy value.
You can see the categories with v.category input=... option=report.
When you use an attribute, you have to have an table connected with the
vector layer.
Achim
___
grass-user mailing list
grass
Hi,
how do you use r.thin? Can you post the complete command?
Am 09.11.2010 06:25, schrieb charlie:
Hi. I am currently using Grass 6.4 on linux and beginning the the
process of vectorizing an old NOAA raster map with r.thin. The map is
approximately 12,000 pixels by 8,000 pixels wide.
Hi Nikhil,
e.g. you can do it with:
v.surf.rst
(http://grass.osgeo.org/gdp/html_grass64/v.surf.rst.html)
Does this help?
Achim
Am 09.11.2010 12:28, schrieb nikmor:
I am trying to import an contour file and trying to come up with an DEM map
from it. I tried it my way but end up with wrong
vector-map - layers - categories - database-connection: thats a bit
confusing, but manageble:
Am 09.11.2010 13:27, schrieb Luisa Peña:
1- How can I know that this particular vector data does/doesn't have a table?
v.db.connect -p map=your_map
2- Since this is a vector, created from a
Hi,
v.db.addcol should be there by default.
Maybe there has to be a table connected to the vector map!
This can bo done with 'v.db.addtable'.
Achim
Am 09.11.2010 12:55, schrieb Luisa Peña:
v.colors map=teste01_l...@regional column=cat color=aspect
I get this error:
Creating column
to improve the manual yet.
Regards,
Achim
___
grass-user mailing list
grass-user@lists.osgeo.org
http://lists.osgeo.org/mailman/listinfo/grass-user
Am 09.11.2010 15:51, schrieb Luisa Peña:
Ok I implented what is in 4r44284
1- v.db.connect -p map=teste_09_LULC
layers 1 table teste09_LULC in database E:\v3/LOCAL/Regional/dbf
through driver dbf with key cat
CVector map teste09_LULC is connected by:
Vector map should be connected by
, the contour
raster and the resulting dem.
Like Achim said, first you must convert the HEIGHT column to integer (I
did this with Openoffice Calc).
--%--
here he means the dbf-file
or you change the table after imported
(using SQL)
I think the output you got is correct. But when i try to change
Did you see:
http://grass.fbk.eu/statsgrass/grass_geostats.html
?
Am 10.11.2010 16:23, schrieb Steve Wangen:
Hello list -
There may be a simple solution to my problem, but I'm a bit new to
grass and haven't been able to track it down yet, so I was wondering
if anyone had suggestions. I'm
: 348.770 Number of faces: 0
|
| Number of centroids:213.666 Number of kernels:0
|
...that´ll be a HUGE vector file!
So it´s normal, that it takes some minutes.
Achim
___
grass-user mailing list
grass
with d.histogram.
Regards,
Achim
___
grass-user mailing list
grass-user@lists.osgeo.org
http://lists.osgeo.org/mailman/listinfo/grass-user
of the vector map. If the
vector map is not yet linked to any table, new a database link is
established based on the MAPSET database settings (see db.connect).
--%--).
Is that what you are looking for?
Regards,
Achim
___
grass-user mailing list
grass-user
Am 15.11.2010 10:47, schrieb Markus Metz:
Since the vector to be exported is a converted raster, apparently
without further modification, why not exporting the raster?
...I thought the same. What do you want to do?
Maybe you can tile the raster and create some or many vector maps?
Achim
had decimal numbers or strings instead of numbers?
Thanks
Gabriele
Strings in a raster map?
With decimal numbers it works the same way. Or you multiply the map with
number of digits and handle it as an integer map (maybe that does not
work because of limited bit information per pixel).
Achim
Am 15.11.2010 23:36, schrieb Gabriele N.:
Hi Achim.
Returning to the problem of the three maps in some cases remain the problem
(for example for the map C), because some numbers are 2 digits and another
with 1 digit.
For example, with values A = 12 B = 3323 C = 1 instead of 1233231 I have
Hi,
I did something similar with a river network.
1. Make a network out of it.
2. Count lines connected to each node. (v.net option=nreport)
3. Check the lines: if a related node has only one connected line,
delete the line.
Cheers,
Achim
Am 17.11.2010 13:49, schrieb Patrick_schirmer:
Hi
Hello Luisa,
(see below)
Am 17.11.2010 16:21, schrieb Luisa Peña:
Hello achim
1- I converted from raster to vector (r.to.vect)
r.to.vect -v input=teste01_l...@regional output=vector01 feature=area
Extracting areas...
Building topology for vector map vector01...
Registering primitives...
833
Well:
Q-GIS?
Achim
Am 25.11.2010 09:03, schrieb Giacomo Piva:
Ciao a tutti,
Hi all,
Does someone know a good GeoTIFF (and other image format) free and
multi-platform viewer software?
Thank you
___
grass-user mailing list
grass-user
Hi,
My guess:
7GB is a lot.
334611 primitives are many, but not that much,
but:
you got 18 millions vertices! (18459167)
That seems a lot.
Depending on how tangled your lines are, a potential mass of features
can be generated.
Just a guess,
best,
Achim
Am 26.11.2010 03:14, schrieb Pierre
Hi Sab,
v.stats _IS_ very slow (for many primitives).
Maybe you could analyse the database, where you can calculate the
statistics you like. For spatial statistics you can analyse uploaded
coordinates of centroids or sizes of area.
Hope it helps a bit,
Achim
Am 26.11.2010 12:32, schrieb
101 - 195 of 195 matches
Mail list logo