Opened 13 years ago
Closed 12 years ago
#3709 closed defect (fixed)
Georeference problem with Envisat Images
Reported by: | timatkinson | Owned by: | chaitanya |
---|---|---|---|
Priority: | high | Milestone: | 1.8.1 |
Component: | GDAL_Raster | Version: | unspecified |
Severity: | critical | Keywords: | Envisat GeoReference |
Cc: | warmerdam, Even Rouault, antonio |
Description
Gdal 1.6.2 (as featured in Fedora 13) and gdal 1.7.2 (as featured in Ubuntu 10.4) both have a problem with Georeference of Envisat Images
It appears they are reading both the Geolocation grid ADS (correct) and the Doppler Centroid Grid (incorrect) as GCPs.
The numbers in the grids are similar but actually different and hence you get two sets of ground coordinates per one GCP Line,Column
Run the command line below (*.N1 any envisat file will do)
gdalinfo *.N1 | grep \(0.5,0.5\) and you'll instantly see the problem
Change History (10)
comment:1 by , 13 years ago
Cc: | added |
---|---|
Component: | default → GDAL_Raster |
comment:2 by , 13 years ago
Tim,
Can you provide some sample data? You can email it to me (chaitanya.ch * gmail.com) if you want it to be confidential.
comment:3 by , 13 years ago
Cc: | added; removed |
---|---|
Owner: | changed from | to
Status: | new → assigned |
I was able to reproduce the problem with a sample dataset from Anders Berg. I'm working on this.
comment:4 by , 13 years ago
Severity: | normal → critical |
---|
For Envisat ASAR, the mistake is that the parameter "line_num" of Geolocation Grid is used as the Y-coordinate (tested with gdalinfo v1.7). In the ESA documentation for ASAR Level 1 products, there is a warning that "line_num" is the line number relative to the "granule within the slice", and not of the image itself. http://envisat.esa.int/handbooks/asar/CNTR6-6-9.htm#eph.asar.asardf.asarrec.ASAR_Geo_Grid_ADSR
Instead, the Y-coordinate should be calculated as a cumulative sum from the "num_lines" parameter: [1 cumsum(num_lines)] in Matlab-notation.
Note that Ticket #3160 refers to the same issue. A patch is posted there, but I am not sure if that solves the issue.
I propose to increase the priority of this ticket, as Envisat is one of the major satellites, and the data are not very useful without georeferencing.
comment:5 by , 13 years ago
Priority: | normal → high |
---|
comment:6 by , 12 years ago
Cc: | added |
---|
comment:7 by , 12 years ago
(I see you've cc'ed me. Just to mention I have no particular competence/interest on this file format and don't intend to take any particular action on this)
comment:8 by , 12 years ago
Cc: | added |
---|
comment:9 by , 12 years ago
I have recompiled GDAL with the patch which was posted by Antonio in Ticket #3160, and can confirm that it solves the problem.
Are there any developers who would like to apply this patch to the official repository, and then close these two tickets (#3160 and #3709)?
http://trac.osgeo.org/gdal/attachment/ticket/3160/gdal-envisat-gcp.patch
comment:10 by , 12 years ago
Milestone: | → 1.8.1 |
---|---|
Resolution: | → fixed |
Status: | assigned → closed |
Chaitanya,
Could you see if you can find envisat data to replicate this? It may be hard to find some.