<br><br><div class="gmail_quote">On Thu, Apr 18, 2013 at 12:10 PM, Rutger <span dir="ltr"><<a href="mailto:kassies@gmail.com" target="_blank">kassies@gmail.com</a>></span> wrote:<br><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">
Hey,<br>
<br>
<br>
Etienne Tourigny-3 wrote<br>
<div>> I use a similar command, without the -a_ullr argument<br>
<br>
</div>And does it then capture the correct corner coordinates? If i run your<br>
command it does assign the sinusoidal projection, but takes the pixel<br>
coordinates (eg 0-2048) as the corner coordinates instead of the sinusoidal<br>
ones.<br></blockquote><div><br></div>Actually, I used this with the level-3 tiled products such as MCD12Q1 or MCD45A1. Never tried with level-2 or swath products.<br>Yes the corners are fine, I can convert a few of them to WGS84 and then stitch them together with gdalwarp.<br>
<br>Etienne</div><div class="gmail_quote"><br><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">
<br>
@Laura,<br>
Indeed we are talking about different products, sorry for the confusion. You<br>
mentioned in your first email the products MYD14.A2 and MOD14.A2. I thought<br>
you meant the "MYD14 A2" product, but you attached a "MYD14" product, which<br>
is different. Its obvious where the confusion comes from. :)<br>
<br>
These are all the fire products:<br>
MYD14 Aqua Thermal Anomalies & Fire Swath 1000m 5 min<br>
MOD14A2 Terra Thermal Anomalies & Fire Tile 1000m 8 day<br>
MYD14A2 Aqua Thermal Anomalies & Fire Tile 1000m 8 day<br>
MOD14A1 Terra Thermal Anomalies & Fire Tile 1000m Daily<br>
MYD14A1 Aqua Thermal Anomalies & Fire Tile 1000m Daily<br>
MOD14 Terra Thermal Anomalies & Fire Swath 1000m 5 min<br>
<br>
The "MOD14" & "MYD14" products are indeed swath, and should be used together<br>
with the MOD03 with the same timestamp.<br>
<br>
I think the best way to process them is by using "MRT Swath" which you can<br>
download at:<br>
<a href="https://lpdaac.usgs.gov/tools" target="_blank">https://lpdaac.usgs.gov/tools</a><br>
<br>
MRT Swath is way beyond the scope of this mailing list, but feel free to<br>
send me a direct email about it if you have questions.<br>
<br>
I have used Pyresample as an alternative which works very well, but requires<br>
the use of Python:<br>
<a href="http://code.google.com/p/pyresample/" target="_blank">http://code.google.com/p/pyresample/</a><br>
<br>
To end with an actual GDAL related solution, there is an option to warp an<br>
image using gdalwarp and so-called "geolocation arrays". I can get it to<br>
work with the MODIS swaths, but i dont like the results at all. It is<br>
previously discussed at:<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>
The MODIS swath products arent equally gridded and still contain what is<br>
called the 'bow-tie' effect. (overlapping scans at high zenith angles). I<br>
dont know if the GDAL geoloc option is supposed to handle this properly. The<br>
image a posted earlier in the thread mentioned above clearly shows artifacts<br>
which i cannot explain or resolve.<br>
<br>
<br>
Regards,<br>
Rutger<br>
<br>
<br>
<br>
<br>
<br>
<br>
<br>
<br>
--<br>
View this message in context: <a href="http://osgeo-org.1560.x6.nabble.com/MODIS-reprojection-tp5047700p5048019.html" target="_blank">http://osgeo-org.1560.x6.nabble.com/MODIS-reprojection-tp5047700p5048019.html</a><br>
<div><div>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>
</div></div></blockquote></div><br>