Mapinfo Utility for gdal_translate.
I made a small perl utility to automate the gcp’s from Mapinfo Raster TABS to gdal_translate command line. Currently your tabfile would have:
!table !version 300 !charset WindowsLatin1 Definition Table File "beijing_6th_1.jpg" Type "RASTER" (116.522865,40.016316) (347,184) Label "Pt 1", (116.681215,40.015286) (7729,243) Label "Pt 2", (116.679777,39.777904) (7666,14674) Label "Pt 3", (116.523827,39.779108) (397,14606) Label "Pt 4" CoordSys Earth Projection 1, 104 Units "degree" RasterStyle 4 1 RasterStyle 7 1677695 |
How to use
gdal_mapinfo
-
ls *.TAB > init.sh
-
vi init.sh to reflect the ff:
perl gdal_mapinfo.pl Beijing_6th_1.TAB >> final.bat perl gdal_mapinfo.pl Beijing_6th_10.TAB >> final.bat perl gdal_mapinfo.pl Beijing_6th_11.TAB >> final.bat perl gdal_mapinfo.pl Beijing_6th_12.TAB >> final.bat perl gdal_mapinfo.pl Beijing_6th_13.TAB >> final.bat perl gdal_mapinfo.pl Beijing_6th_2.TAB >> final.bat perl gdal_mapinfo.pl Beijing_6th_3.TAB >> final.bat perl gdal_mapinfo.pl Beijing_6th_4.TAB >> final.bat perl gdal_mapinfo.pl Beijing_6th_5.TAB >> final.bat perl gdal_mapinfo.pl Beijing_6th_6.TAB >> final.bat perl gdal_mapinfo.pl Beijing_6th_7.TAB >> final.bat perl gdal_mapinfo.pl Beijing_6th_8.TAB >> final.bat perl gdal_mapinfo.pl Beijing_6th_9.TAB >> final.bat |
- The resulting final.bat should have the ff:
gdal_translate -gcp 347 184 116.522865 40.016316 -gcp 7729 243 116.681215 40.015286 -gcp 7666 14674 116.679777 39.777904 -gcp 397 14606 116.523827 39.779108 -of GTiff Beijing_6th_1.jpg I:\satimages\translated\Beijing_6th_1_translated.tif
gdalwarp -s_srs epsg:4326 -t_srs epsg:4326 I:\satimages\translated\Beijing_6th_1_translated.tif
I:\satimages\warped\Beijing_6th_1.tif </pre>