I have some shapefiles of the world from Natural Earth at WGS84. I would like to convert them to the Miller Projection. I have ogr2ogr, but I’m not sure what the parameters for this conversion should be, or if this conversion can be done with it. Test if target srs epsg:54003 gives correct result. You can check with some coordinates with gdaltransform first: Because the eastern part of Russia crosses the 180° meridian, you have to cut your shapefile there. I have written a tutorial for cutting at any meridian to use the Natural Earth shapefiles here: QGIS display world country shape files centered on pacific ocean using Robinson, Miller Cylindrical or other projection The problem is that straight lines in one CRS won’t necessarily be straight in another – unless they’re both conformal – so transforming a line segment will position its end points accurately, but no new points will be added in between, so the line as a whole won’t be accurate. In your case, if you transform then rasterize, and the features are smallish in size, generally covering less that a square kilometre or so, there will be no appreciable error at any reasonable resolution of raster. But as you head towards the poles, or as your features get bigger, you’ll notice greater error. So there will come a point at which the accuracy of your source data is better than the accuracy of your projected data – your allowable error. Computationally it’s generally quicker to transform vector data before rasterizing it, so you would have to weigh that up if you have any time constraints. If you have time for some experimentation, I would generate two rasters using a subset of your data – one using the transform/rasterize process, and one the other way round. Then combine the two bitmaps with an exclusive-or operation, which will highlight areas where there is greater than one pixel error, and you can make a judgement on whether any error is allowable. You can also compare time taken for processing to further inform your decision. Source.