[Qgis-user] problems to cut images - black result

[Jan Tappenbeck]

Hi !

i have a little big problem. my job is to cut image in userdefine tiles.

in the first command you can see a call with success - in the secound 
the result is only black images.

here are some informations.

command and log of image-cut with success

Bild 1
5681411.1 - 5685564.3 | 362500.00 - 5683500.00

gdalwarp -te 360638.4 5683500.00 362500.00 5685564.3 -of GTiff 
E:/datenlager_yyyyy/20130218_xxxxx_RasterbildTransformation_GK2ETRS/ETRS/1988/1988.tif 
E:/datenlager_yyyyy/20130218_xxxxx_RasterbildTransformation_GK2ETRS/ETRS/1988/1988_1.tif

Creating output file that is 2975P x 3299L.
Processing input file 
E:/datenlager_yyyyy/20130218_xxxxx_RasterbildTransformation_GK2ETRS/ETRS/1988/1988.tif.
Using band 4 of source image as alpha.
Using band 4 of destination image as alpha.
0...10...20...30...40...50...60...70...80...90...100 - done.
Input file size is 2975, 3299

now i create a new command to cut parts of a large image

new cut-definiton for a 4GB-Image

;;Z7-S6-Bild: xxxxx_365500_5673000
;;Koord - oben links: 365500 - 5673000
;;Koord - unten rechts: 366500 - 5672000

gdalwarp -te 365500 5672000 366500 5673000 -of GTiff 
E:/datenlager_yyyyy/20130218_xxxxx_RasterbildTransformation_GK2ETRS/ETRS/2002/2002.tif 
E:/datenlager_yyyyy/20130218_xxxxx_RasterbildTransformation_GK2ETRS/ETRS/2002/2002__365500_5673000.tif

Creating output file that is 3201P x 3201L.
Processing input file 
E:/datenlager_yyyyy/20130218_xxxxx_RasterbildTransformation_GK2ETRS/ETRS/2002/2002.tif.
Using band 4 of source image as alpha.
Using band 4 of destination image as alpha.
0...10...20...30...40...50...60...70...80...90...100 - done.
;;Z7-S7-Bild: xxxxx_365500_5672000

data of the large image

Driver: GTiff/GeoTIFF
Files: 
E:/datenlager_yyyyy/20130218_xxxxx_RasterbildTransformation_GK2ETRS/ETRS/2002/2002.tif
Size is 33281, 33281
Coordinate System is:
PROJCS["ETRS89 / UTM zone 32N",
     GEOGCS["ETRS89",
         DATUM["European_Terrestrial_Reference_System_1989",
             SPHEROID["GRS 1980",6378137,298.2572221010002,
                 AUTHORITY["EPSG","7019"]],
             AUTHORITY["EPSG","6258"]],
         PRIMEM["Greenwich",0],
         UNIT["degree",0.0174532925199433],
         AUTHORITY["EPSG","4258"]],
     PROJECTION["Transverse_Mercator"],
     PARAMETER["latitude_of_origin",0],
     PARAMETER["central_meridian",9],
     PARAMETER["scale_factor",0.9996],
     PARAMETER["false_easting",500000],
     PARAMETER["false_northing",0],
     UNIT["metre",1,
         AUTHORITY["EPSG","9001"]],
     AUTHORITY["EPSG","25832"]]
Origin = (358564.810486490660000,5689657.740791935500000)
Pixel Size = (0.312428079160783,-0.312428079160783)
Metadata:
   AREA_OR_POINT=Area
Image Structure Metadata:
   INTERLEAVE=PIXEL
Corner Coordinates:
Upper Left  (  358564.810, 5689657.741) (  6d58' 9.79"E, 51d20'26.10"N)
Lower Left  (  358564.810, 5679259.822) (  6d58'24.62"E, 51d14'49.72"N)
Upper Right (  368962.729, 5689657.741) (  7d 7' 6.89"E, 51d20'35.08"N)
Lower Right (  368962.729, 5679259.822) (  7d 7'20.63"E, 51d14'58.66"N)
Center      (  363763.770, 5684458.781) (  7d 2'45.48"E, 51d17'42.48"N)
Band 1 Block=33281x1 Type=Byte, ColorInterp=Red
   Mask Flags: PER_DATASET ALPHA
Band 2 Block=33281x1 Type=Byte, ColorInterp=Green
   Mask Flags: PER_DATASET ALPHA
Band 3 Block=33281x1 Type=Byte, ColorInterp=Blue
   Mask Flags: PER_DATASET ALPHA
Band 4 Block=33281x1 Type=Byte, ColorInterp=Alpha

the results are every time image wiht the same size and all are black - 
info:
Driver: GTiff/GeoTIFF
Files: 
E:/datenlager_yyyyy/20130218_xxxxx_RasterbildTransformation_GK2ETRS/ETRS/2002/2002__358500_5678000.tif
Size is 3201, 3201
Coordinate System is:
PROJCS["ETRS89 / UTM zone 32N",
     GEOGCS["ETRS89",
         DATUM["European_Terrestrial_Reference_System_1989",
             SPHEROID["GRS 1980",6378137,298.2572221010002,
                 AUTHORITY["EPSG","7019"]],
             AUTHORITY["EPSG","6258"]],
         PRIMEM["Greenwich",0],
         UNIT["degree",0.0174532925199433],
         AUTHORITY["EPSG","4258"]],
     PROJECTION["Transverse_Mercator"],
     PARAMETER["latitude_of_origin",0],
     PARAMETER["central_meridian",9],
     PARAMETER["scale_factor",0.9996],
     PARAMETER["false_easting",500000],
     PARAMETER["false_northing",0],
     UNIT["metre",1,
         AUTHORITY["EPSG","9001"]],
     AUTHORITY["EPSG","25832"]]
Origin = (358500.000000000000000,5678000.000000000000000)
Pixel Size = (0.312402374258044,-0.312402374258044)
Metadata:
   AREA_OR_POINT=Area
Image Structure Metadata:
   INTERLEAVE=PIXEL
Corner Coordinates:
Upper Left  (  358500.000, 5678000.000) (  6d58'23.07"E, 51d14' 8.91"N)
Lower Left  (  358500.000, 5677000.000) (  6d58'24.49"E, 51d13'36.55"N)
Upper Right (  359500.000, 5678000.000) (  6d59'14.61"E, 51d14' 9.80"N)
Lower Right (  359500.000, 5677000.000) (  6d59'16.02"E, 51d13'37.44"N)
Center      (  359000.000, 5677500.000) (  6d58'49.55"E, 51d13'53.18"N)
Band 1 Block=3201x1 Type=Byte, ColorInterp=Red
   Mask Flags: PER_DATASET ALPHA
Band 2 Block=3201x1 Type=Byte, ColorInterp=Green
   Mask Flags: PER_DATASET ALPHA
Band 3 Block=3201x1 Type=Byte, ColorInterp=Blue
   Mask Flags: PER_DATASET ALPHA
Band 4 Block=3201x1 Type=Byte, ColorInterp=Alpha

did anyone have a idea why the result images are only black. if the 
regtang-defintion was wrong normally i get a message that the border is 
outside of the image.

regards Jan :-)