I have had similar experiences using the gdal_translate utility. However if you reproject the data using gdalwarp, it will recognize the projection info and the -o flag does not need to be used.
Casey
grass-user-request@lists.osgeo.org wrote:
Send grass-user mailing list submissions to
grass-user@lists.osgeo.orgTo 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
grass-user-request@lists.osgeo.orgYou can reach the person managing the list at
grass-user-owner@lists.osgeo.orgWhen 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 (michael.barton@asu.edu)----------------------------------------------------------------------
Message: 1
Date: Wed, 25 Mar 2009 17:12:59 +0100
From: Moritz Lennert <mlennert@club.worldonline.be>
Subject: Re: [GRASS-user] extracting probabilities from i.maxlik
To: Georg Kaspar <georg@geofs.de>
Cc: grass-user@lists.osgeo.org
Message-ID: <49CA580B.7090006@club.worldonline.be>
Content-Type: text/plain; charset=UTF-8; format=flowedOn 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 <mlennert@club.worldonline.be>
Subject: Re: [GRASS-user] Vector map attributes
To: Craig Leat <craig.leat@gmail.com>
Cc: GRASS users list <grass-user@lists.osgeo.org>
Message-ID: <49CA5A5D.4010801@club.worldonline.be>
Content-Type: text/plain; charset=ISO-8859-1; format=flowedOn 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 <maris.gis@gmail.com>
Subject: Re: [GRASS-user] EPSG 3035
To: Laura Poggio <laura.poggio@gmail.com>
Cc: grass-user@lists.osgeo.org
Message-ID:
<2b3d8c1c0903250923s1141327p1b7ef2656c021420@mail.gmail.com>
Content-Type: text/plain; charset=UTF-8My 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 <mlennert@club.worldonline.be>
Subject: Re: [GRASS-user] EPSG 3035
To: Laura Poggio <laura.poggio@gmail.com>
Cc: grass-user@lists.osgeo.org
Message-ID: <49CA5AC2.50501@club.worldonline.be>
Content-Type: text/plain; charset=ISO-8859-1; format=flowedOn 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 <laura.poggio@gmail.com>
Subject: Re: [GRASS-user] EPSG 3035
To: Moritz Lennert <mlennert@club.worldonline.be>
Cc: grass-user@lists.osgeo.org
Message-ID:
<e027b8750903251005wa86d524p88484ea817c3d5d0@mail.gmail.com>
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=athematicThe parameters seems to be correct. I will then use the -o option
Thank you very much
Laura
2009/3/25 Moritz Lennert <mlennert@club.worldonline.be>
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: michael.barton@asu.edu
Subject: [GRASS-user] Mail System Error - Returned Mail
To: grassuser@grass.itc.it
Message-ID: <200903251828.n2PIS84v021273@grass.itc.it>
Content-Type: text/plain; charset="us-ascii"³dbS¡ö
ÓhCFw½Ôßà
5û·zb}iú¬éCÎö2çh·OÝsN¨_ýÌÀéöèln}qTÈ1zèxPz«"¼sÐ*©R»ÙË´`<ÔfÄGyMåÄumæfÃ)ú¬Âb2G?úè5Ϥhml¿È÷ªÙmNsGjG¾m¨;Ç} ÉG.§yÈéÄø_RÖÞB"¾cæèÂüMZ1´Yå3hÑÙ¨Xy'&f
¼5ßû¾
*M])4³ªÄòIL1zXß¾,̧íìU>&Dd½a1ÕG°ü÷Ïöú1mr¿F<2,JÑ{/÷%ÚSgºöæ
1h9ÌÚ¥òܧÍÑÝ
¨¹qãÊè(É}WH[æ#ñxcMG"NuÌÅyñ¦:±ð`aOÉE~
«OÑ°× ËUE,!òwg¬ÑÄÐ6þVçÅä;ûÄ~Â"þCþª2ÞQ:BbB¶ ü#Ñ7ÒéñlU»q®ç6hBê¾ùØÃ4KÌ-
`]¶ûÈÊåü¹~ØëÛÏæj:$çõ¤ÚNDWqça!fú3ë6è(â¹rënD
7&NÇH>¥²§Jk,þ¾½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
grass-user@lists.osgeo.org
http://lists.osgeo.org/mailman/listinfo/grass-userEnd of grass-user Digest, Vol 35, Issue 55
******************************************