/////////////////////////////////////////////////////////////////////////////////////////////////
//                                                                                             //
//                                  Wildtrack Project - Lowland Tapir                          //
//                                                                                             //
// Authors: 
// Amendments: Cross hair mod by JP + point 26 adjusted                                                                                //
// Name: Sky     Date: Apr 10      Office: JMP Cary   JMP Version: 11                       //
// 1) Embedded crosshair function into Wildtrack JSL to improve accuracy of landmark points    //
// 2) Improve the funcatioality of the "New Picture" button to speed up process of collecting  //
//    landmark points from images                                                              //
// 3) Calculate specific areas of the foot (using area of polygons):                           //
//    Area of polygon = 1/2 ( x1(y2 − yn) + xn(y1 − yn−1) + x2(y3 − y1) + x3(y4 − y2) + ...)   //
// 4) Increase font size                                                                       //
// 5) Add Explorer functionality for saving images and table for the buttons "Save As",        //
//    "Save Table" and "New Picture"                                                           //
//                                                                                             //
/////////////////////////////////////////////////////////////////////////////////////////////////

/////////////////////////////////////////////////////////////////////////////////////////////////
//                                                                                             //
//                               Instructions for Operation                     	           //             
//											                                                   //
// 1. Run this script. The result will be a window with buttons and an empty graph             //
//											                                                   //
// 2. Drag a rhino spoor picture into the graph. On the Mac, you can just drag in a            //
//    jpeg file from the Finder. On Windows the only format that seems to work well            //
//    is BMP, so you may have to convert the JPEG to BMP before you drag it in.                //
//    The one application that holds a BMP picture that seems to drag well right               //
//    into the JMP graphics frame is Microsoft Word. Others don't work so well,                //
//    and the only thing that seems to work with those is to copy and paste into               //
//    the area below the frame, then use the fat-plus cursor to select, then drag              //
//    it into the frame. Soon we will provide a "Paste" right-click command in                 //
//    the Graphics frame to make this work.                                                    //
//                                                                                             //
// 3. If the graph is not the right size, then shift-drag from the corner to make              //
//    it the right size. IMPORTANT: the graph must be square, which you can double-check       //
//    by right-clicking in the graph but not on the picture and doing the Size/Scale           //
//    submenu command 'FrameSize' and verifying the coordinates are the same for horiz         //
//    and vert.                                                                                //
//																							   //
// 4. Click and Enter into the edit fields. Currently there are three edit fields,             //
//    labeled: Rhino, Track, Location.													       //
//																							   //
// 5. Click the landmark points in, in the order prompted by the message in the lower left.    //
//    If you click a point by mistake, you can click 'Delete Last Point' and click again.      //
//    Also, if the points do not line up perfectly, you can carefully click on a point and     //
//    drag it to a better position. You may wish to drag it away to get a better look at the   //
//    picture, then drag it back in. The messsage will say "Done" when all landmark points     //
//    are in.																				   //
//																							   //
// 6. Click 'Derived Points'. Now all the derived points will be shown with red labels.        //
//    If they don't look right, you can click "Clear Derived" to remove them, and refine       //
//    the original landmark points by dragging them. Then click "Derived Points" again.        //
//																							   //
// 7. When you are happy with the points, now you can add them to a data table.                //
//    If you don't already have the data table open, click 'Open' to open an old               //
//    data table, or 'New' to start a new table. Then to add the new row, click                //
//    the 'Append Row' button.                                                                 //
//                                                                                             //
// 8. Periodically, click Save or use save in the File menu to make sure the new               //
//    data makes it to disk.                                                                   //
//																							   //
// 9. For the next picture, click the 'New Picture' button to clear the landmarks,             //
//    right-click in the picture and do the "Delete" command. Then go to step 2.               //
//                                                                                             //
/////////////////////////////////////////////////////////////////////////////////////////////////
namesdefaulttohere(1);
// Instruction Image
//instructionImagePath = "C:\Users\JMP\Desktop\script images\am tiger footprint temp c.jpg";
instructionImagePath ="C:\Users\skyal\AppData\Local\SAS\JMP\Addins\com.sas.fitv1b\lowland tapir footprint b.jpg";
// Initial Point Set
// xx and yy are the placed markers in the graph
// xxd and yyd are the properly numbered points that include derived points
// xxa and yya are after normalizing the data to centimeters

getPictSize = function({pseg},										//Joseph ammend 14-01-15
	/*{i=0,tmp,blob=chartoblob(pseg<<getjournal),szblob=length(blob)},
	while(blobtochar(blobpeek(blob,i,7))!="picture"&i<szblob,i++);
	write("\!NgetPictSize: offset = ",i,", size blob = ",szblob);
	if(i>=szblob,
		empty(),
		tmp=blobtochar(blobpeek(blob,i+8,100));
		parse(substr(tmp,1,contains(tmp,")")))
	)*/	
	substitute(pseg<<getsize,{},expr(size()))
);
//
getRotationAngle = function({xx,yy},
	{angl},
	if(nrow(xx)==2,
		angl=calcAngle(xx[2],yy[1],xx[1],yy[1],xx[2],yy[2]); 
		if((xx[1]>xx[2]&yy[1]<yy[2])|(xx[1]<xx[2]&yy[1]>yy[2]),-1,1)*angl,
		0
	)
);

