एक अन्य विकल्प केवल ओआर 2 आर 2 में रिप्रोडक्शन विकल्पों का उपयोग करना होगा, निश्चित रूप से जेईक्यूएल, फियोना या जियोस्क्रिप्ट दृष्टिकोण की तुलना में एक हैकर दृष्टिकोण लेकिन प्रभावी कोई नहीं-कम। ध्यान दें कि से और अनुमानों को वास्तव में मूल आकार-प्रकार का वास्तविक प्रक्षेपण होने की आवश्यकता नहीं है, जब तक कि एकमात्र चीज जो s_srs और t_srs में उपयोग किए जाने वाले अनुमानों के बीच बदल रही है, झूठे पूर्वलेखन और नॉर्थिंग हैं। इस उदाहरण में मैं सिर्फ Google Mercator का उपयोग कर रहा हूं। मुझे यकीन है कि आधार के रूप में उपयोग करने के लिए बहुत सरल समन्वय प्रणाली है, लेकिन यह मेरे सामने कॉपी / पेस्ट करने के लिए सही था।
ogr2ogr -s_srs EPSG:900913 -t_srs 'PROJCS["Google Mercator",GEOGCS["WGS 84",DATUM["World Geodetic System 1984",SPHEROID["WGS 84",6378137.0,298.257223563,AUTHORITY["EPSG","7030"]],AUTHORITY["EPSG","6326"]],PRIMEM["Greenwich",0.0,AUTHORITY["EPSG","8901"]],UNIT["degree",0.017453292519943295],AXIS["Geodetic latitude",NORTH],AXIS["Geodetic longitude",EAST],AUTHORITY["EPSG","4326"]],PROJECTION["Mercator_1SP"],PARAMETER["semi_minor",6378137.0],PARAMETER["latitude_of_origin",0.0],PARAMETER["central_meridian",0.0],PARAMETER["scale_factor",1.0],PARAMETER["false_easting",1000.0],PARAMETER["false_northing",1000.0],UNIT["m",1.0],AXIS["Easting",EAST],AXIS["Northing",NORTH],AUTHORITY["EPSG","900913"]]' -f "ESRI Shapefile" shift.shp original.shp
या टाइपिंग / पेस्टिंग को बचाने के लिए, निम्न को projcs.txt
(ऊपर के समान, लेकिन एकल उद्धरण को हटाते हुए) सहेजें :
-s_srs EPSG:900913
-t_srs PROJCS["Google Mercator",GEOGCS["WGS 84",DATUM["World Geodetic System 1984",SPHEROID["WGS 84",6378137.0,298.257223563,AUTHORITY["EPSG","7030"]],AUTHORITY["EPSG","6326"]],PRIMEM["Greenwich",0.0,AUTHORITY["EPSG","8901"]],UNIT["degree",0.017453292519943295],AXIS["Geodetic latitude",NORTH],AXIS["Geodetic longitude",EAST],AUTHORITY["EPSG","4326"]],PROJECTION["Mercator_1SP"],PARAMETER["semi_minor",6378137.0],PARAMETER["latitude_of_origin",0.0],PARAMETER["central_meridian",0.0],PARAMETER["scale_factor",1.0],PARAMETER["false_easting",1000.0],PARAMETER["false_northing",1000.0],UNIT["m",1.0],AXIS["Easting",EAST],AXIS["Northing",NORTH],AUTHORITY["EPSG","900913"]]
और फिर चलाएं:
ogr2ogr --optfile projcs.txt shifted.shp input.shp