[gdal-dev] reproject python numpy binary swath/lat/lon

[Etienne Tourigny]

To use geolocation arrays with gdalwarp you probably need to add the
-geoloc argument

Etienne

On Tue, Jun 12, 2012 at 6:08 AM, Rutger <kassies at gmail.com> wrote:
> Hello David,
>
> A while back i did some tests with GDAL geolocation arrays using MODIS swath
> data. I found that specifying the "s_srs" tag when running gdalwarp gets rid
> of the 'too many points' error, even if you already have a SRS in your VRT.
> If you use gdalwarp, your target grid can be completely described by the
> -t_srs -tr and -t switches. The Python equivalent of these are the
> dataset.SetProjection() and dataset.SetGeoTransform(). Getting the input SRS
> from only the coordinate values seems impossible to me, the SRS should be
> either read from some metadata or specified by the user.
>
> For me the geolocation arrays worked technically well, and very convenient.
> However. i found the resulting image quality very disappointing (for my
> specific MODIS case) compared to MRTswath and Pyresample. This might be due
> to the irregular coordinate grids as a result of the swiping image sensor or
> something like that. Maybe in different cases it would give good quality
> results.
>
> I ended up using the Pyresample module which is part from the Pytroll
> project. See:
> http://pytroll.org/index.html http://pytroll.org/index.html
> This worked great for my case.
>
>
> Regards,
> Rutger
>
>
> --
> View this message in context: http://osgeo-org.1560.n6.nabble.com/gdal-dev-reproject-python-numpy-binary-swath-lat-lon-tp4978609p4980686.html
> Sent from the GDAL - Dev mailing list archive at Nabble.com.
> _______________________________________________
> gdal-dev mailing list
> gdal-dev at lists.osgeo.org
> http://lists.osgeo.org/mailman/listinfo/gdal-dev