Regards,
Guillaume
On 03/16/2015 04:23 PM, Jordi Inglada wrote:
Guillaume Pasero <[email protected]> wrote:
Hi Jordi,
gdal is saying that :
- your image has no coordinate system, which is okay for an image in
sensor geomtry
I understand.
- the GCPs in your image are defined in the EPSG:4326
coordinate
system, which means that the ground coordinated of your GCPs are
expressed in plain lon/lat .
OK.
OTB needs a sensor model to reproject your image. Check with
"otbcli_ReadImageInfo -keywordlist 1" if there is a RPC model
detected.
I dont't think there is an RPC model:
-----------
Image OSSIM keywordlist (optional):
Ossim Keyword list:
Output parameters value:
indexx: 0
indexy: 0
sizex: 25397
sizey: 20713
spacingx: 1
spacingy: 1
originx: 0.5
originy: 0.5
estimatedgroundspacingx: 900.0325928
estimatedgroundspacingy: 328.9576721
numberbands: 1
sensor:
id:
time:
ullat: 0
ullon: 0
urlat: 0
urlon: 0
lrlat: 0
lrlon: 0
lllat: 0
lllon: 0
town:
country:
rgb.r: 0
rgb.g: 1
rgb.b: 2
projectionref:
keyword: Ossim Keyword list:
--------------
Otherwise, you could try to estimate an RPC model from the
GCPs
(application GenerateRPCSensorModel).
I have >250 GCPs in the image, so that should do it.
What is the format of the tie points file for this application?
Thank you.
Jordi
Regards,
Guillaume
On 03/13/2015 03:06 PM, Jordi Inglada wrote:
Hi all,
I have been trying to reproject Sentinel-1 imagery with OTB and I am having
difficulties.
I am using IW-GRD (non-complex) images which are in geotiff format and contain
GCPs. It's been a long time since I worked with image projections, so I may be
misunderstanding things.
First of all, OTB tells me that:
Estimated ground spacing (in meters): [900.033,328.958]
which is definitely wrong. Having a look at the meta-data (with both
ReadImageInfo and gdalinfo) I get coherent results in terms of GCPs,
projections and sizes, but there are things which seem suspicious like the fact
that there is no coordinate system defined and that there are 2 EPSG codes
(gdalinfo):
----------------
Size is 25397, 20713
Coordinate System is `'
GCP Projection =
GEOGCS["WGS 84",
DATUM["WGS_1984",
SPHEROID["unnamed",6378137,298.257223560493,
AUTHORITY["EPSG","4326"]],
AUTHORITY["EPSG","6326"]],
PRIMEM["Greenwich",0],
UNIT["degree",0.0174532925199433],
AUTHORITY["EPSG","4326"]]
GCP[ 0]: Id=1, Info=
(0,0) -> (3.35838036356165,44.0387571612464,573.503620014526)
etc.
-----------------
and ReadImageInformation says the same but gives me funny origin and spacing:
-----------------
Origin : [0.5,0.5]
Spacing : [1,1]
Estimated ground spacing (in meters): [900.033,328.958]
-----------------
For geometry and EPSG-codes experts, do you think that there are problems with
the meta-data?
I could maybe fix the problem with gdal_edit?
Has anyone already worked with Sentinel-1 geometry corrections using
gdal/ossim/OTB?
Thank you.
Jordi
--
* Guillaume PASERO
Ingénieur d'études et développement
Business Unit E-SPACE & Geo Information - Département
APPLICATIONS
CS Systèmes d'Information
Parc de la Grande Plaine - 5, Rue Brindejonc des Moulinais -
BP 15872
31506 Toulouse Cedex 05 - FRANCE
+33 561 17 64 21 - [email protected]
*
--
--
* Guillaume PASERO
Ingénieur d'études et développement
Business Unit E-SPACE & Geo Information - Département
APPLICATIONS
CS Systèmes d'Information
Parc de la Grande Plaine - 5, Rue Brindejonc des Moulinais -
BP 15872
31506 Toulouse Cedex 05 - FRANCE
+33 561 17 64 21 - [email protected]
*
--
