I have 1-meter resolution imagery thatI downloaded from this site:
http://viewer.nationalmap.gov/viewer/. Due to the size of the raster data,
the results are broken into chucks and the downloaded format is JP2. I used
the following scripts to convert the JP2 to GeoTiff:
for file in *.jp2;do
filename=`echo "$file" |cut -d'.' -f1`
gdal_translate -of GTiff $file $filename.tif
gdalwarp -t_srs EPSG:3857 $filename.tif $filename-warped.tif
gdal_translate -co "TILED=YES" -co "BLOCKXSIZE=512" -co "BLOCKYSIZE=512"
$filename-warped.tif $filename.tif
gdaladdo -r cubic $filename.tif 2 4 8 16 32 64 128
rm $filename-warped.tif
done
Below is the output of gdalinfo of one of the GeoTiff file:
Files: m_3711909_ne_11_1_20120628_20120904.tif
Size is 6562, 7948
Coordinate System is:
PROJCS["WGS 84 / Pseudo-Mercator",
GEOGCS["WGS 84",
DATUM["WGS_1984",
SPHEROID["WGS 84",6378137,298.257223563,
AUTHORITY["EPSG","7030"]],
AUTHORITY["EPSG","6326"]],
PRIMEM["Greenwich",0],
UNIT["degree",0.0174532925199433],
AUTHORITY["EPSG","4326"]],
PROJECTION["Mercator_1SP"],
PARAMETER["central_meridian",0],
PARAMETER["scale_factor",1],
PARAMETER["false_easting",0],
PARAMETER["false_northing",0],
UNIT["metre",1,
AUTHORITY["EPSG","9001"]],
EXTENSION["PROJ4","+proj=merc +a=6378137 +b=6378137 +lat_ts=0.0
+lon_0=0.0 +x_0=0.0 +y_0=0 +k=1.0 +units=m +nadgrids=@null +wktext
+no_defs"],
AUTHORITY["EPSG","3857"]]
Origin = (-13352064.088617129251361,4562412.907068637199700)
Pixel Size = (1.267235134041296,-1.267235134041296)
Metadata:
AREA_OR_POINT=Area
Image Structure Metadata:
INTERLEAVE=PIXEL
Corner Coordinates:
Upper Left (-13352064.089, 4562412.907) (119d56'37.08"W, 37d52'46.09"N)
Lower Left (-13352064.089, 4552340.922) (119d56'37.08"W, 37d48'28.87"N)
Upper Right (-13343748.492, 4562412.907) (119d52' 8.16"W, 37d52'46.09"N)
Lower Right (-13343748.492, 4552340.922) (119d52' 8.16"W, 37d48'28.87"N)
Center (-13347906.290, 4557376.915) (119d54'22.62"W, 37d50'37.51"N)
Band 1 Block=512x512 Type=Byte, ColorInterp=Red
Overviews: 3281x3974, 1641x1987, 821x994, 411x497, 206x249, 103x125, 52x63
Mask Flags: PER_DATASET ALPHA
Overviews of mask band: 3281x3974, 1641x1987, 821x994, 411x497, 206x249,
103x125, 52x63
Band 2 Block=512x512 Type=Byte, ColorInterp=Green
Overviews: 3281x3974, 1641x1987, 821x994, 411x497, 206x249, 103x125, 52x63
Mask Flags: PER_DATASET ALPHA
Overviews of mask band: 3281x3974, 1641x1987, 821x994, 411x497, 206x249,
103x125, 52x63
Band 3 Block=512x512 Type=Byte, ColorInterp=Blue
Overviews: 3281x3974, 1641x1987, 821x994, 411x497, 206x249, 103x125, 52x63
Mask Flags: PER_DATASET ALPHA
Overviews of mask band: 3281x3974, 1641x1987, 821x994, 411x497, 206x249,
103x125, 52x63
Band 4 Block=512x512 Type=Byte, ColorInterp=Alpha
Overviews: 3281x3974, 1641x1987, 821x994, 411x497, 206x249, 103x125, 52x63
Depending on the size of the selected area, the number of files can be in
the 50s or even 100s. I tested with a sample imagery that has 24 JP2 files
and create a store using ImageMosaic. The performance was not good. Can
someone tell me what is the best way to perpare this kind of imagery for
best performance?
Many thanks in advance,
Regards,
Tam
--
View this message in context:
http://osgeo-org.1560.x6.nabble.com/Prepare-GeoTiff-data-for-optimal-performance-tp5136963.html
Sent from the GeoServer - User mailing list archive at Nabble.com.
------------------------------------------------------------------------------
"Accelerate Dev Cycles with Automated Cross-Browser Testing - For FREE
Instantly run your Selenium tests across 300+ browser/OS combos. Get
unparalleled scalability from the best Selenium testing platform available.
Simple to use. Nothing to install. Get started now for free."
http://p.sf.net/sfu/SauceLabs
_______________________________________________
Geoserver-users mailing list
[email protected]
https://lists.sourceforge.net/lists/listinfo/geoserver-users
