[gdal-dev] reprojecting coastcolour (/meris) using python via GCPs

[Ivan Price]

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


-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://lists.osgeo.org/pipermail/gdal-dev/attachments/20140217/34aaacf3/attachment-0001.html>