[GRASS-user] i.rectify: how to make it works?!?!?

[Enrico g]

 Thank you John for the answer.
 That's a choiche (I will try it after lunch time) but I'm studying a
more user friendly approach as I'm looking to find a collection of
tools to use in our monkey... sorry, office environment....
 (So basically now I'm focused on GRASS while later maybe I'll try
QGIS et  similia to find out something more "usable")

> --------------------------------------------------
> From: "John Stevenson" <john.stevenson at manchester.ac.uk>
> Sent: Tuesday, April 07, 2009 1:26 PM
> To: "Enrico Graziani" <mrenrich84 at gmail.com>
> Cc: <grass-user at lists.osgeo.org>
> Subject: Re: [GRASS-user] i.rectify: how to make it works?!?!?
>
>> Enrico Graziani wrote:
>>>
>>> Hello!
>>> I'm trying to use imagery tools to georeference/rectify an image taken
>>> from google earth...
>>> Here's the procedure I've used without success:
>>>  - From GEarth saved image to a local file (jpeg) ( I've tried to convert
>>> this jpeg to png but this didn't helped me. In both case I've keeped the
>>> image as 24 bit...)
>>> - Created a location wgs84 utm29 (North!) with region coordinate taken
>>> from GEarth (in local format)
>>> - Created a location xy in grass (6.2.3) with region as big as the
>>> picture
>>> - Imported the image with g.in.gdal: this creates 3 raster maps, blue
>>> green and red
>>> - r.composite the 3 rasters to one "gearth" raster
>>> - i.group to create group "group" containing the gearth raster
>>> - i.target giving target the wgs84 utm29 location (mapset user)
>>> - i.points to give 3 GCPs
>>> (- tried also to run i.group to remove the gearth raster from r.composite
>>> and add the 3 rgb rasters map)
>>> - i.rectify -a group=group input=group extension=_rect order=1
>>>  Everything goes fine till here, except for the fact that if I add the -c
>>> option to i.rectify it will always fail with "error while writing to temp
>>> file".
>>> However, when I go to other location the map doesn't display.
>>> If I've rectified the 3 rgb rasters I can try to display the map with an
>>> RGB layer or creating a single raster map with r.composite but, first of
>>> all, it takes forever (running "top" while doing this displayed a 99% of CPU
>>> use by d.rgb or r.composite) and, last but not last, the map doesn't
>>> display. If I select on the map display "zoom to selected map" it will zoom
>>> to default region (the selected region in this location) even if I used
>>> i.rectify without the -c flag!
>>> Worster if I try to rectify the "r.composite"d raster (the one containing
>>> the 3 imported channel)  first and then rectify it. Everything goes ok but
>>> if I try to display it grass just stuck and I need to restart the x server
>>> to quit the grass windows...
>>>  I had a look at the data in the wgs84_utm29 location and noticed that
>>> this last "r.composite"d raster hasn't been "converted". Data in the cell
>>> directory is around 24 kb and there is a 148 kb file in the
>>> cell_misc/$RASTER_NAME directory ..   starting jpeg was a couple of MBs..
>>> Stranger in the other case (rectifying the 3 channels and then displaying
>>> rgb layer or r.composite them): the files in the cell directory are around
>>> 15 kb and 148 kb in the cell_misc folders, but after r.composite, the raster
>>> coming from this has 6.9 MB in the cell folder and 107.3 MB in the cell_misc
>>> folder of data.
>>> All the files I'm talking about in the cell misc folder are called "null"
>>> and dolphin recognized them by the MIME type as jpeg2000, but gwenview isn't
>>> able to display them...
>>>  My configuration is the installation of the GRASS Live DVD coming from
>>> the Trento University in a VMWare machine running with VMPlayer 2.0.5 on
>>> Debian Etch.
>>> This live should be Kubuntu Hardy ( ?) based.
>>> ------------------------------------------------------------------------
>>>
>>> _______________________________________________
>>> grass-user mailing list
>>> grass-user at lists.osgeo.org
>>> http://lists.osgeo.org/mailman/listinfo/grass-user
>>>
>> For a quick georeferencing, with just 3 gcps, you can just use gdal.
>>
>> gdal_translate -a_srs "EPSG:4326" -of GTiff -gcp x1 y1 X1 Y1 -gcp x2 y2 X2
>> Y2 -gcp x3 y3 X3 Y3 gearth.jpg gearth.tif
>>
>> gdalwarp -t_srs "EPSG:4326" -rc  gearth.tif gearth_rectified.tif
>>
>> Then load the image with r.in.gdal.
>> You can get the gcps by hovering the mouse over landmarks in The Gimp
>> (x,y) and Google Earth (X,Y).
>> You can find the EPSG code for your utm region here:
>> http://www.epsg-registry.org/
>>
>> Cheers
>>
>> John
>>
>> --
>>
>>
>> Dr John Stevenson
>> Postdoctoral Research Associate
>> School of Earth, Atmospheric and Environmental Sciences
>> Williamson Building (Room 2.42)
>> University of Manchester
>> Manchester M13 9PL, UK
>> tel. +44(0)161 306 6585; fax. +44(0)161 306 9361;
>> john.stevenson at manchester.ac.uk
>