[Geoserver-users] MrSID to GeoTiff problem?

Hello everybody,

I'm struggeling with this problem for days now. I used SID files with
the GeoServer and the GDAL Plugin, but unfortunately geoserver needs a
lot of more time to display the sid files than a for instance a
geotiff. Therefor I want to used the GeoTiff files (I know they are
much bigger), but converting SID to geotiff with MrSID Geoviewer does
a wrong or bad geo convertion. Maybe some of u can help me with that.
The interestion thing is, geoserver takes the sid files as they are, I
just need to tell him its EPGS:4326 and it works.

So, my question is: How can I convert the SID Files to GeoTiff witht
the correct coordinates, so Geoserver can used it?

Original Metadata From SID:

                GEOTIFF_CHAR::
                               GTModelTypeGeoKey: ModelTypeProjected
                               GTRasterTypeGeoKey: RasterPixelIsArea
                               GeographicTypeGeoKey: GCS_WGS_84
                               ProjCoordTransGeoKey: CT_TransverseMercator
                               ProjLinearUnitsGeoKey: Linear_Meter
                               ProjectedCSTypeGeoKey: User-Defined
                               ProjectionGeoKey: User-Defined

                GEOTIFF_NUM::
                               1024:: GTModelTypeGeoKey: 1
                               1025:: GTRasterTypeGeoKey: 1
                               2048::GeographicTypeGeoKey: 4326
                               3072::ProjectedCSTypeGeoKey: 32767
                               3073::PCSCitationGeoKey:
Projection=TransverseMercator
                               3074::ProjectionGeoKey: 32767
                               3075:: ProjCoordTransGeoKey: 1
                               3076::ProjLinearUnitsGeoKey: 9001
                               3080:: ProjNatOriginLongGeoKey: 31
                               3081::ProjNatOriginLatGeoKey: 0
                               3082:: ProjFalseEastingGeoKey: 0
                               3083:: ProjFalseNorthingGeoKey: 0
                               3092:: ProjScaleAtNatOriginGeoKey: 1

                IMAGE::
                               BITS_PER_SAMPLE: 8
                               COLOR_SCHEME: 3
                               COMPRESSION_BLOCK_SIZE: 512
                               COMPRESSION_GAMMA: 2
                               COMPRESSION_NLEV: 7
                               COMPRESSION_VERSION[3]:
                                               2 0 0
                               COMPRESSION_WEIGHT: 2
                               CREATION_DATE: ThuOct2711:11:262005
                               DATA_TYPE: 0
                               DYNAMIC_RANGE_LEVEL: 127,5
                               DYNAMIC_RANGE_WINDOW: 256
                               HEIGHT: 7797
                               INPUT_FILE_SIZE: 55717908
                               INPUT_FORMAT: TIFFw/WorldFile
                               INPUT_NAME:
E:\Rectified_Images\Tif\WGS84\498_437_Phelindaba\498_437_01_0022_rect.tif
                               TARGET_COMPRESSION_RATIO: 10
                               WIDTH: 7145
                               XY_ORIGIN[2]:^
                                               -25479 -2710534
                               X_RESOLUTION: 0,75
                               Y_RESOLUTION: 0,75

GeoServer Information of SID File:
...
  <styles default = "raster" />
  <envelope crs = "GEOGCS['WGS 84',
  DATUM['World Geodetic System 1984',
    SPHEROID['WGS 84', 6378137.0, 298.257223563, AUTHORITY['EPSG','7030']],
    AUTHORITY['EPSG','6326']],
  PRIMEM['Greenwich', 0.0, AUTHORITY['EPSG','8901']],
  UNIT['degree', 0.017453292519943295],
  AXIS['Geodetic longitude', EAST],
  AXIS['Geodetic latitude', NORTH],
  AUTHORITY['EPSG','4326']]" srsName = "EPSG:4326" >
    <pos>30.748510284876335 -24.551435586381977</pos>
    <pos>30.801485617549876 -24.498563349392</pos>
  </envelope>
  <grid dimension = "2" >
    <low>0 0 </low>
    <high>7145 7797 </high>
    <axisName>GRAY</axisName>
    <geoTransform >
      <scaleX>0.75</scaleX>
      <scaleY>-0.75</scaleY>
      <shearX>0.0</shearX>
      <shearY>0.0</shearY>
      <translateX>-25479.75</translateX>
      <translateY>-2710533.25</translateY>
    </geoTransform>
  </grid>
...

SID to GeoTiff:
0,75
0
0
-0,75
-25479,00
-2710534,00

Thanks,
Sebastian

Ciao Sebastian,
quick question, when you configured your SID files, did you customized
the tiling settings? Did you also enabled multithreading read
operations?
Some MrSID files are encoded in stripes (line by line), hence unless
you ask geoserver to retile them on the fly it gets pretty slow at
serving them.

Simone.
-------------------------------------------------------
Eng. Simone Giannecchini
GeoSolutions S.A.S.
Owner - Software Engineer
Via Carignoni 51
55041 Camaiore (LU)
Italy

phone: +39 0584983027
fax: +39 0584983027
mob: +39 333 8128928

http://www.geo-solutions.it
http://simboss.blogspot.com/
http://www.linkedin.com/in/simonegiannecchini

-------------------------------------------------------

