Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1185 1688 (N38W33 L=286)  1190    250        2 (1076/1659-1315/1706) -22.4 W25-W42 N34-N41
      2 1648 1645 (N22W68 L=321)  2567    259        4 (1502/1593-1749/1700) -27.6 W55-W79 N15-N28
      3 1554 1535 (N19W50 L=303)  1810    362        2 (1408/1465-1675/1606) -49.1 W42-W56 N11-N26
      4 1436 1478 (N19W38 L=291)  1485    398        1 (1335/1346-1584/1602) -79.2 W34-W43 N07-N28
      5  375 1300 (N35E39 L=214)  2707    617        1 ( 127/1270- 677/1323) -11.8 E82-E15 N26-N41
      6  890 1103 (N11E06 L=248)   417    102        2 ( 851/1058- 921/1154) -90.0 E07-E04 N08-N14
      7 1683 1105 (S08W43 L=296)   283     58        1 (1658/1083-1711/1128)  14.0 W40-W45 S09-S07
      8  452  996 (N15E36 L=217)   262     55        0 ( 416/ 975- 474/1012) -36.9 E38-E34 N13-N16
      9 1484  931 (S13W25 L=278)   204     55        0 (1458/ 915-1516/ 950) -45.0 W23-W26 S13-S10
     10 1905  890 (S28W64 L=318)  1333    100        0 (1854/ 866-1955/ 917)  -5.4 W56-W77 S29-S25
     11 1741  813 (S27W43 L=296)   552    152        0 (1708/ 748-1775/ 871)  40.6 W39-W46 S30-S24
     12 1594  652 (S33W28 L=282)   667    145        1 (1552/ 615-1648/ 719) -48.4 W25-W33 S36-S27
     13  550  638 (S07E38 L=215)   611    192        0 ( 496/ 554- 580/ 714)  84.8 E40-E36 S12-S01
     14  472  581 (S09E46 L=207)   297     48        1 ( 454/ 554- 491/ 604)  71.6 E47-E44 S10-S07
     15  819  507 (S21E26 L=227)   585    136        1 ( 793/ 448- 843/ 577)  41.2 E30-E22 S24-S17
     16 1312  411 (S41W02 L=255)   547    103        0 (1271/ 374-1344/ 458)  33.7 E01-W05 S43-S38
     17  997  196 (S49E34 L=219)  1918    371        4 ( 857/ 123-1162/ 293)   6.6 E57-E13 S52-S45
  

    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