<div dir="ltr">Thank you Rutger, it was me who was missing things - I've been confused about to what level the GDAL warping is appropriate, and every part of this discussion has been very helpful. It seems very much like the orthorectification process where you need a DEM as well as details on camera geometry, but my understanding on the relationship here is very dim. <div><br><div>I'd love for lower level techniques like the MODIS Swath to grid and Seadas L1A to L2 chain to be available via GDAL, but I have a lot to learn about all the steps. I have no idea how appropriate it is for GDAL to aspire to these things, but I'm very keen to learn.  </div><div><br></div><div>Cheers, Mike<br><br><div class="gmail_quote">On Thu, 21 May 2015 at 17:48 Rutger <<a href="mailto:kassies@gmail.com">kassies@gmail.com</a>> wrote:<br><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">Hey,<br>
<br>
Maybe i'm missing something, but for as far as i know, the MODIS MOD021KM<br>
product is still a 'swath', and therefore still containing the bow-tie<br>
effect as the result of the scanning sensor. Wouldn't this mean that<br>
gdalwarp will never work? The lat/lon arrays aren't continuous, so any<br>
fitting (GCP/TPS) wont be accurate. And upsampling the arrays with cubic,<br>
like Even did, will result in interpolating pixels which are next each other<br>
in the array/file, but aren't next to each other 'on the ground', wont work<br>
well.<br>
<br>
I think you need to either use a dedicated algorithm (like gradient search),<br>
or something very generic (like the KDTree of Pyresample) to overcome the<br>
bow-tie effect.<br>
<br>
I have spend alot of time experimenting with this in GDAL, since i try to do<br>
as much with GDAL as possible. But projecting a swath file from a scanning<br>
sensor (like MODIS or VIIRS) is something i have never gotten to work. For<br>
some historic discussion about this see for example:<br>
<a href="http://osgeo-org.1560.x6.nabble.com/gdal-dev-reproject-python-numpy-binary-swath-lat-lon-td4978609.html" target="_blank">http://osgeo-org.1560.x6.nabble.com/gdal-dev-reproject-python-numpy-binary-swath-lat-lon-td4978609.html</a><br>
<br>
<br>
Regards,<br>
Rutger<br>
<br>
<br>
<br>
--<br>
View this message in context: <a href="http://osgeo-org.1560.x6.nabble.com/gdal-dev-gdalwarp-problems-with-MODIS-L1B-tp5206583p5206825.html" target="_blank">http://osgeo-org.1560.x6.nabble.com/gdal-dev-gdalwarp-problems-with-MODIS-L1B-tp5206583p5206825.html</a><br>
Sent from the GDAL - Dev mailing list archive at Nabble.com.<br>
_______________________________________________<br>
gdal-dev mailing list<br>
<a href="mailto:gdal-dev@lists.osgeo.org" target="_blank">gdal-dev@lists.osgeo.org</a><br>
<a href="http://lists.osgeo.org/mailman/listinfo/gdal-dev" target="_blank">http://lists.osgeo.org/mailman/listinfo/gdal-dev</a><br>
</blockquote></div></div></div></div>