Mapinfo Utility for gdal_translate.

rupert — Wed, 20 Jun 2007 17:44:36 +0000

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

1. ls *.TAB > init.sh

2. 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

3. 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

mobile and gis dev notes » perl

Mapinfo Utility for gdal_translate.