On Tue, Nov 18, 2008 at 3:40 PM, Sebastian Mueller <basti2nd@anonymised.com> wrote:

Hello everybody,

I'm struggeling with this problem for days now. I used SID files with
the GeoServer and the GDAL Plugin, but unfortunately geoserver needs a
lot of more time to display the sid files than a for instance a
geotiff. Therefor I want to used the GeoTiff files (I know they are
much bigger), but converting SID to geotiff with MrSID Geoviewer does
a wrong or bad geo convertion. Maybe some of u can help me with that.
The interestion thing is, geoserver takes the sid files as they are, I
just need to tell him its EPGS:4326 and it works.

So, my question is: How can I convert the SID Files to GeoTiff witht
the correct coordinates, so Geoserver can used it?

Original Metadata From SID:

               GEOTIFF_CHAR::
                              GTModelTypeGeoKey: ModelTypeProjected
                              GTRasterTypeGeoKey: RasterPixelIsArea
                              GeographicTypeGeoKey: GCS_WGS_84
                              ProjCoordTransGeoKey: CT_TransverseMercator
                              ProjLinearUnitsGeoKey: Linear_Meter
                              ProjectedCSTypeGeoKey: User-Defined
                              ProjectionGeoKey: User-Defined

               GEOTIFF_NUM::
                              1024:: GTModelTypeGeoKey: 1
                              1025:: GTRasterTypeGeoKey: 1
                              2048::GeographicTypeGeoKey: 4326
                              3072::ProjectedCSTypeGeoKey: 32767
                              3073::PCSCitationGeoKey:
Projection=TransverseMercator
                              3074::ProjectionGeoKey: 32767
                              3075:: ProjCoordTransGeoKey: 1
                              3076::ProjLinearUnitsGeoKey: 9001
                              3080:: ProjNatOriginLongGeoKey: 31
                              3081::ProjNatOriginLatGeoKey: 0
                              3082:: ProjFalseEastingGeoKey: 0
                              3083:: ProjFalseNorthingGeoKey: 0
                              3092:: ProjScaleAtNatOriginGeoKey: 1

               IMAGE::
                              BITS_PER_SAMPLE: 8
                              COLOR_SCHEME: 3
                              COMPRESSION_BLOCK_SIZE: 512
                              COMPRESSION_GAMMA: 2
                              COMPRESSION_NLEV: 7
                              COMPRESSION_VERSION[3]:
                                              2 0 0
                              COMPRESSION_WEIGHT: 2
                              CREATION_DATE: ThuOct2711:11:262005
                              DATA_TYPE: 0
                              DYNAMIC_RANGE_LEVEL: 127,5
                              DYNAMIC_RANGE_WINDOW: 256
                              HEIGHT: 7797
                              INPUT_FILE_SIZE: 55717908
                              INPUT_FORMAT: TIFFw/WorldFile
                              INPUT_NAME:
E:\Rectified_Images\Tif\WGS84\498_437_Phelindaba\498_437_01_0022_rect.tif
                              TARGET_COMPRESSION_RATIO: 10
                              WIDTH: 7145
                              XY_ORIGIN[2]:^
                                              -25479 -2710534
                              X_RESOLUTION: 0,75
                              Y_RESOLUTION: 0,75

GeoServer Information of SID File:
...
<styles default = "raster" />
<envelope crs = "GEOGCS['WGS 84',
DATUM['World Geodetic System 1984',
   SPHEROID['WGS 84', 6378137.0, 298.257223563, AUTHORITY['EPSG','7030']],
   AUTHORITY['EPSG','6326']],
PRIMEM['Greenwich', 0.0, AUTHORITY['EPSG','8901']],
UNIT['degree', 0.017453292519943295],
AXIS['Geodetic longitude', EAST],
AXIS['Geodetic latitude', NORTH],
AUTHORITY['EPSG','4326']]" srsName = "EPSG:4326" >
   <pos>30.748510284876335 -24.551435586381977</pos>
   <pos>30.801485617549876 -24.498563349392</pos>
</envelope>
<grid dimension = "2" >
   <low>0 0 </low>
   <high>7145 7797 </high>
   <axisName>GRAY</axisName>
   <geoTransform >
     <scaleX>0.75</scaleX>
     <scaleY>-0.75</scaleY>
     <shearX>0.0</shearX>
     <shearY>0.0</shearY>
     <translateX>-25479.75</translateX>
     <translateY>-2710533.25</translateY>
   </geoTransform>
</grid>
...

SID to GeoTiff:
0,75
0
0
-0,75
-25479,00
-2710534,00

Thanks,
Sebastian

-------------------------------------------------------------------------
This SF.Net email is sponsored by the Moblin Your Move Developer's challenge
Build the coolest Linux based applications with Moblin SDK & win great prizes
Grand prize is a trip for two to an Open Source event anywhere in the world
http://moblin-contest.org/redirect.php?banner_id=100&url=/
_______________________________________________
Geoserver-users mailing list
Geoserver-users@lists.sourceforge.net
https://lists.sourceforge.net/lists/listinfo/geoserver-users

Simone Giannecchini ha scritto:

Ciao Sebastian,
quick question, when you configured your SID files, did you customized
the tiling settings? Did you also enabled multithreading read
operations?
Some MrSID files are encoded in stripes (line by line), hence unless
you ask geoserver to retile them on the fly it gets pretty slow at
serving them.

