[gdal-dev] NTF/JPEG2000 Question
Logan Greenlee
lgreenlee at digitalresultsgroup.com
Sun Jul 17 11:16:10 EDT 2011
Sorry - here are the GCP's:
Size is 9216, 7168
Coordinate System is `'
GCP Projection =
GEOGCS["WGS 84",
DATUM["WGS_1984",
SPHEROID["WGS 84",6378137,298.257223563,
AUTHORITY["EPSG","7030"]],
TOWGS84[0,0,0,0,0,0,0],
AUTHORITY["EPSG","6326"]],
PRIMEM["Greenwich",0,
AUTHORITY["EPSG","8901"]],
UNIT["degree",0.0174532925199433,
AUTHORITY["EPSG","9108"]],
AUTHORITY["EPSG","4326"]]
GCP[ 0]: Id=UpperLeft, Info=
(0.5,0.5) -> (-114.341666666667,33.3408333333333,0)
GCP[ 1]: Id=UpperRight, Info=
(9215.5,0.5) -> (-114.146111111111,33.3516666666667,0)
GCP[ 2]: Id=LowerRight, Info=
(9215.5,7167.5) -> (-114.146666666667,33.2147222222222,0)
GCP[ 3]: Id=LowerLeft, Info=
(0.5,7167.5) -> (-114.341388888889,33.2038888888889,0)
Is there any particular reason mapserver can't use the GCP's for display? Is it that they are not necessarily corner points of the image?
Thanks!
Logan
-----Original Message-----
From: Even Rouault [mailto:even.rouault at mines-paris.org]
Sent: Sunday, July 17, 2011 05:34
To: Logan Greenlee
Cc: gdal-dev at lists.osgeo.org
Subject: Re: [gdal-dev] NTF/JPEG2000 Question
Selon Logan Greenlee <lgreenlee at digitalresultsgroup.com>:
I think that in older GDAL versions the NITF driver was computing an approximate geotransform even if it is was not always accurate. Now the driver must be more careful and return GCPs when returning the geotransform wouldn't be appropriate.
One option to use this image in mapserver is to use gdalwarp to produce an image with a geotransform.
You didn't include the GCPs in your output of gdalinfo. It would have helped to evaluate if the behaviour is consistant or not.
> I'm not sure if this is the right venue - please direct me to another list.
>
> I am having an issue processing DigitalGlobe NTF/JPEG2000 imagery with
> GDAL and MapServer.
>
> I am currently using Version 1.8.1 MSVC2008 (Win64) -release with
> mapserver
> 6.0.1 found here:
>
> http://www.gisinternals.com/sdk/
>
> I also have an older GDAL version 1.6.0dev provided by an old install
> of FWTools.
>
> 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.
>
> 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.
>
> Here is the gdalinfo output from 1.8.1
>
> >gdalinfo img_1.NTF
> GDAL: Auto register C:\GDALMS\bin\gdal\plugins\gdal_BAG.dll using
> GDALRegister_B AG.
> GDAL: Auto register C:\GDALMS\bin\gdal\plugins\gdal_FITS.dll using
> GDALRegister_ FITS.
> GDAL: Auto register C:\GDALMS\bin\gdal\plugins\gdal_GMT.dll using
> GDALRegister_G MT.
> GDAL: Auto register C:\GDALMS\bin\gdal\plugins\gdal_HDF5.dll using
> GDALRegister_ HDF5.
> GDAL: Auto register C:\GDALMS\bin\gdal\plugins\gdal_HDF5Image.dll
> using GDALRegi ster_HDF5Image.
> GDAL: Auto register C:\GDALMS\bin\gdal\plugins\gdal_MrSID.dll using
> GDALRegister _MrSID.
> GDAL: Auto register C:\GDALMS\bin\gdal\plugins\gdal_netCDF.dll using
> GDALRegiste r_netCDF.
> GDAL: GDALOpen(/vsisubfile/3958_401263470,img_1.NTF,
> this=00000000037E0C40) succeeds as JP2OpenJPEG.
> GDAL: NITFDataset::Open() wasn't able to derive a first order
> geotransform. It will be returned as GCPs.
> GDAL: GDALOpen(img_1.NTF, this=0000000000C2C23
> 0) succeeds as NITF.
> Driver: NITF/National Imagery Transmission Format
> GDAL: GDALDefaultOverviews::OverviewScan()
>
> ... stuff ...
>
> Resulting corner coordinates:
>
> Corner Coordinates:
> Upper Left ( 0.0, 0.0)
> Lower Left ( 0.0, 7168.0)
> Upper Right ( 9216.0, 0.0)
> Lower Right ( 9216.0, 7168.0)
> Center ( 4608.0, 3584.0)
> Band 1 Block=1024x1024 Type=UInt16, ColorInterp=Undefined
> GDAL: GDALDefaultOverviews::OverviewScan()
> Band 2 Block=1024x1024 Type=UInt16, ColorInterp=Blue Band 3
> Block=1024x1024 Type=UInt16, ColorInterp=Green Band 4 Block=1024x1024
> Type=UInt16, ColorInterp=Undefined Band 5 Block=1024x1024 Type=UInt16,
> ColorInterp=Red Band 6 Block=1024x1024 Type=UInt16,
> ColorInterp=Undefined Band 7 Block=1024x1024 Type=UInt16,
> ColorInterp=Undefined Band 8 Block=1024x1024 Type=UInt16,
> ColorInterp=Undefined
>
>
> The GCP's appear correct.
>
>
> >gdalinfo --version
> GDAL 1.6.0dev, FWTools 2.2.8, released 2008/10/29
>
> gdalinfo img_1.NTF
> Warning 1: Image with 11 bits per sample will not be interpreted properly.
> Driver: NITF/National Imagery Transmission Format
> Files: img_1.NTF
> Size is 9216, 7168
> Coordinate System is:
> GEOGCS["WGS 84",
> DATUM["WGS_1984",
> SPHEROID["WGS 84",6378137,298.257223563,
> AUTHORITY["EPSG","7030"]],
> TOWGS84[0,0,0,0,0,0,0],
> AUTHORITY["EPSG","6326"]],
> PRIMEM["Greenwich",0,
> AUTHORITY["EPSG","8901"]],
> UNIT["degree",0.0174532925199433,
> AUTHORITY["EPSG","9108"]],
> AXIS["Lat",NORTH],
> AXIS["Long",EAST],
> AUTHORITY["EPSG","4326"]]
> GeoTransform =
> -114.3414689117551, 2.117622234280029e-005, -1.937894360779166e-008
> 33.34084229934283, 1.175619461023539e-006, -1.910763840441486e-005
> Metadata:
> NITF_FHDR=NITF02.10
>
> ... stuff ....
>
> Image Structure Metadata:
> COMPRESSION=JPEG2000
> RPC Metadata:
> LINE_OFF=3425
> LINE_SCALE=3525
> SAMP_OFF=4412
> SAMP_SCALE=4450
> LONG_OFF=-114.2441
> LONG_SCALE=0.0979
> LAT_OFF=33.2772
> LAT_SCALE=0.0742
> HEIGHT_OFF=518
> HEIGHT_SCALE=500
> MIN_LONG=-114.29305
> MAX_LONG=-114.19515
> MIN_LAT=33.2401
> MAX_LAT=33.3143
> Corner Coordinates:
> Upper Left (-114.3414689, 33.3408423) (114d20'29.29"W,
> 33d20'27.03"N) Lower Left (-114.3416078, 33.2038787)
> (114d20'29.79"W, 33d12'13.96"N) Upper Right (-114.1463088,
> 33.3516768) (114d 8'46.71"W, 33d21'6.04"N) Lower Right (-114.1464478, 33.2147133) (114d 8'47.21"W, 33d12'52.97"N)
> Center (-114.2439583, 33.2777778) (114d14'38.25"W, 33d16'40.00"N)
> Band 1 Block=9216x1 Type=UInt16, ColorInterp=Undefined
> Overviews: arbitrary
> Band 2 Block=9216x1 Type=UInt16, ColorInterp=Blue
> Overviews: arbitrary
> Band 3 Block=9216x1 Type=UInt16, ColorInterp=Green
> Overviews: arbitrary
> Band 4 Block=9216x1 Type=UInt16, ColorInterp=Undefined
> Overviews: arbitrary
> Band 5 Block=9216x1 Type=UInt16, ColorInterp=Red
> Overviews: arbitrary
> Band 6 Block=9216x1 Type=UInt16, ColorInterp=Undefined
> Overviews: arbitrary
> Band 7 Block=9216x1 Type=UInt16, ColorInterp=Undefined
> Overviews: arbitrary
> Band 8 Block=9216x1 Type=UInt16, ColorInterp=Undefined
> Overviews: arbitrary
>
> Any help would be appreciated!
>
More information about the gdal-dev
mailing list