#!/bin/bash
<br>#
<br># Developed by Fred Weinhaus 6/20/2009 .......... revised 1/23/2013
<br>#
<br># USAGE: spherize [-d diameter(s)] [-a amp] [-z zoom] [-b bcolor ] [-s] [-t] infile outfile
<br># USAGE: spherize [-h or -help]
<br>#
<br># OPTIONS:
<br>#
<br># -d        diameter(s)         diameters; comma separate x and y diameters;
<br>#                               integers>0 and less than or equal to image
<br>#                               dimensions; if only one value provided,
<br>#                               it will be used for both; default is the
<br>#                               image "width,height"
<br># -a        amp                 distortion amplification factor; float>=0;
<br>#                               default=0
<br># -z        zoom                zoom factor; scales the input image on the
<br>#                               sphere; float>=0; default=1
<br># -b        bcolor              background color; any valid IM color or
<br>#                               "none" for transparent or "image" for input
<br>#                               image background; default=black
<br># -s                            make result with both diameters equal to the
<br>#                               smaller of the x or y diameters
<br># -t                            trim image to diameters
<br>#
<br>###
<br>#
<br># NAME: SPHERIZE
<br>#
<br># PURPOSE: To warp an image onto a (hemi-)sphere.
<br>#
<br># DESCRIPTION: UNROTATE warps an image onto a (hemi-)sphere and views it
<br># orthographically. The x and y diameters can be specified along with the
<br># distortion amplification and input image zoom. The background around
<br># the sphere is also user controlled. This is a limited version of my
<br># bubblewarp script that is not limited by the use of the IM -fx function
<br># and therefore is very fast.
<br>#
<br>#
<br># Arguments:
<br>#
<br># -d diameter(s) --- DIAMETER(S) is a comma separate list of the desired
<br># x and y diameters of the sphere. If only one value is provided, then it
<br># will be used for both diameters. Values are integers greater than zero
<br># and equal to or smaller than the dimensions of the image. The default
<br># values are the size of the image. Thus if the image is square, then a
<br># spherical effect will be produced. But if the image is not square, then
<br># elliptical effect will be produced.
<br>#
<br># -a amp --- AMP is the distortion amplification factor. Values are floats
<br># greater than or equal to zero. The default=1. Values larger than one will
<br># make more distortion and values less than one will make less distortion.
<br>#
<br># -z zoom --- ZOOM is the zoom factor for the input image before warping it
<br># onto the sphere. Values are floats greater than or equal to zero. The
<br># default=1. Values larger than one will magnify the input image on the sphere.
<br># Values less than one will minify the input image on the sphere.
<br>#
<br># -b bcolor ... BCOLOR is the color for the background outside the sphere.
<br># Any valid IM color is allowed or one may specify "none" for a transparent
<br># background or "image" to have the input image as the background. The default
<br># is black.
<br>#
<br># -s ... Make the result with both diameters equal to the smaller of the
<br># x or y diameters.
<br>#
<br># -t ... Trim the output image to the diameters.
<br>#
<br># NOTE: This script is a limited version of my bubblewarp script, but is
<br># much faster, since it does not use -fx.
<br>#
<br># NOTE: This script requires IM 6.5.3-9 or higher due to the use of
<br># -compose distort. Many thanks to Anthony Thyssen for his hard work
<br># developing it.
<br>#
<br># CAVEAT: No guarantee that this script will work on all platforms,
<br># nor that trapping of inconsistent parameters is complete and
<br># foolproof. Use At Your Own Risk.
<br>#
<br>######
<br>#
<br># set default values;
<br>dx=""					# x diameter
<br>dy=""					# y diameter
<br>amp=1					# distortion amplification
<br>zoom=1					# input image zoom factor
<br>bcolor="black"			# background color; none for transparency; or "image"
<br>symmetry="no"			# sphere with both diameters equal to min of x and y diameters
<br>trim="no"				# trim output to diameters
<br>
<br># set directory for temporary files
<br>dir="."    # suggestions are dir="." or dir="/tmp"
<br>
<br># set up functions to report Usage and Usage with Description
<br>PROGNAME=`type $0 | awk '{print $3}'`  # search for executable on path
<br>PROGDIR=`dirname $PROGNAME`            # extract directory of program
<br>PROGNAME=`basename $PROGNAME`          # base name of program
<br>usage1()
<br>	{
<br>	echo >&2 ""
<br>	echo >&2 "$PROGNAME:" "$@"
<br>	sed >&2 -n '/^###/q;  /^#/!q;  s/^#//;  s/^ //;  4,$p' "$PROGDIR/$PROGNAME"
<br>	}
<br>usage2()
<br>	{
<br>	echo >&2 ""
<br>	echo >&2 "$PROGNAME:" "$@"
<br>	sed >&2 -n '/^######/q;  /^#/!q;  s/^#*//;  s/^ //;  4,$p' "$PROGDIR/$PROGNAME"
<br>	}
<br>
<br># function to report error messages
<br>errMsg()
<br>	{
<br>	echo ""
<br>	echo $1
<br>	echo ""
<br>	usage1
<br>	exit 1
<br>	}
<br>
<br># function to test for minus at start of value of second part of option 1 or 2
<br>checkMinus()
<br>	{
<br>	test=`echo "$1" | grep -c '^-.*$'`   # returns 1 if match; 0 otherwise
<br>    [ $test -eq 1 ] && errMsg "$errorMsg"
<br>	}
<br>
<br># test for correct number of arguments and get values
<br>if [ $# -eq 0 ]
<br>	then
<br>	# help information
<br>	echo ""
<br>	usage2
<br>	exit 0
<br>elif [ $# -gt 12 ]
<br>	then
<br>	errMsg "--- TOO MANY ARGUMENTS WERE PROVIDED ---"
<br>else
<br>	while [ $# -gt 0 ]
<br>		do
<br>		# get parameters
<br>		case "$1" in
<br>	  -h|-help)    # help information
<br>				   echo ""
<br>				   usage2
<br>				   ;;
<br>			-d)    # diameters
<br>				   shift  # to get the next parameter
<br>				   # test if parameter starts with minus sign
<br>				   errorMsg="--- INVALID DIAMETERS SPECIFICATION ---"
<br>				   checkMinus "$1"
<br>				   diameters=$1
<br>				   test=`echo "$1" | tr "," " " | wc -w`
<br>				   [ $test -eq 1 -o $test -gt 2 ] && errMsg "--- INCORRECT NUMBER OF DIAMETERS SUPPLIED ---"
<br>				   diameters=`expr "$1" : '\([0-9]*,[0-9]*\)'`
<br>				   [ "$diameters" = "" ] && errMsg "--- DIAMETERS=$diameters MUST BE A PAIR OF NON-NEGATIVE INTEGERS SEPARATED BY A COMMA ---"
<br>				   diameters="$1,"
<br>				   dx=`echo "$diameters" | cut -d, -f1`
<br>				   dy=`echo "$diameters" | cut -d, -f2`
<br>				   # further testing done later
<br>				   ;;
<br>			-a)    # amp
<br>				   shift  # to get the next parameter
<br>				   # test if parameter starts with minus sign
<br>				   errorMsg="--- INVALID AMPLIFICATION SPECIFICATION ---"
<br>				   checkMinus "$1"
<br>				   amp=`expr "$1" : '\([.0-9]*\)'`
<br>				   [ "$amp" = "" ] && errMsg "--- AMPLIFICATION=$amp MUST BE A NON-NEGATIVE FLOATING POINT VALUE (with no sign) ---"
<br>				   ;;
<br>			-z)    # zoom
<br>				   shift  # to get the next parameter
<br>				   # test if parameter starts with minus sign
<br>				   errorMsg="--- INVALID ZOOM SPECIFICATION ---"
<br>				   checkMinus "$1"
<br>				   zoom=`expr "$1" : '\([.0-9]*\)'`
<br>				   [ "$zoom" = "" ] && errMsg "--- ZOOM=$zoom MUST BE A NON-NEGATIVE FLOATING POINT VALUE (with no sign) ---"
<br>				   ;;
<br>			-b)    # get bcolor
<br>				   shift  # to get the next parameter
<br>				   # test if parameter starts with minus sign
<br>				   errorMsg="--- INVALID BACKGROUND COLOR SPECIFICATION ---"
<br>				   checkMinus "$1"
<br>				   bcolor="$1"
<br>				   ;;
<br>			-s)    # set symmetry
<br>				   symmetry="yes"
<br>				   ;;
<br>			-t)    # set trim
<br>				   trim="yes"
<br>				   ;;
<br>			 -)    # STDIN and end of arguments
<br>				   break
<br>				   ;;
<br>			-*)    # any other - argument
<br>				   errMsg "--- UNKNOWN OPTION ---"
<br>				   ;;
<br>			*)     # end of arguments
<br>				   break
<br>				   ;;
<br>		esac
<br>		shift   # next option
<br>	done
<br>	#
<br>	# get infile and outfile
<br>	infile=$1
<br>	outfile=$2
<br>fi
<br>
<br># test that infile provided
<br>[ "$infile" = "" ] && errMsg "NO INPUT FILE SPECIFIED"
<br>
<br># test that outfile provided
<br>[ "$outfile" = "" ] && errMsg "NO OUTPUT FILE SPECIFIED"
<br>
<br># setup temporaries
<br>tmpA1="$dir/spherize_A_$$.mpc"
<br>tmpA2="$dir/spherize_A_$$.cache"
<br>tmpB1="$dir/spherize_B_$$.mpc"
<br>tmpB2="$dir/spherize_B_$$.cache"
<br>tmpF1="$dir/spherize_T_$$.mpc"
<br>tmpF2="$dir/spherize_T_$$.cache"
<br>tmpM1="$dir/spherize_M_$$.mpc"
<br>tmpM2="$dir/spherize_M_$$.cache"
<br>tmpX1="$dir/spherize_X_$$.mpc"
<br>tmpX2="$dir/spherize_X_$$.cache"
<br>tmpY1="$dir/spherize_Y_$$.mpc"
<br>tmpY2="$dir/spherize_Y_$$.cache"
<br>trap "rm -f $tmpA1 $tmpA2 $tmpB1 $tmpB2 $tmpF1 $tmpF2 $tmpM1 $tmpM2 $tmpX1 $tmpX2 $tmpY1 $tmpY2; exit 0" 0
<br>trap "rm -f $tmpA1 $tmpA2 $tmpB1 $tmpB2 $tmpM1 $tmpM2 $tmpX1 $tmpX2 $tmpY1 $tmpY2; exit 1" 1 2 3 15
<br>
<br>
<br># read the input image into the TMP cached image.
<br>/bin/convert -quiet -regard-warnings "$infile" +repage "$tmpA1" ||
<br>	errMsg "--- FILE $infile NOT READABLE OR HAS ZERO SIZE ---"
<br>
<br># test for minimum IM version required
<br># IM 6.5.3.9 or higher to conform to new -compose distort
<br>im_version=`/bin/convert -list configure | \
<br>	sed '/^LIB_VERSION_NUMBER /!d; s//,/;  s/,/,0/g;  s/,0*\([0-9][0-9]\)/\1/g' | head -n 1`
<br>[ "$im_version" -lt "06050309" ] && errMsg "--- REQUIRES IM VERSION 6.5.3-9 OR HIGHER ---"
<br>
<br># get image dimensions
<br>ww=`/bin/convert $tmpA1 -ping -format "%w" info:`
<br>hh=`/bin/convert $tmpA1 -ping -format "%h" info:`
<br>
<br># set default diameters
<br>if [ "$dx" = "" ]; then
<br>	dx=$ww
<br>fi
<br>if [ "$dy" = "" ]; then
<br>	dy=$hh
<br>fi
<br>
<br># compute min diameter
<br>dm=`/bin/convert xc: -format "%[fx:min($dx,$dy)]" info:`
<br>if [ "$symmetry" = "yes" ]; then
<br>	dx=$dm
<br>	dy=$dm
<br>fi
<br>
<br># test diameters
<br>[ $dx -eq 0 -o $dx -gt $ww ] && errMsg "--- X DIAMETER MUST BE GREATER THAN 0 AND LESS OR EQUAL TO IMAGE WIDTH ---"
<br>[ $dy -eq 0 -o $dy -gt $hh ] && errMsg "--- Y DIAMETER MUST BE GREATER THAN 0 AND LESS OR EQUAL TO IMAGE HEIGHT ---"
<br>
<br>
<br># scale radius so that that image fits circle or ellipse at any diameter with zoom=1
<br># invert zoom so increases with larger values
<br># use absolute coord displacement
<br># xamount=(rx/zoom)*(ww/dx)=ww/(2*zoom); where rx = x radius
<br>xamount=`/bin/convert xc: -format "%[fx:$ww/(2*$zoom+quantumscale)]" info:`
<br>yamount=`/bin/convert xc: -format "%[fx:$hh/(2*$zoom+quantumscale)]" info:`
<br>
<br>
<br># create radial gradient
<br>/bin/convert -size ${dm}x${dm} radial-gradient: -negate $tmpF1
<br>
<br># convert circular gradient to elliptical
<br>if [ $dx -ne $dy ]; then
<br>	/bin/convert $tmpF1 -resize ${dx}x${dy}! $tmpF1
<br>fi
<br>
<br># create mask
<br>/bin/convert $tmpF1 -negate -gravity center -background black -extent ${ww}x${hh} -threshold 0 $tmpM1
<br>
<br># set up amplify for distortion
<br>if [ "$amp" = "1" ]; then
<br>	amplify=""
<br>else
<br>	amplify="-evaluate pow $amp"
<br>fi
<br>
<br># compute spherize function F=(2/pi)*asin(rd)/rd
<br>/bin/convert $tmpF1 \( +clone -function arcsin '2,0,2,0' $amplify \) \
<br>	-compose divide -composite $tmpF1
<br>
<br># set up x displacement
<br># compute x gradient increasing with i
<br>/bin/convert -size ${dy}x${dx} gradient: -rotate 90 $tmpX1
<br>
<br># compute u*v=xd*asin(rd)/rd; u=gradient 0,1; v=asin(rd)/rd
<br># but where want xd in range -1,1; so must use (2*u-1)
<br># then modify to 0.5*u*v+0.5 to recenter at mid gray
<br># old fx method: 0.5*((2*u-1)*v)+0.5 = u*v - 0.5*v + 0.5
<br># convert $tmpX1 $tmpF1 -monitor -fx "u*v-0.5*v+0.5" $tmpX1
<br>
<br>if [ "$im_version" -ge "06050403" ]; then
<br>	/bin/convert $tmpX1 $tmpF1 +swap \
<br>	-compose Mathematics -set option:compose:args "1,0,-0.5,0.5" \
<br>	-composite $tmpX1
<br>else
<br>
<br>	# multiply arcsin(r)/(r) by abs of gradient relative to midgray as zero
<br>	# abs is computed as -solarize 50% -level 50,0%
<br>	# abs has V shape and ranges 0 to 1
<br>	# then scale so that range is 0.5 to 1
<br>	#
<br>	# convert gradient to sign of gradient
<br>	# sign is computed as -threshold 50%
<br>	# thus zero (for neg) on left and one (for pos) on right
<br>	# use sign as mask to separate the right (pos) side,
<br>	# then negate mask and arcsin-absgrad product
<br>	# to get the left (neg) side and then combine with right side
<br>	/bin/convert $tmpX1 $tmpF1 \
<br>		\( -clone 0 -solarize 50% -level 50,0% \) \
<br>		\( -clone 1,2 -compose multiply -composite -function polynomial "0.5,0.5" \) \
<br>		\( -clone 0 -threshold 50% \) \
<br>		\( -clone 3,4 -compose multiply -composite \) \
<br>		\( -clone 3,4 -negate -compose multiply -composite \) \
<br>		-delete 0-4 -compose plus -composite \
<br>		$tmpX1
<br>fi
<br>
<br># y displacement
<br>if [ $dx -eq $dy ]; then
<br>	# map is symmetric so just rotate x displacement
<br>	/bin/convert $tmpX1 -rotate 90 $tmpY1
<br>else
<br>	# compute y gradient increasing with j
<br>	/bin/convert -size ${dx}x${dy} gradient: -negate $tmpY1
<br>
<br>	if [ "$im_version" -ge "06050403" ]; then
<br>		/bin/convert $tmpY1 $tmpF1 +swap \
<br>		-compose Mathematics -set option:compose:args "1,0,-0.5,0.5" \
<br>		-composite $tmpY1
<br>	else
<br>		/bin/convert $tmpY1 $tmpF1 \
<br>			\( -clone 0 -solarize 50% -level 50,0% \) \
<br>			\( -clone 1,2 -compose multiply -composite -function polynomial "0.5,0.5" \) \
<br>			\( -clone 0 -threshold 50% \) \
<br>			\( -clone 3,4 -compose multiply -composite \) \
<br>			\( -clone 3,4 -negate -compose multiply -composite \) \
<br>			-delete 0-4 -compose plus -composite \
<br>			$tmpY1
<br>	fi
<br>fi
<br>
<br># displace image (center relative)
<br># absolute coord displacement
<br>/bin/convert $tmpA1 $tmpX1 $tmpY1 \
<br>	-gravity center -compose distort \
<br>	-set option:compose:args ${xamount}x${yamount} -composite \
<br>	$tmpB1
<br>
<br># set up trim
<br>if [ "$trim" = "yes" ]; then
<br>	trimming="-gravity center -crop ${dx}x${dy}+0+0 +repage"
<br>else
<br>	trimming=""
<br>fi
<br>
<br># composite mask
<br>if [ "$bcolor" = "none" ]; then
<br>    # echo "xparent background"
<br>    /bin/convert $tmpB1 $tmpM1 \
<br>        -alpha off -compose copy_opacity -composite $trimming $outfile
<br>elif [ "$bcolor" = "image" -a "$symmetry" = "no" ]; then
<br>    # echo "img background w/sym"
<br>	/bin/convert $tmpB1 $trimming $outfile
<br>elif [ "$bcolor" = "image" -a "$symmetry" = "yes" ]; then
<br>    # echo "img background w/o sym"
<br>	/bin/convert $tmpA1 $tmpB1 $tmpM1 -compose over -composite $trimming $outfile
<br>else
<br>    # echo "background color $bcolor"
<br>	/bin/convert \( $tmpB1 $tmpM1 \
<br>		-gravity center -alpha off -compose copy_opacity -composite \) \
<br>		\( -size ${ww}x${hh} xc:"$bcolor" \) +swap \
<br>		-gravity center -compose over -composite $trimming $outfile
<br>fi
<br>
<br>exit 0
<br>