From 791fabebf688d4d83fed3c28e3040ea5145ec097 Mon Sep 17 00:00:00 2001
From: Aslak Grinsted <ag+github@glaciology.net>
Date: Mon, 19 Sep 2016 18:41:38 +0200
Subject: [PATCH] Changed demo batura to use OC

Fix arrow orientation on batura example.
---
 demobatura.m | 5 +++--
 1 file changed, 3 insertions(+), 2 deletions(-)

diff --git a/demobatura.m b/demobatura.m
index 375ff0b..b75a6c7 100644
--- a/demobatura.m
+++ b/demobatura.m
@@ -17,6 +17,7 @@
 [A,x,y,Ia]=geoimread(fullfile(datafolder,'batura_2001.tif')); 
 [B,xb,yb,Ib]=geoimread(fullfile(datafolder,'batura_2002.tif')); 
 deltax=x(2)-x(1);%m/pixel
+deltay=y(2)-y(1);%m/pixel
 
 
 %make regular grid of points to track:
@@ -33,7 +34,7 @@
 mask=inpolygon(pu,pv,roi(:,1),roi(:,2));
 pu(~mask)=nan; %inserting nans at some locations will tell template match to skip these locations
 
-[du,dv,C,Cnoise,pu,pv]=templatematch(A,B,pu,pv,'showprogress',{'2001' '2002'});
+[du,dv,C,Cnoise,pu,pv]=templatematch(A,B,pu,pv,'showprogress',true,'method','oc');
 close all
 
 %visualize the results
@@ -44,7 +45,7 @@
 hold on
 signal2noise=C./Cnoise;
 keep=(signal2noise>2)&(C>.6);
-V=(du+dv*1i)*deltax; %m/yr
+V=(du*deltax)+(dv*1i)*deltay; %m/yr
 Vn=abs(V);
 alphawarp(px,py,Vn,.2+keep*.5)
 quiver(px(keep),py(keep),real(V(keep))./Vn(keep),imag(V(keep))./Vn(keep),0.2,'k') %arrows show direction.