[gdal-dev] Display ADRG data with multiple subdatasets

Sun Dec 27 09:39:04 EST 2009

Ayhan,

this obviously confirms that ADRG subdataset was never tested before.  
I'd be willing to improve the support for them in GDAL itself rather 
than forcing each user to use their little own hacks.

The ADRG driver cannot currently open directly the .IMG files, but I 
suspect this would be propably the right thing to do for subdataset 
access. This might be probably something similar to what is suggested 
for the ASRP/USRP driver in http://trac.osgeo.org/gdal/ticket/3194

If I could have access to your data, it will be of course the most 
convenient solution. Otherwise could you list all the filenames of your 
ADRG product (from what you say I suspect there's only one GEN file but 
several IMG files ?). What would also greatly help is that you provide 
the output of the '8211dump' utility on both the .THF and .GEN files. 
The source for this utility is in the frmts/iso8211 subdirectory of GDAL 
source distribution. If you cd into this directory and do 'make 
8211dump' (on Linux), it should get compiled (but you might need to try 
this on gdal-trunk, as I'm not sure the GNUmakefile in the 1.5 branch 
will compile it without a few changes).

Even

Ayhan TEKGÜL a écrit :
> 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" 
> <even.rouault at mines-paris.org>
> To: "Ayhan TEKGÜL" <ayhan.tekgul at netcad.com.tr>
> Cc: <gdal-dev at lists.osgeo.org>
> 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
>>> gdal-dev at lists.osgeo.org
>>> http://lists.osgeo.org/mailman/listinfo/gdal-dev
>>
>>
>
>
>
>