[Qgis-user] Re: RasterCalc not writing projection info

[Agustin Lobo]

Alexander,
I've tried with 0.1.18 (1.4.0 from binaries for ubuntu 9.04), same result
Are you sure that your plugin is writing all fields for the projection? Have you 
tested
with gdalinfo? Note the output for test2.tif (input) and test2_0.tif (output).
Note that you can get this output from within qgis using Information
in Raster (plugin gdalutilities)

Agus

alobo at alobo-laptop:/media/Transcend/MASTER_ICTA2007_2008/miniproj_GISInP2009/MARICEL2$ 
gdalinfo test2.tif
Driver: GTiff/GeoTIFF
Files: test2.tif
Size is 256, 256
Coordinate System is:
PROJCS["ED50 / UTM zone 31N",
     GEOGCS["ED50",
         DATUM["European_Datum_1950",
             SPHEROID["International 1924",6378388,297.000000000005,
                 AUTHORITY["EPSG","7022"]],
             AUTHORITY["EPSG","6230"]],
         PRIMEM["Greenwich",0],
         UNIT["degree",0.0174532925199433],
         AUTHORITY["EPSG","4230"]],
     PROJECTION["Transverse_Mercator"],
     PARAMETER["latitude_of_origin",0],
     PARAMETER["central_meridian",3],
     PARAMETER["scale_factor",0.9996],
     PARAMETER["false_easting",500000],
     PARAMETER["false_northing",0],
     UNIT["metre",1,
         AUTHORITY["EPSG","9001"]],
     AUTHORITY["EPSG","23031"]]
Origin = (424389.000000000000000,4635822.000000000000000)
Pixel Size = (345.007812500000000,194.019531250000000)
Metadata:
   AREA_OR_POINT=Area
Image Structure Metadata:
   INTERLEAVE=BAND
Corner Coordinates:
Upper Left  (  424389.000, 4635822.000) (  2d 5'20.10"E, 41d52'11.88"N)
Lower Left  (  424389.000, 4685491.000) (  2d 4'56.97"E, 42d19'2.06"N)
Upper Right (  512711.000, 4635822.000) (  3d 9'11.42"E, 41d52'24.52"N)
Lower Right (  512711.000, 4685491.000) (  3d 9'15.31"E, 42d19'14.90"N)
Center      (  468550.000, 4660656.500) (  2d37'10.89"E, 42d 5'47.83"N)
Band 1 Block=256x4 Type=Float64, ColorInterp=Gray


alobo at alobo-laptop:/media/Transcend/MASTER_ICTA2007_2008/miniproj_GISInP2009/MARICEL2$ 
gdalinfo test2_0.tif
Driver: GTiff/GeoTIFF
Files: test2_0.tif
Size is 256, 256
Coordinate System is `'
Origin = (424389.000000000000000,4635822.000000000000000)
Pixel Size = (345.007812500000000,194.019531250000000)
Image Structure Metadata:
   INTERLEAVE=BAND
Corner Coordinates:
Upper Left  (  424389.000, 4635822.000)
Lower Left  (  424389.000, 4685491.000)
Upper Right (  512711.000, 4635822.000)
Lower Right (  512711.000, 4685491.000)
Center      (  468550.000, 4660656.500)
Band 1 Block=256x32 Type=Byte, ColorInterp=Gray
alobo at alobo-laptop:/media/Transcend/MASTER_ICTA2007_2008/miniproj_GISInP2009/MARICEL2$ 





Alexander Bruy wrote:
> On Wed, 10 Feb 2010 19:08:15 +0100
> Agustin Lobo <alobolistas at gmail.com> wrote:
> 
>> Alexander,
>> As far as I'm able to see, RasterCalc is not writing projection information
>> to geotif files:
>> I just multiplied by 0 a layer coming form a tif file (test2.tif) and saved
>> as test2_0.tif This is what gdalinfo reports (note field "projection"):
> 
> Thanks, Agustin
> 
> which version of RasterCalc do you use? I've test 0.1.17 and 0.1.18
> versions and projection info is writed correctly. Can you try with 0.1.18?
>