JPEG 2000 format instead of TIFF for a RASTER layer

Ed McNierney ed at TOPOZONE.COM
Thu Aug 26 05:09:53 PDT 2004


Brant -
 
You've posted two LAYER definitions - what's the second one for?  The
default field name used to look up the name of the image file is not
"PATH" - that's why there's a TILEITEM statement in the second LAYER
definition you've posted.  But there is no such TILEITEM statement in
the LAYER where you're trying to use the JPEG2000 files.
 
    - Ed
Ed McNierney
President and Chief Mapmaker
TopoZone.com / Maps a la carte, Inc.
73 Princeton Street, Suite 305
North Chelmsford, MA  01863
Phone: +1 978 251-4242   Fax: +1 978 251-1396 


________________________________

From: UMN MapServer Users List [mailto:MAPSERVER-USERS at LISTS.UMN.EDU] On
Behalf Of Brant Anderson
Sent: Thursday, August 26, 2004 3:34 AM
To: MAPSERVER-USERS at LISTS.UMN.EDU
Subject: Re: [UMN_MAPSERVER-USERS] JPEG 2000 format instead of TIFF for
a RASTER layer


Here is a line from the apache web server log file:

[Thu Aug 26 02:29:45 2004] [error] [client 192.168.1.101] malformed
header from
script. Bad header=GeoJp2 info found...: /usr/lib/cgi-bin/mapserv421


Maybe that has something to do with it. If so, any ideas how to fix
this?

Brant



