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

[Maxim Dubinin]

Agus, it would be nice if you can share a raster and formula, that
you're testing.

Maxim

Вы писали 12 февраля 2010 г., 0:34:43:

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

AL> Agus

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


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





AL> 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?
>> 

AL> _______________________________________________
AL> Qgis-user mailing list
AL> Qgis-user at lists.osgeo.org
AL> http://lists.osgeo.org/mailman/listinfo/qgis-user