Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  877 1619 (N47W10 L= 71)   343     43        0 ( 857/1584- 894/1644)  59.0 W08-W12 N43-N48
      2  958 1619 (N44W16 L= 77)   324     43        0 ( 934/1591- 979/1642)  63.4 W14-W18 N41-N45
      3 1583 1448 (N14W48 L=110)   689    145        2 (1565/1379-1608/1527)  56.3 W45-W52 N09-N18
      4 1463 1372 (N14W36 L= 98)   462     98        1 (1404/1338-1499/1395) -15.9 W32-W39 N11-N15
      5  859 1280 (N26E03 L= 58)   195     87        0 ( 812/1269- 904/1292) -18.4 E06-W00 N24-N26
      6  788 1197 (N23E09 L= 52)   308     64        1 ( 756/1163- 812/1220) -76.0 E11-E07 N20-N24
      7 1758 1110 (S10W49 L=110)   414    108        0 (1720/1067-1793/1143)   6.3 W44-W53 S10-S08
      8 1670 1042 (S10W40 L=101)   429     96        1 (1640/1004-1697/1092) -74.1 W37-W41 S13-S06
      9  477  987 (N18E35 L= 26)   570    146        1 ( 433/ 943- 533/1038)  18.4 E38-E29 N15-N19
     10  230  969 (N21E58 L=  3)  1066    219        1 ( 156/ 911- 299/1029) -19.4 E68-E51 N16-N25
     11  793  934 (N08E15 L= 46)   253     83        1 ( 777/ 897- 814/ 980)  59.0 E17-E13 N05-N10
     12  159  883 (N17E70 L=351)   562     78        1 ( 141/ 843- 176/ 921) -90.0 E72-E67 N14-N19
     13 1703  766 (S28W38 L= 99)   574    164        0 (1667/ 715-1730/ 859)  53.1 W33-W41 S30-S22
     14 1332  834 (S12W12 L= 74)   264     65        0 (1307/ 802-1354/ 853)  33.7 W10-W13 S13-S10
     15 1609  783 (S24W30 L= 91)   670    171        1 (1562/ 714-1659/ 847)  29.1 W25-W34 S26-S21
     16 1198  663 (S18W00 L= 61)   921    393        1 (1038/ 637-1379/ 697) -18.4 E09-W12 S22-S13
     17 1715  595 (S41W39 L=100)   499     56        1 (1688/ 574-1737/ 626)  21.8 W35-W42 S42-S38
     18  409  462 (S12E60 L=  1)  1238    234        1 ( 384/ 373- 432/ 568)  37.4 E69-E51 S18-S05
     19  696  427 (S20E38 L= 23)   384     64        1 ( 678/ 394- 716/ 460)  68.2 E40-E35 S22-S17
  

    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