xx = [];
yy = [];
xxd = [];
yyd = [];

// xx=[144,138,362,686,314,666,232,350,608,736,238,670,354,448];
// yy=[682,214,676,636,586,514,546,490,442,452,368,314,146,170];

// Labels beside landmark points
pointLabel = {"Scale1","Scale2", // rule
   "1","2","3","4","5" ,"6","7","8","9","10","11","12","13","14","15",// placed
  "D16","D17","D18","D19","D20","D21","D22"}; // derived
numberOfLandmarksAndRules = 17;
numberOfDerivedPoints = 7;

LandmarkCircleColor = "white";
LandmarkTextColor = "white";
DerivedCircleColor = "yellow";
DerivedTextColor = "yellow";

// Labels for prompting
promptLabel = {"on ruler                                                           ","on ruler                                                           ",
  "Landmark 1=Highest point of front toe                               ",
  "Landmark 2=Highest point of toe 4                                   ",
  "Landmark 3=Highest point of toe 2                                   ",
  "Landmark 4=Intersection of L2-5 with left edge of front toe         ",
  "Landmark 5=Intersection of L2-5 with right edge of front toe        ",
  "Landmark 6=Most lateral point of left toe                           ",
  "Landmark 7=Most lateral point of right toe                          ",
  "Landmark 8=Intersection of L6-7 with left edge of front toe         ",
  "Landmark 9=Intersection of L6-7 with right edge of front toe        ",
  "Landmark 10=Most medial point of left toe                           ",
  "Landmark 11=Most medial point of right toe                          ",
  "Landmark 12=Lowest point of left toe                                ",
  "Landmark 13=Lowest point of right toe                               ",
  "Landmark 14=Lowest point of heel                                    ",
  "Landmark 15=Centre point of heel                                    ",
  "","","","","","","","","","","","","","",""
};

nLabel = nItems(pointLabel);

scaleFactor = 10;
depth = 0;

// Normalization of Data
CalcNormalize = expr(
   np = nrow(xx);

   //--rescale to centimeters, distance between 1 and 2 is 10 cm
   scale = scaleFactor/sqrt((yy[2]-yy[1])^2+(xx[2]-xx[1])^2);
   xxa = xxd*scale;
   yya = yyd*scale;

   //--center on 32nd point
   xxa = xxa-xxa[19];
   yya = yya-yya[19];
);




// function to calculate the intersection of two lines
// see http://astronomy.swin.edu.au/~pbourke/geometry/lineline2d/
CalcIntersection = function({p1x,p1y,p2x,p2y,p3x,p3y,p4x,p4y}, 
   tempDenom= (p4y-p3y)*(p2x-p1x)-(p4x-p3x)*(p2y-p1y);
   tempNumerA = (p4x-p3x)*(p1y-p3y)-(p4y-p3y)*(p1x-p3x);
   tempNumerB = (p2x-p1x)*(p1y-p3y)-(p2y-p1y)*(p1x-p3x);
   result = 
   if (tempDenom==0, {0,0},
     uA = tempNumerA / tempDenom;
     uB = tempNumerB / tempDenom;
     evalList({ p1x+uA*(p2x-p1x),
                p1y+uA*(p2y-p1y)}));
);


// test: CalcIntersection(1,0,0,1,0,0,1,1); should be {.5,.5}

// function to calculate distance between two points
calcDist = function({x1,y1,x2,y2},
  sqrt((x2-x1)^2+(y2-y1)^2);
);

// function to calculate the angle between three points
CalcAngle = function({x1,y1,x2,y2,x3,y3},
 L12 = calcDist(x1,y1,x2,y2);
 L13 = calcDist(x1,y1,x3,y3);
 L23 = calcDist(x2,y2,x3,y3);
 abs(arcCosine((L12^2+L23^2-L13^2)/(2*L12*L23)))*360/(2*pi());
);
// test: CalcAngle(0,0,1,0,0,1); // should be 45 degrees


CalcAngleIntersection = function({p1x,p1y,p2x,p2y,p3x,p3y,p4x,p4y}, 
 {tx0,ty0} = CalcIntersection(p1x,p1y,p2x,p2y,p3x,p3y,p4x,p4y);
 CalcAngle(p1x,p1y,tx0,ty0,p4x,p4y);
);


// The data on the graph is in xx and yy
// Here we reorganize those points into the numbered points,
// deriving new points into the vectors xxd and yyd

