Changeset 922 for trunk/src/phys_ima.m

Timestamp:

Nov 24, 2015, 10:48:46 AM (8 years ago)

Author:

sommeria

Message:

corrections quentin

File:

• trunk/src/phys_ima.m (modified) (2 diffs)

trunk/src/phys_ima.m

@@ -40 +40 @@
 Rangy(1)=max(ycorner);
 test_multi=(max(npx)~=min(npx)) || (max(npy)~=min(npy)); %different image lengths
-npX=1+round((Rangx(2)-Rangx(1))/min(dx));% nbre of pixels in the new image (use the finest resolution min(dx) in the set of images)
-npY=1+round((Rangy(1)-Rangy(2))/min(dy));
+% npX=1+round((Rangx(2)-Rangx(1))/min(dx));% nbre of pixels in the new image (use the finest resolution min(dx) in the set of images)
+% npY=1+round((Rangy(1)-Rangy(2))/min(dy));
+
+npX=1+round((Rangx(2)-Rangx(1))/max(dx));% nbre of pixels in the new image (use the largest resolution max(dx) in the set of images)
+npY=1+round((Rangy(1)-Rangy(2))/max(dy));
+
+
 x=linspace(Rangx(1),Rangx(2),npX);
 y=linspace(Rangy(1),Rangy(2),npY);
@@ -69 +74 @@
 %            end
         end
-        xima=0.5:npx-0.5;%image coordinates of corners
-        yima=npy-0.5:-1:0.5;
+        xima=0.5:npx(icell)-0.5;%image coordinates of corners
+        yima=npy(icell)-0.5:-1:0.5;
         [XIMA_init,YIMA_init]=meshgrid(xima,yima);%grid of initial image in px coordinates
         [XIMA,YIMA]=px_XYZ(XmlData{icell}.GeometryCalib,X,Y,zphys);% image coordinates for each point in the real