Opened 10 years ago

Closed 10 years ago

Last modified 10 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 10 years ago.
00000023.GNJ (124.8 KB) - added by gaopeng 10 years ago.

Download all attachments as: .zip

Change History (10)

Changed 10 years ago by gaopeng

Attachment: 00000023.GN9 added

comment:1 Changed 10 years ago by Even Rouault

I think the fix is in trunk. See #2940

comment:2 Changed 10 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 10 years ago by Even Rouault

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

comment:4 Changed 10 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 10 years ago by warmerdam

Resolution: fixed
Status: assignedclosed

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

Changed 10 years ago by gaopeng

Attachment: 00000023.GNJ added

comment:6 Changed 10 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 10 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 10 years ago by warmerdam

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