DerivePoints = expr(
 // Adjust for depth, if necessary
 if (depth > 0.0, 
  depthFactor = (19.469172 + 0.7021119 * depth) / 20;
  
  // Don't adjust first two points
  xx1 = xx[1]; yy1 = yy[1];
  xx2 = xx[2]; yy2 = yy[2];
  xx -= 500;
  yy -= 500;
  xx *= depthFactor;
  yy *= depthFactor;
  xx += 500;
  yy += 500;
  xx[1] = xx1; yy[1] = yy1;
  xx[2] = xx2; yy[2] = yy2;
 );
  
 nrx = nrow(xx);
 if (nrow(xx)<15 | nrow(yy)<15,print("NEED MORE POINTS",nrow(xx),nrow(yy)));
 xxx = xx[3::nrow(xx),0]; 
 yyy = yy[3::nrow(yy),0];
 xxd = j(22,1,0);
 xxd[1] = xxx[1];
 xxd[2] = xxx[2];
 xxd[3] = xxx[3];
 xxd[4] = xxx[4];
 xxd[5] = xxx[5];
 xxd[6] = xxx[6];
 xxd[7] = xxx[7];
 xxd[8] = xxx[8];
 xxd[9] = xxx[9];
 xxd[10] = xxx[10];
 xxd[11] = xxx[11];
 xxd[12] = xxx[12];
 xxd[13] = xxx[13];
 xxd[14] = xxx[14];
 xxd[15] = xxx[15];
 yyd = j(22,1,0);
 yyd[1] = yyy[1];
 yyd[2] = yyy[2];
 yyd[3] = yyy[3];
 yyd[4] = yyy[4];
 yyd[5] = yyy[5];
 yyd[6] = yyy[6];
 yyd[7] = yyy[7];
 yyd[8] = yyy[8];
 yyd[9] = yyy[9];
 yyd[10] = yyy[10];
 yyd[11] = yyy[11];
 yyd[12] = yyy[12];
 yyd[13] = yyy[13];
 yyd[14] = yyy[14];
 yyd[15] = yyy[15];
 temp = (yyd[2]-yyd[ 3])/(xxd[2]-xxd[ 3]);
 xTemp = xxd[1] + 1;
 yTemp = yyd[1] + 1*temp;
 {xxd[16],yyd[16]} = CalcIntersection(xxd[ 1],yyd[ 1],xTemp  ,yTemp  ,xxd[15],yyd[15],xxd[ 2],yyd[ 2]);
 temp = (yyd[2]-yyd[ 3])/(xxd[2]-xxd[ 3]);
 xTemp = xxd[1] + 1;
 yTemp = yyd[1] + 1*temp;
 {xxd[17],yyd[17]} = CalcIntersection(xxd[1],yyd[1],xTemp  ,yTemp  ,xxd[15],yyd[15],xxd[3],yyd[3]);
 {xxd[18],yyd[18]} = CalcIntersection(xxd[ 2],yyd[ 2],xxd[ 3],yyd[ 3],xxd[15],yyd[15],xxd[ 1],yyd[ 1]);
 {xxd[19],yyd[19]} = CalcIntersection(xxd[10],yyd[10],xxd[11],yyd[11],xxd[15],yyd[15],xxd[ 1],yyd[ 1]);
 {xxd[20],yyd[20]} = CalcIntersection(xxd[12],yyd[12],xxd[10],yyd[10],xxd[15],yyd[15],xxd[ 2],yyd[ 2]);
 {xxd[21],yyd[21]} = CalcIntersection(xxd[13],yyd[13],xxd[11],yyd[11],xxd[15],yyd[15],xxd[ 3],yyd[ 3]);
 {xxd[22],yyd[22]} = CalcIntersection(xxd[12],yyd[12],xxd[13],yyd[13],xxd[15],yyd[15],xxd[ 1],yyd[ 1]);
);
// Here are all the numeric column names for the output data table
colList = {
		"X1",	"X2",	"X3",	"X4",	"X5",
		"X6",	"X7",	"X8",	"X9",	"X10",
		"X11",	"X12",	"X13",	"X14",	"X15",
		"X16",	"X17",	"X18",	"X19",	"X20",
		"X21",	"X22",	
		"Y1",	"Y2",	"Y3",	"Y4",	"Y5",
		"Y6",	"Y7",	"Y8",	"Y9",	"Y10",
		"Y11",	"Y12",	"Y13",	"Y14",	"Y15",
		"Y16",	"Y17",	"Y18",	"Y19",	"Y20",
		"Y21",	"Y22",
		"V1" ,	"V2" ,	"V3" ,	"V4" ,	"V5" ,
		"V6" ,	"V7" ,	"V8" ,	"V9" ,	"V10",
		"V11",	"V12",	"V13",	"V14",	"V15",
		"V16",	"V17",	"V18",	"V19",	"V20",
		"V21",	"V22",	"V23",	"V24",	"V25",
		"V26",	"V27",	"V28",	"V29",	"V30",
		"V31",	"V32",	"V33",	"V34",	"V35",
		"V36",	"V37",	"V38",	"V39",	"V40",
		"V41",	"V42",	"V43",	"V44",	"V45",
		"V46",	"V47",	"V48",	"V49",	"V50",
		"V51",	"V52",	"V53",	"V54",	"V55",
		"V56",	"V57",	"V58",	"V59",	"V60",
		"V61",	"V62",	"V63",	"V64",	"V65",
		"V66",	"V67",	"V68",	"V69",	"V70",
		"V71",	"V72",	"V73",	"V74",	"V75",
		"V76",	"V77","V78",	"V79","V80",
		"V81",	"V82",	"V83",	"V84",	"V85",
		"V86",	"V87",	"V88",	"V89",	"V90","V91",
		"Area1","Area2","Area3","Area4","Area5","Area6"
		};

