[GRASS-user] Re: EPSG 3035

Casey Vandenberg casey.vandenberg at sjgeophysics.com
Wed Mar 25 15:51:00 EDT 2009 

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 at lists.osgeo.org wrote:
> Send grass-user mailing list submissions to
> 	grass-user at lists.osgeo.org
>
> 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
> 	grass-user-request at lists.osgeo.org
>
> You can reach the person managing the list at
> 	grass-user-owner at lists.osgeo.org
>
> 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 (michael.barton at asu.edu)
>
>
> ----------------------------------------------------------------------
>
> Message: 1
> Date: Wed, 25 Mar 2009 17:12:59 +0100
> From: Moritz Lennert <mlennert at club.worldonline.be>
> Subject: Re: [GRASS-user] extracting probabilities from i.maxlik
> To: Georg Kaspar <georg at geofs.de>
> Cc: grass-user at lists.osgeo.org
> Message-ID: <49CA580B.7090006 at club.worldonline.be>
> 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 <mlennert at club.worldonline.be>
> Subject: Re: [GRASS-user] Vector map attributes
> To: Craig Leat <craig.leat at gmail.com>
> Cc: GRASS users list <grass-user at lists.osgeo.org>
> Message-ID: <49CA5A5D.4010801 at club.worldonline.be>
> 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 <maris.gis at gmail.com>
> Subject: Re: [GRASS-user] EPSG 3035
> To: Laura Poggio <laura.poggio at gmail.com>
> Cc: grass-user at lists.osgeo.org
> Message-ID:
> 	<2b3d8c1c0903250923s1141327p1b7ef2656c021420 at mail.gmail.com>
> 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 <mlennert at club.worldonline.be>
> Subject: Re: [GRASS-user] EPSG 3035
> To: Laura Poggio <laura.poggio at gmail.com>
> Cc: grass-user at lists.osgeo.org
> Message-ID: <49CA5AC2.50501 at club.worldonline.be>
> 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 <laura.poggio at gmail.com>
> Subject: Re: [GRASS-user] EPSG 3035
> To: Moritz Lennert <mlennert at club.worldonline.be>
> Cc: grass-user at lists.osgeo.org
> Message-ID:
> 	<e027b8750903251005wa86d524p88484ea817c3d5d0 at 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=athematic
>
> The parameters seems to be correct. I will then use the -o option
>
> Thank you very much
>
> Laura
>
> 2009/3/25 Moritz Lennert <mlennert at 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 at asu.edu
> Subject: [GRASS-user] Mail System Error - Returned Mail
> To: grassuser at grass.itc.it
> Message-ID: <200903251828.n2PIS84v021273 at grass.itc.it>
> 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
> grass-user at lists.osgeo.org
> http://lists.osgeo.org/mailman/listinfo/grass-user
>
>
> End of grass-user Digest, Vol 35, Issue 55
> ******************************************
>

More information about the grass-user mailing list