Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  372 1114 (S14E43 L=208)  1343    348        3 ( 286/1011- 480/1205)  51.8 E49-E36 S21-S07
      2  255 1133 (S16E54 L=197)   376     59        1 ( 239/1098- 275/1164)  90.0 E56-E53 S17-S13
      3 1868 1050 (N21W58 L=310)   493    138        0 (1832/1030-1889/1094)  51.3 W53-W60 N19-N24
      4  181  815 (S37E58 L=194)  1207    275        0 ( 133/ 777- 221/ 861)  32.5 E64-E53 S39-S32
      5  408  779 (S34E35 L=216)   670    144        0 ( 349/ 741- 470/ 823) -10.3 E40-E29 S35-S32
      6  805  761 (S25E06 L=246)  1735    462        0 ( 566/ 679- 999/ 812)   9.8 E22-W07 S28-S21
      7 1371  664 (S15W31 L=282)  1788    467        3 (1168/ 631-1582/ 715)  17.8 W16-W44 S18-S07
      8  601  535 (S45E16 L=236)  1131    308        3 ( 523/ 502- 731/ 584)  14.7 E23-E04 S48-S40
      9 1246  381 (S35W36 L=288)   310     42        1 (1218/ 365-1272/ 401)  56.3 W34-W37 S36-S33
     10 1410  279 (S36W55 L=307)   504     79        0 (1370/ 266-1450/ 295)  38.7 W51-W57 S37-S33
  

    OPENR, u, binfile, /COMP, /GET_LUN
    head = FLTARR(6) ;[xsize,ysize,radius,hh,mm,ss]
    READU, u, head
    data = INTARR(FIX(head[0]),FIX(head[1]))
    READU, u, data
    FREE_LUN, u
    

• each filament has the value of Fil_ID
• all values are in pixel

  c = WHERE(data EQ Fil_ID, area)
  f = '(2(A,I6))'
  PRINT, FORM=f, "Area of Fil_ID", Fil_ID," =", area 

  data = ROTATE(data,7)
  Solar_y = 2048-Solar_y