// add all the column names to the output table
addAllColumns = expr(
      idCol    = dtD<<NewColumn("Tapir",Character(25));
      trackCol = dtD<<NewColumn("Track",Character(24));
      lrfhCol  = dtD<<NewColumn("Trail",Character(24));
	  dateCol  = dtD<<NewColumn("Date",Character(24));
	  timeCol  = dtD<<NewColumn("Time",Character(24));
	  GPSCol  = dtD<<NewColumn("GPS",Character(25));
      show("---A---");
	nv = nitems(colList);
	for(i=1,i<=nv,i++, dtD<<NewColumn(colList[i],Numeric));
      show("----B---");
);

// This adds one row to the data table, setting many columns
addTheRow = expr(
	dtD<<AddRow({
   // ID Fields
   Tapir  = tb1<<GetText,
   Track  = tb2<<GetText,
   Trail   = tb3<<GetText,
   Date   = tb4<<GetText,
   Time   = tb5<<GetText,
   GPS   = tb6<<GetText,
   // Landmark and Derived Points
		X1=xxa[1],
	 	X2=xxa[2],
	 	X3=xxa[3],
	 	X4=xxa[4],
	 	X5=xxa[5],
	 	X6=xxa[6],
	 	X7=xxa[7],
	 	X8=xxa[8],
	 	X9=xxa[9],
	 	X10=xxa[10],
	 	X11=xxa[11],
	 	X12=xxa[12],
	 	X13=xxa[13],
	 	X14=xxa[14],
	 	X15=xxa[15],
	 	X16=xxa[16],
	 	X17=xxa[17],
	 	X18=xxa[18],
	 	X19=xxa[19],
	 	X20=xxa[20],
	 	X21=xxa[21],
	 	X22=xxa[22],
	 	Y1=yya[1],
	 	Y2=yya[2],
	 	Y3=yya[3],
	 	Y4=yya[4],
	 	Y5=yya[5],
	 	Y6=yya[6],
	 	Y7=yya[7],
	 	Y8=yya[8],
	 	Y9=yya[9],
	 	Y10=yya[10],
	 	Y11=yya[11],
	 	Y12=yya[12],
	 	Y13=yya[13],
	 	Y14=yya[14],
	 	Y15=yya[15],
	 	Y16=yya[16],
	 	Y17=yya[17],
	 	Y18=yya[18],
	 	Y19=yya[19],
	 	Y20=yya[20],
	 	Y21=yya[21],
  	 	Y22=yya[22],
	 	
   // Derived Values
		V1 = sqrt((:X2 -:X16 )^2+(:Y2 -:Y16 )^2),
		V2 = sqrt((:X2 -:X20)^2+(:Y2 -:Y20)^2),
		V3 = sqrt((:X20-:X15)^2+(:Y20-:Y15)^2),
		V4 = sqrt((:X1-:X18)^2+(:Y1-:Y18)^2),
		V5 = sqrt((:X18 -:X19 )^2+(:Y18 -:Y19 )^2),
		V6 = sqrt((:X19 -:X22 )^2+(:Y19 -:Y22 )^2),
		V7 = sqrt((:X22 -:X15)^2+(:Y22 -:Y15)^2),
		V8 = sqrt((:X17-:X3)^2+(:Y17-:Y3)^2),
		V9 = sqrt((:X3-:X21)^2+(:Y3-:Y21)^2),
		V10 = sqrt((:X21-:X15)^2+(:Y21-:Y15)^2),
		V11 = sqrt((:X2 -:X8 )^2+(:Y2 -:Y8 )^2),
		V12 = sqrt((:X3 -:X9)^2+(:Y3 -:Y9)^2),
		V13 = sqrt((:X2-:X10)^2+(:Y2-:Y10)^2),
		V14 = sqrt((:X3-:X11)^2+(:Y3-:Y11)^2),
		V15 = sqrt((:X10 -:X15 )^2+(:Y10 -:Y15 )^2),
		V16 = sqrt((:X11 -:X15 )^2+(:Y11 -:Y15 )^2),
		V17 = sqrt((:X10 -:X6 )^2+(:Y10 -:Y6 )^2),
		V18 = sqrt((:X11 -:X7 )^2+(:Y11 -:Y7 )^2),
		V19 = sqrt((:X2 -:X1 )^2+(:Y2 -:Y1 )^2),
		V20 = sqrt((:X1 -:X3 )^2+(:Y1 -:Y3 )^2),
		V21 = sqrt((:X6-:X1)^2+(:Y6-:Y1)^2),
		V22 = sqrt((:X1-:X7)^2+(:Y1-:Y7)^2),
		V23 = sqrt((:X12-:X1)^2+(:Y12-:Y1)^2),
		V24 = sqrt((:X1-:X13)^2+(:Y1-:Y13)^2),
		V25 = sqrt((:X12-:X20)^2+(:Y12-:Y20)^2),
		V26 = sqrt((:X20-:X13)^2+(:Y20-:Y13)^2),
		V27 = sqrt((:X13-:X21)^2+(:Y13-:Y21)^2),
		V28 = sqrt((:X21-:X20)^2+(:Y21-:Y20)^2),
		V29 = sqrt((:X21-:X12)^2+(:Y21-:Y12)^2),
		V30 = sqrt((:X20-:X10)^2+(:Y20-:Y10)^2),
		V31 = sqrt((:X21 -:X11)^2+(:Y21 -:Y11)^2),
		V32 = sqrt((:X2 -:X19)^2+(:Y2 -:Y19 )^2),
		V33 = sqrt((:X19 -:X3)^2+(:Y19 -:Y3 )^2),
		V34 = sqrt((:X8 -:X15)^2+(:Y8 -:Y15 )^2),
		V35 = sqrt((:X9 -:X15)^2+(:Y9 -:Y15 )^2),
		V36 = sqrt((:X14 -:X15)^2+(:Y14 -:Y15 )^2),
		V37 = sqrt((:X8 -:X18)^2+(:Y8 -:Y18 )^2),
		V38 = sqrt((:X9 -:X18)^2+(:Y9 -:Y18 )^2),
		V39 = sqrt((:X20 -:X22)^2+(:Y20 -:Y22 )^2),
		V40 = sqrt((:X21 -:X22)^2+(:Y21 -:Y22 )^2),
		V41 = sqrt((:X16 -:X1)^2+(:Y16 -:Y1 )^2),
		V42 = sqrt((:X1 -:X17)^2+(:Y1 -:Y17 )^2),
		V43 = sqrt((:X2 -:X4)^2+(:Y2 -:Y4 )^2),
		V44 = sqrt((:X4 -:X18)^2+(:Y4 -:Y18 )^2),
		V45 = sqrt((:X18 -:X5)^2+(:Y18 -:Y5)^2),
		V46 = sqrt((:X5 -:X3)^2+(:Y5 -:Y3 )^2),
		V47 = sqrt((:X6 -:X8)^2+(:Y6 -:Y8 )^2),
		V48 = sqrt((:X8 -:X9)^2+(:Y8 -:Y9 )^2),
		V49 = sqrt((:X9 -:X7)^2+(:Y9 -:Y7 )^2),
		V50 = sqrt((:X10 -:X19)^2+(:Y10 -:Y19 )^2),
		V51 = sqrt((:X19 -:X11)^2+(:Y19 -:Y11 )^2),
		V52 = sqrt((:X12 -:X22)^2+(:Y12 -:Y22)^2),
		V53 = sqrt((:X22 -:X13)^2+(:Y22 -:Y13 )^2),
		V54 = sqrt((:X14 -:X15)^2+(:Y14 -:Y15 )^2),
		V55 = CalcAngle(:X16,:Y16,:X15,:Y15,:X1,:Y1),
		V56 = CalcAngle(:X17,:Y17,:X15,:Y15,:X1,:Y1),
		V57 = CalcAngle(:X6,:Y6,:X15,:Y15,:X1,:Y1),
		V58 = CalcAngle(:X7,:Y7,:X15,:Y15,:X1,:Y1),
		V59 = CalcAngle(:X12,:Y12,:X15,:Y15,:X1,:Y1),
		V60 = CalcAngle(:X13,:Y13,:X15,:Y15,:X1,:Y1),
		V61 = CalcAngle(:X10,:Y10,:X15,:Y15,:X1,:Y1),
		V62 = CalcAngle(:X11,:Y11,:X15,:Y15,:X1,:Y1),
		V63 = CalcAngle(:X16,:Y16,:X22,:Y22,:X1,:Y1),
		V64 = CalcAngle(:X17,:Y17,:X22,:Y22,:X1,:Y1),
		V65 = CalcAngle(:X6,:Y6,:X22,:Y22,:X1,:Y1),
		V66 = CalcAngle(:X7,:Y7,:X22,:Y22,:X1,:Y1),
		V67 = CalcAngle(:X12,:Y12,:X22,:Y22,:X1,:Y1),
		V68 = CalcAngle(:X13,:Y13,:X22,:Y22,:X1,:Y1),
		V69 = CalcAngle(:X10,:Y10,:X22,:Y22,:X1,:Y1),
		V70 = CalcAngle(:X11,:Y11,:X22,:Y22,:X1,:Y1),
		V71 = CalcAngle(:X2,:Y2,:X1,:Y1,:X22,:Y22),
		V72 = CalcAngle(:X3,:Y3,:X1,:Y1,:X22,:Y22),
		V73 = CalcAngle(:X4,:Y4,:X1,:Y1,:X22,:Y22),
		V74 = CalcAngle(:X5,:Y5,:X1,:Y1,:X22,:Y22),
		V75 = CalcAngle(:X6,:Y6,:X1,:Y1,:X22,:Y22),
		V76 = CalcAngle(:X7,:Y7,:X1,:Y1,:X22,:Y22),
		V77 = CalcAngle(:X12,:Y12,:X1,:Y1,:X22,:Y22),
		V78 = CalcAngle(:X13,:Y13,:X1,:Y1,:X22,:Y22),
		V79 = CalcAngle(:X2,:Y2,:X10,:Y10,:X4,:Y4),
		V80 = CalcAngle(:X3,:Y3,:X11,:Y11,:X5,:Y5),
		V81 = CalcAngle(:X4,:Y4,:X1,:Y1,:X18,:Y18),
		V82 = CalcAngle(:X5,:Y5,:X1,:Y1,:X18,:Y18),
		V83 = CalcAngle(:X8,:Y8,:X1,:Y1,:X22,:Y22),
		V84 = CalcAngle(:X9,:Y9,:X1,:Y1,:X22,:Y22),
		V85 = CalcAngle(:X12,:Y12,:X10,:Y10,:X22,:Y22),
		V86 = CalcAngle(:X13,:Y13,:X11,:Y11,:X22,:Y22),
		V87 = CalcAngle(:X1,:Y1,:X16,:Y16,:X2,:Y2),
		V88 = CalcAngle(:X1,:Y1,:X17,:Y17,:X3,:Y3),
		V89 = CalcAngle(:X1,:Y1,:X8,:Y8,:X7,:Y7),
 		V90 = CalcAngle(:X1,:Y1,:X9,:Y9,:X6,:Y6),
 		V91 = sqrt((:X6 -:X7)^2+(:Y6 -:Y7 )^2),
 		Area1=(-1/2)*( xxa[ 1]*(yya[ 3] − yya[ 2]) + xxa[ 3]*(yya[ 7] − yya[ 1]) + xxa[ 7]*(yya[13] − yya[ 3]) 
		             + xxa[13]*(yya[15] − yya[ 7]) + xxa[15]*(yya[14] − yya[13]) + xxa[14]*(yya[12] − yya[15])
		             + xxa[12]*(yya[ 6] − yya[14]) + xxa[ 6]*(yya[ 2] − yya[12]) + xxa[ 2]*(yya[ 1] − yya[ 6])
		             ),
		             
		Area2=(-1/2)*( xxa[ 1]*(yya[ 3] − yya[ 2]) + xxa[ 3]*(yya[ 7] − yya[ 1]) + xxa[ 7]*(yya[13] − yya[ 3]) 
		             + xxa[13]*(yya[12] − yya[ 7])
		             + xxa[12]*(yya[ 6] − yya[13]) + xxa[ 6]*(yya[ 2] − yya[12]) + xxa[ 2]*(yya[ 1] − yya[ 6])
		             ),
		             
		Area3=(-1/2)*( xxa[ 1]*(yya[ 3] − yya[ 2]) + xxa[ 3]*(yya[ 7] − yya[ 1]) + xxa[ 7]*(yya[6] − yya[ 3]) 
		             + xxa[ 6]*(yya[ 2] − yya[7]) + xxa[ 2]*(yya[ 1] − yya[ 6])
		             ),
		             
		Area4=(-1/2)*( xxa[ 1]*(yya[ 5] − yya[ 4]) + xxa[ 5]*(yya[ 9] − yya[ 1]) + xxa[ 9]*(yya[19] − yya[ 5]) 
		             + xxa[19]*(yya[8] − yya[ 9]) + xxa[8]*(yya[4] − yya[19]) + xxa[4]*(yya[1] − yya[8])
		             ),
		             
		Area5=(-1/2)*( xxa[ 2]*(yya[ 10] − yya[ 6]) + xxa[ 10]*(yya[ 12] − yya[ 2]) + xxa[ 12]*(yya[6] − yya[ 10]) 
		             + xxa[6]*(yya[2] − yya[12])
		             ),
		             
		Area6=(-1/2)*( xxa[ 3]*(yya[ 7] − yya[11]) + xxa[ 7]*(yya[13] − yya[ 3]) + xxa[ 13]*(yya[11] − yya[ 7]) 
		             + xxa[11]*(yya[ 3] − yya[13])
		             )              
		                                       
		                                       			
	});
);