Can we make those params the default? All the SID files I stumbled
upon so far need them to get decent performance.

Cheers
Andrea

--
Andrea Aime
OpenGeo - http://opengeo.org
Expert service straight from the developers.

I use the standard configuration of the geoserver, do i need to set
that in a special config file? Or just during setting up the layer? Or
do I have to set something special with those parameters:

USE_MULTITHREADING:
SUGGESTED_TILE_SIZE:
USE_JAI_IMAGEREAD:

Nevertheless, I would still like to know how to convert the SID Files
to GeoTiff in the correct way?

Thanks again and regards,
Sebastian

2008/11/18 Andrea Aime <aaime@anonymised.com>:

Simone Giannecchini ha scritto:

Ciao Sebastian,
quick question, when you configured your SID files, did you customized
the tiling settings? Did you also enabled multithreading read
operations?
Some MrSID files are encoded in stripes (line by line), hence unless
you ask geoserver to retile them on the fly it gets pretty slow at
serving them.

Can we make those params the default? All the SID files I stumbled
upon so far need them to get decent performance.

Cheers
Andrea

--
Andrea Aime
OpenGeo - http://opengeo.org
Expert service straight from the developers.

I'd suggest something like the following:

USE_MULTITHREADING: true
SUGGESTED_TILE_SIZE: 512,512
USE_JAI_IMAGEREAD: true

Let me Know how it goes with these params.

As far as converting is concerned, I'd use gdal, or you can do that
with geotools as well.
Look for gdal translate utility and make sure you give it the right options.
Simone.

-------------------------------------------------------
Eng. Simone Giannecchini
GeoSolutions S.A.S.
Owner - Software Engineer
Via Carignoni 51
55041 Camaiore (LU)
Italy

phone: +39 0584983027
fax: +39 0584983027
mob: +39 333 8128928

http://www.geo-solutions.it
http://simboss.blogspot.com/
http://www.linkedin.com/in/simonegiannecchini

-------------------------------------------------------

On Tue, Nov 18, 2008 at 4:31 PM, Sebastian Mueller <basti2nd@anonymised.com> wrote:

I use the standard configuration of the geoserver, do i need to set
that in a special config file? Or just during setting up the layer? Or
do I have to set something special with those parameters:

USE_MULTITHREADING:
SUGGESTED_TILE_SIZE:
USE_JAI_IMAGEREAD:

Nevertheless, I would still like to know how to convert the SID Files
to GeoTiff in the correct way?

Thanks again and regards,
Sebastian

2008/11/18 Andrea Aime <aaime@anonymised.com>:

Simone Giannecchini ha scritto:

Ciao Sebastian,
quick question, when you configured your SID files, did you customized
the tiling settings? Did you also enabled multithreading read
operations?
Some MrSID files are encoded in stripes (line by line), hence unless
you ask geoserver to retile them on the fly it gets pretty slow at
serving them.

Can we make those params the default? All the SID files I stumbled
upon so far need them to get decent performance.

Cheers
Andrea

--
Andrea Aime
OpenGeo - http://opengeo.org
Expert service straight from the developers.

It would surely be feasible, even though I am not a fan of automagic
configuration.
Andrea, could you opne up a jira for this?

Simone.
-------------------------------------------------------
Eng. Simone Giannecchini
GeoSolutions S.A.S.
Owner - Software Engineer
Via Carignoni 51
55041 Camaiore (LU)
Italy

phone: +39 0584983027
fax: +39 0584983027
mob: +39 333 8128928

http://www.geo-solutions.it
http://simboss.blogspot.com/
http://www.linkedin.com/in/simonegiannecchini

-------------------------------------------------------

On Tue, Nov 18, 2008 at 3:53 PM, Andrea Aime <aaime@anonymised.com> wrote:

Simone Giannecchini ha scritto:

Ciao Sebastian,
quick question, when you configured your SID files, did you customized
the tiling settings? Did you also enabled multithreading read
operations?
Some MrSID files are encoded in stripes (line by line), hence unless
you ask geoserver to retile them on the fly it gets pretty slow at
serving them.

Can we make those params the default? All the SID files I stumbled
upon so far need them to get decent performance.

Cheers
Andrea

--
Andrea Aime
OpenGeo - http://opengeo.org
Expert service straight from the developers.

Simone Giannecchini ha scritto:

It would surely be feasible, even though I am not a fan of automagic
configuration.
Andrea, could you opne up a jira for this?

There:
http://jira.codehaus.org/browse/GEOS-2395

Cheers
Andrea

--
Andrea Aime
OpenGeo - http://opengeo.org
Expert service straight from the developers.

thx
-------------------------------------------------------
Eng. Simone Giannecchini
GeoSolutions S.A.S.
Owner - Software Engineer
Via Carignoni 51
55041 Camaiore (LU)
Italy

phone: +39 0584983027
fax: +39 0584983027
mob: +39 333 8128928

http://www.geo-solutions.it
http://simboss.blogspot.com/
http://www.linkedin.com/in/simonegiannecchini

-------------------------------------------------------

On Tue, Nov 18, 2008 at 5:12 PM, Andrea Aime <aaime@anonymised.com> wrote:

Simone Giannecchini ha scritto:

It would surely be feasible, even though I am not a fan of automagic
configuration.
Andrea, could you opne up a jira for this?

