Opened 10 years ago

Closed 10 years ago

Last modified 10 years ago

#5395 closed defect (fixed)

RPC GeoTransform Approximation Wrong

Reported by: warmerdam Owned by: warmerdam
Priority: normal Milestone: 1.11.0
Component: Algorithms Version: svn-trunk
Severity: normal Keywords: RPC
Cc:

Description (last modified by warmerdam)

The logic for computing the linear approximation of the RPC (adfGTFromLL) in GDALCreateRPCTransformer() is mixed up. This leads to situations where pixel/line to lat/long conversions do not converge as can be seen in the following example with the attached RPC in a VRT. 

$ gdaltransform  rpc_5395.vrt -t_srs WGS84 
0 0
130.85009097174 -7.59034547771687 0
warmerdam@gdal.cmo:~/pl/planet_common/t965
$ gdaltransform  -i rpc_5395.vrt -t_srs WGS84 
28.26163232 -27.79853245
1.24353157389123e-05 -2.78099425941036e-05 0

The mapping from 0,0 should be 28.26163232 -27.79853245 but comes out with a weird result which on debugging is due to a lack of convergence.

Attachments (1)

rpc_5395.vrt​ (2.3 KB ) - added by warmerdam 10 years ago.
RPC Demonstrating Issue

Download all attachments as: .zip

Change History (4)

by warmerdam, 10 years ago

Attachment: rpc_5395.vrt​ added

RPC Demonstrating Issue

comment:1 by warmerdam, 10 years ago

Description: modified (diff)
Status: new → assigned

comment:2 by warmerdam, 10 years ago

Resolution: → fixed
Status: assigned → closed

Fixed in trunk (r26963) with tests (r26962).

comment:3 by warmerdam, 10 years ago

This *ought* to be back ported.

Note: See TracTickets for help on using tickets.