Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  900 1562 (N20E23 L=281)   351    111        0 ( 860/1516- 920/1601) -23.2 E27-E20 N17-N22
      2 1700 1246 (N27W37 L=341)  1773    433        3 (1664/1106-1759/1442) -71.6 W31-W43 N18-N39
      3 1581 1303 (N26W26 L=330)   312    104        0 (1563/1253-1608/1345)  90.0 W24-W27 N22-N29
      4 1878 1245 (N35W58 L=  3)   886    145        0 (1833/1178-1909/1297) -24.4 W51-W64 N31-N36
      5  157  992 (S29E67 L=238)  1538    167        2 (  87/ 975- 244/1015)   5.4 E86-E54 S30-S27
      6  860  825 (S22E04 L=301)   271    108        0 ( 840/ 777- 883/ 869) -59.0 E05-E01 S24-S19
      7 1729  839 (N05W49 L=354)   664    111        1 (1678/ 820-1791/ 854)  15.3 W44-W55 N02-N07
      8 1140  636 (S25W19 L=323)   553    247        0 (1028/ 613-1259/ 666)  36.9 W12-W24 S28-S19
      9 1446  482 (S23W45 L=350)  1114    276        1 (1373/ 377-1505/ 577) -70.3 W42-W49 S30-S16
     10 1568  494 (S18W54 L=358)   420     69        1 (1535/ 474-1597/ 515)  45.0 W51-W55 S19-S16
     11 1544  416 (S23W58 L=  2)  1304    229        4 (1438/ 365-1636/ 460)  42.3 W52-W65 S28-S17
     12 1417  286 (S34W61 L=  5)   393     66        0 (1400/ 263-1442/ 318) -12.5 W56-W65 S35-S32
  

    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