Hi Even,
When I check gdalinfo for transh01.thf file, the result is,
-------------------------------------------------------
Driver: ADRG/ARC Digitized Raster Graphics
Files: transh01.thf
Size is 13696, 2944
Coordinate System is:
GEOGCS["WGS 84", DATUM["WGS_1984", SPHEROID["WGS
84",6378137,298.257223563]],
PRIMEM["Greenwich",0], UNIT["degree",0.0174532925199433],
AUTHORITY["EPSG","4326"]]
Origin = (-10.588233333333333,52.738852777777780)
Pixel Size = (0.005908613445378,-0.004478503184713)
Band number: 3
Metadata:
ADRG_SCA=5000000
ADRG_NAM=GNUR0102
Corner Coordinates:
Upper Left ( -10.5882333, 52.7388528) ( 10d35'17.64"W, 52d44'19.87"N)
Lower Left ( -10.5882333, 39.5541394) ( 10d35'17.64"W, 39d33'14.90"N)
Upper Right ( 70.3361364, 52.7388528) ( 70d20'10.09"E, 52d44'19.87"N)
Lower Right ( 70.3361364, 39.5541394) ( 70d20'10.09"E, 39d33'14.90"N)
Center ( 29.8739515, 46.1464961) ( 29d52'26.23"E, 46d 8'47.39"N)
Band 1 Block=128x128 Type=Byte, ColorInterp=Red
NoData Value=0
Band 2 Block=128x128 Type=Byte, ColorInterp=Green
NoData Value=0
Band 3 Block=128x128 Type=Byte, ColorInterp=Blue
NoData Value=0
------------------------------------------------------------------------------------
When I check gdalinfo for gnur0101.gen file, the result is,
-----------------------------------------------------------------------------------
Driver: ADRG/ARC Digitized Raster Graphics
Files: gnur0101.gen
Size is 13696, 2944
Coordinate System is:
GEOGCS["WGS 84", DATUM["WGS_1984", SPHEROID["WGS
84",6378137,298.257223563]],
PRIMEM["Greenwich",0], UNIT["degree",0.0174532925199433],
AUTHORITY["EPSG","4326"]]
Origin = (-10.588233333333333,52.738852777777780)
Pixel Size = (0.005908613445378,-0.004478503184713)
Band number: 3
Metadata:
ADRG_SCA=5000000
ADRG_NAM=GNUR0102
Corner Coordinates:
Upper Left ( -10.5882333, 52.7388528) ( 10d35'17.64"W, 52d44'19.87"N)
Lower Left ( -10.5882333, 39.5541394) ( 10d35'17.64"W, 39d33'14.90"N)
Upper Right ( 70.3361364, 52.7388528) ( 70d20'10.09"E, 52d44'19.87"N)
Lower Right ( 70.3361364, 39.5541394) ( 70d20'10.09"E, 39d33'14.90"N)
Center ( 29.8739515, 46.1464961) ( 29d52'26.23"E, 46d 8'47.39"N)
Band 1 Block=128x128 Type=Byte, ColorInterp=Red
NoData Value=0
Band 2 Block=128x128 Type=Byte, ColorInterp=Green
NoData Value=0
Band 3 Block=128x128 Type=Byte, ColorInterp=Blue
NoData Value=0
--------------------------------------------------------------------------------------
There are 4 subdatasets, but gdalinfo shows just the first one in metadata
as ADRG_NAM=GNUR0102 with size 13696, 2944
In ADRGDataset::GetGENListFromTHF function of adrgdataset.cpp file, I tried
to catch IMG files and add to fileNames. Then gdalinfo showed all the img
files in subdatasets metadata.
But this time GDALGetRasterCount gives 0 band. Because in ADRGDataset::Open
function, when count of fileNames is bigger than 1, all fileNames[i] are
being added as subdatasets to poDS and poDS returns without nBands are set.
As I understand that, when subdatasets exist, subdataset files should be
handled separately. But gdalinfo does not recognize the subdataset files
(xxx.img).
So maybe I should handle the adrg data, which has multiple subdatasets, by
my own code. But this is not a convenient way of reading adrg data.
Best regards...
Ayhan TEKGÜL
----- Original Message -----
From: "Even Rouault" <[email protected]>
To: "Ayhan TEKGÜL" <[email protected]>
Cc: <[email protected]>
Sent: Saturday, December 26, 2009 8:10 PM
Subject: Re: [gdal-dev] Display ADRG data with multiple subdatasets
Ayhan,
the driver has support for ADRG subdatasets (although I'm not 100%
positive this has been tested). Theoretically, if you do a gdalinfo on the
.THF file and it has several subdatasets, they should be listed. But
basically the list is just the list of related .GEN files, so you can
directly do gdalinfo on the .GEN files.
Best regards,
Even
Ayhan TEKGÜL a écrit :
Hi, all;
I am using gdal-1.5.0 to display ADRG data. But when ADRG data has
multiple subdatasets, I can display just one image.
I have tried to search the wiki and gdal-dev archives, but I cannot get
any similar problem with this issue or solution to this problem. (Maybe I
have missed.)
Is there anybody who knows to display adrg data with multiple
subdatasets?
Best regards...
Ayhan TEKGÜL
------------------------------------------------------------------------
_______________________________________________
gdal-dev mailing list
[email protected]
http://lists.osgeo.org/mailman/listinfo/gdal-dev
_______________________________________________
gdal-dev mailing list
[email protected]
http://lists.osgeo.org/mailman/listinfo/gdal-dev
