Ciao Jonathan, please, read below... Regards, Simone Giannecchini == Our support, Your Success! Visit http://opensdi.geo-solutions.it for more information. ==
Ing. Simone Giannecchini @simogeo Founder/Director GeoSolutions S.A.S. Via Poggio alle Viti 1187 55054 Massarosa (LU) Italy phone: +39 0584 962313 fax: +39 0584 1660272 mob: +39 333 8128928 http://www.geo-solutions.it http://twitter.com/geosolutions_it ------------------------------------------------------- On Tue, Nov 12, 2013 at 10:20 AM, Jonathan Moules <jonathanmou...@warwickshire.gov.uk> wrote: > Hi Both, > Thanks for the quick and helpful replies. I'll do some experimentation with > these things. And there I was thinking I'd figured out GDAL. :-) You might > want to put some of this stuff in your "GeoServer on Steroids" presentation. > actually that is a good idea, we have some more things coming down the pipe > I'd thought that overviews were compressed - useful to know they're not; > that'll be a saving of a couple of GB. > that's true on very large imagery > I could have sworn someone told me JP2 couldn't do four bands. Hmmm. For now > I'll stick with GeoTIFFs - I don't have any 4 band JP2's and our vendor > wasn't able to create them (nor can FME), so I can't provide a sample I'm > afraid. > > > Will the PHOTOMETRIC= YCBCR work with a four band RGBI? > It should not. > > The GDALINFO for a source Geotiff is: > >> Driver: GTiff/GeoTIFF >> Files: SP4044.tif >> SP4044.tfw >> Size is 8000, 8000 >> Coordinate System is: >> PROJCS["OSGB 1936 / British National Grid", >> GEOGCS["OSGB 1936", >> DATUM["OSGB_1936", >> SPHEROID["Airy 1830",6377563.396,299.3249753150316, >> AUTHORITY["EPSG","7001"]], >> AUTHORITY["EPSG","6277"]], >> PRIMEM["Greenwich",0], >> UNIT["degree",0.0174532925199433], >> AUTHORITY["EPSG","4277"]], >> PROJECTION["Transverse_Mercator"], >> PARAMETER["latitude_of_origin",49], >> PARAMETER["central_meridian",-2], >> PARAMETER["scale_factor",0.9996012717], >> PARAMETER["false_easting",400000], >> PARAMETER["false_northing",-100000], >> UNIT["metre",1, >> AUTHORITY["EPSG","9001"]], >> AUTHORITY["EPSG","27700"]] >> Origin = (440000.000000000000000,245000.000000000000000) >> Pixel Size = (0.125000000000000,-0.125000000000000) >> Metadata: >> TIFFTAG_XRESOLUTION=72 >> TIFFTAG_YRESOLUTION=72 >> TIFFTAG_RESOLUTIONUNIT=2 (pixels/inch) >> AREA_OR_POINT=Area >> Image Structure Metadata: >> COMPRESSION=LZW >> INTERLEAVE=PIXEL >> Corner Coordinates: >> Upper Left ( 440000.000, 245000.000) ( 1d24'57.45"W, 52d 6'5.31"N) >> Lower Left ( 440000.000, 244000.000) ( 1d24'57.87"W, 52d 5'32.94"N) >> Upper Right ( 441000.000, 245000.000) ( 1d24'4.89"W, 52d 6'5.04"N) >> Lower Right ( 441000.000, 244000.000) ( 1d24'5.32"W, 52d 5'32.68"N) >> Center ( 440500.000, 244500.000) ( 1d24'31.39"W, 52d 5'48.99"N) >> Band 1 Block=8000x32 Type=Byte, ColorInterp=Red >> Band 2 Block=8000x32 Type=Byte, ColorInterp=Green >> Band 3 Block=8000x32 Type=Byte, ColorInterp=Blue >> Band 4 Block=8000x32 Type=Byte, ColorInterp=Undefined > > > Thanks again, > Jonathan > > On 12 November 2013 08:50, Simone Giannecchini > <simone.giannecch...@geo-solutions.it> wrote: >> >> Ciao Jonathan, >> please read my comments inline below... >> >> Regards, >> Simone Giannecchini >> == >> Our support, Your Success! Visit http://opensdi.geo-solutions.it for >> more information. >> == >> >> Ing. Simone Giannecchini >> @simogeo >> Founder/Director >> >> GeoSolutions S.A.S. >> Via Poggio alle Viti 1187 >> 55054 Massarosa (LU) >> Italy >> phone: +39 0584 962313 >> fax: +39 0584 1660272 >> mob: +39 333 8128928 >> >> http://www.geo-solutions.it >> http://twitter.com/geosolutions_it >> >> ------------------------------------------------------- >> >> >> On Tue, Nov 12, 2013 at 9:27 AM, Jonathan Moules >> <jonathanmou...@warwickshire.gov.uk> wrote: >> > Hi List, >> > We've recently taken deliver of some aerial photography. >> > >> > As lossless GeoTIFFs, the RedGreenBlue data is 278GB >> > As a compressed JP2 file the data is 26GB - the JP2 is almost as good as >> > the >> > GeoTIFFs and comes with pyramids built in. >> >> Unless you have an ECW or Kakadu license JP2 reading won't be as fast >> as GeoTiff, I mean, nowhere near to it :) >> >> > >> > But for speed purposes with GeoServer I've created an imageMosaic of 13 >> > Jpeg >> > compressed GeoTIFFs with overviews/tiles etc. This comes out at a >> > massive >> > 106GB and the quality is worse than the JP2 file despite being four >> > times >> > the size. >> > >> > Can anyone advise how I can improve these files to make them smaller >> > (ideally with better quality too)? I created each of the GeoTIFFs with >> > the >> > following two commands >> > >> > Merge, tile, and compress all at once. >> >> >> >> gdal_merge -q -o file_name.tif -of GTiff -co TILED=YES -co BIGTIFF=YES >> >> -co >> >> COMPRESS=JPEG -co JPEG_QUALITY=50 -co BLOCKXSIZE=512 -co BLOCKYSIZE=512 >> >> --optfile tiff_list.txt >> >> I believe quality is too low. This makes the files smaller but also >> makes the jpeg artifacts quite visibile. I would go with the default >> which is 75%. >> Moreover, you should use this switch PHOTOMETRIC= YCBCR to improve >> JPEG compression. >> >> > >> > >> > Pyramids >> >> >> >> gdaladdo -r average %THIS_DIR%.tif 2 4 8 16 32 64 128 256 >> >> This is correct but could be done better. Since the base level is >> compressed I would compress overviews as well. >> gdaladdo -r cubic -config COMPRESS_OVERVIEW JPEG --config >> PHOTOMETRIC_OVERVIEW YCBCR %THIS_DIR%.tif 2 4 8 16 32 64 128 256 >> >> This should provide more quality but also relatively small size. >> >> > >> > >> > >> > Using JP2 with GDAL is tempting (though supposedly much slower) - but >> > we're >> > likely to be putting up four band (RGB + Infrared) data, which I think >> > we >> > can use GeoServer to visualise (still to play with). JP2 can't handle >> > four >> > bands. >> >> Ehm, I am playing with 4 bands JP2 just these days, hence I am not so >> sure this is true (which means I believe it is not :) ) >> >> I am doing some work to improve performance + quality of 16 bits 4 >> bands imagery, if you move on along the JP2 lines I could make good >> use of a sample. >> >> > >> > Am I doing something wrong with the GeoTIFF creation? >> > >> > Thoughts welcome. >> > Thanks, >> > Jonathan >> > >> > This transmission is intended for the named addressee(s) only and may >> > contain sensitive or protectively marked material up to RESTRICTED and >> > should be handled accordingly. Unless you are the named addressee (or >> > authorised to receive it for the addressee) you may not copy or use it, >> > or >> > disclose it to anyone else. If you have received this transmission in >> > error >> > please notify the sender immediately. All email traffic sent to or from >> > us, >> > including without limitation all GCSX traffic, may be subject to >> > recording >> > and/or monitoring in accordance with relevant legislation. >> > >> > ------------------------------------------------------------------------------ >> > November Webinars for C, C++, Fortran Developers >> > Accelerate application performance with scalable programming models. >> > Explore >> > techniques for threading, error checking, porting, and tuning. Get the >> > most >> > from the latest Intel processors and coprocessors. See abstracts and >> > register >> > >> > http://pubads.g.doubleclick.net/gampad/clk?id=60136231&iu=/4140/ostg.clktrk >> > _______________________________________________ >> > Geoserver-users mailing list >> > Geoserver-users@lists.sourceforge.net >> > https://lists.sourceforge.net/lists/listinfo/geoserver-users >> > > > > > This transmission is intended for the named addressee(s) only and may > contain sensitive or protectively marked material up to RESTRICTED and > should be handled accordingly. Unless you are the named addressee (or > authorised to receive it for the addressee) you may not copy or use it, or > disclose it to anyone else. If you have received this transmission in error > please notify the sender immediately. All email traffic sent to or from us, > including without limitation all GCSX traffic, may be subject to recording > and/or monitoring in accordance with relevant legislation. ------------------------------------------------------------------------------ November Webinars for C, C++, Fortran Developers Accelerate application performance with scalable programming models. Explore techniques for threading, error checking, porting, and tuning. Get the most from the latest Intel processors and coprocessors. See abstracts and register http://pubads.g.doubleclick.net/gampad/clk?id=60136231&iu=/4140/ostg.clktrk _______________________________________________ Geoserver-users mailing list Geoserver-users@lists.sourceforge.net https://lists.sourceforge.net/lists/listinfo/geoserver-users
