Opened 9 years ago

Closed 9 years ago

Last modified 9 years ago

#3337 closed defect (fixed)

CADRG data in the polar regions does not display properly.

Reported by: gaopeng Owned by: warmerdam
Priority: normal Milestone: 1.7.2
Component: GDAL_Raster Version: 1.6.2
Severity: normal Keywords: NITF CADRG
Cc:

Description

Their projection isinterpreted as GCS WGS84. It should be North Pole Stereographic.

Attachments (2)

00000023.GN9 (124.8 KB) - added by gaopeng 9 years ago.
00000023.GNJ (124.8 KB) - added by gaopeng 9 years ago.

Download all attachments as: .zip

Change History (10)

Changed 9 years ago by gaopeng

Attachment: 00000023.GN9 added

comment:1 Changed 9 years ago by Even Rouault

I think the fix is in trunk. See #2940

comment:2 Changed 9 years ago by warmerdam

Status: newassigned

The relavent patch was r16750. Unfortunately it depends on being able to instantiate a coordinate transformer which is likely not possible in the ArcGIS environment (no PROJ.4). Consulting with Gao on how to proceed.

comment:3 Changed 9 years ago by Even Rouault

Frank, in r18640 (trunk), I've fixed a leak of pszGCPProjection that was never freed

comment:4 Changed 9 years ago by warmerdam

Milestone: 1.6.41.7.1

The transformation into AEQD projection was already in the 1.6-esri branch, however as noted it doesn't do anything in the context of ArcGIS since the transformer cannot be instantiated (nor PROJ.DLL).

As discussed with Gao, we will fallback to publishing the GCPPROJECTIONX metadata in the IMAGE_STRUCTURE domain similarly to what was done in the BSB driver. This essentially tells the projection the GCPs should be converted to in order to better represent the natural coordinate system of the image.

I have implemented this in trunk (r18637), 1.7 (r18638,r18643) and 1.6-esri (r18639, r18642).

In the process of working on this I also identified and corrected (#3347) which ultimately turned out not to be related.

comment:5 Changed 9 years ago by warmerdam

Resolution: fixed
Status: assignedclosed

Added a autotest of the polarsteregraphic case with and without PROJ (r18646) in trunk.

Changed 9 years ago by gaopeng

Attachment: 00000023.GNJ added

comment:6 Changed 9 years ago by gaopeng

Resolution: fixed
Status: closedreopened

I am reopening this ticket with a new test case. It works for the north pole data. However, for the south pole the data shows up as WGS1984 rather than South_Pole_Azimuthal_Equidistant. Attached is a test image.

comment:7 Changed 9 years ago by warmerdam

Keywords: CADRG added
Resolution: fixed
Status: reopenedclosed

The code was supporting "zone 9" files which are north pole. I added "zone J" to the filter per the specification and now southern polar files seems to work properly too. The fix is in trunk (r18837), 1.7 (r18838) and 1.6-esri (r18839).

comment:8 Changed 9 years ago by warmerdam

Milestone: 1.7.11.7.2
Note: See TracTickets for help on using tickets.