Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1301 1808 (N53W17 L=262)   632     74        1 (1267/1770-1325/1841)  58.0 W12-W21 N49-N57
      2  579 1679 (N34E42 L=203)  1123    129        0 ( 553/1595- 607/1732) -54.0 E46-E37 N27-N38
      3 1534 1486 (N29W32 L=277)   358     62        0 (1487/1466-1565/1502)   8.1 W27-W34 N27-N30
      4  262 1427 (N15E61 L=184)   649    115        0 ( 239/1390- 289/1483) -37.9 E65-E56 N12-N19
      5 1752 1177 (N12W48 L=293)   529     54        0 (1719/1145-1778/1203) -11.3 W45-W50 N09-N13
      6  137  789 (S23E74 L=171)  1169    120        1 ( 112/ 741- 157/ 857) -74.1 E79-E70 S26-S19
      7 1726  643 (S20W56 L=301)   475     50        1 (1704/ 614-1758/ 668)  63.4 W53-W58 S21-S17
      8 1374  634 (S25W27 L=272)   292     63        1 (1349/ 605-1399/ 660)  56.3 W25-W28 S27-S23
      9 1553  358 (S41W58 L=303)  1153     86        1 (1509/ 315-1579/ 412) -26.6 W52-W65 S43-S36
     10  623  316 (S56E34 L=211)   396     62        0 ( 598/ 293- 639/ 348)  24.0 E39-E30 S58-S53
     11  560  301 (S58E45 L=200)   833     57        1 ( 534/ 268- 588/ 339)  26.6 E54-E38 S61-S54
  

    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