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
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
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
of many programs to create
an empty one (eg. there is an firefox add-on that can do that).
Store the db (eg. myname.sqlite) in you gis-folder and connect (in
this example case) to:
$GISDBASE/myname.sqlite
Best regards,
Achim
Any other suggestion?
2010/11/16 Achim Kisseler a...@jupiter.uni
Am 15.11.2010 16:41, schrieb Luisa Peña:
In this email,
http://lists.osgeo.org/pipermail/grass-user/2010-November/058644.html,
you suggested me to use sqlite as database. I'm willing to try that. How
can I do that with a vector created from r.to.vect with categories
values in cat column
Hi
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
Am 15.11.2010 23:10, schrieb Gabriele N.:
Excellent Alex:). The maps are as follows:
A - value 0, 1, 2, ... 19
B - value 0, 1001, 1002, 1003 3351
C - value from 0, 1, 2, 3, 4 .. 15
So, I did as you suggested:
r.mapcalc output = (( A *100)+( B *100)+ C )
right?
And if I
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
of the DEM after exagerating it.Now I am looking forward to drape an aerial
image of the same area on the DEM that I have got. In order to do this I
loaded the DEM map and then draped the imported aerial map (imported using
r.in.gdal, but that single jpeg image(satimage.jpg) got imported into 3
Hi,
Am 11.11.2010 16:08, schrieb Luis Lisboa:
| Type of Map: vector (level: 2)
|
|
|
| Number of points: 0 Number of areas:
213.680 |
| Number of lines:0 Number of islands:
73.649 |
| Number of boundaries:
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
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
Am 09.11.2010 14:07, schrieb Martin Landa:
why confusing?
Vector object categories and attribute management
http://grass.osgeo.org/grass64/manuals/html64_user/vectorintro.html
yes, everything is written there, but:
as I started, when it came to these
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
Hi Nikhil,
And one more thing is
that how have you converted the height to integer please give me
details of
this too. I am really very thankful for your time and effort.
There are many ways. For example like Micha said:
--%--
I'm including a mapset bombay with the contour vecto, the
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
Am 05.11.2010 15:31, schrieb Luisa Peña:
v.to.rast --overwrite input=teste01_centr...@brazil
output=teste01centr...@brazil use=val value=5
But I do have 1 (last) question:
When I use the Query pixel value feature, I get this:
568209.941176|9248338.294118||5|Value 5 Does this means that each
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 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!
As
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 Dave,
what I did in this case, was, as you said, to convert the centroids to
points and then run v.what.rast.
One possibility to get the information to your table:
1- Convert the polygons (the categories) to raster with r.to.vect.
2- Do a v.what.rast with this raster-area map.
3- Make a
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
Dear list,
I know that it is not possible in grass65, but can one export the
graphics display to svg or eps in grass7 (from gui-display)?
Alternatively I would like to load my workspace settings to a xmon in
order to export from there. Is that possible?
(I really dont want to rebuild the
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
On 01.05.2010 07:46, Jarek Jasiewicz wrote:
Interested question, but I think the easiest answer is carving or stream
burining (r.stream.carve). Next reproduce r.watershed or
r.stream.extract you shall receive (almost) ideally confidential network
with supplied
HTH
Jarek
Hi Nikos,Jarek,
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
Hi Randy,
I need to create 6 regional maps of the continental U.S., north and
south by west, central, and east. My first thought was to create 7
locations (one for the entire continental US) but my second thought was
to create just one location with the default region as the entire
continental
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
Hi Randy,
I need to create 6 regional maps of the continental U.S., north and
south by west, central, and east. My first thought was to create 7
locations (one for the entire continental US) but my second thought was
to create just one location with the default region as the entire
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
--%--
I already checked the folders, but I cannot find any abnormalities.
What could it be?
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
Hi,
Error:
Assertion failed: nrows * ncols == str-stream_len(), file
grass2str.h,line 144
This application has requested the Runtime to terminate it in an unusual
way. Please contact the application's support team for more information.
about the error: I dont know, but:
can you give
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
Hi Gary,
so
r.mapcalc result=if( isnull(Bones_AVE),0,outputmap)
Is that right, no need for A=Bones_AVE?
could the output map be the original? Or does it have to be a new raster?
r.mapcalc result=if( isnull(Bones_AVE),0,Bones_AVE)
and yes, it can be the original
r.mapcalc Bones_AVE=if(
On 18.03.2010 17:35, sarah moore wrote:
To fellow grass gis users:
I am having trouble with vector works data going into grass gis. Has
anyone else ran into this problem? I have not tried using syntax,
instead I used the pull down files: File-Import Vector map-DXF import-
resulting in illegal
Hello,
I wonder if anyone else has the same problem or if the problem is known:
Using wxGUI the monitor has problems with showing all selected layers in
the layer list. Very often some are missing.
That happens the more often, the more layers are in the list.
To handle this, I can change
Hello Hamish,
now it works. Thank you!
Explanation:
did you set attrcol= to the name of the column containing the text from
the Labels tab?
yes
does disp=cat work?
Yes, that works.
BUT:
I didnt set disp=attr
AND
attrcol at the same time.
Thats what I missed!
(What does not work:
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:
That erosion processes are highly non-linear. Running physical based
erosion models in backwards-mode is not possible I would bet.
I guess one have to work with ground layers, geological land rising and
lowering, climatic changes and so on.
But well-wrought concept ideas are always welcome.
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
Hi Hamish,
thanks a lot for reply!
display=attr ?
thats what I did, but I dont get even simple text
Building a label map with v.label and loading it, I get an
error: map not found.
displaying fancy paint labels is done by d.labels not by d.vect.
(third button from right in wxGUI layer
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
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
the
Hello grass-users,
I have a problem running a r.mapcalc procedure:
given two maps
dir_x,
dir_y
with 0,1,-1 entries
I want to run:
r.mapcalc the_map=if(isnull(the_map),the_map[dir_x,dir_y],null())
in order to make a r.water.outlet on many many areas at the same time.
The error I get is:
doesn't seem to accept functions (you also
could try rand(), ...)
Marco
Achim Kisseler schrieb:
Hello grass-users,
I have a problem running a r.mapcalc procedure:
given two maps
dir_x,
dir_y
with 0,1,-1 entries
I want to run:
r.mapcalc the_map=if(isnull(the_map),the_map[dir_x,dir_y],null
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
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
. I need to be able to run
recursive queries in PostGIS to pick up connected boundaries.
Cheers
Will
2010/1/15 Achim Kisseler a...@jupiter.uni-freiburg.de
mailto:a...@jupiter.uni-freiburg.de
Hi Will,
for lines you can use
v.to.db -p map=... option=start
-p to print your
Hi,
What you are trying to do is a on-the-fly projection. GRASS does not
do this, because that doesn't work properly (as I understand, but I
don't know why).
What you have to do is to make another location with the different
coordinate system and import your data there.
Then go to your
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
beware it will only be valid for small ranges of latitude.
otherwise the weights of the longitudes at high latitudes will
have more weight in the average than they should.
Hi Hamish,
yes, you are right.
In this case the area had to be projected in a projection, where every
raster-cell has
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
Hi all,
I found a solution for the problem of getting the cog for each area with
weight of a value map. Its simply r.stats and one! awk-command. And its
quite fast!
-%--
#
## get stats into file
r.stats -c -n -g input=$1_areas,$1_values output=temp_stats
#
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
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
Hi,
I am using
v.to.db map=af_grid option=area columns=area type=centroid units=kilometers
(in GRASS65 opensuse 11.1 64)
to calculate the area size. Some of the areas are very wrong. Its a grid
in lat-lon, where areas should be nearly equal in size.
Is this a known bug? I don't guess so.
Do I
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
Hello again,
I am still searching a method to generate a flow accumulation map from a
flow direction map.
I got the direction map from an arcgis user. In arcgis the module
FLOWACCUMULATION(FLOW_DIR) would do this, but I want to use an open
source solution.
If anyone could give me a hint, where
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 with a
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
Patrice Dumas schrieb:
On Thu, Nov 19, 2009 at 01:17:13PM +0100, Achim Kisseler wrote:
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
Problem
Hi,
use r.to.vect to map every cell to a polygon, project it to the lambert
map, do the calculation on the polygon, reproject it back to the
you dont need to project it to lambert:
just run v.to.db -p option=area map=...
then you get the size for every lines cell.
To get the raster-map:
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
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
Dear grass users,
has anyone an idea how to handle v.to.rast, if vector-file is empty? -
it produces an error and no map is created:
Loading data...
ERROR: Database connection not defined for layer 2
I want an null-raster-map instead.
I created the empty vector map with v.extract (in an
:
Achim Kisseler wrote:
Dear grass users,
has anyone an idea how to handle v.to.rast, if vector-file is empty? -
it produces an error and no map is created:
Loading data...
ERROR: Database connection not defined for layer 2
I want an null-raster-map instead.
r.mapcalc nullmap = null
I'm not sure if I understood you right:
the border of a watershed is unique, if you see the basin starting from
outlet of a river into the ocean.
But there are unlimited subbasins (up to the resolution of the raster)!
So you have to define a characteristic for the watersheds, eg. a basin
size
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
Hi Dylan,
v.to.rast has the option:
--%--
labelcolumn=name
Name of column used as raster category labels
--%--
And the raster values have to be numeric:
--%--
column=name
Name of column for attr parameter (data type must be numeric)
--%--
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
Hi Hamish,
thanks! I did like you said and it works fine.
Because the database (sqlite) has its own ids, it can handle the grass
internal cats even when they are negative.
Multiplying forward and backward let me identify old numbers.
I was a bit afraid, that the numbers result are to hight,
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.
What would be a way? The only I can imagine so far is:
- collecting coordinates
- building lines from them
- update length of lines
-
Hi,
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 problem: v.to.db does not update negative cats! Is there any way to
handle this?
I dont want to loose the identification, nor miss the -v Flag:
It seems, that there are line-segments and a buffer is made for each. Do
you mean that?
Otherwise check:
http://grass.itc.it/grass64/manuals/html64_user/v.buffer.html
eg.
v.buffer input=map output=buffer type=line distance=100
Flags:
-s
Make outside corners straight
-c
Don't make caps
...
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
89 matches
Mail list logo