I am using the command: gdalwarp -t_srs WGS84 -te 144 -36 150 -32 -crop_to_cutline AustTopo250_SI55.tif
I need to change the pixel resolution to 50m x 50m but the option is given as -tr xres yres: set output file resolution (in target georeferenced units).
The target unit is degrees. So the question is: How do I generate the output file as 50m x 50m pixels?
Feel free to explain relationship between metres and degrees.
Target File:
E:\>gdalinfo AustTopo250_SI55.tif Driver: GTiff/GeoTIFF Files: AustTopo250_SI55.tif Size is 18750, 12500 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 = (144.000000000000000,-32.000000000000000) Pixel Size = (0.000320000000000,-0.000320000000000) Metadata: AREA_OR_POINT=Area COLORSPACE=RGB COMPRESSION_RATE_TARGET=9 VERSION=2 Image Structure Metadata: INTERLEAVE=PIXEL Corner Coordinates: Upper Left ( 144.0000000, -32.0000000) (144d 0' 0.00"E, 32d 0' 0.00"S) Lower Left ( 144.0000000, -36.0000000) (144d 0' 0.00"E, 36d 0' 0.00"S) Upper Right ( 150.0000000, -32.0000000) (150d 0' 0.00"E, 32d 0' 0.00"S) Lower Right ( 150.0000000, -36.0000000) (150d 0' 0.00"E, 36d 0' 0.00"S) Center ( 147.0000000, -34.0000000) (147d 0' 0.00"E, 34d 0' 0.00"S) Band 1 Block=18750x1 Type=Byte, ColorInterp=Red Description = Red Band 2 Block=18750x1 Type=Byte, ColorInterp=Green Description = Green Band 3 Block=18750x1 Type=Byte, ColorInterp=Blue Description = Blue Source File:
E:\>gdalinfo 76943_Aust_gda94.ecw Driver: ECW/ERDAS Compressed Wavelets (SDK 3.x) Files: 76943_Aust_gda94.ecw Size is 128490, 109375 Coordinate System is: GEOGCS["GEOCENTRIC DATUM of AUSTRALIA", DATUM["GDA94", SPHEROID["GRS80",6378137,298.257222101]], PRIMEM["Greenwich",0], UNIT["degree",0.0174532925199433]] Origin = (112.883360000000040,-8.999920000000039) Pixel Size = (0.000320000000000,-0.000320000000000) Metadata: COLORSPACE=RGB COMPRESSION_RATE_TARGET=9 VERSION=2 Corner Coordinates: Upper Left ( 112.8833600, -8.9999200) (112d53' 0.10"E, 8d59'59.71"S) Lower Left ( 112.8833600, -43.9999200) (112d53' 0.10"E, 43d59'59.71"S) Upper Right ( 154.0001600, -8.9999200) (154d 0' 0.58"E, 8d59'59.71"S) Lower Right ( 154.0001600, -43.9999200) (154d 0' 0.58"E, 43d59'59.71"S) Center ( 133.4417600, -26.4999200) (133d26'30.34"E, 26d29'59.71"S) Band 1 Block=256x256 Type=Byte, ColorInterp=Red Description = Red Overviews: 64245x54687, 32122x27343, 16061x13671, 8030x6835, 4015x3417, 2007x1708, 1003x854, 501x427, 250x213 Band 2 Block=256x256 Type=Byte, ColorInterp=Green Description = Green Overviews: 64245x54687, 32122x27343, 16061x13671, 8030x6835, 4015x3417, 2007x1708, 1003x854, 501x427, 250x213 Band 3 Block=256x256 Type=Byte, ColorInterp=Blue Description = Blue Overviews: 64245x54687, 32122x27343, 16061x13671, 8030x6835, 4015x3417, 2007x1708, 1003x854, 501x427, 250x213
أكثر...
I need to change the pixel resolution to 50m x 50m but the option is given as -tr xres yres: set output file resolution (in target georeferenced units).
The target unit is degrees. So the question is: How do I generate the output file as 50m x 50m pixels?
Feel free to explain relationship between metres and degrees.
Target File:
E:\>gdalinfo AustTopo250_SI55.tif Driver: GTiff/GeoTIFF Files: AustTopo250_SI55.tif Size is 18750, 12500 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 = (144.000000000000000,-32.000000000000000) Pixel Size = (0.000320000000000,-0.000320000000000) Metadata: AREA_OR_POINT=Area COLORSPACE=RGB COMPRESSION_RATE_TARGET=9 VERSION=2 Image Structure Metadata: INTERLEAVE=PIXEL Corner Coordinates: Upper Left ( 144.0000000, -32.0000000) (144d 0' 0.00"E, 32d 0' 0.00"S) Lower Left ( 144.0000000, -36.0000000) (144d 0' 0.00"E, 36d 0' 0.00"S) Upper Right ( 150.0000000, -32.0000000) (150d 0' 0.00"E, 32d 0' 0.00"S) Lower Right ( 150.0000000, -36.0000000) (150d 0' 0.00"E, 36d 0' 0.00"S) Center ( 147.0000000, -34.0000000) (147d 0' 0.00"E, 34d 0' 0.00"S) Band 1 Block=18750x1 Type=Byte, ColorInterp=Red Description = Red Band 2 Block=18750x1 Type=Byte, ColorInterp=Green Description = Green Band 3 Block=18750x1 Type=Byte, ColorInterp=Blue Description = Blue Source File:
E:\>gdalinfo 76943_Aust_gda94.ecw Driver: ECW/ERDAS Compressed Wavelets (SDK 3.x) Files: 76943_Aust_gda94.ecw Size is 128490, 109375 Coordinate System is: GEOGCS["GEOCENTRIC DATUM of AUSTRALIA", DATUM["GDA94", SPHEROID["GRS80",6378137,298.257222101]], PRIMEM["Greenwich",0], UNIT["degree",0.0174532925199433]] Origin = (112.883360000000040,-8.999920000000039) Pixel Size = (0.000320000000000,-0.000320000000000) Metadata: COLORSPACE=RGB COMPRESSION_RATE_TARGET=9 VERSION=2 Corner Coordinates: Upper Left ( 112.8833600, -8.9999200) (112d53' 0.10"E, 8d59'59.71"S) Lower Left ( 112.8833600, -43.9999200) (112d53' 0.10"E, 43d59'59.71"S) Upper Right ( 154.0001600, -8.9999200) (154d 0' 0.58"E, 8d59'59.71"S) Lower Right ( 154.0001600, -43.9999200) (154d 0' 0.58"E, 43d59'59.71"S) Center ( 133.4417600, -26.4999200) (133d26'30.34"E, 26d29'59.71"S) Band 1 Block=256x256 Type=Byte, ColorInterp=Red Description = Red Overviews: 64245x54687, 32122x27343, 16061x13671, 8030x6835, 4015x3417, 2007x1708, 1003x854, 501x427, 250x213 Band 2 Block=256x256 Type=Byte, ColorInterp=Green Description = Green Overviews: 64245x54687, 32122x27343, 16061x13671, 8030x6835, 4015x3417, 2007x1708, 1003x854, 501x427, 250x213 Band 3 Block=256x256 Type=Byte, ColorInterp=Blue Description = Blue Overviews: 64245x54687, 32122x27343, 16061x13671, 8030x6835, 4015x3417, 2007x1708, 1003x854, 501x427, 250x213
أكثر...