Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1141 1779 (N44E19 L= 17)  2813    443        3 (1018/1641-1293/1927) -19.2 E45-E02 N37-N54
      2  525 1762 (N25E64 L=332)   925     95        3 ( 487/1730- 568/1791)  38.7 E67-E60 N22-N27
      3  513 1692 (N20E57 L=339)  2542    244        4 ( 435/1617- 618/1775)  90.0 E63-E54 N13-N27
      4  883 1719 (N30E32 L=  5)   389     70        1 ( 865/1678- 902/1759) -56.3 E34-E29 N27-N33
      5  601 1684 (N21E49 L=348)  1019    161        2 ( 542/1627- 653/1735)  77.5 E52-E46 N16-N25
      6  593 1512 (N10E40 L=357)  1548    420        3 ( 464/1374- 723/1663)  74.1 E44-E37 N00-N21
      7 1177 1483 (N23E04 L= 32)   216     72        0 (1157/1454-1202/1523) -36.9 E06-E01 N21-N25
      8 1728 1441 (N41W39 L= 75)   478     46        1 (1700/1416-1751/1467) -36.9 W35-W41 N39-N42
      9  399 1069 (S20E41 L=356)   362     77        0 ( 369/1040- 422/1103)  53.1 E42-E38 S21-S17
     10  334 1036 (S23E46 L=351)   666    171        1 ( 299/ 944- 359/1098) -48.8 E49-E42 S28-S19
     11  485 1005 (S22E33 L=  4)   305     57        0 ( 447/ 989- 517/1025)  31.0 E35-E30 S23-S19
     12 1253  237 (S44W54 L= 91)   516     98        0 (1216/ 210-1301/ 267)   4.4 W48-W61 S45-S43
  

    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