As far as I know, the gdal warp api is not exposed to python. I have no idea on reading the MERIS data with the gdal netcdf driver - what I the problem?
But you might be able to use this workaround: instead of creating a source dataset using the memory driver, create a file on disk (in gtiff format) with GCPS and write data to disk, and then use gdalwarp on that. Etienne On Mon, Feb 17, 2014 at 6:46 AM, Ivan Price <[email protected]> wrote: > > > Hello, > > > > I am trying to reproject a window inside a coastcolour (=MERIS) image. As > far as I can see GDAL cannot read the coastcolour data directly, so i am > reading the coastcolour netcdf in python, building a source dataset using > the memory driver, adding GCPS (1 for every 10th pixel) and writing the > data to it, then reprojecting the source dataset to a destination dataset > which is a spatial subset of the original in wgs84 lat/long. > > > > This works fine and is relatively fast, but the reprojection is not > accurate, the results are out by about 6-10 pixels (in various directions). > On reading the forums it seems if i was using gdalwarp i would be using > -tps, however the ReprojectImage() function does not seem to offer this > parameter ? And i don't have the option of using the commandline tool as > even gdal 1.10 cannot recognise the coastcolour data. > > > > So i guess i have 2 questions.. has anyone had any success reading > coastcolour data with the gdal command line tools, and secondly: > > > > how can i get ReprojectImage() to be more accurate, given i have a GCP for > every pixel ? > > > > thanks and regards, > > > > -ivan > > > > > > _______________________________________________ > gdal-dev mailing list > [email protected] > http://lists.osgeo.org/mailman/listinfo/gdal-dev >
_______________________________________________ gdal-dev mailing list [email protected] http://lists.osgeo.org/mailman/listinfo/gdal-dev
