Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1513 1647 (N25W52 L=210)  4697    748        6 (1268/1578-1722/1691) -24.0 W34-W72 N16-N33
      2 1599 1494 (N14W50 L=208)   693    145        0 (1540/1428-1661/1541) -60.9 W47-W52 N09-N18
      3  281 1361 (N43E47 L=111)   426     80        0 ( 265/1334- 298/1393) -51.3 E49-E45 N40-N45
      4  206  974 (N20E57 L=101)   855    157        0 ( 183/ 900- 235/1052)  52.1 E60-E52 N15-N24
      5 1836  823 (S30W50 L=208)  3347    452        7 (1701/ 668-1901/1006)  34.3 W36-W58 S36-S20
      6  871  975 (N06E10 L=148)   509    169        0 ( 790/ 930- 952/ 999)  -5.2 E15-E04 N03-N07
      7 1708  826 (S26W38 L=196)   680    177        0 (1653/ 751-1763/ 903)  28.6 W32-W44 S28-S22
      8 1085  804 (S09E02 L=156)   187     60        0 (1050/ 789-1122/ 818) -14.0 E04-W00 S09-S08
      9 1059  533 (S24E11 L=147)  1389    530        1 ( 916/ 356-1215/ 736)  28.1 E27-W03 S32-S16
     10  771  661 (S10E24 L=134)   367    123        1 ( 757/ 593- 790/ 711)  66.8 E25-E21 S13-S06
     11  482  483 (S13E50 L=108)   553    166        4 ( 450/ 438- 522/ 524)  24.0 E54-E45 S15-S11
     12  521  316 (S25E57 L=101)   553     88        0 ( 500/ 285- 539/ 353)  68.2 E59-E55 S27-S22
  

    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