There:
http://jira.codehaus.org/browse/GEOS-2395

Cheers
Andrea

--
Andrea Aime
OpenGeo - http://opengeo.org
Expert service straight from the developers.

Thanks for the info... i seems that it might got a little bit faster.
But still as performant as GeoTiff or JPG.

I'm still wondering how to convert the MrSID into JPG+WorldFile or
GeoTiff. Problem is still the same, coordinates are totaly wrong when
saved out of MrSID viewer. (Maybe coordinates aren't wrong, but just
dont fit into worldfile or geotiff format)

Thanks,
Sebastian

2008/11/18 Simone Giannecchini <simboss1@anonymised.com>:

thx
-------------------------------------------------------
Eng. Simone Giannecchini
GeoSolutions S.A.S.
Owner - Software Engineer
Via Carignoni 51
55041 Camaiore (LU)
Italy

phone: +39 0584983027
fax: +39 0584983027
mob: +39 333 8128928

http://www.geo-solutions.it
http://simboss.blogspot.com/
http://www.linkedin.com/in/simonegiannecchini

-------------------------------------------------------

On Tue, Nov 18, 2008 at 5:12 PM, Andrea Aime <aaime@anonymised.com> wrote:

Simone Giannecchini ha scritto:

It would surely be feasible, even though I am not a fan of automagic
configuration.
Andrea, could you opne up a jira for this?

There:
http://jira.codehaus.org/browse/GEOS-2395

Cheers
Andrea

--
Andrea Aime
OpenGeo - http://opengeo.org
Expert service straight from the developers.

Sebastian Mueller ha scritto:

Thanks for the info... i seems that it might got a little bit faster.
But still as performant as GeoTiff or JPG.

I'm still wondering how to convert the MrSID into JPG+WorldFile or
GeoTiff. Problem is still the same, coordinates are totaly wrong when
saved out of MrSID viewer. (Maybe coordinates aren't wrong, but just
dont fit into worldfile or geotiff format)

I'd suggest to use gdal_translate to perform the conversion. The issue
is finding a build that supports MrSid out of the box.
My FWTools windows install does support MrSid but I don't remember
whether I had to tweak it somehow or it just supports it out of the
box...

Cheers
Andrea

--
Andrea Aime
OpenGeo - http://opengeo.org
Expert service straight from the developers.

MrSID support in fwtools should be configured out of the box. FW is
subscribed to the list so if I am wrong he will probably correct me.

Simone.
-------------------------------------------------------
Eng. Simone Giannecchini
GeoSolutions S.A.S.
Owner - Software Engineer
Via Carignoni 51
55041 Camaiore (LU)
Italy

phone: +39 0584983027
fax: +39 0584983027
mob: +39 333 8128928

http://www.geo-solutions.it
http://simboss.blogspot.com/
http://www.linkedin.com/in/simonegiannecchini

-------------------------------------------------------

On Tue, Nov 25, 2008 at 10:31 AM, Andrea Aime <aaime@anonymised.com> wrote:

Sebastian Mueller ha scritto:

Thanks for the info... i seems that it might got a little bit faster.
But still as performant as GeoTiff or JPG.

I'm still wondering how to convert the MrSID into JPG+WorldFile or
GeoTiff. Problem is still the same, coordinates are totaly wrong when
saved out of MrSID viewer. (Maybe coordinates aren't wrong, but just
dont fit into worldfile or geotiff format)

I'd suggest to use gdal_translate to perform the conversion. The issue
is finding a build that supports MrSid out of the box.
My FWTools windows install does support MrSid but I don't remember
whether I had to tweak it somehow or it just supports it out of the
box...

Cheers
Andrea

--
Andrea Aime
OpenGeo - http://opengeo.org
Expert service straight from the developers.

Ciao Sebastian,
can you please tell me how you have set the parameters for mrsid ?

Simone.

-------------------------------------------------------
Eng. Simone Giannecchini
GeoSolutions S.A.S.
Owner - Software Engineer
Via Carignoni 51
55041 Camaiore (LU)
Italy

phone: +39 0584983027
fax: +39 0584983027
mob: +39 333 8128928

http://www.geo-solutions.it
http://simboss.blogspot.com/
http://www.linkedin.com/in/simonegiannecchini

-------------------------------------------------------

On Tue, Nov 25, 2008 at 10:22 AM, Sebastian Mueller <basti2nd@anonymised.com> wrote:

Thanks for the info... i seems that it might got a little bit faster.
But still as performant as GeoTiff or JPG.

I'm still wondering how to convert the MrSID into JPG+WorldFile or
GeoTiff. Problem is still the same, coordinates are totaly wrong when
saved out of MrSID viewer. (Maybe coordinates aren't wrong, but just
dont fit into worldfile or geotiff format)

Thanks,
Sebastian

2008/11/18 Simone Giannecchini <simboss1@anonymised.com>:

thx
-------------------------------------------------------
Eng. Simone Giannecchini
GeoSolutions S.A.S.
Owner - Software Engineer
Via Carignoni 51
55041 Camaiore (LU)
Italy

phone: +39 0584983027
fax: +39 0584983027
mob: +39 333 8128928

http://www.geo-solutions.it
http://simboss.blogspot.com/
http://www.linkedin.com/in/simonegiannecchini

-------------------------------------------------------