if( isEmpty(initData),initData = 0);

// function to output data to a row in a data table
AddRowToTable = expr(

   // the first time through, we need to begin the data table initData=0
   if (initData==0,
       dtD = newTable("track");
       addAllColumns;
       initData=1);

   // add the row and populate with data
   addTheRow;
);

// dtD = Open();

// This clears data to start a new picture
doNewPicture = expr(
     xx = []; yy = []; xxd=[]; yyd=[];
);

//Define landmark points using crosshair
dragScript = Expr(
				xPos = x;
				yPos = y	
				);
mouseUpScript = Expr(
					Print("At mouse up: x = "||Char(Round(x,2))||"  y = "||Char(Round(y,2)));
					xx = xx |/ Matrix(xPos);
                    yy = yy |/ Matrix(yPos);
					// Looks like it's possible to drag the handle outside the range defined
					// by XScale and YScale?? If this happens, just reset the handle position 
					// to (0,0) rather than leave it "outside" the GraphBox
					if(xPos > 1000 |xPos < 0 |yPos > 1000 |yPos < 0, xPos = 500; yPos = 500);

					);
savepic=EXPR(

selectedDir = PickDirectory();

result = dialog(vList("Enter Path",pathName=EditText(selectedDir||"track1.jpg",width(400)),Button("OK")));
path = result["pathName"];
wind[FrameBox(1)]<<Save Picture(path,JPEG);
	
); // savepic close


