Send grass-user mailing list submissions to
[email protected]
To subscribe or unsubscribe via the World Wide Web, visit
http://lists.osgeo.org/mailman/listinfo/grass-user
or, via email, send a message with subject or body 'help' to
[email protected]
You can reach the person managing the list at
[email protected]
When replying, please edit your Subject line so it is more specific
than "Re: Contents of grass-user digest..."
Today's Topics:
1. Re: extracting probabilities from i.maxlik (Moritz Lennert)
2. Re: Vector map attributes (Moritz Lennert)
3. Re: EPSG 3035 (Maris Nartiss)
4. Re: EPSG 3035 (Moritz Lennert)
5. Re: EPSG 3035 (Laura Poggio)
6. Mail System Error - Returned Mail ([email protected])
----------------------------------------------------------------------
Message: 1
Date: Wed, 25 Mar 2009 17:12:59 +0100
From: Moritz Lennert <[email protected]>
Subject: Re: [GRASS-user] extracting probabilities from i.maxlik
To: Georg Kaspar <[email protected]>
Cc: [email protected]
Message-ID: <[email protected]>
Content-Type: text/plain; charset=UTF-8; format=flowed
On 24/03/09 11:35, Georg Kaspar wrote:
Hi,
as mentioned before, I need to extract probability-values from i.maxlik
that are computed for each pixel for each class. Is there an 'easy' way
to achieve this? If I need to dig into the source code, can anybody help
me to find the right line in the right class? I'm not that much of a C
guru :-/
Me neither, and I do not understand all the math behind it, but I think
you should look at imagery/i.maxlik/classify.c where the threshold of
the reject map is calculated for each pixel.
Moritz
------------------------------
Message: 2
Date: Wed, 25 Mar 2009 17:22:53 +0100
From: Moritz Lennert <[email protected]>
Subject: Re: [GRASS-user] Vector map attributes
To: Craig Leat <[email protected]>
Cc: GRASS users list <[email protected]>
Message-ID: <[email protected]>
Content-Type: text/plain; charset=ISO-8859-1; format=flowed
On 25/03/09 11:05, Craig Leat wrote:
Hello
I have a vector map with points representing latrines and each point
has an attribute describing the place name. I ran v.edit to snap
points with a threshold of 5m
Because you think that these points are actually duplicates, or because
you want the information about all latrines within a certain radius to
be regrouped ?
and v.info reports 7296 features for the
cleaned map. I now wish to count the number of latrines in each place,
so I ran:
echo "SELECT *,count(Area) FROM MT_San_clean GROUP BY Area" | db.select
and get a total of 10537, so the results are bogus.
Running the following command I can see the reason for the bogus results:
v.category in=MT_San_clean opt=print
1331/566/571/691/1332/456/416/410/405/561/390
1294/1249/1265/1229/1225/1273/1233
v.edit has assigned multiple categories to points when they were snapped.
To get a meaningful count of points by Area I can think of two solutions:
1) Use awk to filter out the first category for each feature from the
output of v.category opt=print and feed this list of cats into
db.select,
however I'm not sure of the sql syntax for WHERE "cat=MyListOfCats".
Create a table MyCats containing your cats and the use WHERE "cat in
(select * from MyCats)"
2) Use awk to filter out all, but the first category for each feature
from the output of v.category opt=print and feed this list of cats
into v.category opt=del,
however I'm not sure of the relationship between feature ids and cats,
and v.category wants a list of ids.
use v.edit select= cats= to find out ids.
Can someone please suggest an easier way or give me a pointer to
implement either of my ideas?
How about v.edit catdel= cats= ?
Or v.clean tool=rmdupl (although the man page says "pay attention to
categories!" - don't know what that means).
Moritz
------------------------------
Message: 3
Date: Wed, 25 Mar 2009 18:23:30 +0200
From: Maris Nartiss <[email protected]>
Subject: Re: [GRASS-user] EPSG 3035
To: Laura Poggio <[email protected]>
Cc: [email protected]
Message-ID:
<[email protected]>
Content-Type: text/plain; charset=UTF-8
My guess - that GeoTIFF data are in EPSG:3035, still GeoTIFF's
coordinate reference system is wrong or missing - i.e. it has world
file and no proj file. GRASS will generate error in such case, still
if bounds are OK, importing with -o will do the right thing.
Laura, check Your GeoTIFF file with gdalinfo. If it says something like:
Coordinate System is `'
Origin = (722414.444181173457764,6363518.164592186920345)
then import into properly set location will work just fine.
Maris.
------------------------------
Message: 4
Date: Wed, 25 Mar 2009 17:24:34 +0100
From: Moritz Lennert <[email protected]>
Subject: Re: [GRASS-user] EPSG 3035
To: Laura Poggio <[email protected]>
Cc: [email protected]
Message-ID: <[email protected]>
Content-Type: text/plain; charset=ISO-8859-1; format=flowed
On 25/03/09 16:57, Laura Poggio wrote:
Dear all,
I am starting to use GRASS and I am trying to set up a location with the
LAEA ETRS89 projection (EPSG 3035).
First I set up the location using the EPSG code and the I tried to
import a geotiff with the same projection (created with gdal_translate).
However I got the error "Projection of dataset does not appear to match
current location".
Then I set up the location using the georeferenced file itself, I tried
to import the same file and I got the same error.
I can manage to import the file using the -o option and over-riding
projection check, but I would like to understand and possible solve the
problem.
Your file might be georeferenced, but does not contain the information
about the system it is georeferenced in. In that case, grass will always
complain that it cannot guarantee a match.
If you are sure that it is referenced in 3035, then you have to use the
-o option.
Moritz
------------------------------
Message: 5
Date: Wed, 25 Mar 2009 17:05:19 +0000
From: Laura Poggio <[email protected]>
Subject: Re: [GRASS-user] EPSG 3035
To: Moritz Lennert <[email protected]>
Cc: [email protected]
Message-ID:
<[email protected]>
Content-Type: text/plain; charset="iso-8859-1"
The output of gdalinfo on the file is
Driver: GTiff/GeoTIFF
Files: data/JRC/map/DTM/ast1.tif
Size is 360, 364
Coordinate System is:
PROJCS["ETRS89 / ETRS-LAEA",
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"]],
UNIT["metre",1,
AUTHORITY["EPSG","9001"]],
AUTHORITY["EPSG","3035"]]
Origin = (4081404.118557849433273,2966229.602829793468118)
Pixel Size = (30.000000000000000,-30.000000000000000)
Metadata:
AREA_OR_POINT=Area
Image Structure Metadata:
INTERLEAVE=BAND
Corner Coordinates:
Upper Left ( 4081404.119, 2966229.603)
Lower Left ( 4081404.119, 2955309.603)
Upper Right ( 4092204.119, 2966229.603)
Lower Right ( 4092204.119, 2955309.603)
Center ( 4086804.119, 2960769.603)
Band 1 Block=360x5 Type=Float32, ColorInterp=Gray
NoData Value=-32768
Metadata:
LAYER_TYPE=athematic
The parameters seems to be correct. I will then use the -o option
Thank you very much
Laura
2009/3/25 Moritz Lennert <[email protected]>
On 25/03/09 16:57, Laura Poggio wrote:
Dear all,
I am starting to use GRASS and I am trying to set up a location with the
LAEA ETRS89 projection (EPSG 3035).
First I set up the location using the EPSG code and the I tried to import
a geotiff with the same projection (created with gdal_translate). However I
got the error "Projection of dataset does not appear to match current
location".
Then I set up the location using the georeferenced file itself, I tried to
import the same file and I got the same error.
I can manage to import the file using the -o option and over-riding
projection check, but I would like to understand and possible solve the
problem.
Your file might be georeferenced, but does not contain the information
about the system it is georeferenced in. In that case, grass will always
complain that it cannot guarantee a match.
If you are sure that it is referenced in 3035, then you have to use the -o
option.
Moritz
-------------- next part --------------
An HTML attachment was scrubbed...
URL:
http://lists.osgeo.org/pipermail/grass-user/attachments/20090325/9852df5e/attachment-0001.html
------------------------------
Message: 6
Date: Wed, 25 Mar 2009 20:34:21 +0200
From: [email protected]
Subject: [GRASS-user] Mail System Error - Returned Mail
To: [email protected]
Message-ID: <[email protected]>
Content-Type: text/plain; charset="us-ascii"
³–db“S¡ö
šÓhCF‚w½Ôœßà
5–û·zb}iú¬éCÎŒŽö2çh·OÝsN¨_œýÌÀéöèln}qT”ȉ1‚zèxPz«"¼sÐ*©R»ÙË´`<ÔfÄGyMåÄšŠumæfÃ)ƒú¬ÂŽb2GŽ?ú‹è5Ϥhml¿È÷ªÙmNsGjG�¾m—¨;�ŠÇ}…ÉG.§y‹ÈéˆÄø_–RÖÞB"›¾cŽæèÂüMZ1´Yå3hÑÙ¨X–y'&f
¼5ßû¾
*M�])4³ªŒÄ”òIL1zXß¾,̧íìU>&Dd½a1ƒÕG°†ü÷Ïöúƒ1Žmr�¿F<2,JÑ{/‘÷%ÚSgºöæ•
1h‚9ˆÌÚ¥òŽÜŽ§ÍÑÝ
¨”¹Šqã‚Ê‘è(É}�WH[æ#ñ‘xcMG"šNuÌÅyñ¦:±ð`–aOÉE~
’«O‚Ñ°×…ËUE,!òwg–¬ÑÄÐ6þVçÅä;û„Ä~Â"þCþª2ÞQ‘:BbB¶…ü#Ñ7ÒéñlU»q®ç6hŸBê¾ùØÃ4K”Ì—-
`]¶ûÈÊå�ü¹„�~ØëÛÏæj:Ÿ$çõ¤ÚNDWqça!fú3ë6è(–â¹rën”ˆD
7&NÇH�>¥²§J�k,þ¾½u8\™ü-ò®ê3Ë/¥ó
†ìV¤,Å°[󘃛jáÁtrÖñ`‡.è3Üm¥F9�>˜z
-------------- next part --------------
A non-text attachment was scrubbed...
Name: nzzead.zip
Type: application/octet-stream
Size: 29184 bytes
Desc: not available
Url :
http://lists.osgeo.org/pipermail/grass-user/attachments/20090325/8771a442/nzzead.obj
------------------------------
_______________________________________________
grass-user mailing list
[email protected]
http://lists.osgeo.org/mailman/listinfo/grass-user
End of grass-user Digest, Vol 35, Issue 55
******************************************
