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>
<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
<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
<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 <http://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
<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
<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, visithttps://groups.google.com/d/optout
<https://groups.google.com/d/optout>.
>
>
--
Julien MICHEL
CNES - DSO/SI/2A
--
--
Check the OTB FAQ at
http://www.orfeo-toolbox.org/FAQ.html
<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
<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
<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]
<mailto:[email protected]>.
For more options, visit https://groups.google.com/d/optout.
