Thank both of you Dominik and Loïc

alignExtent() did the job after regressing the wright points vs bad points and predicting a new correct extent. After setting the new extent, the points match correctly known features of the image.

I would give a try to gdalwarp to add GCPs to the dataset.


On 17/08/15 14:51, Dominik Schneider wrote:
Marcelo,
did you try alignExtent() in the raster package?


Dominik Schneider
o 303.735.6296 | c 518.956.3978


On Mon, Aug 17, 2015 at 8:34 AM, Loïc Dutrieux <[email protected]> wrote:

Hi Marcelo,

I don't see anything in raster or rgdal that would do that directly.
However, this can be done I believe using gdal_translate / gdal_edit (for
adding the GCPs to the dataset) + gdalwarp (for warping). Both
gdal_translate and gdalwarp utilities have wrappers in the gdalUtils
package.

Cheers,
Loïc Dutrieux


On 08/17/2015 11:43 AM, Marcelo Kittlein wrote:

Hi all

I want to change a MSS landsat image which has bad georeferencing.

I have a grid of points with correct georeferences and want to use them
to reproject the image to a correct reference system.

Is there some way to this in R using "raster" or other package?

in idrisi I used to resample the image using a correspondence file
(.cor) with xold yold xnew ynew values.

I wonder if there is something similar I can do in R? Thanks in advance

Marcelo Kittlein

_______________________________________________
R-sig-Geo mailing list
[email protected]
https://stat.ethz.ch/mailman/listinfo/r-sig-geo

_______________________________________________
R-sig-Geo mailing list
[email protected]
https://stat.ethz.ch/mailman/listinfo/r-sig-geo

        [[alternative HTML version deleted]]

_______________________________________________
R-sig-Geo mailing list
[email protected]
https://stat.ethz.ch/mailman/listinfo/r-sig-geo

_______________________________________________
R-sig-Geo mailing list
[email protected]
https://stat.ethz.ch/mailman/listinfo/r-sig-geo

Reply via email to