Hi Devs, I have a geotiff file which contains several datasets, with this I just need to perform one gdalwarp operation on this geotiff file to apply transformation on 6 datasets inside it and the operation produces just one geotiff output file.
If I translate this geotiff file into netcdf file, now I need to perform 6
operations to apply transformation on each of the 6 datasets inside it and
produces 6 separate netcdf files which is quite cumbersome.
I am wondering if it is possible to modify the current netCDF driver so that it
behaves like geotiff driver when operating on multiple bands (datasets) netcdf
file.
Any opinion on this is really appreciated.
This might be a potential issue though:
When the several datasets have the same data type, it is not an issue but
when the several datasets have different data type (some datasets are float
while some other metadata are integer), then it might be an issue as when we do
gdalwarp operation selecting one method of interpolation, that chosen
interpolation method will be applied across all data types for all datasets in
that single file.
=============================
gdalinfo on the geotiff file: and netcdf file:
=============================
Driver: GTiff/GeoTIFF
Files: multipledataset.tif
Size is 4, 3
Coordinate System is:
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"]]
Origin = (141.000867233935992,-34.060737220013806)
Pixel Size = (0.107999999999998,-0.107999999999998)
Metadata:
AREA_OR_POINT=Area
Image Structure Metadata:
INTERLEAVE=PIXEL
Corner Coordinates:
Upper Left ( 141.0008672, -34.0607372) (141d 0' 3.12"E, 34d 3'38.65"S)
Lower Left ( 141.0008672, -34.3847372) (141d 0' 3.12"E, 34d23' 5.05"S)
Upper Right ( 141.4328672, -34.0607372) (141d25'58.32"E, 34d 3'38.65"S)
Lower Right ( 141.4328672, -34.3847372) (141d25'58.32"E, 34d23' 5.05"S)
Center ( 141.2168672, -34.2227372) (141d13' 0.72"E, 34d13'21.85"S)
Band 1 Block=4x3 Type=Float64, ColorInterp=Gray
Band 2 Block=4x3 Type=Float64, ColorInterp=Undefined
Band 3 Block=4x3 Type=Float64, ColorInterp=Undefined
Band 4 Block=4x3 Type=Float64, ColorInterp=Undefined
Band 5 Block=4x3 Type=Float64, ColorInterp=Undefined
Band 6 Block=4x3 Type=Float64, ColorInterp=Undefined
=============================
gdalinfo on the netcdf file:
=============================
Driver: netCDF/Network Common Data Format
Files: multipledataset.nc
Size is 512, 512
Coordinate System is `'
Metadata:
NC_GLOBAL#Conventions=CF-1.5
NC_GLOBAL#GDAL=GDAL 1.9.0, released 2011/12/29
NC_GLOBAL#history=Thu Mar 08 14:03:01 2012: GDAL Create( multipledataset.nc,
... )
Subdatasets:
SUBDATASET_1_NAME=NETCDF:"multipledataset.nc":Band1
SUBDATASET_1_DESC=[3x4] Band1 (64-bit floating-point)
SUBDATASET_2_NAME=NETCDF:"multipledataset.nc":Band2
SUBDATASET_2_DESC=[3x4] Band2 (64-bit floating-point)
SUBDATASET_3_NAME=NETCDF:"multipledataset.nc":Band3
SUBDATASET_3_DESC=[3x4] Band3 (64-bit floating-point)
SUBDATASET_4_NAME=NETCDF:"multipledataset.nc":Band4
SUBDATASET_4_DESC=[3x4] Band4 (64-bit floating-point)
SUBDATASET_5_NAME=NETCDF:"multipledataset.nc":Band5
SUBDATASET_5_DESC=[3x4] Band5 (64-bit floating-point)
SUBDATASET_6_NAME=NETCDF:"multipledataset.nc":Band6
SUBDATASET_6_DESC=[3x4] Band6 (64-bit floating-point)
Corner Coordinates:
Upper Left ( 0.0, 0.0)
Lower Left ( 0.0, 512.0)
Upper Right ( 512.0, 0.0)
Lower Right ( 512.0, 512.0)
Center ( 256.0, 256.0)
Regards,
Hendy<<attachment: multipledataset.tif>>
multipledataset.nc
Description: Cdf file
_______________________________________________ gdal-dev mailing list [email protected] http://lists.osgeo.org/mailman/listinfo/gdal-dev
