Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  639 1796 (N73E25 L= 74)   505     51        0 ( 621/1767- 659/1822)  90.0 E29-E21 N69-N76
      2  817 1609 (N49E01 L= 98)   729    107        2 ( 768/1556- 867/1638) -29.7 E04-W03 N44-N51
      3  823 1325 (N29E07 L= 92)   427    126        0 ( 789/1278- 879/1375)  26.6 E09-E01 N26-N31
      4  428 1152 (N25E38 L= 61)   929    239        1 ( 330/1076- 492/1215) -35.5 E47-E33 N19-N30
      5  245 1100 (N24E57 L= 42)  1333    289        2 ( 180/ 998- 318/1208)  36.9 E65-E48 N18-N30
      6 1029  466 (S28E12 L= 87)   293     58        0 (1002/ 441-1059/ 493)  21.8 E14-E09 S29-S25
      7  494  431 (S20E53 L= 47)   547     70        1 ( 472/ 401- 527/ 454)   0.0 E55-E49 S22-S18
      8  629  363 (S27E45 L= 54)  1264    233        1 ( 573/ 273- 709/ 445)  82.4 E47-E41 S35-S20
  

    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