// Graphics Window with Script

wdw=expr(
wind = NewWindow("WILDTRACK - FIT : Feature Extraction",
//image1 = new image( "C:\Users\JMP\Desktop\script images\amur tiger zoe h.jpg" );
image1 = try(new image( "C:\Users\skyal\AppData\Local\SAS\JMP\Addins\com.sas.fitv1b\tapir wide.jpg" ),newimage(j(200,200,1)));
//image2 = new image( "C:\Users\JMP\Desktop\script images\amur tiger zoe a.jpg" );
image2 = try(new image( "C:\Users\skyal\AppData\Local\SAS\JMP\Addins\com.sas.fitv1b\amur tiger y.jpg" ),newimage(j(200,200,1)));
hsplitterbox(
//hlistBox(
  vlistBox(
  	Lineup Box( N Col( 3 ),	spacing( 5 ),
  image1),
  //image2),
  Lineup Box( N Col( 2 ),	spacing( 22 ),
	Text box ("                Lowland tapir",<<Set Wrap( 500 ), <<SETFONTSTYLE("bold"),<<setfontsize(14),<<Justify Text("center")),
	Text box ("(Tapirus terrestris) ",<<Set Wrap( 500 ), <<SETFONTSTYLE("italic"),<<setfontsize(14),<<Justify Text("center"))),
	Lineup Box( N Col( 1 ),	spacing( 20 ),
	Text box ("       Feature extraction template for left hind footrint image",<<Set Wrap( 500 ),<<SETFONTSTYLE(""),<<setfontsize(11),<<Justify Text("center"))),
     spacerbox(size(1,5)),
     Lineup Box( N Col( 1 ),	spacing( 20 ),
  	 try(newImage(instructionImagePath),newimage(j(200,200,1)))),
  ),
 vlistBox(
 // Fields for Entering ID data
	lineupbox(ncol(1),
									lineupbox(ncol(8),spacing(0),
										buttonbox("Relaunch Analysis",<<setfunction(function({t},include(_fitdir||"Lowland tapir JPscript.jsl");((t<<topparent)<<closewindow)))),
										buttonbox("Pairwise Data Analysis",<<setfunction(function({t},include(_fitdir||"Lowland tapir JMautoscript.jsl");((t<<topparent)<<closewindow)))),
										//buttonbox("Mapping",<<setfunction(function({t},mainmenu("Mapping");_fitnw<<closewindow))),
										buttonbox("Validated Discriminant Analysis",<<setfunction(function({t},include(_fitdir||"DiscValid Lowland tapir.jsl");((t<<topparent)<<closewindow)))),
										buttonbox("Mapping",<<setfunction(function({t},include(_fitdir||"DanHomeRange Lowland tapir.jsl");((t<<topparent)<<closewindow)))),
										buttonbox("Main Menu",<<setfunction(function({t},include(_fitdir||"FIT LOWLAND TAPIR.jsl");((t)<<topparent)<<closewindow))),
										spacerbox(size(70,1)),	
										buttonbox("Help",web( "file://" || _addinHelpPath || "Add-in Help.pdf")),
	),
	spacerbox(size(1,2)),
									spacerbox(size(1,2),color(gray))
	),								
	lineupbox(ncol(1),
	Lineup Box( N Col( 130 ), spacing( 1 ),
	text box ("                                        Image drag and drop window for feature extraction",<<Set Wrap(700),<<SETFONTSTYLE("bold"),<<setfontsize(12),<<Justify Text(center))),
	spacerbox(size(1,2)),
	spacerbox(size(1,2),color(gray))
	),
	
 hlistBox(
   
   TextBox("  Tapir "), tb1=TextEditBox("                        "),
   TextBox("  Track "), tb2=TextEditBox("                       "),
   TextBox("  Trail "), tb3=TextEditBox("                       "),
   TextBox("  Date "), tb4=TextEditBox("                       "),
   TextBox("  Time "), tb5=TextEditBox("                       "),
   TextBox("  GPS "), tb6=TextEditBox("                        "),
),
 // Buttons for the landmarking work
 Lineup Box( N Col( 1 ), spacing( 5 ),
 hlistBox(
 
 //hlistbox(
	//vlistbox(
		buttonbox("Resize",
			<<setfunction(
				function({t},
					{scl,tmp, width, height, ref=(t<<topparent)[framebox(1)], pseg=ref<<findseg("Pict Seg")},
					_tm=tickseconds();tmp=getPictSize(pseg);write("\!t time = ",tickseconds()-_tm);
					tmp=substitute(nameexpr(tmp),expr(size),{});
					scl=max(tmp[1],tmp[2])/1000;
					eval(substitute(expr(pseg<<specifysize(0,t,r,0)),expr(t),tmp[2]/scl,expr(r),tmp[1]/scl)); 
					ref<<inval;	pseg<<lock(1)
				)
			)
		),
		buttonbox("Rotate",
			<<setfunction(
				function({t},
					{angl,seg,ref=t<<topparent},
					seg=ref[framebox(1)]<<findseg("Pict Seg");
					angl=getRotationAngle(xx,yy);
					write("\!NIn Rotate!!!!\!Nangl = ",angl);
					eval(substitute(expr(seg<<rotate(angl)),expr(angl),angl));
					ref<<inval;
					xx=yy=[];
					ref[framebox(1)]<<removegraphicsscript(1);
					ref[framebox(1)]<<addgraphicsscript(
						PenColor("Green");
						VLine(xPos); HLine(yPos);
						Handle( xPos, yPos, dragScript, mouseUpScript );
						textsize(8); TextColor("black"); Backcolor("white");
						Text(Erased, {750, 980}, "xPos ", Round(xPos,2), "   yPos ", Round(yPos,2) );
						//--- Drawing Script ---
						// show derived points
						nDerived = nrow(yyd);
						if(nDerived,
							textColor(DerivedTextColor);pencolor(DerivedCircleColor); // 0 is black, 3 is red, 2 is white
							for(i=numberOfLandmarksAndRules-2+1,i<=nDerived,i++,
								circle({xxd[i],yyd[i]},6);text({xxd[i]+10,yyd[i]+5},char(i))
							)
						);
						// show landmark points
						pencolor(5); // the tiny draggable landmark point formed by DragMarker.  Any pencolor seems to make it black; no pencolor seems to make it white.
						DragMarker(xx,yy);
						nLandmark = nrow(xx);
						textColor(LandmarkTextColor);pencolor(LandmarkCircleColor);
						for(i=1,i<=nLandmark,i++,   // draw circles around landmark points
							circle({xx[i],yy[i]},10); text({xx[i]+10,yy[i]+5},pointLabel[i])
						);
						textsize(8); textColor("black");
						MouseTrap({},xx=xx|/matrix({x});yy=yy|/matrix({y}));
						msg1 = if(nrow(xx)<numberOfLandmarksAndRules,pointLabel[nrow(xx)+1]||" "||promptLabel[nrow(xx)+1],"DONE");
						Text(Erased,{5,980},msg1);
					)
				)
			)
		),
	//),
	
//)  // End of hlistbox
	
   spacerbox(size(55,1)),
   ButtonBox("Delete Last Point", xx = xx[1::nrow(xx)-1,0]; yy = yy[1::nrow(yy)-1,0]; wind<<reshow; ),
   spacerbox(size(35,1)),
   ButtonBox("Derived Points",  DerivePoints;  wind<<reshow;  ),
   ButtonBox("Clear Derived",  xxd=[]; yyd=[]; wind<<reshow;  ),
   spacerbox(size(55,1)),
   ButtonBox("Open Table",dtD = Open();initData=1;),
   ButtonBox("New Table", dtD = newTable("trackABC.jmp"); addAllColumns; initData=1; ),
   spacerbox(size(50,1)),
   
),
),
// Buttons for Data Table Setup
Lineup Box( N Col( 1 ), spacing( 5 ),
 hlistBox(
  // ButtonBox("Open Table",dtD = Open();initData=1;),
  // ButtonBox("New Table", dtD = newTable("trackABC.jmp"); addAllColumns; initData=1; ),


  Lineup Box( N Col( 3 ),	spacing( 0 ),
 hlistbox(
 	TextBox("SCALE FACTOR"), scaleFactorText = TextEditBox (scaleFactor, << setwidth(20),<<SETFONTSTYLE("bold"),<< setScript(scaleFactor = num(scaleFactorText << getText))), Textbox("cm"), 
 	TextBox("   SUBSTRATE DEPTH"), substrateDepthText = TextEditBox(depth, << setwidth(20),<<SETFONTSTYLE("bold"),<< setScript(depth = num(substrateDepthText << getText))), Textbox("cm"))
	),
  spacerbox(size(32,1)),
  ButtonBox("Append Row",  
		CalcNormalize; 
		local({angl,rot,tmp},
			angl=getRotationAngle(xxa,yya);
			rot=(cos(angl)||-sin(angl))|/(sin(angl)||cos(angl));
			show(xxa,yya,angl);
			tmp=rot*(xxa`|/yya`);
			xxa=tmp[1,0]`;yya=tmp[2,0]`;
			show(xxa,yya)
		);
		AddRowToTable; 
	),
   spacerbox(size(32,1)),
   ButtonBox("Save Table",  
   
   selectedDir = PickDirectory();

result = dialog(vList("Enter Path",pathName=EditText(selectedDir||"table.jmp",width(400)),Button("OK")));
path = result["pathName"];
dtD<<Save(path);
 
   dtD<<Save; ),
   ButtonBox("Save Picture", 
   savepic;
   ),
   
   spacerbox(size(32,1)),
   ButtonBox("New Picture",  
   
   savepic;
   doNewPicture; wind<<reshow;
   wind<<close window;
   wdw;
             ),
   
 
	//buttonbox("Help",web( "file://" || _addinHelpPath || "Add-in Help.pdf"))																									
 ),
 ),
xPos = 500;
yPos = 500;
xx = [];
yy = [];
//hlistbox(

	GraphBox(
		FrameSize(600,600),
		X Scale(0,1000),
		y Scale(0,1000),
		DoubleBuffer,
		textsize(8); textColor("black"); Backcolor("white");
		//Text(Erased,{5,980},"Rotation points");
		Text(Erased,{5,980},"Drag & Drop, Resize, Rotate with Rotation Points");
		pen color( "Red" );markersize(3);
		dragMarker(xx,yy);
		for(i=1,i<=nrow(xx),i++,circle({xx[i],yy[i]},20));
		mouseTrap({},if(nrow(xx)<2,xx=xx|/matrix({x});yy=yy|/matrix({y}))), <<setbackgroundcolor ("gray")
	)  // End of GraphBox
//)  // End of hlistbox
)  // End of vlistbox
)  // End of splitterbox
); // End of window

);

wdw;