Opened 3 years ago

Closed 3 years ago

Last modified 3 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

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
$ 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 3 years ago.
RPC Demonstrating Issue

Download all attachments as: .zip

Change History (4)

Changed 3 years ago by warmerdam

Attachment: rpc_5395.vrt added

RPC Demonstrating Issue

comment:1 Changed 3 years ago by warmerdam

Description: modified (diff)
Status: newassigned

comment:2 Changed 3 years ago by warmerdam

Resolution: fixed
Status: assignedclosed

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

comment:3 Changed 3 years ago by warmerdam

This *ought* to be back ported.

Note: See TracTickets for help on using tickets.