Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  492 1685 (N24E54 L= 40)  1109    173        3 ( 458/1625- 529/1754) -63.4 E57-E50 N19-N29
      2  233 1297 (S04E61 L= 33)  3933    524        9 ( 156/1086- 316/1507) -78.7 E71-E56 S16-N09
      3 1239 1300 (N13W07 L=101)   453    194        0 (1194/1233-1314/1384) -29.1 W02-W12 N10-N17
      4 1708 1308 (N26W39 L=133)   390     71        0 (1681/1262-1726/1343) -53.1 W36-W40 N23-N27
      5  740 1167 (S04E19 L= 75)   403    141        0 ( 684/1107- 786/1230)  69.4 E21-E17 S08-N00
      6 1802 1051 (N13W50 L=145)   643    136        0 (1775/ 972-1818/1117) -58.0 W46-W53 N08-N16
      7  283 1037 (S19E50 L= 44)   729    144        2 ( 258/ 956- 320/1099) -53.1 E53-E45 S23-S15
      8  722  842 (S23E15 L= 80)   231     86        0 ( 680/ 830- 772/ 858)  18.4 E17-E11 S24-S21
      9 1815  637 (S09W67 L=161)   646    247        4 (1775/ 607-1860/ 672)   4.1 W60-W74 S10-S05
     10 1543  369 (S31W56 L=150)   730    173        1 (1481/ 352-1609/ 390)  21.8 W50-W60 S33-S28
  

    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