[Qgis-user] Clipped tif not displaying

Steve Steve stevetatman.research at live.co.uk
Thu Feb 25 12:00:09 PST 2016


I have noticed that although Gdalinfo -stats shows the output tif as having data:

    STATISTICS_MAXIMUM=223
    STATISTICS_MEAN=-1.4013838186591e+38
    STATISTICS_MINIMUM=-3.3999999521444e+38
    STATISTICS_STDDEV=1.6735675394745e+38

When I open it in Qgis it shows 
Min -3.4e+38
Max -3.4e+35

Is this a problem with my file or with Qgis?
Cheers, Steve

From: stevetatman.research at live.co.uk
To: qgis-user at lists.osgeo.org
Date: Thu, 25 Feb 2016 19:49:09 +0000
Subject: [Qgis-user] Clipped tif not displaying




I have used this base file:

bio9_49.tif
Driver: GTiff/GeoTIFF
Files: bio9_49.tif
Size is 3600, 3600
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 = (90.000000000000000,-30.000000000000000)
Pixel Size = (0.008333333333333,-0.008333333333333)
Metadata:
  AREA_OR_POINT=Area
Image Structure Metadata:
  COMPRESSION=LZW
  INTERLEAVE=BAND
Corner Coordinates:
Upper Left  (  90.0000000, -30.0000000) ( 90d 0' 0.00"E, 30d 0' 0.00"S)
Lower Left  (  90.0000000, -60.0000000) ( 90d 0' 0.00"E, 60d 0' 0.00"S)
Upper Right ( 120.0000000, -30.0000000) (120d 0' 0.00"E, 30d 0' 0.00"S)
Lower Right ( 120.0000000, -60.0000000) (120d 0' 0.00"E, 60d 0' 0.00"S)
Center      ( 105.0000000, -45.0000000) (105d 0' 0.00"E, 45d 0' 0.00"S)
Band 1 Block=3600x1 Type=Float32, ColorInterp=Gray
  Min=127.000 Max=250.000 
  Minimum=127.000, Maximum=250.000, Mean=0.000, StdDev=0.000
  NoData Value=-3.39999999999999996e+38
  Metadata:
    STATISTICS_MAXIMUM=250
    STATISTICS_MEAN=0
    STATISTICS_MINIMUM=127
    STATISTICS_STDDEV=0

I then used gdal_translate -projwin 114.5 -33.0 120.0 -35.5 -of GTiff bio9_49.tif clip_9.tif

clip_9.tif
Driver: GTiff/GeoTIFF
Files: clip_9.tif
       clip_9.tif.aux.xml
Size is 660, 300
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 = (114.500000000000000,-33.000000000000000)
Pixel Size = (0.008333333333333,-0.008333333333333)
Metadata:
  AREA_OR_POINT=Area
Image Structure Metadata:
  INTERLEAVE=BAND
Corner Coordinates:
Upper Left  ( 114.5000000, -33.0000000) (114d30' 0.00"E, 33d 0' 0.00"S)
Lower Left  ( 114.5000000, -35.5000000) (114d30' 0.00"E, 35d30' 0.00"S)
Upper Right ( 120.0000000, -33.0000000) (120d 0' 0.00"E, 33d 0' 0.00"S)
Lower Right ( 120.0000000, -35.5000000) (120d 0' 0.00"E, 35d30' 0.00"S)
Center      ( 117.2500000, -34.2500000) (117d15' 0.00"E, 34d15' 0.00"S)
Band 1 Block=660x3 Type=Float32, ColorInterp=Gray
  Min=-339999995214440013906259394729726705664.000 Max=223.000 
  Minimum=-339999995214440013906259394729726705664.000, Maximum=223.000, Mean=-140138381865910009253764042088668725248.000, StdDev=167356753947449990662275439509815427072.000
  NoData Value=-3.39999999999999996e+38
  Metadata:
    STATISTICS_MAXIMUM=223
    STATISTICS_MEAN=-1.4013838186591e+38
    STATISTICS_MINIMUM=-3.3999999521444e+38
    STATISTICS_STDDEV=1.6735675394745e+38

Qgis will only display the clipped file as monochrome over the whole land surface.

Can anyone help solve this, it's driving me mad!
Cheers,
Steve
 		 	   		  

_______________________________________________
Qgis-user mailing list
Qgis-user at lists.osgeo.org
List info: http://lists.osgeo.org/mailman/listinfo/qgis-user
Unsubscribe: http://lists.osgeo.org/mailman/listinfo/qgis-user 		 	   		  
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://lists.osgeo.org/pipermail/qgis-user/attachments/20160225/0a7f187a/attachment.html>


More information about the Qgis-user mailing list