Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  997 1628 (N36W01 L=184)   395     53        1 ( 963/1607-1043/1650) -14.0 E01-W03 N34-N38
      2 1713 1548 (N29W55 L=238)   409     72        1 (1696/1516-1727/1573)  53.1 W53-W57 N26-N30
      3  256 1287 (N17E54 L=129)  1552    231        1 ( 203/1203- 305/1363) -42.5 E60-E48 N11-N22
      4  932 1327 (N16E05 L=179)   296     92        0 ( 885/1314- 985/1343) -18.4 E07-E01 N15-N17
      5 1591 1188 (N06W37 L=220)   312     48        1 (1564/1164-1621/1212) -33.7 W34-W38 N04-N07
      6 1262 1020 (S04W14 L=197)   444     70        1 (1233/ 978-1283/1063)  78.7 W12-W15 S05-N00
      7  851  530 (S32E14 L=169)   313     54        0 ( 825/ 504- 876/ 559) -53.1 E15-E12 S34-S30
      8  724  534 (S32E23 L=160)   406     70        1 ( 678/ 514- 766/ 553)   0.0 E26-E20 S32-S30
  

    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