[PATCH] NOAA level 1b geolocation interpolation
|Reported by:||arb||Owned by:||dron|
|Cc:||neteler, jurien, spareeth, mmetz|
The current noaa level 1b geolocation relies on sparse GCPs. This might be fine when reprojecting a small area but when reprojecting a whole swath width it produces significant errors.
It might be better if the GCPs were interpolated inside the driver and returned as a geolocation array, or if a new Lagrangian interpolation method was added to GDAL.
See ticket 2388: http://trac.osgeo.org/gdal/ticket/2388 although the sample code there doesn't handle the case where longitude wraps around (360 back to 0), or switches sign (180 to -180). When copying the 51 longitudes into a 2048-element array (prior to interpolation) it would be advisable to check for a difference > 180 compared to the previous longitude and just continue with the same magnitude/sign, then perform some modulo function after interpolation. Hope that helps!
Change History (15)
comment:3 Changed 8 years ago by
|Component:||default → GDAL_Raster|
|Version:||unspecified → svn-trunk|
comment:6 follow-up: 7 Changed 6 years ago by
|Summary:||NOAA level 1b geolocation interpolation → [PATCH] NOAA level 1b geolocation interpolation|