Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1478 1435 (N28W36 L=301)  1307    265        1 (1422/1362-1576/1535) -63.4 W31-W41 N22-N34
      2  272 1429 (N32E66 L=199)   772     92        1 ( 236/1399- 306/1457)  -7.1 E75-E59 N29-N33
      3  220 1330 (N25E68 L=197)   533     73        1 ( 189/1318- 258/1349)   4.8 E74-E61 N24-N25
      4 1048 1306 (N22W03 L=268)   216     77        1 (1011/1294-1088/1324) -11.3 W00-W05 N21-N23
      5  311 1211 (N18E52 L=213)   320     70        1 ( 294/1172- 323/1252) -51.3 E54-E50 N15-N20
      6  736 1158 (N14E18 L=247)   289    115        0 ( 704/1109- 779/1182)  59.0 E19-E14 N10-N15
      7 1492  855 (S09W30 L=295)   509    142        1 (1470/ 801-1526/ 920)  60.3 W27-W32 S11-S04
      8  800  895 (S02E15 L=250)   265     79        0 ( 758/ 878- 834/ 919)  11.3 E17-E12 S03-S01
      9 1754  860 (S11W52 L=317)   583     55        2 (1730/ 831-1777/ 890)  63.4 W49-W54 S12-S08
     10 1756  786 (S16W53 L=318)   440    115        0 (1714/ 739-1787/ 833)  50.2 W48-W55 S18-S12
     11  433  579 (S22E46 L=219)   730    146        1 ( 412/ 518- 460/ 652)  71.6 E49-E43 S26-S17
     12  293  565 (S24E64 L=201)   786    101        1 ( 261/ 528- 329/ 598)  12.1 E70-E56 S26-S21
  

    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