On Tue, Nov 18, 2008 at 5:12 PM, Andrea Aime <aaime@anonymised.com> wrote:

Simone Giannecchini ha scritto:

It would surely be feasible, even though I am not a fan of automagic
configuration.
Andrea, could you opne up a jira for this?

There:
http://jira.codehaus.org/browse/GEOS-2395

Cheers
Andrea

--
Andrea Aime
OpenGeo - http://opengeo.org
Expert service straight from the developers.

Ciao Simone,

I took the seetings u suggested.

USE_MULTITHREADING: true
SUGGESTED_TILE_SIZE: 512,512
USE_JAI_IMAGEREAD: true

Regards,
Sebastian

Ciao sebastian,
joìust out of curiosity, can you please copy and paste a gdalinfo for
one of your files.
In generale MrSID != more speed but I am curious about the size of
those files. It couold even be impossible to convert them to geotiff
due to size limitations.

SImone.
-------------------------------------------------------
Eng. Simone Giannecchini
GeoSolutions S.A.S.
Owner - Software Engineer
Via Carignoni 51
55041 Camaiore (LU)
Italy

phone: +39 0584983027
fax: +39 0584983027
mob: +39 333 8128928

http://www.geo-solutions.it
http://simboss.blogspot.com/
http://www.linkedin.com/in/simonegiannecchini

-------------------------------------------------------

On Tue, Nov 25, 2008 at 10:45 AM, Sebastian Mueller <basti2nd@anonymised.com> wrote:

Ciao Simone,

I took the seetings u suggested.

USE_MULTITHREADING: true
SUGGESTED_TILE_SIZE: 512,512
USE_JAI_IMAGEREAD: true

Regards,
Sebastian

Sure... here's the info:

