[GIS] Rasterize shapefile to a new gtiff file with gdal_rasterize utility,however,I cannot open the tif file

gdal

Now I'm tring to rasterize a shapefile with a utility "gdal_rasterize" and want to generate a new gtiff file,but when I open the rasterized gtiff file,it prompts "Windows Photo Viewer can't open this picture because the file appears to be damaged,corrupted,or is too large.

……………………………………………………………………………..
……………………………………………………………………………..
Command Line:
gdal_rasterize -burn 255 -burn 0 -burn 0 -ts 800 800 -l Countries02 "D:\Map Data\Countries02.shp" "D:\Result\Countries02.tif"

Map Data:
https://skydrive.live.com/redir?resid=7EF4FFBCD523D8E0!431&authkey=!AKiLt1a0bpOGi7w

The output gdalinfo returns is as follow:
Driver: GTiff/GeoTIFF
Files: d:\Result\World.tif
Size is 800, 800
Coordinate System is `'
Origin = (-180.000000000000000,83.623596191406250)
Pixel Size = (0.450000000000000,-0.217029495239258)
Image Structure Metadata:
INTERLEAVE=PIXEL
Corner Coordinates:
Upper Left (-180.0000000, 83.6235962)
Lower Left (-180.0000000, -90.0000000)
Upper Right ( 180.0000000, 83.6235962)
Lower Right ( 180.0000000, -90.0000000)
Center ( 0.0000000, -3.1882019)
Band 1 Block=800×1 Type=Float64, ColorInterp=Gray
Minimum=0.000, Maximum=255.000, Mean=87.551, StdDev=121.080
Metadata:
STATISTICS_MAXIMUM=255
STATISTICS_MEAN=87.5514609375
STATISTICS_MINIMUM=0
STATISTICS_STDDEV=121.07999102565
Band 2 Block=800×1 Type=Float64, ColorInterp=Undefined
Minimum=0.000, Maximum=0.000, Mean=0.000, StdDev=0.000
Metadata:
STATISTICS_MAXIMUM=0
STATISTICS_MEAN=0
STATISTICS_MINIMUM=0
STATISTICS_STDDEV=0
Band 3 Block=800×1 Type=Float64, ColorInterp=Undefined
Minimum=0.000, Maximum=0.000, Mean=0.000, StdDev=0.000
Metadata:
STATISTICS_MAXIMUM=0
STATISTICS_MEAN=0
STATISTICS_MINIMUM=0
STATISTICS_STDDEV=0

Best Answer

Three suggestions:

Try outputting to another format. Maybe Jpeg or PNG. These will likely work better with Windows Photo Viewer.
Download IrfanView. This will display your tiff quite well.
If neither work - could you run 'gdalinfo -stats Countries02.tif' and post the output.

Edit:

The output is float64. This is likely the issue. Try 'gdal_translate -ot Byte -scale Countries02.tif Countires02_scaled.tif'

I'm not sure why rasterize is creating 64-bit output, but not much reads 64-bit. Also, -scale should not be necessary, but is precautionary.

Related Solutions

[GIS] Georeferencing png using gdal

Here's a world file for you (as per Gene's suggestion). Cut and paste the following into a text editor and name the file the same as your full-sized map but give it the extension '.pgw' instead on '.png':

1.0
0.0
0.0
-1.0
0.0
3231.0

I've made the following assumptions:

Your map has a resolution (pixel size) of 1m square (i.e. each pixel is a reasonable stride length for a fantasy game humanoid character). If I am wrong, edit lines 1 and 4 accordingly.
I have assumed that your map is not rotated.
I have assumed that your LL corner is at 0,0 in your world. If your pixel size is bigger than 1m square then you will need to multiply your Y coordinate (3231.0) accordingly.
Technically the UL coordinate (last two lines) is the center of your pixel. I have been lazy and set it to the corner. This in theory means that your coordinates would be consistently half a meter out but it's a fantasy world so nobody will know; it is consistent so won't screw up any distant calculations; and it saves a negative value as one of your origin coordinates. Adjust if you are sufficiently worried :)

I would then re-tile your map from the georeferenced full-sized one (this saves you making a world file for each tile)

Now, if you load your fantasy map into a GIS or similar system, you will probably be asked for the Spatial Reference System (SRS). The World file doesn't encapsulate that. It just locates and orients the image file in space relative to a nominal 0,0 origin. Since this is a fantasy map, all you need to do is pretend that it is in the real world somewhere and pick any coordinate system that gives you a flat projection with units in meters (or feet if you prefer but you will may have to adjust the values for resolution in the world file). A couple of suggestions for SRS could be:

EPSG: 3857 (aka 900913)
EPSG: 27700 (UK National Grid - the advantage of a national grid system UK or otherwise, is that you mostly stay clear of the poles and the international dateline and can hopefully avoid any weirdness which may result from such a location)

Your choice might not matter at all... on the other hand, you choice of SRS might be influenced by what you intend to do with your KML/CSV file (like you might be making a game with a post-apocalyptic theme but set in Quebec ...with zombies...). If you are using Google Earth or similar as a 'platform' for your game/fantasy world, so long as all your other data are in the same coordinate system, all should be well.

However, I do wonder if you may be better sticking with your PNG file as a text format like KML/CSV will be big. A PNG may be better for an internet-based game (reducing bandwidth). If you have a stand-alone scenario, most programming languages understand image formats well (especially if you call the GDAL libraries) and so I'm not convinced that a text-based format is even necessary here - but that's your choice of course.

[GIS] Can gdalwarp reproject from ESPG:4326 (WGS84) to UTM

Try with other zone it has no problem, I guess GDAL cannot wrap it outside.

gdalwarp -t_srs '+proj=utm +zone=12 +datum=WGS84 +units=m +no_defs' test.tif test_utm1.tif

Best Answer

Related Solutions

[GIS] Georeferencing png using gdal

[GIS] Can gdalwarp reproject from ESPG:4326 (WGS84) to UTM

Related Question