[Gdal-dev] Gdal_translate Error

[David W. Lowther]

List,

After some more investigation I have found that the images that do not
work correctly have been compressed with Encoder Version 3.0.0.6:

Input filename: ortho_e1-1_ok109 [TIFF]
Output filename:
F:\1300511\DOQQ\OKLAHOMA_OK_40109_88_2003\ortho_e1-1_ok109.sid
Input image size: 10460.2 MB (10968298530 bytes)
Output file size: 216.9 MB (227481179 bytes)
Target Encoding Ratio: 50.00:1
Actual Encoding Ratio: 48.22:1
File type: GeoTIFF
Georeferencing: GeoTIFF
Encoder Version: 3.0.0.6
Output Format: MrSID Generation 2
Start Time: Thu Sep 25 15:21:05 2003 
Completion Time: Thu Sep 25 16:31:14 2003

I've got a 1.1.9 build of GDAL with the following MrSid options:

MRSID_DIR =	e:\build\win6_dsdk
MRSID_LIB =	/LIBPATH:$(MRSID_DIR)\bin\lib_md \
		MrSIDDecode_md.lib lt_xTrans_md.lib advapi32.lib
user32.lib

Do I need to recompile with a newer version to make this work?

Thanks,
Dave

-----Original Message-----
From: gdal-dev-bounces at xserve.flids.com
[mailto:gdal-dev-bounces at xserve.flids.com] On Behalf Of David W. Lowther
Sent: Wednesday, August 18, 2004 9:34 AM
To: gdal-dev at xserve.flids.com
Subject: RE: [Gdal-dev] Gdal_translate Error

Frank,

Here's the complete command and results:

E:\Dave\GDAL\apps>gdal_translate -projwin 635544.7539 3930207.3199
648110.83864 3915163.12478 h:\deqdata\oklahoma\Oklahoma.sid
h:\deqdata\oklahoma\x1_35097D4.tif -co "TFW=YES"
Input file size is 63929, 57190
Computed -srcwin 22976 35037 12566 15044 from projected window.
0.ERROR 1: Internalize failed
ERROR 1: GetBlockRef failed at X block offset 0, Y block offset 0


Here's the SDW:

             1.000000000000000
             0.000000000000000
             0.000000000000000
            -1.000000000000000
        612568.500000000000000
       3965244.500000000000000

Can you tell if my -projwin is outside? I think it's OK.

BTW, I'd be happy to send someone a copy of this image if it would help
in debugging this situation.

Thanks,

Dave

-----Original Message-----
From: Frank Warmerdam [mailto:warmerdam at pobox.com] 
Sent: Monday, August 16, 2004 10:33 AM
To: David W. Lowther
Cc: gdal-dev at xserve.flids.com
Subject: Re: [Gdal-dev] Gdal_translate Error

David W. Lowther wrote:
> Hi List,
> 
>  
> 
> We're processing some large SID images back to TIF to project. Some 
> images work just fine, others (most) have issues like this:
> 
>  
> 
> gdal_translate -projwin 635544.7539 3930207.3199 648110.83864 
> 3915163.12478 h:\deqdata\oklahoma\Oklahoma.sid 
> h:\deqdata\oklahoma\x1_35097D4.tif -co "TFW=YES"
> 
> Input file size is 63929, 57190
> 
> Computed -srcwin 22976 35037 12566 15044 from projected window.
> 
> 0.ERROR 1: Internalize failed
> 
> ERROR 1: GetBlockRef failed at X block offset 0, Y block offset

David,

My guess would be that your -projwin definition touches an edge and for
some
reason related to rounding or whatever the requested window ends up
being a
pixel or more off the image.  This might result in an error similar to
what
you see.  Can you verify if that might be it?

Altneratively, is it possible these are color images?  If so, the
resulting
TIFF would be larger than 4GB and would run into some sort of problem.

PS. the error message should have included a y block offset which might
have
been helpful.


Best regards,

-- 
---------------------------------------+--------------------------------
------
I set the clouds in motion - turn up   | Frank Warmerdam,
warmerdam at pobox.com
light and sound - activate the windows | http://pobox.com/~warmerdam
and watch the world go round - Rush    | Geospatial Programmer for Rent

_______________________________________________
Gdal-dev mailing list
Gdal-dev at xserve.flids.com
http://xserve.flids.com/mailman/listinfo/gdal-dev