Driver: MrSID/Multi-resolution Seamless Image Database (MrSID)
Files: 2430AC_25_498_396_05_5389_PAN_RECT.SID
Size is 7118, 7824
Coordinate System is:
PROJCS["IMAGINE GeoTIFF Support
Copyright 1991 - 2001 by ERDAS, Inc. All Rights Reserved
@(#)$RCSfile: egtf.c $ $Revision: 1.7 $ $Date: 2001/09/27 14:59:29EDT $
Projection Name = Transverse Mercator
Units = ",
    GEOGCS["WGS 84",
        DATUM["WGS_1984",
            SPHEROID["WGS 84",6378137,298.2572235630016,
                AUTHORITY["EPSG","7030"]],
            AUTHORITY["EPSG","6326"]],
        PRIMEM["Greenwich",0],
        UNIT["degree",0.0174532925199433],
        AUTHORITY["EPSG","4326"]],
    PROJECTION["Transverse_Mercator"],
    PARAMETER["latitude_of_origin",0],
    PARAMETER["central_meridian",31],
    PARAMETER["scale_factor",1],
    PARAMETER["false_easting",0],
    PARAMETER["false_northing",0],
    UNIT["metre",1,
        AUTHORITY["EPSG","9001"]]]
Origin = (-81236.375000000000000,-2705186.625000000000000)
Pixel Size = (0.750000000000000,-0.750000000000000)
Metadata:
  GEOTIFF_CHAR__GeographicTypeGeoKey=GCS_WGS_84
  GEOTIFF_CHAR__GTModelTypeGeoKey=ModelTypeProjected
  GEOTIFF_CHAR__GTRasterTypeGeoKey=RasterPixelIsArea
  GEOTIFF_CHAR__ProjCoordTransGeoKey=CT_TransverseMercator
  GEOTIFF_CHAR__ProjectedCSTypeGeoKey=User-Defined
  GEOTIFF_CHAR__ProjectionGeoKey=User-Defined
  GEOTIFF_CHAR__ProjLinearUnitsGeoKey=Linear_Meter
  GEOTIFF_NUM__1024__GTModelTypeGeoKey=1
  GEOTIFF_NUM__1025__GTRasterTypeGeoKey=1
  GEOTIFF_NUM__1026__GTCitationGeoKey=IMAGINE GeoTIFF Support
Copyright 1991 - 2001 by ERDAS, Inc. All Rights Reserved
@(#)$RCSfile: egtf.c $ $Revision: 1.7 $ $Date: 2001/09/27 14:59:29EDT $
Projection Name = Transverse Mercator
Units = meters
GeoTIFF Units = meters
  GEOTIFF_NUM__2048__GeographicTypeGeoKey=4326
  GEOTIFF_NUM__3072__ProjectedCSTypeGeoKey=32767
  GEOTIFF_NUM__3073__PCSCitationGeoKey=IMAGINE GeoTIFF Support
Copyright 1991 - 2001 by ERDAS, Inc. All Rights Reserved
@(#)$RCSfile: egtf.c $ $Revision: 1.7 $ $Date: 2001/09/27 14:59:29EDT $
Projection = Transverse Mercator
  GEOTIFF_NUM__3074__ProjectionGeoKey=32767
  GEOTIFF_NUM__3075__ProjCoordTransGeoKey=1
  GEOTIFF_NUM__3076__ProjLinearUnitsGeoKey=9001
  GEOTIFF_NUM__3080__ProjNatOriginLongGeoKey=31.000000
  GEOTIFF_NUM__3081__ProjNatOriginLatGeoKey=0.000000
  GEOTIFF_NUM__3082__ProjFalseEastingGeoKey=0.000000
  GEOTIFF_NUM__3083__ProjFalseNorthingGeoKey=0.000000
  GEOTIFF_NUM__3092__ProjScaleAtNatOriginGeoKey=1.000000
  IMAGE__COMPRESSION_BLOCK_SIZE=512
  IMAGE__COMPRESSION_GAMMA=2.000000
  IMAGE__COMPRESSION_NLEV=7
  IMAGE__COMPRESSION_VERSION=2,0,0
  IMAGE__COMPRESSION_WEIGHT=2.000000
  IMAGE__CREATION_DATE=Wed Apr 07 10:22:21 2004

  IMAGE__DYNAMIC_RANGE_LEVEL=127.500000
  IMAGE__DYNAMIC_RANGE_WINDOW=256.000000
  IMAGE__INPUT_FILE_SIZE=55699988.000000
  IMAGE__INPUT_FORMAT=TIFF w/ World File
  IMAGE__INPUT_NAME=D:\Rectified_Images\Tif\Wgs84\498_396_Ga-Malwane\498_396_05_5389_rect.tif
  IMAGE__TARGET_COMPRESSION_RATIO=10.000000
  IMAGE__WIDTH=7118
  IMAGE__HEIGHT=7824
  IMAGE__COLOR_SCHEME=3
  IMAGE__DATA_TYPE=0
  IMAGE__BITS_PER_SAMPLE=8
  IMAGE__XY_ORIGIN=-81236.000000,-2705187.000000
  IMAGE__X_RESOLUTION=0.750000
  IMAGE__Y_RESOLUTION=0.750000
  IMAGE__WKT=PROJCS["IMAGINE GeoTIFF Support
Copyright 1991 - 2001 by ERDAS, Inc. All Rights Reserved
@(#)$RCSfile: egtf.c $ $Revision: 1.7 $ $Date: 2001/09/27 14:59:29EDT $
Projection Name = Transverse Mercator
Units = ",GEOGCS["WGS 84",DATUM["WGS_1984",SPHEROID["WGS
84",6378137,298.2572235630016,AUTHORITY["EPSG","7030"]],AUTHORITY["EPSG","6326"]],PRIMEM["Greenwich",0],UNIT["degree",0.0174532925199433],AUTHORITY["EPSG","4326"]],PROJECTION["Transverse_Mercator"],PARAMETER["latitude_of_origin",0],PARAMETER["central_meridian",31],PARAMETER["scale_factor",1],PARAMETER["false_easting",0],PARAMETER["false_northing",0],UNIT["metre",1,AUTHORITY["EPSG","9001"]]]
Corner Coordinates:
Upper Left ( -81236.375,-2705186.625) ( 30d11'55.82"E, 24d26'54.16"S)
Lower Left ( -81236.375,-2711054.625) ( 30d11'54.61"E, 24d30'4.86"S)
Upper Right ( -75897.875,-2705186.625) ( 30d15'5.34"E, 24d26'55.13"S)
Lower Right ( -75897.875,-2711054.625) ( 30d15'4.21"E, 24d30'5.83"S)
Center ( -78567.125,-2708120.625) ( 30d13'30.00"E, 24d28'30.00"S)
Band 1 Block=1024x128 Type=Byte, ColorInterp=Gray
  Minimum=37.000, Maximum=211.000, Mean=143.704, StdDev=26.567
  Overviews: 3559x3912, 1780x1956, 890x978, 445x489, 223x245, 112x123, 56x62

Best regards,
Sebastian

Ciao Sebastian,
the file does not seem to be so big, so yeah, a geotiff with internal
overviews and tiling would probably be faster than mrsid with high
compression ratio.

Simone.

-------------------------------------------------------
Eng. Simone Giannecchini
GeoSolutions S.A.S.
Owner - Software Engineer
Via Carignoni 51
55041 Camaiore (LU)
Italy

phone: +39 0584983027
fax: +39 0584983027
mob: +39 333 8128928

http://www.geo-solutions.it
http://simboss.blogspot.com/
http://www.linkedin.com/in/simonegiannecchini

-------------------------------------------------------

On Tue, Nov 25, 2008 at 10:55 AM, Sebastian Mueller <basti2nd@anonymised.com> wrote:

Sure... here's the info:

Driver: MrSID/Multi-resolution Seamless Image Database (MrSID)
Files: 2430AC_25_498_396_05_5389_PAN_RECT.SID
Size is 7118, 7824
Coordinate System is:
PROJCS["IMAGINE GeoTIFF Support
Copyright 1991 - 2001 by ERDAS, Inc. All Rights Reserved
@(#)$RCSfile: egtf.c $ $Revision: 1.7 $ $Date: 2001/09/27 14:59:29EDT $
Projection Name = Transverse Mercator
Units = ",
   GEOGCS["WGS 84",
       DATUM["WGS_1984",
           SPHEROID["WGS 84",6378137,298.2572235630016,
               AUTHORITY["EPSG","7030"]],
           AUTHORITY["EPSG","6326"]],
       PRIMEM["Greenwich",0],
       UNIT["degree",0.0174532925199433],
       AUTHORITY["EPSG","4326"]],
   PROJECTION["Transverse_Mercator"],
   PARAMETER["latitude_of_origin",0],
   PARAMETER["central_meridian",31],
   PARAMETER["scale_factor",1],
   PARAMETER["false_easting",0],
   PARAMETER["false_northing",0],
   UNIT["metre",1,
       AUTHORITY["EPSG","9001"]]]
Origin = (-81236.375000000000000,-2705186.625000000000000)
Pixel Size = (0.750000000000000,-0.750000000000000)
Metadata:
GEOTIFF_CHAR__GeographicTypeGeoKey=GCS_WGS_84
GEOTIFF_CHAR__GTModelTypeGeoKey=ModelTypeProjected
GEOTIFF_CHAR__GTRasterTypeGeoKey=RasterPixelIsArea
GEOTIFF_CHAR__ProjCoordTransGeoKey=CT_TransverseMercator
GEOTIFF_CHAR__ProjectedCSTypeGeoKey=User-Defined
GEOTIFF_CHAR__ProjectionGeoKey=User-Defined
GEOTIFF_CHAR__ProjLinearUnitsGeoKey=Linear_Meter
GEOTIFF_NUM__1024__GTModelTypeGeoKey=1
GEOTIFF_NUM__1025__GTRasterTypeGeoKey=1
GEOTIFF_NUM__1026__GTCitationGeoKey=IMAGINE GeoTIFF Support
Copyright 1991 - 2001 by ERDAS, Inc. All Rights Reserved
@(#)$RCSfile: egtf.c $ $Revision: 1.7 $ $Date: 2001/09/27 14:59:29EDT $
Projection Name = Transverse Mercator
Units = meters
GeoTIFF Units = meters
GEOTIFF_NUM__2048__GeographicTypeGeoKey=4326
GEOTIFF_NUM__3072__ProjectedCSTypeGeoKey=32767
GEOTIFF_NUM__3073__PCSCitationGeoKey=IMAGINE GeoTIFF Support
Copyright 1991 - 2001 by ERDAS, Inc. All Rights Reserved
@(#)$RCSfile: egtf.c $ $Revision: 1.7 $ $Date: 2001/09/27 14:59:29EDT $
Projection = Transverse Mercator
GEOTIFF_NUM__3074__ProjectionGeoKey=32767
GEOTIFF_NUM__3075__ProjCoordTransGeoKey=1
GEOTIFF_NUM__3076__ProjLinearUnitsGeoKey=9001
GEOTIFF_NUM__3080__ProjNatOriginLongGeoKey=31.000000
GEOTIFF_NUM__3081__ProjNatOriginLatGeoKey=0.000000
GEOTIFF_NUM__3082__ProjFalseEastingGeoKey=0.000000
GEOTIFF_NUM__3083__ProjFalseNorthingGeoKey=0.000000
GEOTIFF_NUM__3092__ProjScaleAtNatOriginGeoKey=1.000000
IMAGE__COMPRESSION_BLOCK_SIZE=512
IMAGE__COMPRESSION_GAMMA=2.000000
IMAGE__COMPRESSION_NLEV=7
IMAGE__COMPRESSION_VERSION=2,0,0
IMAGE__COMPRESSION_WEIGHT=2.000000
IMAGE__CREATION_DATE=Wed Apr 07 10:22:21 2004

IMAGE__DYNAMIC_RANGE_LEVEL=127.500000
IMAGE__DYNAMIC_RANGE_WINDOW=256.000000
IMAGE__INPUT_FILE_SIZE=55699988.000000
IMAGE__INPUT_FORMAT=TIFF w/ World File
IMAGE__INPUT_NAME=D:\Rectified_Images\Tif\Wgs84\498_396_Ga-Malwane\498_396_05_5389_rect.tif
IMAGE__TARGET_COMPRESSION_RATIO=10.000000
IMAGE__WIDTH=7118
IMAGE__HEIGHT=7824
IMAGE__COLOR_SCHEME=3
IMAGE__DATA_TYPE=0
IMAGE__BITS_PER_SAMPLE=8
IMAGE__XY_ORIGIN=-81236.000000,-2705187.000000
IMAGE__X_RESOLUTION=0.750000
IMAGE__Y_RESOLUTION=0.750000
IMAGE__WKT=PROJCS["IMAGINE GeoTIFF Support
Copyright 1991 - 2001 by ERDAS, Inc. All Rights Reserved
@(#)$RCSfile: egtf.c $ $Revision: 1.7 $ $Date: 2001/09/27 14:59:29EDT $
Projection Name = Transverse Mercator
Units = ",GEOGCS["WGS 84",DATUM["WGS_1984",SPHEROID["WGS
84",6378137,298.2572235630016,AUTHORITY["EPSG","7030"]],AUTHORITY["EPSG","6326"]],PRIMEM["Greenwich",0],UNIT["degree",0.0174532925199433],AUTHORITY["EPSG","4326"]],PROJECTION["Transverse_Mercator"],PARAMETER["latitude_of_origin",0],PARAMETER["central_meridian",31],PARAMETER["scale_factor",1],PARAMETER["false_easting",0],PARAMETER["false_northing",0],UNIT["metre",1,AUTHORITY["EPSG","9001"]]]
Corner Coordinates:
Upper Left ( -81236.375,-2705186.625) ( 30d11'55.82"E, 24d26'54.16"S)
Lower Left ( -81236.375,-2711054.625) ( 30d11'54.61"E, 24d30'4.86"S)
Upper Right ( -75897.875,-2705186.625) ( 30d15'5.34"E, 24d26'55.13"S)
Lower Right ( -75897.875,-2711054.625) ( 30d15'4.21"E, 24d30'5.83"S)
Center ( -78567.125,-2708120.625) ( 30d13'30.00"E, 24d28'30.00"S)
Band 1 Block=1024x128 Type=Byte, ColorInterp=Gray
Minimum=37.000, Maximum=211.000, Mean=143.704, StdDev=26.567
Overviews: 3559x3912, 1780x1956, 890x978, 445x489, 223x245, 112x123, 56x62

Best regards,
Sebastian

-------------------------------------------------------------------------
This SF.Net email is sponsored by the Moblin Your Move Developer's challenge
Build the coolest Linux based applications with Moblin SDK & win great prizes
Grand prize is a trip for two to an Open Source event anywhere in the world
http://moblin-contest.org/redirect.php?banner_id=100&url=/
_______________________________________________
Geoserver-users mailing list
Geoserver-users@lists.sourceforge.net
https://lists.sourceforge.net/lists/listinfo/geoserver-users

Ciao Simone,

Thanks for that info, so at least I know it could be faster... but now
is the question, who to do achive that?

Regards and many thanks,
Sebastian

- download fwtools
- gdal_translate -co "TILED=YES" -co "BLOCKXSIZE=512" -co
"BLOCKYSIZE=512" input.sid output.tif

-------------------------------------------------------
Eng. Simone Giannecchini
GeoSolutions S.A.S.
Owner - Software Engineer
Via Carignoni 51
55041 Camaiore (LU)
Italy

phone: +39 0584983027
fax: +39 0584983027
mob: +39 333 8128928

http://www.geo-solutions.it
http://simboss.blogspot.com/
http://www.linkedin.com/in/simonegiannecchini

-------------------------------------------------------

On Tue, Nov 25, 2008 at 12:23 PM, Sebastian Mueller <basti2nd@anonymised.com> wrote:

Ciao Simone,

Thanks for that info, so at least I know it could be faster... but now
is the question, who to do achive that?

Regards and many thanks,
Sebastian

Thanks for the help and all the info! Thank u very much.

Did u or someone else played around with other tilesizes? Like 256 or
1024? Does is "dramatically" change the performance, or is 512 "the
best" tilesize?

Regards,
Sebastian

2008/11/25 Simone Giannecchini <simone.giannecchini@anonymised.com>:

- download fwtools
- gdal_translate -co "TILED=YES" -co "BLOCKXSIZE=512" -co
"BLOCKYSIZE=512" input.sid output.tif

-------------------------------------------------------
Eng. Simone Giannecchini
GeoSolutions S.A.S.
Owner - Software Engineer
Via Carignoni 51
55041 Camaiore (LU)
Italy

phone: +39 0584983027
fax: +39 0584983027
mob: +39 333 8128928

http://www.geo-solutions.it
http://simboss.blogspot.com/
http://www.linkedin.com/in/simonegiannecchini

-------------------------------------------------------

On Tue, Nov 25, 2008 at 12:23 PM, Sebastian Mueller <basti2nd@anonymised.com> wrote:

Ciao Simone,

Thanks for that info, so at least I know it could be faster... but now
is the question, who to do achive that?

Regards and many thanks,
Sebastian

512 is more or less a de-facto standard.
The correc tile size depdends pretty much on what you are going to do
with your data, hence probably you could capture some real use cases
into some test cases and then run them towards the same data but
converted with different tile size.

Simone.
-------------------------------------------------------
Eng. Simone Giannecchini
GeoSolutions S.A.S.
Owner - Software Engineer
Via Carignoni 51
55041 Camaiore (LU)
Italy

phone: +39 0584983027
fax: +39 0584983027
mob: +39 333 8128928

http://www.geo-solutions.it
http://simboss.blogspot.com/
http://www.linkedin.com/in/simonegiannecchini

-------------------------------------------------------

On Tue, Nov 25, 2008 at 1:29 PM, Sebastian Mueller <basti2nd@anonymised.com> wrote:

Thanks for the help and all the info! Thank u very much.

Did u or someone else played around with other tilesizes? Like 256 or
1024? Does is "dramatically" change the performance, or is 512 "the
best" tilesize?

Regards,
Sebastian

2008/11/25 Simone Giannecchini <simone.giannecchini@anonymised.com>:

- download fwtools
- gdal_translate -co "TILED=YES" -co "BLOCKXSIZE=512" -co
"BLOCKYSIZE=512" input.sid output.tif

-------------------------------------------------------
Eng. Simone Giannecchini
GeoSolutions S.A.S.
Owner - Software Engineer
Via Carignoni 51
55041 Camaiore (LU)
Italy

phone: +39 0584983027
fax: +39 0584983027
mob: +39 333 8128928

http://www.geo-solutions.it
http://simboss.blogspot.com/
http://www.linkedin.com/in/simonegiannecchini

-------------------------------------------------------

On Tue, Nov 25, 2008 at 12:23 PM, Sebastian Mueller <basti2nd@anonymised.com> wrote:

Ciao Simone,

Thanks for that info, so at least I know it could be faster... but now
is the question, who to do achive that?

Regards and many thanks,
Sebastian