[Qgis-user] Re: RasterCalc not writing projection info
Maxim Dubinin
sim at gis-lab.info
Thu Feb 11 22:45:17 PST 2010
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
More information about the Qgis-user
mailing list