Thanks Manuel.
I could download the file from address you have provided.
This is the final code that I used to generate an orthorectified SAR image
with the 2157 projection.
*otbcli_OrthoRectification -io.in SARimage.tif -io.out SARoutput.tif
-map.epsg.code 2157 -elev.dem EUDEM_folders/EUDEM_2157 -elev.geoid
Input_DEM_egm96.grd -opt.gridspacing 60.0*
Again it seems to take the default projection '32629' - information as
follows
Driver: GTiff/GeoTIFF
Files: vh_trial1_tc.tif
Size is 28319, 23641
Coordinate System is:
PROJCS["WGS 84 / UTM zone 29N",
GEOGCS["WGS 84",
DATUM["WGS_1984",
SPHEROID["WGS 84",6378137,298.257223563,
AUTHORITY["EPSG","7030"]],
AUTHORITY["EPSG","6326"]],
PRIMEM["Greenwich",0],
UNIT["degree",0.0174532925199433],
AUTHORITY["EPSG","4326"]],
PROJECTION["Transverse_Mercator"],
PARAMETER["latitude_of_origin",0],
PARAMETER["central_meridian",-9],
PARAMETER["scale_factor",0.9996],
PARAMETER["false_easting",500000],
PARAMETER["false_northing",0],
UNIT["metre",1,
AUTHORITY["EPSG","9001"]],
AUTHORITY["EPSG","32629"]]
Origin = (368425.093750000000000,5943843.000000000000000)
Pixel Size = (10.013348579406738,-10.013348579406738)
Metadata:
AREA_OR_POINT=Area
TIFFTAG_DATETIME=2014:11:28 00:38:00
TIFFTAG_IMAGEDESCRIPTION=Sentinel-1A IW GRD HR L1
TIFFTAG_SOFTWARE=Sentinel-1 IPF 002.36
Image Structure Metadata:
INTERLEAVE=BAND
Corner Coordinates:
Upper Left ( 368425.094, 5943843.000) ( 10d59'22.66"W, 53d37'36.57"N)
Lower Left ( 368425.094, 5707117.426) ( 10d53'44.22"W, 51d29'59.50"N)
Upper Right ( 651993.112, 5943843.000) ( 6d42' 6.74"W, 53d37'16.67"N)
Lower Right ( 651993.112, 5707117.426) ( 6d48'37.56"W, 51d29'41.07"N)
Center ( 510209.103, 5825480.213) ( 8d50'57.61"W, 52d34'45.38"N)
Band 1 Block=28319x1 Type=Float32, ColorInterp=Gray
NoData Value=0
I am not sure what is causing the image to take the default values and not
the input projection information.
Also, the DEM cuts the water portion of the image - this is also not
evident in the output image.
Could somebody please help me with this?
I really appreciate any help/suggestions.
Thanks!
On Tuesday, September 12, 2017 at 3:55:21 PM UTC+1, Manuel Grizonnet wrote:
>
> the geoid file can be download here:
>
> https://git.orfeo-toolbox.org/otb-data.git/blob/HEAD:/Input/DEM/egm96.grd
>
> <https://git.orfeo-toolbox.org/otb-data.git/blob/HEAD:/Input/DEM/egm96.grd>
> You need also the .hdr file:
>
>
> https://git.orfeo-toolbox.org/otb-data.git/blob/HEAD:/Input/DEM/egm96.grd.hdr
>
> The URL you've got point to the mercurial repository which is not
> available anymore. Let me know if you find this adresse in the OTB
> documentation or website, in this case it should be updated.
>
> Regards,
>
>
> 2017-09-12 12:56 GMT+02:00 preet.balaji20 <[email protected]
> <javascript:>>:
>
>> Hi Julien,
>>
>> Thanks very much for your suggestions. I am trying to download geoid
>> files, do you know from where I can download the files? I tried this -
>> http://hg.orfeo-toolbox.org/OTB-Data/raw-file/404aa6e4b3e0/Input/DEM/egm96.grd,
>>
>> but the page is not working.
>>
>> Thanks and Regards
>> Preethi
>>
>>
>> On Tuesday, September 12, 2017 at 10:17:33 AM UTC+1, Julien Michel wrote:
>>
>>> Also note that you should change the default -opt.gridspacing value (in
>>> meters) to something like twice the resolution of your DEM or so
>>> (otherwise ortho-recitification will take forever to complete).
>>>
>>> Regards,
>>>
>>> Julien
>>>
>>> Le 12/09/2017 à 11:08, Julien Michel a écrit :
>>> > Hi Preethi,
>>> >
>>> > First of all, you do not need to use the DEMConvert application
>>> > anymore. Simply put DEM.tif in a directory and set it as -elev.dem
>>> > parameter. Second, depending on your DEM convention, you also need to
>>> > set a geoid file to get accurate results.
>>> >
>>> > Last, the output image is not supposed to be in the same SRS as the
>>> > DEM. In fact the output SRS is tuned with the -map parameter. By
>>> > default it will use the closest UTM zone. If you want to use a
>>> > specific SRS, you can set it by its EPSG code : -map epsg -map.epsg
>>> 2157
>>> >
>>> > Hope that helps,
>>> >
>>> > Regards,
>>> >
>>> > Julien
>>> >
>>> >
>>> > Le 12/09/2017 à 11:08, preet.balaji20 a écrit :
>>> >>
>>> >>
>>> >> Dear all,
>>> >>
>>> >> I am trying to orthorectify Sentinel-1 GRD Radiometric calibrated
>>> >> C-band SAR image using EU digital elevation model of 30m resolution
>>> >> using the command in Linux terminal
>>> >>
>>> >> *otbcli_OrthoRectification -io.in SARimage.tif -io.out
>>> >> outputimage.tif -elev.dem /DEMdirectory. *
>>> >>
>>> >> The DEM used was reprojected to Irish Transverse Mercator –
>>> epsg:2157.
>>> >>
>>> >> The DEM.tif file was converted to general raster file compatible with
>>> >> the OTB DEM handling using the DEMconvert option. Three files with
>>> >> the extensions .ras, .geom and .omd were saved in the dem directory -
>>> >> this directory was given as input for orthorectification
>>> >>
>>> >>
>>> >> *DEM file information:*
>>> >>
>>> >> Size is 9663, 7541
>>> >>
>>> >> Coordinate System is:
>>> >>
>>> >> PROJCS["IRENET95 / Irish Transverse Mercator",
>>> >>
>>> >> GEOGCS["IRENET95",
>>> >>
>>> >> DATUM["IRENET95",
>>> >>
>>> >> SPHEROID["GRS 1980",6378137,298.2572221010002,
>>> >>
>>> >> AUTHORITY["EPSG","7019"]],
>>> >>
>>> >> TOWGS84[0,0,0,0,0,0,0],
>>> >>
>>> >> AUTHORITY["EPSG","6173"]],
>>> >>
>>> >> PRIMEM["Greenwich",0],
>>> >>
>>> >> UNIT["degree",0.0174532925199433],
>>> >>
>>> >> AUTHORITY["EPSG","4173"]],
>>> >>
>>> >> PROJECTION["Transverse_Mercator"],
>>> >>
>>> >> PARAMETER["latitude_of_origin",53.5],
>>> >>
>>> >> PARAMETER["central_meridian",-8],
>>> >>
>>> >> PARAMETER["scale_factor",0.99982],
>>> >>
>>> >> PARAMETER["false_easting",600000],
>>> >>
>>> >> PARAMETER["false_northing",750000],
>>> >>
>>> >> UNIT["metre",1,
>>> >>
>>> >> AUTHORITY["EPSG","9001"]],
>>> >>
>>> >> AUTHORITY["EPSG","2157"]]
>>> >>
>>> >> Origin = (419504.356750835664570,709074.848424524534494)
>>> >>
>>> >> Pixel Size = (25.288322889600000,-25.288322889600000)
>>> >>
>>> >> Metadata:
>>> >>
>>> >> AREA_OR_POINT=Area
>>> >>
>>> >> Image Structure Metadata:
>>> >>
>>> >> INTERLEAVE=BAND
>>> >>
>>> >> Corner Coordinates:
>>> >>
>>> >> Upper Left ( 419504.357, 709074.848) ( 10d41'44.00"W, 53d 6' 6.05"N)
>>> >>
>>> >> Lower Left ( 419504.357, 518375.606) ( 10d35'37.23"W, 51d23'21.80"N)
>>> >>
>>> >> Upper Right ( 663865.421, 709074.848) ( 7d 2'44.58"W, 53d 7'42.18"N)
>>> >>
>>> >> Lower Right ( 663865.421, 518375.606) ( 7d 4'54.58"W, 51d24'52.20"N)
>>> >>
>>> >> Center ( 541684.889, 613725.227) ( 8d51'16.06"W, 52d16'19.64"N)
>>> >>
>>> >> Band 1 Block=9663x1 Type=Float32, ColorInterp=Gray
>>> >>
>>> >> NoData Value=0
>>> >>
>>> >>
>>> >> It seems that the DEM is not being considered as an input here. The
>>> >> ouput image has the default projection – UTM zone 29. But it is
>>> >> supposed to consider the projection of the input DEM which is ITM
>>> >> (2157). Please find the attached orthorectified image obtained using
>>> >> the DEM files.
>>> >>
>>> >> The Orthorectification application in OTB version 6.0.0 was also
>>> >> tried which did not give the desired result.
>>> >>
>>> >> Can somebody please help me in understanding what is it that I'm
>>> >> missing/doing wrong?
>>> >>
>>> >> I really appreciate the help.
>>> >>
>>> >> Thanks in advance
>>> >> Preethi
>>> >>
>>> >> --
>>> >> --
>>> >> 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]
>>> >> <mailto:[email protected]>.
>>> >> For more options, visit https://groups.google.com/d/optout.
>>> >
>>> >
>>>
>>>
>>> --
>>> Julien MICHEL
>>> CNES - DSO/SI/2A
>>>
>>> --
>> --
>> 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]
>> <javascript:>
>> To unsubscribe from this group, send email to
>> [email protected] <javascript:>
>> 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] <javascript:>.
>> For more options, visit https://groups.google.com/d/optout.
>>
>
>
>
> --
> Manuel Grizonnet
>
--
--
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.
