Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  688 1767 (N30E48 L=197)  1327    244        4 ( 631/1726- 760/1811)  -7.6 E55-E40 N27-N32
      2  928 1558 (N23E21 L=224)   852    219        3 ( 899/1465- 950/1653) -47.7 E27-E17 N17-N28
      3  217 1340 (S06E67 L=179)   568     99        0 ( 194/1298- 252/1379)  78.7 E69-E63 S08-S03
      4 1549 1307 (N26W24 L=269)   579    105        2 (1525/1250-1572/1365) -56.3 W21-W26 N22-N29
      5  858 1255 (N03E15 L=230)   712    123        3 ( 821/1198- 898/1313)  81.9 E17-E13 N00-N07
      6 1661 1145 (N20W35 L=281)   591    100        1 (1623/1092-1691/1192) -39.8 W32-W38 N17-N22
      7 1798  959 (N14W52 L=297)   701    139        2 (1778/ 901-1819/1045) -58.0 W48-W54 N10-N18
      8 1323  325 (S35W46 L=291)  1014    161        3 (1262/ 304-1393/ 351)  15.3 W40-W51 S38-S32
  

    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