Brant Anderson wrote: 

        Here are some more details:
        Mapfile: ( Full mapfile can be found here:
http://mudlabs.homelinux.org:8181/mapserver_dev/Douglas.map )
        ---->snip<----
          LAYER
            NAME 'index'
            TYPE TILEINDEX
            DATA 'dougIndex_jp2.shp' # JPEG2000 Images in same folder as
Shapefiles (TESTING)
            #DATA 'dougIndex_ortho.shp' # TIFF Images in 'ortho'
sub-folder (WORKING)
            #DATA 'dougIndex.shp' # TIFF Images in same folder as
Shapefiles (WORKING)
          END
	
          LAYER
            NAME 'test_raster'
            STATUS OFF
            TYPE RASTER
            TILEINDEX 'index'
            TILEITEM 'PATH'
          END
        ---->snip<----
	
        The 'tile index' layer was created via an ArcGIS extension
written by Mike Juvrud of Mud Labs. What it does is takes the
Orthophoto/raster image from a MrSID format, and chops it up into tiles,
outputting the tiles as TIFF formatted images (max 50MB per image). At
the same time, it generates the 'tile index' layer, filling in the
corresponding filename in an attribute named 'PATH'.
	
        After verifying this method worked using TIFF files, the 'PATH'
attribute was changed from a .tif extension to a .jp2 extension. So if I
had an image named '104.tif' the new name would be '104.jp2'
	
        I used the 'jasper' utility to convert the TIFF image to a
JPEG2000 using the command: 'jasper -f 104.tif -t tif -F 104.jp2 -T jp2'
	
        Here is the output from ogrinfo of a TIFF version, and another
of the JPEG2000 version.
	
        ==TIFF: (command used: 'ogrinfo dougIndex.shp dougIndex')
        INFO: Open of `dougIndex.shp'
        using driver `ESRI Shapefile' successful.
         
        Layer name: dougIndex
        Geometry: Polygon
        Feature Count: 7
        Extent: (285931.949063, 5109160.000042) - (289100.000009,
5109580.000043)
        Layer SRS WKT:
        (unknown)
        PATH: String (15.0)
        STATUS: Integer (4.0)
        IndexID: Integer (9.0)
        OGRFeature(dougIndex):0
          PATH (String) = 104.tif
          STATUS (Integer) = 1
          IndexID (Integer) = 104
          POLYGON ((286040.000 5109160.000,285931.949
5109160.000,285945.481 5109580.000,286040.000 5109580.000,286040.000
5109160.000))
         
        OGRFeature(dougIndex):1
          PATH (String) = 105.tif
          STATUS (Integer) = 1
          IndexID (Integer) = 105
          POLYGON ((286550.000 5109160.000,286040.000
5109160.000,286040.000 5109580.000,286550.000 5109580.000,286550.000
5109160.000))
         
        OGRFeature(dougIndex):2
          PATH (String) = 106.tif
          STATUS (Integer) = 1
          IndexID (Integer) = 106
          POLYGON ((287060.000 5109160.000,286550.000
5109160.000,286550.000 5109580.000,287060.000 5109580.000,287060.000
5109160.000))
         
        OGRFeature(dougIndex):3
          PATH (String) = 107.tif
          STATUS (Integer) = 1
          IndexID (Integer) = 107
          POLYGON ((287570.000 5109160.000,287060.000
5109160.000,287060.000 5109580.000,287570.000 5109580.000,287570.000
5109160.000))
         
        OGRFeature(dougIndex):4
          PATH (String) = 108.tif
          STATUS (Integer) = 1
          IndexID (Integer) = 108
          POLYGON ((288080.000 5109160.000,287570.000
5109160.000,287570.000 5109580.000,288080.000 5109580.000,288080.000
5109160.000))
         
        OGRFeature(dougIndex):5
          PATH (String) = 109.tif
          STATUS (Integer) = 1
          IndexID (Integer) = 109
          POLYGON ((288590.000 5109160.000,288080.000
5109160.000,288080.000 5109580.000,288590.000 5109580.000,288590.000
5109160.000))
         
        OGRFeature(dougIndex):6
          PATH (String) = 110.tif
          STATUS (Integer) = 1
          IndexID (Integer) = 110
          POLYGON ((289100.000 5109160.000,288590.000
5109160.000,288590.000 5109580.000,289100.000 5109580.000,289100.000
5109160.000))
        ==========
	
        ==JPEG2000: (command used: 'ogrinfo dougIndex_jp2.shp
dougIndex_jp2')
        INFO: Open of `dougIndex_jp2.shp'
        using driver `ESRI Shapefile' successful.
         
        Layer name: dougIndex_jp2
        Geometry: Polygon
        Feature Count: 7
        Extent: (286040.000007, 5109160.000042) - (289100.000009,
5109580.000043)
        Layer SRS WKT:
        (unknown)
        PATH: String (15.0)
        STATUS: Integer (4.0)
        IndexID: Integer (9.0)
        OGRFeature(dougIndex):0
          PATH (String) = 104.jp2
          STATUS (Integer) = 1
          IndexID (Integer) = 104
          POLYGON ((286040.000 5109160.000,285931.949
5109160.000,285945.481 5109580.000,286040.000 5109580.000,286040.000
5109160.000))
	
        OGRFeature(dougIndex_jp2):1
          PATH (String) = 105.jp2
          STATUS (Integer) = 1
          IndexID (Integer) = 105
          POLYGON ((286550.000 5109160.000,286040.000
5109160.000,286040.000 5109580.000,286550.000 5109580.000,286550.000
5109160.000))
         
        OGRFeature(dougIndex_jp2):2
          PATH (String) = 106.jp2
          STATUS (Integer) = 1
          IndexID (Integer) = 106
          POLYGON ((287060.000 5109160.000,286550.000
5109160.000,286550.000 5109580.000,287060.000 5109580.000,287060.000
5109160.000))
         
        OGRFeature(dougIndex_jp2):3
          PATH (String) = 107.jp2
          STATUS (Integer) = 1
          IndexID (Integer) = 107
          POLYGON ((287570.000 5109160.000,287060.000
5109160.000,287060.000 5109580.000,287570.000 5109580.000,287570.000
5109160.000))
         
        OGRFeature(dougIndex_jp2):4
          PATH (String) = 108.jp2
          STATUS (Integer) = 1
          IndexID (Integer) = 108
          POLYGON ((288080.000 5109160.000,287570.000
5109160.000,287570.000 5109580.000,288080.000 5109580.000,288080.000
5109160.000))
         
        OGRFeature(dougIndex_jp2):5
          PATH (String) = 109.jp2
          STATUS (Integer) = 1
          IndexID (Integer) = 109
          POLYGON ((288590.000 5109160.000,288080.000
5109160.000,288080.000 5109580.000,288590.000 5109580.000,288590.000
5109160.000))
         
        OGRFeature(dougIndex_jp2):6
          PATH (String) = 110.jp2
          STATUS (Integer) = 1
          IndexID (Integer) = 110
          POLYGON ((289100.000 5109160.000,288590.000
5109160.000,288590.000 5109580.000,289100.000 5109580.000,289100.000
5109160.000))
        ==========
	
	
        I hope I covered all the bases. If you need any more info, sling
your questions my way.
	
	
        Thank you,
	
        Brant Anderson
	
	
	

                Ed McNierney wrote: 

                        Brant -
                         
                        You mentioned "tile", so I assume you're using a
TILEINDEX.  You mention changing the PATH attribute of the shapefile -
could you be a little more specific?  Could you post the LAYER
definition in your MAP file for this layer?  I would expect that you
would use gdaltindex to create a new TILEINDEX shapefile with the new
filenames in it (remember that gdaltindex will add to an existing file,
so you'll want to create a new one).  Then you would change the
TILEINDEX statement in your MAP file to refer to this new shapefile.
                         
                            - Ed
                        Ed McNierney
                        President and Chief Mapmaker
                        TopoZone.com / Maps a la carte, Inc.
                        73 Princeton Street, Suite 305
                        North Chelmsford, MA  01863
                        Phone: +1 978 251-4242   Fax: +1 978 251-1396 


________________________________

                        From: UMN MapServer Users List
[mailto:MAPSERVER-USERS at LISTS.UMN.EDU] On Behalf Of Brant Anderson
                        Sent: Thursday, August 26, 2004 1:21 AM
                        To: MAPSERVER-USERS at LISTS.UMN.EDU
                        Subject: [UMN_MAPSERVER-USERS] JPEG 2000 format
instead of TIFF for a RASTER layer
                	
                	
                        I am trying to tile ortho photos using the JPEG
2000 (jp2) format. I have created a working tif version.
                	
                        I converted the TIF images to JP2 files using
jasper. Then I updated the PATH attribute of the shapefile to reflect
the new file extension. This, however, created unexpected results. I am
getting an Internal Server Error.
                	
                        Here is my MapServer output: (mapserv -v)
                        MapServer version 4.2.2 OUTPUT=PNG OUTPUT=JPEG
OUTPUT=WBMP OUTPUT=PDF SUPPORTS=PROJ SUPPORTS=FREETYPE
SUPPORTS=WMS_SERVER INPUT=EPPL7 INPUT=OGR INPUT=GDAL INPUT=SHAPEFILE
                	
                        Here is my gdal-config output: (gdal-config
--formats)
                        gxf gtiff hfa aigrid aaigrid ceos ceos2 iso8211
xpm sdts raw dted mem jdem envisat elas fit vrt usgsdem l1b nitf bmp
pcidsk bsb jpeg2000 gif jpeg png grass
                	
                	
                        Does anybody know if it is possible to use JPEG
2000 instead of TIFF for this?
                	
                        I would like to use this format if possible, as
the file sizes are 2-3x smaller than with the .tif, and the quality
differences are unnoticeable.
                	
                        Any help/suggestions welcome.
                	
                	
                        Thank you,
                	
                        Brant Anderson
                	




-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://lists.osgeo.org/pipermail/mapserver-users/attachments/20040826/f1ea2810/attachment.htm>


More information about the MapServer-users mailing list