[mapserver-users] Processing JP2000 files
EVANS, JAMES R GS-13 USAF ACC 84 RADES/SCZE
james.evans.4 at us.af.mil
Thu Mar 20 08:20:40 PDT 2014
Hi Jukka,
Here's the info I get on one of the files using the 64Bit build of
GISInternals that seems to completely fail to read the file:
Driver: JP2OpenJPEG/JPEG-2000 driver based on OpenJPEG library
Files: m_4111149_sw_12_1_20110810_20111011.jp2
Size is 5950, 7630
Coordinate System is:
PROJCS["NAD83 / UTM zone 12N",
GEOGCS["NAD83",
DATUM["North_American_Datum_1983",
SPHEROID["GRS 1980",6378137,298.2572221010002,
AUTHORITY["EPSG","7019"]],
AUTHORITY["EPSG","6269"]],
PRIMEM["Greenwich",0],
UNIT["degree",0.0174532925199433],
AUTHORITY["EPSG","4269"]],
PROJECTION["Transverse_Mercator"],
PARAMETER["latitude_of_origin",0],
PARAMETER["central_meridian",-111],
PARAMETER["scale_factor",0.9996],
PARAMETER["false_easting",500000],
PARAMETER["false_northing",0],
UNIT["metre",1,
AUTHORITY["EPSG","9001"]],
AUTHORITY["EPSG","26912"]]
Origin = (415740.000000000000000,4560370.000000000000000)
Pixel Size = (1.000000000000000,-1.000000000000000)
Image Structure Metadata:
INTERLEAVE=PIXEL
Corner Coordinates:
Upper Left ( 415740.000, 4560370.000) (112d 0'17.18"W, 41d11'25.11"N)
Lower Left ( 415740.000, 4552740.000) (112d 0'13.41"W, 41d 7'17.72"N)
Upper Right ( 421690.000, 4560370.000) (111d56' 1.79"W, 41d11'27.26"N)
Lower Right ( 421690.000, 4552740.000) (111d55'58.28"W, 41d 7'19.86"N)
Center ( 418715.000, 4556555.000) (111d58' 7.66"W, 41d 9'22.51"N)
Band 1 Block=4096x4096 Type=Byte, ColorInterp=Red
Overviews: 2975x3815, 1487x1907, 743x953, 371x476, 185x238
Overviews: arbitrary
Mask Flags: PER_DATASET ALPHA
Overviews of mask band: 2975x3815, 1487x1907, 743x953, 371x476, 185x238
Band 2 Block=4096x4096 Type=Byte, ColorInterp=Green
Overviews: 2975x3815, 1487x1907, 743x953, 371x476, 185x238
Overviews: arbitrary
Mask Flags: PER_DATASET ALPHA
Overviews of mask band: 2975x3815, 1487x1907, 743x953, 371x476, 185x238
Band 3 Block=4096x4096 Type=Byte, ColorInterp=Blue
Overviews: 2975x3815, 1487x1907, 743x953, 371x476, 185x238
Overviews: arbitrary
Mask Flags: PER_DATASET ALPHA
Overviews of mask band: 2975x3815, 1487x1907, 743x953, 371x476, 185x238
Band 4 Block=4096x4096 Type=Byte, ColorInterp=Alpha
Overviews: 2975x3815, 1487x1907, 743x953, 371x476, 185x238
Overviews: arbitrary
Here's the info from the MS4W released version that has been working
flawlessly for months:
Driver: JP2MrSID/MrSID JPEG2000
Files: m_3609464_sw_15_1_20100622_20101021.jp2
Size is 6311, 7598
Coordinate System is:
PROJCS["NAD83 / UTM zone 15N",
GEOGCS["NAD83",
DATUM["North_American_Datum_1983",
SPHEROID["GRS 1980",6378137,298.2572221010002,
AUTHORITY["EPSG","7019"]],
AUTHORITY["EPSG","6269"]],
PRIMEM["Greenwich",0],
UNIT["degree",0.0174532925199433],
AUTHORITY["EPSG","4269"]],
PROJECTION["Transverse_Mercator"],
PARAMETER["latitude_of_origin",0],
PARAMETER["central_meridian",-93],
PARAMETER["scale_factor",0.9996],
PARAMETER["false_easting",500000],
PARAMETER["false_northing",0],
UNIT["metre",1,
AUTHORITY["EPSG","9001"]],
AUTHORITY["EPSG","26915"]]
Origin = (398304.000000000000000,3991767.000000000000000)
Pixel Size = (1.000000000000000,-1.000000000000000)
Metadata:
GEOTIFF_CHAR__GTModelTypeGeoKey=ModelTypeProjected
GEOTIFF_CHAR__GTRasterTypeGeoKey=RasterPixelIsArea
GEOTIFF_CHAR__ProjectedCSTypeGeoKey=PCS_NAD83_UTM_zone_15N
GEOTIFF_CHAR__ProjLinearUnitsGeoKey=Linear_Meter
GEOTIFF_NUM__1024__GTModelTypeGeoKey=1
GEOTIFF_NUM__1025__GTRasterTypeGeoKey=1
GEOTIFF_NUM__1026__GTCitationGeoKey=NAIP 2010 m_3609464_sw_15_1_20100622
GEOTIFF_NUM__3072__ProjectedCSTypeGeoKey=26915
GEOTIFF_NUM__3073__PCSCitationGeoKey=NAD83 / UTM zone 15N
GEOTIFF_NUM__3076__ProjLinearUnitsGeoKey=9001
IMAGE__ENCODING_APPLICATION=GeoExpress 7.0.0.2052
IMAGE__INPUT_FILE_SIZE=191865212.000000
IMAGE__INPUT_FORMAT=GeoTIFF
IMAGE__INPUT_NAME=/hsm/naip/36094/64/m_3609464_sw_15_1_20100622_20101021.tif
Image Structure Metadata:
INTERLEAVE=PIXEL
Corner Coordinates:
Upper Left ( 398304.000, 3991767.000) ( 94d 7'45.35"W, 36d 3'54.62"N)
Lower Left ( 398304.000, 3984169.000) ( 94d 7'41.83"W, 35d59'48.06"N)
Upper Right ( 404615.000, 3991767.000) ( 94d 3'33.10"W, 36d 3'56.92"N)
Lower Right ( 404615.000, 3984169.000) ( 94d 3'29.80"W, 35d59'50.35"N)
Center ( 401459.500, 3987968.000) ( 94d 5'37.52"W, 36d 1'52.51"N)
Band 1 Block=1024x128 Type=Byte, ColorInterp=Undefined
Minimum=49.000, Maximum=232.000, Mean=121.048, StdDev=29.936
Overviews: 3156x3799, 1578x1900, 789x950, 395x475, 198x238, 99x119, 50x60
Band 2 Block=1024x128 Type=Byte, ColorInterp=Undefined
Minimum=65.000, Maximum=219.000, Mean=134.718, StdDev=24.551
Overviews: 3156x3799, 1578x1900, 789x950, 395x475, 198x238, 99x119, 50x60
Band 3 Block=1024x128 Type=Byte, ColorInterp=Undefined
Minimum=82.000, Maximum=219.000, Mean=132.334, StdDev=17.415
Overviews: 3156x3799, 1578x1900, 789x950, 395x475, 198x238, 99x119, 50x60
Band 4 Block=1024x128 Type=Byte, ColorInterp=Undefined
Minimum=39.000, Maximum=255.000, Mean=186.412, StdDev=22.709
Overviews: 3156x3799, 1578x1900, 789x950, 395x475, 198x238, 99x119, 50x60
I'll have to check the version I'm using at home tonight. It came with the
32bit build of GISInternals, and can read the files, but they are washed
out, and boundaries of the tiles show some blended overlap.
James
-----Original Message-----
From: Rahkonen Jukka (Tike) [mailto:jukka.rahkonen at mmmtike.fi]
Sent: Thursday, March 20, 2014 7:45 AM
To: EVANS, JAMES R GS-13 USAF ACC 84 RADES/SCZE;
mapserver-users at lists.osgeo.org
Subject: VS: [mapserver-users] Processing JP2000 files
Hi,
Are your JPEG2000 images the same? SOT marker is JPEG2000 stuff and it
means "start of tilepart". If images are the same then the message means
that something has changed in GDAL part. GDAL has several JPEG200 drivers
and they all behave in a bit different way. Check which driver you have in
your old system and which one you have now with "gdalinfo --formats".
I warn that from the free alternatives only OpenJPEG driver is something
that can be considered for real use.
-Jukka Rahkonen-
James_in_Utah wrote:
>
> There must have been something wrong with my shapefile, because I
> tried this at home and I was able to get the same data to work. Now
> the problem is that the image I see on my client is all washed out,
> not a variant RGB image but mostly brownish with weak greens. None of
> the reds or blues are making it through. Plus, where the tiles
> overlap I can see the edges are blended from the overlapping tiles.
> This wasn't the case with the older GDAL from MS4W. Has something
changed? Is there some directive in the layer that I'm missing?
> Thanks,
> James
>
>
>
>
> --
> View this message in context: http://osgeo-
> org.1560.x6.nabble.com/Processing-JP2000-files-tp5129997p5130136.html
> Sent from the Mapserver - User mailing list archive at Nabble.com.
> _______________________________________________
> mapserver-users mailing list
> mapserver-users at lists.osgeo.org
> http://lists.osgeo.org/mailman/listinfo/mapserver-users
-------------- next part --------------
A non-text attachment was scrubbed...
Name: smime.p7s
Type: application/pkcs7-signature
Size: 6912 bytes
Desc: not available
URL: <http://lists.osgeo.org/pipermail/mapserver-users/attachments/20140320/7fbd0132/attachment-0001.bin>
More information about the mapserver-users
mailing list