<html xmlns:v="urn:schemas-microsoft-com:vml" xmlns:o="urn:schemas-microsoft-com:office:office" xmlns:w="urn:schemas-microsoft-com:office:word" xmlns:m="http://schemas.microsoft.com/office/2004/12/omml" xmlns="http://www.w3.org/TR/REC-html40"><head><meta http-equiv=Content-Type content="text/html; charset=us-ascii"><meta name=Generator content="Microsoft Word 14 (filtered medium)"><style><!--
/* Font Definitions */
@font-face
{font-family:"Cambria Math";
panose-1:2 4 5 3 5 4 6 3 2 4;}
@font-face
{font-family:Calibri;
panose-1:2 15 5 2 2 2 4 3 2 4;}
/* Style Definitions */
p.MsoNormal, li.MsoNormal, div.MsoNormal
{margin:0in;
margin-bottom:.0001pt;
font-size:11.0pt;
font-family:"Calibri","sans-serif";}
a:link, span.MsoHyperlink
{mso-style-priority:99;
color:blue;
text-decoration:underline;}
a:visited, span.MsoHyperlinkFollowed
{mso-style-priority:99;
color:purple;
text-decoration:underline;}
span.EmailStyle17
{mso-style-type:personal-compose;
font-family:"Calibri","sans-serif";
color:windowtext;}
.MsoChpDefault
{mso-style-type:export-only;
font-family:"Calibri","sans-serif";}
@page WordSection1
{size:8.5in 11.0in;
margin:1.0in 1.0in 1.0in 1.0in;}
div.WordSection1
{page:WordSection1;}
--></style><!--[if gte mso 9]><xml>
<o:shapedefaults v:ext="edit" spidmax="1026" />
</xml><![endif]--><!--[if gte mso 9]><xml>
<o:shapelayout v:ext="edit">
<o:idmap v:ext="edit" data="1" />
</o:shapelayout></xml><![endif]--></head><body lang=EN-US link=blue vlink=purple><div class=WordSection1><p class=MsoNormal>I’m not sure if this is the right venue – please direct me to another list.<o:p></o:p></p><p class=MsoNormal><o:p> </o:p></p><p class=MsoNormal>I am having an issue processing DigitalGlobe NTF/JPEG2000 imagery with GDAL and MapServer.<o:p></o:p></p><p class=MsoNormal><o:p> </o:p></p><p class=MsoNormal>I am currently using Version 1.8.1 MSVC2008 (Win64) –release with mapserver 6.0.1 found here:<o:p></o:p></p><p class=MsoNormal><o:p> </o:p></p><p class=MsoNormal><a href="http://www.gisinternals.com/sdk/">http://www.gisinternals.com/sdk/</a><o:p></o:p></p><p class=MsoNormal><o:p> </o:p></p><p class=MsoNormal>I also have an older GDAL version 1.6.0dev provided by an old install of FWTools.<o:p></o:p></p><p class=MsoNormal><o:p> </o:p></p><p class=MsoNormal>I have found that there is a fairly significant difference between the two versions in the way that they handle NTF/JPEG2000 imagery from digital globe.<o:p></o:p></p><p class=MsoNormal><o:p> </o:p></p><p class=MsoNormal>The issue that I am having is that 1.6.0 properly displays the four corner points of the image. Version 1.8.1 does not dispalay, and instead displays the GCP points. The GCP points look correct but the image does not render correctly using shp2img from map server. It uses the corrdinate system that is shown below which is pixel relative.<o:p></o:p></p><p class=MsoNormal><o:p> </o:p></p><p class=MsoNormal>Here is the gdalinfo output from 1.8.1<o:p></o:p></p><p class=MsoNormal><o:p> </o:p></p><p class=MsoNormal>>gdalinfo img_1.NTF<o:p></o:p></p><p class=MsoNormal>GDAL: Auto register C:\GDALMS\bin\gdal\plugins\gdal_BAG.dll using GDALRegister_B<o:p></o:p></p><p class=MsoNormal>AG.<o:p></o:p></p><p class=MsoNormal>GDAL: Auto register C:\GDALMS\bin\gdal\plugins\gdal_FITS.dll using GDALRegister_<o:p></o:p></p><p class=MsoNormal>FITS.<o:p></o:p></p><p class=MsoNormal>GDAL: Auto register C:\GDALMS\bin\gdal\plugins\gdal_GMT.dll using GDALRegister_G<o:p></o:p></p><p class=MsoNormal>MT.<o:p></o:p></p><p class=MsoNormal>GDAL: Auto register C:\GDALMS\bin\gdal\plugins\gdal_HDF5.dll using GDALRegister_<o:p></o:p></p><p class=MsoNormal>HDF5.<o:p></o:p></p><p class=MsoNormal>GDAL: Auto register C:\GDALMS\bin\gdal\plugins\gdal_HDF5Image.dll using GDALRegi<o:p></o:p></p><p class=MsoNormal>ster_HDF5Image.<o:p></o:p></p><p class=MsoNormal>GDAL: Auto register C:\GDALMS\bin\gdal\plugins\gdal_MrSID.dll using GDALRegister<o:p></o:p></p><p class=MsoNormal>_MrSID.<o:p></o:p></p><p class=MsoNormal>GDAL: Auto register C:\GDALMS\bin\gdal\plugins\gdal_netCDF.dll using GDALRegiste<o:p></o:p></p><p class=MsoNormal>r_netCDF.<o:p></o:p></p><p class=MsoNormal>GDAL: GDALOpen(/vsisubfile/3958_401263470,img_1.NTF, this=00000000037E0C40) succeeds as JP2OpenJPEG.<o:p></o:p></p><p class=MsoNormal>GDAL: NITFDataset::Open() wasn't able to derive a first order<o:p></o:p></p><p class=MsoNormal>geotransform. It will be returned as GCPs.<o:p></o:p></p><p class=MsoNormal>GDAL: GDALOpen(img_1.NTF, this=0000000000C2C23<o:p></o:p></p><p class=MsoNormal>0) succeeds as NITF.<o:p></o:p></p><p class=MsoNormal>Driver: NITF/National Imagery Transmission Format<o:p></o:p></p><p class=MsoNormal>GDAL: GDALDefaultOverviews::OverviewScan()<o:p></o:p></p><p class=MsoNormal><o:p> </o:p></p><p class=MsoNormal>… stuff …<o:p></o:p></p><p class=MsoNormal><o:p> </o:p></p><p class=MsoNormal>Resulting corner coordinates:<o:p></o:p></p><p class=MsoNormal><o:p> </o:p></p><p class=MsoNormal>Corner Coordinates:<o:p></o:p></p><p class=MsoNormal>Upper Left ( 0.0, 0.0)<o:p></o:p></p><p class=MsoNormal>Lower Left ( 0.0, 7168.0)<o:p></o:p></p><p class=MsoNormal>Upper Right ( 9216.0, 0.0)<o:p></o:p></p><p class=MsoNormal>Lower Right ( 9216.0, 7168.0)<o:p></o:p></p><p class=MsoNormal>Center ( 4608.0, 3584.0)<o:p></o:p></p><p class=MsoNormal>Band 1 Block=1024x1024 Type=UInt16, ColorInterp=Undefined<o:p></o:p></p><p class=MsoNormal>GDAL: GDALDefaultOverviews::OverviewScan()<o:p></o:p></p><p class=MsoNormal>Band 2 Block=1024x1024 Type=UInt16, ColorInterp=Blue<o:p></o:p></p><p class=MsoNormal>Band 3 Block=1024x1024 Type=UInt16, ColorInterp=Green<o:p></o:p></p><p class=MsoNormal>Band 4 Block=1024x1024 Type=UInt16, ColorInterp=Undefined<o:p></o:p></p><p class=MsoNormal>Band 5 Block=1024x1024 Type=UInt16, ColorInterp=Red<o:p></o:p></p><p class=MsoNormal>Band 6 Block=1024x1024 Type=UInt16, ColorInterp=Undefined<o:p></o:p></p><p class=MsoNormal>Band 7 Block=1024x1024 Type=UInt16, ColorInterp=Undefined<o:p></o:p></p><p class=MsoNormal>Band 8 Block=1024x1024 Type=UInt16, ColorInterp=Undefined<o:p></o:p></p><p class=MsoNormal><o:p> </o:p></p><p class=MsoNormal><o:p> </o:p></p><p class=MsoNormal>The GCP’s appear correct.<o:p></o:p></p><p class=MsoNormal><o:p> </o:p></p><p class=MsoNormal><o:p> </o:p></p><p class=MsoNormal>>gdalinfo --version<o:p></o:p></p><p class=MsoNormal>GDAL 1.6.0dev, FWTools 2.2.8, released 2008/10/29<o:p></o:p></p><p class=MsoNormal><o:p> </o:p></p><p class=MsoNormal>gdalinfo img_1.NTF<o:p></o:p></p><p class=MsoNormal>Warning 1: Image with 11 bits per sample will not be interpreted properly.<o:p></o:p></p><p class=MsoNormal>Driver: NITF/National Imagery Transmission Format<o:p></o:p></p><p class=MsoNormal>Files: img_1.NTF<o:p></o:p></p><p class=MsoNormal>Size is 9216, 7168<o:p></o:p></p><p class=MsoNormal>Coordinate System is:<o:p></o:p></p><p class=MsoNormal>GEOGCS["WGS 84",<o:p></o:p></p><p class=MsoNormal> DATUM["WGS_1984",<o:p></o:p></p><p class=MsoNormal> SPHEROID["WGS 84",6378137,298.257223563,<o:p></o:p></p><p class=MsoNormal> AUTHORITY["EPSG","7030"]],<o:p></o:p></p><p class=MsoNormal> TOWGS84[0,0,0,0,0,0,0],<o:p></o:p></p><p class=MsoNormal> AUTHORITY["EPSG","6326"]],<o:p></o:p></p><p class=MsoNormal> PRIMEM["Greenwich",0,<o:p></o:p></p><p class=MsoNormal> AUTHORITY["EPSG","8901"]],<o:p></o:p></p><p class=MsoNormal> UNIT["degree",0.0174532925199433,<o:p></o:p></p><p class=MsoNormal> AUTHORITY["EPSG","9108"]],<o:p></o:p></p><p class=MsoNormal> AXIS["Lat",NORTH],<o:p></o:p></p><p class=MsoNormal> AXIS["Long",EAST],<o:p></o:p></p><p class=MsoNormal> AUTHORITY["EPSG","4326"]]<o:p></o:p></p><p class=MsoNormal>GeoTransform =<o:p></o:p></p><p class=MsoNormal> -114.3414689117551, 2.117622234280029e-005, -1.937894360779166e-008<o:p></o:p></p><p class=MsoNormal> 33.34084229934283, 1.175619461023539e-006, -1.910763840441486e-005<o:p></o:p></p><p class=MsoNormal>Metadata:<o:p></o:p></p><p class=MsoNormal> NITF_FHDR=NITF02.10 <o:p></o:p></p><p class=MsoNormal><o:p> </o:p></p><p class=MsoNormal> … stuff ….<o:p></o:p></p><p class=MsoNormal><o:p> </o:p></p><p class=MsoNormal>Image Structure Metadata:<o:p></o:p></p><p class=MsoNormal> COMPRESSION=JPEG2000<o:p></o:p></p><p class=MsoNormal>RPC Metadata:<o:p></o:p></p><p class=MsoNormal> LINE_OFF=3425<o:p></o:p></p><p class=MsoNormal> LINE_SCALE=3525<o:p></o:p></p><p class=MsoNormal> SAMP_OFF=4412<o:p></o:p></p><p class=MsoNormal> SAMP_SCALE=4450<o:p></o:p></p><p class=MsoNormal> LONG_OFF=-114.2441<o:p></o:p></p><p class=MsoNormal> LONG_SCALE=0.0979<o:p></o:p></p><p class=MsoNormal> LAT_OFF=33.2772<o:p></o:p></p><p class=MsoNormal> LAT_SCALE=0.0742<o:p></o:p></p><p class=MsoNormal> HEIGHT_OFF=518<o:p></o:p></p><p class=MsoNormal> HEIGHT_SCALE=500<o:p></o:p></p><p class=MsoNormal> MIN_LONG=-114.29305<o:p></o:p></p><p class=MsoNormal> MAX_LONG=-114.19515<o:p></o:p></p><p class=MsoNormal> MIN_LAT=33.2401<o:p></o:p></p><p class=MsoNormal> MAX_LAT=33.3143<o:p></o:p></p><p class=MsoNormal>Corner Coordinates:<o:p></o:p></p><p class=MsoNormal>Upper Left (-114.3414689, 33.3408423) (114d20'29.29"W, 33d20'27.03"N)<o:p></o:p></p><p class=MsoNormal>Lower Left (-114.3416078, 33.2038787) (114d20'29.79"W, 33d12'13.96"N)<o:p></o:p></p><p class=MsoNormal>Upper Right (-114.1463088, 33.3516768) (114d 8'46.71"W, 33d21'6.04"N)<o:p></o:p></p><p class=MsoNormal>Lower Right (-114.1464478, 33.2147133) (114d 8'47.21"W, 33d12'52.97"N)<o:p></o:p></p><p class=MsoNormal>Center (-114.2439583, 33.2777778) (114d14'38.25"W, 33d16'40.00"N)<o:p></o:p></p><p class=MsoNormal>Band 1 Block=9216x1 Type=UInt16, ColorInterp=Undefined<o:p></o:p></p><p class=MsoNormal> Overviews: arbitrary<o:p></o:p></p><p class=MsoNormal>Band 2 Block=9216x1 Type=UInt16, ColorInterp=Blue<o:p></o:p></p><p class=MsoNormal> Overviews: arbitrary<o:p></o:p></p><p class=MsoNormal>Band 3 Block=9216x1 Type=UInt16, ColorInterp=Green<o:p></o:p></p><p class=MsoNormal> Overviews: arbitrary<o:p></o:p></p><p class=MsoNormal>Band 4 Block=9216x1 Type=UInt16, ColorInterp=Undefined<o:p></o:p></p><p class=MsoNormal> Overviews: arbitrary<o:p></o:p></p><p class=MsoNormal>Band 5 Block=9216x1 Type=UInt16, ColorInterp=Red<o:p></o:p></p><p class=MsoNormal> Overviews: arbitrary<o:p></o:p></p><p class=MsoNormal>Band 6 Block=9216x1 Type=UInt16, ColorInterp=Undefined<o:p></o:p></p><p class=MsoNormal> Overviews: arbitrary<o:p></o:p></p><p class=MsoNormal>Band 7 Block=9216x1 Type=UInt16, ColorInterp=Undefined<o:p></o:p></p><p class=MsoNormal> Overviews: arbitrary<o:p></o:p></p><p class=MsoNormal>Band 8 Block=9216x1 Type=UInt16, ColorInterp=Undefined<o:p></o:p></p><p class=MsoNormal> Overviews: arbitrary<o:p></o:p></p><p class=MsoNormal><o:p> </o:p></p><p class=MsoNormal>Any help would be appreciated!<o:p></o:p></p></div></body></html>