Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1592 1433 (N34W38 L= 33)  1502    277        3 (1534/1310-1681/1533) -47.5 W34-W45 N26-N41
      2 1093 1401 (N25E00 L=355)   618    150        2 (1055/1340-1121/1471)  90.0 E02-W01 N20-N30
      3 1546 1237 (N20W33 L= 27)   285     98        1 (1515/1213-1586/1280) -31.0 W30-W35 N18-N22
      4  918 1228 (N12E09 L=346)   318     95        1 ( 871/1207- 962/1253)  31.0 E11-E06 N10-N14
      5  178  866 (S19E66 L=289)  3243    293        6 ( 140/ 731- 222/ 989)  86.6 E72-E59 S28-S11
      6 1674  883 (S00W45 L= 40)   322     84        1 (1650/ 843-1690/ 922)  68.2 W44-W46 S03-N02
      7 1703  715 (S10W52 L= 47)   922    217        1 (1623/ 659-1757/ 782)  45.0 W46-W56 S14-S05
      8  657  681 (S25E21 L=334)   317     66        0 ( 640/ 643- 678/ 714)  78.7 E21-E19 S27-S22
      9 1307  519 (S28W27 L= 22)   767    158        1 (1236/ 486-1360/ 571)  45.0 W22-W29 S30-S23
     10 1466  496 (S27W40 L= 35)   924    140        0 (1405/ 458-1514/ 545)  49.4 W36-W43 S30-S23
     11  735  362 (S48E15 L=340)  2496    406        5 ( 686/ 200- 787/ 513)  67.0 E27-E09 S68-S35
     12 1613  453 (S28W57 L= 52)  2244    304        6 (1480/ 401-1749/ 483)  10.3 W45-W80 S33-S25
  

    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