Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  692 1738 (N39E36 L=172)   400     30        0 ( 674/1711- 711/1761) -38.7 E39-E33 N36-N40
      2 1451 1543 (N32W25 L=233)   244     49        0 (1433/1512-1467/1570) -63.4 W23-W25 N30-N34
      3  713 1434 (N17E24 L=184)  1142    253        0 ( 670/1330- 802/1531) -45.0 E28-E16 N10-N23
      4 1794 1244 (N18W53 L=261)   401     30        0 (1769/1225-1814/1266)  31.0 W49-W54 N16-N19
      5 1653  919 (S04W42 L=250)   270     38        0 (1636/ 895-1674/ 947) -71.6 W40-W43 S05-S02
      6  315  603 (S37E55 L=153)   703    148        1 ( 267/ 544- 367/ 662) -53.1 E59-E51 S40-S32
      7 1661  444 (S32W63 L=271)   491     72        0 (1627/ 409-1697/ 476)  68.2 W60-W65 S34-S29
      8  816  292 (S57E10 L=199)  2222    367        6 ( 668/ 263- 951/ 325)   3.4 E29-W05 S60-S53
      9 1029  290 (S54W13 L=221)   452     45        1 ( 995/ 270-1059/ 311)  38.7 W09-W15 S56-S52
  

    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