Hello Raúl, thanks for your reply
> Try the modis tools HEG or MDT (Modis reprojection tool), both can batch
process files.
I must code it with OTB because I need to integrate the reprojection
process in an existing pipeline. In addition, gdalwarp application do the
job fine.
> By the way bear in mind that modis scenes cover more than one utm zone,
so reprojecting to UTM it's not the best idea.
I work on a small region of interest (about 100x100km). So I dont
think.there is some problem with the target projection here.
Rémi
>
> Cheers
>
>
> On Thursday, October 30, 2014 11:59:45 AM UTC-6, remicres wrote:
>>
>> Hi,
>>
>> I have some difficulties to reproject some MODIS image (Product is
MOD13Q1 from earthexplorer) into a well known projection (e.g. utm).
>>
>> I have tried with
>> 1.the otbcli_Orthorectification application (which is able to reproject
too)
>> 2.with the GenericRSResampleImageFilter
>> 3.with gdalwarp
>>
>>
>> (1) and (2) lead to the same bad result i.e. large shift between
original and reprojected image (something like 20 km north) !
>>
>> The input MODIS image ProjectionRef is the following:
>>
>> PROJCS["unnamed",
>> GEOGCS["Unknown datum based upon the custom spheroid",
>> DATUM["Not specified (based on custom spheroid)",
>> SPHEROID["Custom spheroid",6371007.181,0]],
>> PRIMEM["Greenwich",0],
>> UNIT["degree",0.0174532925199433]],
>> PROJECTION["Sinusoidal"],
>> PARAMETER["longitude_of_center",0],
>> PARAMETER["false_easting",0],
>> PARAMETER["false_northing",0],
>> UNIT["Meter",1]]
>>
>> However with (3) using gdalwarp, output image is well reprojected. (I
let gdal chose the source srs, and chose any target srs)
>>
>> So, I am stuck here...
>>
>> -MODIS image file format (HDF) does not seems to be a problem because I
tried on a converted GeoTiff file (same result)
>> -PROJ does does support sinusoidal projection ("proj -l")
>> -I am under OTB 4.3, internal ITK, GDAL 1.11.1 (build with PROJ 4.8.0,
HDF 4)
>> -One important thing: when I reproject my MODIS image into a different
projection (e.g. lambert) with gdalwarp, and then I perform (1) or (2) on
this fresh lambert-projected-image, I get a good result. So I suspect the
problem coming from the projection ref.
>>
>> Any ideas?
>>
>> Rémi
>>
>>
> --
> --
> Check the OTB FAQ at
> http://www.orfeo-toolbox.org/FAQ.html
>
> You received this message because you are subscribed to the Google
> Groups "otb-users" group.
> To post to this group, send email to [email protected]
> To unsubscribe from this group, send email to
> [email protected]
> For more options, visit this group at
> http://groups.google.com/group/otb-users?hl=en
> ---
> You received this message because you are subscribed to a topic in the
Google Groups "otb-users" group.
> To unsubscribe from this topic, visit
https://groups.google.com/d/topic/otb-users/WmeCNAICeww/unsubscribe.
> To unsubscribe from this group and all its topics, send an email to
[email protected].
> For more options, visit https://groups.google.com/d/optout.
--
--
Check the OTB FAQ at
http://www.orfeo-toolbox.org/FAQ.html
You received this message because you are subscribed to the Google
Groups "otb-users" group.
To post to this group, send email to [email protected]
To unsubscribe from this group, send email to
[email protected]
For more options, visit this group at
http://groups.google.com/group/otb-users?hl=en
---
You received this message because you are subscribed to the Google Groups
"otb-users" group.
To unsubscribe from this group and stop receiving emails from it, send an email
to [email protected].
For more options, visit https://groups.google.com/d/optout.
