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@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/gdal-dev
