Hi, I use gdal 1.9.0 (trunk), mapserver 6.0.0, ubuntu 11.04
I need to access several gdal subdatasets within gdal dataset as different tileindexes (need to build time-support on top of this). Situation: 1) I have GDAL WCS data source file, with temporal information in it 2) I created shapefile with 2 attributes: location AND time. location cotains the subdataset definition, which can be read by gdal ogrinfo on /tmp/WCS-a7145b1a3d54e8df55f019a136f55f84iYEWBn/20101014T030000_pph.wcs.tindex/20101014T030000_pph.wcs.tindex.shp (which is the generated time index file): ... OGRFeature(20101014T030000_pph.wcs.tindex):18 location (String) = WCS_SDS:time="2010-10-14T21:00:00.000Z",/tmp/WCS-a7145b1a3d54e8df55f019a136f55f84iYEWBn/20101014T030000_pph.wcs time (String) = 2010-10-14T21:00:00.000Z POLYGON ((409900 -11700,12900 -11700,12900 0,409900 0,409900 -11700)) ... gdalinfo on tindex location: gdalinfo 'WCS_SDS:time="2010-10-14T21:00:00.000Z",/tmp/WCS-a7145b1a3d54e8df55f019a136f55f84iYEWBn/20101014T030000_pph.wcs' Driver: WCS/OGC Web Coverage Service Files: none associated Size is 645, 585 Coordinate System is: PROJCS["GGRS87 / Greek Grid", GEOGCS["GGRS87", DATUM["Greek_Geodetic_Reference_System_1987", SPHEROID["GRS 1980",6378137,298.257222101, AUTHORITY["EPSG","7019"]], TOWGS84[-199.87,74.79,246.62,0,0,0,0], AUTHORITY["EPSG","6121"]], PRIMEM["Greenwich",0, AUTHORITY["EPSG","8901"]], UNIT["degree",0.0174532925199433, AUTHORITY["EPSG","9122"]], AUTHORITY["EPSG","4121"]], UNIT["metre",1, AUTHORITY["EPSG","9001"]], PROJECTION["Transverse_Mercator"], PARAMETER["latitude_of_origin",0], PARAMETER["central_meridian",24], PARAMETER["scale_factor",0.9996], PARAMETER["false_easting",500000], PARAMETER["false_northing",0], AUTHORITY["EPSG","2100"], AXIS["Easting",EAST], AXIS["Northing",NORTH]] Origin = (409900.000000000000000,4503200.000000000000000) Pixel Size = (20.000000000000000,-20.000000000000000) Metadata: TIME_POSITION=2010-10-14T21:00:00.000Z Corner Coordinates: Upper Left ( 409900.000, 4503200.000) ( 22d56' 2.07"E, 40d40'29.14"N) Lower Left ( 409900.000, 4491500.000) ( 22d56' 8.10"E, 40d34' 9.76"N) Upper Right ( 422800.000, 4503200.000) ( 23d 5'11.49"E, 40d40'33.86"N) Lower Right ( 422800.000, 4491500.000) ( 23d 5'16.66"E, 40d34'14.45"N) Center ( 416350.000, 4497350.000) ( 23d 0'39.58"E, 40d37'21.89"N) Band 1 Block=645x585 Type=Float64, ColorInterp=Undefined -> works perfectly BUT when I run GetMap request (WMS) on that layer, I got mapscript.MapServerError: msDrawMap(): Image handling error. Failed to draw layer named '20101014T030000_pph'. msDrawRasterLayerLow(): Unable to access file. Corrupt, empty or missing file '/tmp/WCS-a7145b1a3d54e8df55f019a136f55f84iYEWBn/20101014T030000_pph.wcs.tindex/WCS_SDS:time="2010-10-14T03:00:00.000Z",/tmp/WCS-a7145b1a3d54e8df55f019a136f55f84iYEWBn/20101014T030000_pph.wcs' for layer '20101014T030000_pph'. , referer: http://briseide.localhost/map/index2.php WMS Mapfile layer definition looks like LAYER METADATA "wms_timeitem" "time" "wms_timeextent" "2010-10-14T03:00:00.000Z,2010-10-14T04:00:00.000Z,2010-10-14T05:00:00.000Z,2010-10-14T06:00:00.000Z,2010-10-14T07:00:00.000Z,2010-10-14T08:00:00.000Z,2010-10-14T09:00:00.000Z,2010-10-14T10:00:00.000Z,2010-10-14T11:00:00.000Z,2010-10-14T12:00:00.000Z,2010-10-14T13:00:00.000Z,2010-10-14T14:00:00.000Z,2010-10-14T15:00:00.000Z,2010-10-14T16:00:00.000Z,2010-10-14T17:00:00.000Z,2010-10-14T18:00:00.000Z,2010-10-14T19:00:00.000Z,2010-10-14T20:00:00.000Z,2010-10-14T21:00:00.000Z" "wms_title" "20101014T030000_pph" END # METADATA NAME "20101014T030000_pph" PROJECTION "proj=tmerc" "lat_0=0" "lon_0=24" "k=0.9996" "x_0=500000" "y_0=0" "datum=GGRS87" "units=m" "no_defs" END # PROJECTION STATUS OFF TEMPLATE "foo" TILEINDEX "/tmp/WCS-a7145b1a3d54e8df55f019a136f55f84iYEWBn/20101014T030000_pph.wcs.tindex/20101014T030000_pph.wcs.tindex.shp" TILEITEM "location" TYPE RASTER UNITS METERS CLASS STYLE ANGLE 0 OFFSET 0 0 END # STYLE END # CLASS END # LAYER The part "/tmp/WCS-a7145b1a3d54e8df55f019a136f55f84iYEWBn/20101014T030000_pph.wcs.tindex/" of the requested "file name" should *NOT* be there. Questions: 1) What did I wrong with the configuration? 2) Is there other way, how to access temporal-enabled time data, available as GDAL subdataset? Thanks Jachym -- Jachym Cepicky e-mail: jachym.cepicky gmail com URL: http://les-ejk.cz GPG: http://les-ejk.cz/pgp/JachymCepicky.pgp _______________________________________________ mapserver-users mailing list mapserver-users@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/mapserver-users