Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1691 1479 (N11W57 L=258)   541    122        0 (1661/1444-1727/1509)  -7.1 W53-W61 N09-N12
      2 1690 1404 (N07W53 L=254)   337     50        1 (1674/1372-1706/1436)  76.0 W51-W54 N05-N09
      3  996 1206 (N15W03 L=204)   767    132        3 ( 932/1182-1080/1254) -49.4 W00-W07 N11-N18
      4  457 1158 (N25E32 L=169)  1124    251        0 ( 357/1096- 556/1207) -40.2 E39-E26 N19-N30
      5 1019  636 (S17E10 L=191)  1166    230        0 ( 967/ 522-1059/ 735)  69.0 E14-E07 S23-S10
      6 1693  503 (S48W37 L=239)  1253    226        0 (1672/ 434-1714/ 574) -57.3 W34-W43 S55-S41
      7  960  480 (S26E19 L=182)   527     80        1 ( 922/ 446-1006/ 505)   9.5 E21-E15 S27-S23
      8  472  387 (S20E56 L=145)   430     76        1 ( 451/ 369- 501/ 412) -56.3 E58-E54 S21-S18
      9 1326  263 (S55E03 L=198)   718    119        2 (1278/ 244-1373/ 288) -45.0 E06-W01 S58-S51
  

    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