Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1035 1686 (N45W26 L= 75)   944    151        2 ( 968/1650-1113/1711) -20.0 W19-W31 N42-N48
      2 1384 1552 (N26W40 L= 89)   488    108        2 (1344/1527-1427/1569) -31.0 W37-W43 N23-N27
      3 1671 1517 (N14W60 L=109)   419     65        1 (1656/1488-1685/1552)  53.1 W57-W63 N12-N16
      4 1576 1449 (N14W48 L= 97)   778    151        0 (1511/1431-1647/1482) -13.0 W42-W55 N12-N15
      5  929 1274 (N23W01 L= 50)   283     61        1 ( 906/1248- 951/1302) -76.0 W00-W02 N21-N25
      6 1642 1227 (N00W43 L= 92)   302     57        0 (1618/1194-1661/1253) -71.6 W41-W44 S01-N02
      7 1839 1146 (S11W59 L=108)  1438    199        3 (1784/1067-1907/1207)   8.5 W50-W71 S13-S07
      8  635 1060 (N18E22 L= 27)   624    201        0 ( 563/1001- 707/1124)  18.4 E28-E16 N16-N20
      9  173  984 (N23E66 L=343)  1979    167        5 ( 144/ 906- 219/1063)  71.6 E70-E59 N18-N27
     10  340  984 (N20E46 L=  3)   995    257        3 ( 240/ 896- 390/1049) -45.0 E56-E42 N14-N25
     11  971 1010 (N07E03 L= 46)   268     70        1 ( 940/ 967- 996/1036)  45.0 E04-E01 N05-N08
     12  263  870 (N15E56 L=353)   873    166        3 ( 235/ 807- 309/ 956) -84.3 E59-E51 N10-N20
     13 1502  912 (S12W25 L= 74)   320    102        3 (1475/ 874-1537/ 950)  53.1 W23-W26 S14-S10
     14 1747  862 (S24W43 L= 92)   932    208        2 (1684/ 774-1796/ 937)  26.6 W35-W48 S27-S20
     15 1803  811 (S29W49 L= 98)   456     71        0 (1776/ 785-1830/ 842)   9.5 W45-W52 S30-S27
     16 1380  752 (S18W13 L= 62)   988    373        2 (1213/ 727-1545/ 787) -20.0 W03-W25 S22-S13
     17  508  449 (S15E51 L=358)  1605    353        2 ( 439/ 356- 558/ 598)  33.1 E65-E40 S20-S05
     18  855  459 (S22E27 L= 22)   444     78        2 ( 832/ 412- 879/ 494)  51.3 E28-E24 S24-S19
     19  637  287 (S29E53 L=356)   442     73        0 ( 615/ 255- 658/ 317)  39.8 E57-E50 S30-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