Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  583 1563 (N32E40 L=209)  1229    217        3 ( 531/1466- 617/1629) -42.5 E48-E34 N25-N36
      2 1592 1494 (N37W42 L=291)  1185    262        2 (1522/1373-1661/1567) -65.2 W36-W46 N28-N41
      3  627 1429 (N24E31 L=217)   319    122        0 ( 585/1410- 689/1457)   0.0 E35-E26 N22-N25
      4 1495  927 (S01W31 L=280)   644    149        1 (1431/ 906-1577/ 962)   5.2 W26-W37 S02-N01
      5 1821  927 (N01W60 L=308)   565     58        0 (1797/ 900-1842/ 958) -56.3 W57-W62 N00-N03
      6 1749  854 (S04W53 L=302)   594     64        1 (1723/ 820-1776/ 884) -63.4 W50-W55 S05-S01
      7 1821  776 (S08W63 L=312)   546     59        0 (1804/ 743-1840/ 810)  78.7 W60-W65 S10-S05
      8 1420  723 (S14W28 L=277)   942    276        2 (1314/ 656-1522/ 777) -16.4 W20-W37 S17-S10
      9 1227  687 (S18W16 L=265)   498    125        1 (1164/ 658-1275/ 718)  14.0 W11-W19 S19-S15
     10  623  657 (S25E24 L=224)   235     70        0 ( 585/ 647- 662/ 666)   0.0 E27-E21 S25-S24
     11 1376  547 (S26W29 L=278)   254     61        0 (1362/ 516-1391/ 585) -76.0 W28-W30 S27-S23
     12 1443  527 (S27W35 L=284)   259     84        0 (1415/ 493-1464/ 570)  80.5 W33-W36 S29-S23
  

    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