Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  731 1721 (N60E12 L= 23)   339     43        2 ( 708/1697- 750/1745)  68.2 E15-E09 N57-N62
      2  639 1683 (N58E24 L= 11)   768     78        4 ( 610/1642- 669/1715)  90.0 E27-E19 N54-N61
      3  544 1647 (N57E36 L=359)  1208    126        1 ( 498/1591- 596/1688)  54.5 E43-E29 N52-N59
      4  982 1598 (N44W09 L= 44)   651    140        2 ( 914/1573-1049/1621) -26.6 W04-W14 N41-N46
      5 1194 1576 (N38W25 L= 60)   612    118        1 (1144/1555-1256/1594)  -5.7 W20-W30 N36-N39
      6 1476 1451 (N23W40 L= 75)   262     97        0 (1452/1416-1515/1489) -71.6 W39-W41 N20-N26
      7  925 1186 (N18E03 L= 32)   443    147        0 ( 873/1131- 964/1229) -60.3 E06-E01 N14-N21
      8 1718  927 (S13W45 L= 80)   332    102        0 (1691/ 893-1741/ 980)  38.7 W42-W47 S14-S09
      9 1877  845 (S23W64 L= 99)   760    109        0 (1847/ 793-1900/ 889)  35.5 W59-W67 S25-S20
     10 1332  646 (S22W13 L= 48)  1403    301        4 (1270/ 535-1403/ 774) -77.5 W09-W17 S31-S13
     11 1703  651 (S31W44 L= 79)  1210    282        0 (1596/ 580-1793/ 733)  14.0 W33-W53 S33-S27
     12 1557  626 (S28W30 L= 65)   301     48        0 (1532/ 601-1579/ 656)  45.0 W27-W31 S29-S26
  

    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