Legacy ImageMagick Discussions Archive

Code: Select all

#! /bin/bash
# script to tile high res global images


LIST="A1 A2 B1 B2 C1 C2 D1 D2"

for TEXT in $LIST  ; do
  MAP=world.topo.bathy.200407.3x21600x21600.${TEXT}.png
  if [ ! -r tile_${TEXT}.png.0 ] && [ ! -r tile_${TEXT}_00.png ] ; then
    echo "Tiling $MAP"
    convert $MAP -crop 2160x2160! +repage tile_${TEXT}.png
  else
    echo "$MAP already tiled"
  fi

  #rename files if required
  echo "Renaming new tiles if required"
  if [ -r tile_${TEXT}.png.0 ] ; then
    echo "Renaming tiles...."
    rm -f files.txt
    ls tile_${TEXT}.png.* > files.txt
    while read FILE ; do
      TILE_NO=`echo "$FILE" | cut -f3 -d"."`
      TILE_NAME=`echo "$FILE" | cut -f1 -d"."`
      if [ $TILE_NO -le 9 ] ; then
        TILE_NO=0$TILE_NO
      fi
      NEW_NAME=${TILE_NAME}_${TILE_NO}.png
      echo "renaming $FILE to $NEW_NAME"
      mv $FILE $NEW_NAME
    done < files.txt
  fi
done

#write world files....
#A1 origin=180, 90
#B1 origin=270, 90
#C1 origin=0,   90
#D1 origin=90,  90
#all x2 files have Y origin of 0, same x origin

X_PIXEL=0.004166666667
Y_PIXEL=-0.004166666667

echo "generating world files"

for TEXT in $LIST ; do
  rm -f files.txt

  ls tile_${TEXT}*.png > files.txt

  if [ "$TEXT" = "A1" ] ; then
    X_ORIGIN=180
    Y_ORIGIN=90
  elif [ "$TEXT" = "A2" ] ; then
    X_ORIGIN=180
    Y_ORIGIN=0
  elif [ "$TEXT" = "B1" ] ; then
    X_ORIGIN=270
    Y_ORIGIN=90
  elif [ "$TEXT" = "B2" ] ; then
    X_ORIGIN=270
    Y_ORIGIN=0
  elif [ "$TEXT" = "C1" ] ; then
    X_ORIGIN=0
    Y_ORIGIN=90
  elif [ "$TEXT" = "C2" ] ; then
    X_ORIGIN=0
    Y_ORIGIN=0
  elif [ "$TEXT" = "D1" ] ; then
    X_ORIGIN=90
    Y_ORIGIN=90
  elif [ "$TEXT" = "D2" ] ; then
    X_ORIGIN=90
    Y_ORIGIN=0
  fi

  while read FILE ; do
    PREFIX=`echo "$FILE" | cut -f1 -d"."`
    PGW=${PREFIX}.pgw
    echo "Writing $PGW for $FILE"

    ROW=`echo "$FILE" | cut -f3 -d"_" | cut -b1`
    COL=`echo "$FILE" | cut -f3 -d"_" | cut -b2`
    echo " using row $ROW    column $COL"

    # get each new tile origin
    X_TILE=`echo "1.000000000000 * $X_ORIGIN + ( $COL * 9.0 )" | bc`
    Y_TILE=`echo "1.000000000000 * $Y_ORIGIN - ( $ROW * 9.0 )" | bc`
    echo "  Origin: $X_TILE   $Y_TILE"

    rm -f $PGW

    echo "$X_PIXEL"  >  $PGW
    echo "0"         >> $PGW
    echo "0"         >> $PGW
    echo "$Y_PIXEL"  >> $PGW
    echo "$X_TILE" >> $PGW
    echo "$Y_TILE" >> $PGW


  done < files.txt
done