Re: [Qgis-user] Re: RasterCalc not writing projection info
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> al...@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.0005, 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",50], AL> PARAMETER["false_northing",0], AL> UNIT["metre",1, AL> AUTHORITY["EPSG","9001"]], AL> AUTHORITY["EPSG","23031"]] AL> Origin = (424389.000,4635822.000) AL> Pixel Size = (345.0078125,194.01953125000) 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> al...@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.000,4635822.000) AL> Pixel Size = (345.0078125,194.01953125000) 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> al...@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 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@lists.osgeo.org AL> http://lists.osgeo.org/mailman/listinfo/qgis-user ___ Qgis-user mailing list Qgis-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/qgis-user
[Qgis-user] Re: RasterCalc not writing projection info
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 al...@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.0005, 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",50], PARAMETER["false_northing",0], UNIT["metre",1, AUTHORITY["EPSG","9001"]], AUTHORITY["EPSG","23031"]] Origin = (424389.000,4635822.000) Pixel Size = (345.0078125,194.01953125000) 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 al...@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.000,4635822.000) Pixel Size = (345.0078125,194.01953125000) 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 al...@alobo-laptop:/media/Transcend/MASTER_ICTA2007_2008/miniproj_GISInP2009/MARICEL2$ Alexander Bruy wrote: On Wed, 10 Feb 2010 19:08:15 +0100 Agustin Lobo 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? ___ Qgis-user mailing list Qgis-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/qgis-user
[Qgis-user] Re: RasterCalc not writing projection info
On Wed, 10 Feb 2010 19:08:15 +0100 Agustin Lobo 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? -- Alexander Bruy ___ Qgis-user mailing list Qgis-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/qgis-user
